Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Многоцикловая усталость при переменных амплитудах нагружения

..pdf
Скачиваний:
6
Добавлен:
13.11.2023
Размер:
22.76 Mб
Скачать

Рис. 26. К схематизации методом потока дождя. Зависимость о — е

(цифры означают то же, что и на рис. 25).

Рис. 27. К схематизации нагруз­ ки методом полных циклов.

ра записи. Штриховыми ли­ ниями показаны те размахи, которые были исключены на предыдущем просмотре. Оче­ видно, что при схематизации методом полных циклов мож­ но учитывать среднее нап­ ряжение цикла [81].

При сравнении схематиза­ ции по методу полных цик­

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

Для большинства историй нагружения метод «потока дождя» дает результаты, очень близкие к методу полных циклов, хотя сами алгоритмы существенно отличаются. Каж­ дый из этих алгоритмов не лишен недостатков. Так, в изло­ женном алгоритме метода полных циклов необходимо иметь в памяти ЭВМ все экстремумы нагрузки и песколько раз (по числу разрядов) просматривать весь массив экстрему­ мов. Метод потока дождя достаточно сложен для программи­ рования, и, хотя нет необходимости помещать весь массив экстремумов в оперативную память ЭВМ, число вычислений при неоптимальпой реализации этого метода велико.

В данном параграфе излагается алгоритм схематизации нагрузки, который объединяет преимущества (с точки

 

 

 

 

зрения

простоты программи­

 

 

 

 

рования и уменьшения време­

 

 

 

 

ни вычислений)

 

методов

по­

 

 

 

 

тока

дождя

и

полных цик­

 

 

 

 

лов [88]. Этот алгоритм яв­

 

 

 

 

ляется

усовершенствованием

 

 

 

 

алгоритмов, предложенных в

 

 

 

 

работах

[178, 241].

 

 

 

 

 

Основной

принцип оста­

 

 

 

 

ется

прежним — выделение,

 

 

 

 

подсчет

и дальнейшее иск­

 

 

 

 

лючение из

рассмотрения

в

 

 

 

 

первую очередь меньших раз-

 

 

 

 

махов. Алгоритм

схематиза­

Рис. 28.

К алгоритму схематиза­

ции

состоит

из

шагов.

На

ции методом полных циклов;

каждом

шаге

схематизации

а а б — случай гАд >

в и г — слу­

рассматривают

три последо­

чай гАв ^ ТВС‘

 

 

вательных

экстремума

А ,

 

 

 

 

и хс.

 

 

В и С с ординатами х а ,

хв

На первом шаге схематизации

А , В и

 

С — первые

три экстремума.

Точки А , В,

С образуют

размахи гав

=

Jха — хв | и гвс = | хв хс |. Возможны два случая: гАв >

>

гвс и гвс ^ гав- В первом случае

(рис. 28, а) цикл нагру­

жения не выделяется, поскольку следующий по записи размах может быть меньше размаха гвс* Переходим к рассмот­ рению следующей тройки экстремумов по правилу: В обозна­ чаем А, С обозначаем В, за экстремум С принимаем следую­ щий по записи (рис. 28, б). Затем переходим к следующему

шагу схематизации.

Во втором случае (рис. 28, в) считаем, что экстремумы А и В образуют цикл нагружения. Отметим, что размах, предшествующий гав (если он имеется), больше гАв, иначе

он был бы исключен на предыдущем шаге. Амплитуда цик-

1

ХА~Ъ~ХВ

Г\

ла А — В ал = ^

Гав»среднее значение ат = — ^—

. Экстре­

мумы А та. В исключаются из дальнейшего рассмотрения. В качестве экстремумов А и В берем два экстремума, про­ пущенных ранее (рис. 28, г). Если число пропущенных экст­ ремумов меньше двух, в качестве А , В и С берем первые

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

по записи экстремум. В этом случае оставшиеся (уже про­ смотренные и пропущенные) экстремумы образуют последо­

ватеЛьность убывающих размахов (па рис. 29 экстремумы 11—14—15—16). Для учета повреждений и от этих размахов объединяем оставшиеся экстремумы в пары, например 11— 14 и 1516, и соответствующие размади подсчитываем как

циклы. На этом схематизация заканчивается. Преимущество изложенного алгоритма перед алгоритмом

поразрядного исключения метода полных циклов в том, что отпадает необходимость в хранении в оперативной памяти ЭВМ всего массива экстремумов. Достаточно только иметь массив небольшой длины для хранения экстремумов, пропу­ щенных при схематизации. Например, при обработке 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)

 

 

 

 

 

СИ

 

 

К

 

 

 

 

 

 

 

С

 

 

а число превышений уровня о0 за время Т

 

 

ЛГЫ г) =

Т

о с„

 

(°о — Е

[ст])2

(2.49)

2п

о

ехр

 

2о:ск

 

 

 

ск

 

 

 

 

Последняя формула впервые получена Райсом [255]. Ве­ личина

(2.50)

называется эффективным периодом процесса <т(t). Для гар­ монического процесса Те совпадает с его частотой, для уз­

кополосного процесса Те совпадает с центральной частотой

процесса.

Число

выбросов выше среднего уровня Е [а]

обозначим

из

выражения (2.49) следует формула для N 0:

т_

(2.51)

 

Если обозначить моменты спектральной плотности че­

рез 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] аналитическое выражение для

Соседние файлы в папке книги