книги / Многоцикловая усталость при переменных амплитудах нагружения
..pdfРис. 26. К схематизации методом потока дождя. Зависимость о — е
(цифры означают то же, что и на рис. 25).
Рис. 27. К схематизации нагруз ки методом полных циклов.
ра записи. Штриховыми ли ниями показаны те размахи, которые были исключены на предыдущем просмотре. Оче видно, что при схематизации методом полных циклов мож но учитывать среднее нап ряжение цикла [81].
При сравнении схематиза ции по методу полных цик
лов, максимумов и размахов оказывается, что метод полных циклов дает промежуточный результат по отношению к ме тодам максимумов и размахов.
Для большинства историй нагружения метод «потока дождя» дает результаты, очень близкие к методу полных циклов, хотя сами алгоритмы существенно отличаются. Каж дый из этих алгоритмов не лишен недостатков. Так, в изло женном алгоритме метода полных циклов необходимо иметь в памяти ЭВМ все экстремумы нагрузки и песколько раз (по числу разрядов) просматривать весь массив экстрему мов. Метод потока дождя достаточно сложен для программи рования, и, хотя нет необходимости помещать весь массив экстремумов в оперативную память ЭВМ, число вычислений при неоптимальпой реализации этого метода велико.
В данном параграфе излагается алгоритм схематизации нагрузки, который объединяет преимущества (с точки
|
|
|
|
зрения |
простоты программи |
|||||
|
|
|
|
рования и уменьшения време |
||||||
|
|
|
|
ни вычислений) |
|
методов |
по |
|||
|
|
|
|
тока |
дождя |
и |
полных цик |
|||
|
|
|
|
лов [88]. Этот алгоритм яв |
||||||
|
|
|
|
ляется |
усовершенствованием |
|||||
|
|
|
|
алгоритмов, предложенных в |
||||||
|
|
|
|
работах |
[178, 241]. |
|
||||
|
|
|
|
Основной |
принцип оста |
|||||
|
|
|
|
ется |
прежним — выделение, |
|||||
|
|
|
|
подсчет |
и дальнейшее иск |
|||||
|
|
|
|
лючение из |
рассмотрения |
в |
||||
|
|
|
|
первую очередь меньших раз- |
||||||
|
|
|
|
махов. Алгоритм |
схематиза |
|||||
Рис. 28. |
К алгоритму схематиза |
ции |
состоит |
из |
шагов. |
На |
||||
ции методом полных циклов; |
каждом |
шаге |
схематизации |
|||||||
а а б — случай гАд > |
в и г — слу |
рассматривают |
три последо |
|||||||
чай гАв ^ ТВС‘ |
|
|||||||||
|
вательных |
экстремума |
А , |
|||||||
|
|
|
|
|||||||
и хс. |
|
|
В и С с ординатами х а , |
хв |
||||||
На первом шаге схематизации |
А , В и |
|
С — первые |
|||||||
три экстремума. |
Точки А , В, |
С образуют |
размахи гав — |
|||||||
= |
Jха — хв | и гвс = | хв — хс |. Возможны два случая: гАв > |
|||||||||
> |
гвс и гвс ^ гав- В первом случае |
(рис. 28, а) цикл нагру |
жения не выделяется, поскольку следующий по записи размах может быть меньше размаха гвс* Переходим к рассмот рению следующей тройки экстремумов по правилу: В обозна чаем А, С обозначаем В, за экстремум С принимаем следую щий по записи (рис. 28, б). Затем переходим к следующему
шагу схематизации.
Во втором случае (рис. 28, в) считаем, что экстремумы А и В образуют цикл нагружения. Отметим, что размах, предшествующий гав (если он имеется), больше гАв, иначе
он был бы исключен на предыдущем шаге. Амплитуда цик-
1 |
ХА~Ъ~ХВ |
Г\ |
ла А — В ал = ^ |
Гав»среднее значение ат = — ^— |
. Экстре |
мумы А та. В исключаются из дальнейшего рассмотрения. В качестве экстремумов А и В берем два экстремума, про пущенных ранее (рис. 28, г). Если число пропущенных экст ремумов меньше двух, в качестве А , В и С берем первые
экстремумы, оставшиеся в записи. На рис. 29 изображены шаги схематизации нагрузки, изображенной на рис. 27. Анализ алгоритма показывает, что после определенного числа шагов возникает ситуация, когда запись кончится и в качестве экстремума С нельзя будет брать следующий
по записи экстремум. В этом случае оставшиеся (уже про смотренные и пропущенные) экстремумы образуют последо
ватеЛьность убывающих размахов (па рис. 29 экстремумы 11—14—15—16). Для учета повреждений и от этих размахов объединяем оставшиеся экстремумы в пары, например 11— 14 и 15—16, и соответствующие размади подсчитываем как
циклы. На этом схематизация заканчивается. Преимущество изложенного алгоритма перед алгоритмом
поразрядного исключения метода полных циклов в том, что отпадает необходимость в хранении в оперативной памяти ЭВМ всего массива экстремумов. Достаточно только иметь массив небольшой длины для хранения экстремумов, пропу щенных при схематизации. Например, при обработке 200 записей узкополосного процесса нагружения, в каждой из которых содержалось около 800 экстремумов, наибольшее число пропущенных экстремумов (при обработке одной за писи) равнялось 50. Для широкополосных процессов это значение оказалось еще меньше.
Экономия времени вычислений зависит от типа нагрузки. Для процесса нагружения, использованного для сравнения алгоритмов, расчет с помощью двух программ показал умень шение времени счета в два раза. Сравнение с методом потока дождя показывает, что изложенный алгоритм отличается простотой и позволяет сэкономить память ЭВМ и уменьшить
Рис. 29. Последовательные шаги при схематизации нагрузки, изобра женной на рис. 27 (цифры означают выделенные размахи, прочерк — от сутствие размахов, подлежащих выделению).
время обработки записи нагрузки. По сравнению С алгорит мами, приведенными в работах [178, 2411, предложенный алгоритм также обладает преимуществом, которое заклю чается в том, что на каждом шаге схематизации рассматри ваются три, а не четыре, как в этих работах, последователь ных экстремума. В работе [247] приведен алгоритм подсчета петель гистерезиса, аналогичный изложенному.
Последовательность циклов, выделенных при схемати зации, обычно преобразуется в распределение эквивалент ных амплитуд. Выделенные циклы характеризуются амп литудой оа и средним значением ат. По значениям оа и ат находится эквивалентная амплитуда а цикла с нулевым средним значением так, чтобы цикл с амплитудой о был
эквивалентным по повреждающему действию асимметрич ному циклу. При этом могут применяться следующие фор мулы:
а = а0 |
(2.41) |
(здесь среднее значение цикла не учитывается);
о = оа + фа™ |
(2.42) |
(ф — константа материала, ф = 0,1 — 0,3 [81]);
а = |
оа {аа + ат). |
(2.43) |
Асимметрия цикла в ряде случаев учитывается путем при ведения к эквивалентным напряжениям отнулевого цикла по модифицированной формуле Одинга [127]
°о = |
2°п |
(2.44) |
|
° т - ° а |
|||
1 — |
|
+ аа )
где — константа материала, при отсутствии эксперимен тальных данных можно полагать ха — 0,5.
Распределение эквивалентных амплитуд может быть пред ставлено в виде интегральных эмпирических функций рас пределения, в виде гистограмм, а также в виде накопленных частот превышения заданных уровней за определенный пе риод эксплуатации (рис. 30) [267]. При построении гисто грамм диапазон изменения напряжений разбивается на ин тервалы равной длины (<jj, CTi+i) и подсчитывается число максимумов щ, попадающих в этот интервал. Величина
(Ji -f-
пi откладывается по оси ординат в точке о = ------2 -. При
представлении распределения нагрузки в виде накопленных частот по оси ординат откладывается нагрузка Р, а по оси абсцисс — число нагрузок, превышающих этот уровень за
Рис. 30. Обработка записи реального нагружения трансмиссии тракто ра [267]:
а — гистограмма максимумов (j ), минимумов (2) и размахав полных циклов |
(3); |
б — накопление частоты превышения уровней максимумов (J) и минимумов |
(2). |
определенный цикл эксплуатации (например, за определен
ное |
число полетов, за определенный пробег в километрах |
и т. |
п.). |
Для многих конструкций, в частности авиационных и автомобильных, можно допустить существование линейной связи между напряжением и логарифмом накопленного ко личества случаев появления этого напряжения (рис. 31). Такая история деформирования описывается двумя пара метрами: атах и наклоном fe, что аналитически выражается
экспоненциальной |
зависимостью |
|
|
N = Ю<атах_ст)/Л |
(2.45) |
между количеством |
превышений N уровня а |
и самим о. |
В качестве примера на рис. 30 показаны гистограммы максимумов, минимумов нагрузки на трансмиссию трактора
блш/б
пои,
а размахов полных циклов. Кроме того, показаны накоп ленные частоты распределения максимумов и минимумов нагрузки. В большинстве случаев можно принять, что рас пределения максимумов и минимумов нагрузок примерно симметричны относительно горизонтальной линии, про ходящей через точку пересечения этих распределений (см. рис. 31). Нагрузка, соответствующая этой горизонтальной
линии, представляет собой среднюю нагрузку |
процесса. |
В том случае, если необходимо представить |
результаты |
схематизации с учетом среднего значения выделенных цик лов, диапазон изменения амплитуд циклов и средних зна чений разбивается на определенное число разрядов (при мерно 20—30). Результаты схематизации оформляются в виде матрицы, строки которой соответствуют разрядам ампли тудных значений, а столбцы — разрядам средних. Столбец матрицы с номером ; представляет собой гистограмму амп литуд циклов, среднее значение которых попадает в /-й раз ряд соответствующего разбиения. Элемент матрицы озна чает число циклов, амплитуда которых попадает в г-й раз ряд, а среднее значение — в ;*-й разряд. Такая матрица называется корреляционной таблицей двухпараметрической схематизации [128].
3. МЕТОДЫ ТЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ
ДЛЯ ВЕРОЯТНОСТНОГО ОПИСАНИЯ ЭКСТРЕМУМОВ СЛУЧАЙНОЙ НАГРУЗКИ
Применительно к проблеме усталости методы теории слу чайных процессов прилагаются для оценки числа циклов (максимумов) случайной нагрузки, а также для определения теоретических распределений характеристик циклов, вы деляемых в режиме нагружения согласно существующим методам схематизации.
С помощью теории выбросов случайных функций [14] можно оценить среднее в единицу времени число выбросов N (о010 выше уровня а0 и среднее число превышений N (<j01Т) уровня о0 за время Т. Для определения этих величин в об
щем |
случае необходимо знать |
совместную плотность |
р (а, |
||
о, t) |
распределения о |
(t) и о |
(t). |
о (£) и |
а (t) |
Для стационарного |
нормального процесса |
||||
|
|
• |
• |
■ |
|
некоррелированы и р (о, о, t) — р(а) р (а), причем р(о) также
нормально распределена с нулевым средним:
(2.46)
« |
|
|
|
отклонение |
|
• |
||
где aCf<— среднекпадратическое |
процесса o(t), |
|||||||
определяемое по формуле |
|
|
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
& = |
4я» J f W |
(/) df, |
|
|
(2.47) |
|||
|
|
|
О |
|
|
|
|
|
которая следует из выражений (2.23) и (2.17). |
|
|||||||
Согласно теории |
выбросов |
среднее |
число превышений |
|||||
N (а0) уровня ст0 в единицу времени для нормального стацио |
||||||||
нарного процесса определяется |
по формуле |
|
||||||
N (о0) = |
2па, |
|
ехр |
|
(°о — Е fq|)2 |
(2.48) |
||
|
|
2о |
|
|
||||
|
СИ |
|
|
К |
|
|
||
|
|
|
|
|
С |
|
|
|
а число превышений уровня о0 за время Т |
|
|
||||||
ЛГЫ г) = |
Т |
о с„ |
|
(°о — Е |
[ст])2 |
(2.49) |
||
2п |
о |
ехр |
|
2о:ск |
|
|||
|
|
ск |
|
|
|
|
Последняя формула впервые получена Райсом [255]. Ве личина
(2.50)
называется эффективным периодом процесса <т(t). Для гар монического процесса Те совпадает с его частотой, для уз
кополосного процесса Те совпадает с центральной частотой |
||
процесса. |
Число |
выбросов выше среднего уровня Е [а] |
обозначим |
из |
выражения (2.49) следует формула для N 0: |
т_ |
(2.51) |
|
2л |
||
|
||
Если обозначить моменты спектральной плотности че |
||
рез Mh> |
|
|
M h = $ f W (D df, |
(2 52) |
|
0 |
|
то с учетом выражений (2.51) и (2.47) формулу для числа выбросов выше Е [а] можно записать так:
N a = T ] / Ж . |
(2.53) |
Величину N 0 можно интерпретировать как число пересе чений уровня Е [о] с положительной (или отрицательной) производной. Для процесса с нулевым средним Е [о] — 0
и — половина числа нулей за время 7\
«о/3<
0,6
М
Аналогично можно получить формулу для числа максимумов процесса за время Т :
N м = т У |
м л |
(2.54) |
м9 |
||
Л |
|
|
|
10° |
|
J |
/О4 W, циклы |
Для |
широкополосного |
процес |
||||
Рис. 32. |
Связь пересечений |
са нагружения характерно |
нали |
||||||||
чие |
максимумов, |
меньших средне |
|||||||||
N |
(о0 | Т ) и |
уровня о0/о тах |
|||||||||
го уровня, и минимумов, больших |
|||||||||||
согласно |
|
формуле |
(2.55) |
||||||||
|
среднего. Поскольку число макси |
||||||||||
для нормального |
процесса |
||||||||||
(ffmax = |
|
|
|
|
мумов, больших Е [ст], минус число |
||||||
|
|
|
|
минимумов, также больших Е 1а], |
|||||||
равно |
N 0 (что |
следует |
|||||||||
из |
геометрических соображений), |
||||||||||
то |
число |
максимумов, |
больших |
Е [а], |
равно полусумме |
||||||
числа всех максимумов и всех пересечений среднего |
уровня |
||||||||||
с |
положительной |
производной. |
|
|
|
||||||
|
Формулу (2.49) |
для |
числа |
пересечений |
уровня а0 с уче |
том выражения (2.51) можно представить в виде (в случае процесса с нулевым средним Е (о] = 0)
N fao IТ) = N 0exр |
4 |
(2.55) |
|
2D2 |
|||
|
|
||
|
W CK |
|
Соответствующая кривая для N 0= 10е показана на рис. 32.
Отметим, что формула (2.49) представляет собой теорети ческое описание результатов схематизации методом пере сечений заданного уровня для стационарного нормального процесса. С помощью формулы (2.55) можно проверять нор мальность процесса путем сопоставления эксперименталь ной кривой накопленных частот с прямой линией в коор
динатах lg N (а01Т) — оо-
Методы теории выбросов [14] могут быть применены для оценки распределения выбросов за время, равное эффектив ному периоду процесса Те. Тогда распределение выбросов
может быть интерпретировано как распределение амплитуды выброса между соседними пересечениями среднего уровня с положительной производной, и исходный процесс может быть заменен узкополосным процессом с периодом Те и со
ответствующим |
распределением амплитуд. |
|
||
Пусть Р (а > |
а01Те) — вероятность |
превышения уровня |
||
а0 за время Те\ N (о01Те) — среднее число превышений уров- |
||||
н яа0 за время Те, а Р х(а0| Те), |
Р2(°о1 те ) , Р ц(а0| Те) — ве |
|||
роятности одно-, дву-, и ft-кратных превышений |
уровня а0 |
|||
за время Т*. Предполагается, |
что а0 >■ Е [а]. Очевидно, что |
|||
N (JЕ [а]| Те) = 1 [см. формулу |
(2.49)], |
поэтому |
можно счи |
тать, что Рj (а0| Т„) |
Рг(°о I Т,)... Величины |
Р (а >• «т0| Те) |
||||||||
и /V (а0| Тё) |
можно |
следующим |
|
образом |
выразить через |
|||||
P h ( o Q\ Те): |
|
|
|
|
|
|
|
|
|
|
|
|
Р (о > |
а01Тс) = 2 |
|
Ph К 1 Те); |
|
|
|||
|
|
|
|
k = i |
|
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
N (<г01Г.) =■ 2 |
|
(<7„ I Го). |
|
|
||||
Отсюда |
|
|
|
Л=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Р (а > |
а„1Т‘) - |
N (а, | Т,) = |
f |
(1 - |
|
4) Ph(а„1Т,). |
||||
|
|
|
|
|
Л = 2 |
|
|
|
|
|
Поскольку |
Р),(а0\Те) |
(при |
/с ^ |
2) |
можно считать |
малыми |
||||
по сравнению |
с Р, (ст01Те), |
приближенно |
полагаем, |
что |
||||||
Р (а > |
ст01Те) а * N (о01Те) |
|
(при |
ст0> Я [о]). |
|
|||||
Интегральная |
функция распределения |
для |
перегрузок |
|||||||
F (а01Тл) = 1 - |
Р (а > а01Тш) « 1 - |
ЛГ (а01Г*). |
|
Для стационарного нормального процесса из последней формулы можно получить следующую оценку для плотно сти распределения амплитуд выбросов за время, равное
эффективному периоду |
процесса: |
|
|
|
ехр |
(o0-l[q j)2 |
а0> Я (а ] |
|
|
||
Р Ы = |
°ск |
2<4 |
а0 < Е [а] |
|
|
О, |
Таким образом, согласно теории выбросов исходный процесс заменяется приближенно эквивалентным по по вреждающему значению, представляющему собой узкопо лосный процесс с периодом Те (формула (2.50)), средним зна чением Е [о] и распределением Рэлея амплитуд <та5
п
Важной характеристикой случайной нагрузки является коэффициент широкополоспости |5 (или обратная величина
v = -|р называемая коэффициентом нерегулярности),
определяемая как отношение числа максимумов нагрузки к числу пересечений среднего уровня с положительной про изводной:
о |
'ум |
Г м уГ 0 |
(2.56) |
|
р |
N0 |
м % |
||
|
Величина ]3 может изменяться от 1 до бесконечности,
но для реальных процессов (3 редко превышает 4. Например,
для спектра с постоянной |
мощностью в полосе частот от F1 |
||||
до F2 коэффициент широкополосности |
|||||
3 |
Vк4+ А3+ к'1+ к + 1 |
> —— л* 1,34 при к — оо« |
|||
|
|
|
|
|
У 5 |
где к = |
F2IFx. |
|
|
|
|
Для процесса, равного сумме двух узкополосных про |
|||||
цессов с |
центральными |
частотами |
Д и /2 и среднеквадрати |
||
ческими |
отклонениями |
оСк1 и оС1<2 |
[14], |
||
|
о _ |
/ |
(1 + |
V-) (1 + |
|
|
Н |
|
i |
_ J у202 |
Важность определения (3 обусловлена тем, что в ряде работ высказано предположение [57, 59, 134, 166, 186], что усталостная долговечность полностью определяется па раметрами оС1{ и |3.
Величина сгСк не является достаточной для описания слу чайной нагрузки, поскольку, как показывает анализ экспе риментальных данных, число максимумов до разрушения при одинаковых оск зависит от формы спектральной плот ности, т. е. один параметр ос„ не определяет однозначно определенную таким образом долговечность [17, 166, 181, 186, 253, 265 и др.1.
Следует отметить, что влияние среднего значения про цесса нагружения на долговечность не является обычно (за исключением, по-видимому, работы [188]) предметом спе циального исследования. Это связано с тем, что в расчетах долговечности при нерегулярном нагружении можно исполь зовать кривую усталости при гармоническом нагружении со средним значением, равным среднему значению случай ной нагрузки. Тем самым задача упрощается и сводится
только к учету |
нерегулярности |
нагружения |
без дополни |
тельного учета |
асимметрии. |
|
|
Выделение |
параметра (3 как |
одной из |
существенных |
характеристик нагрузки обусловлено тем, что ряд важных теоретических распределений, описывающих статистику экс тремумов нормального стационарного процесса, определя ется величинами ос1< и [3.
В частности, в работах [37, 60] для оценок долговечности предложено применять схематизацию методом максимумов, используя известное [15] аналитическое выражение для