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

книги / Сопротивление усталости и живучесть конструкций при случайных нагрузках

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
13.19 Mб
Скачать

Рис. 13.8. Снижение предела выносливости по'мере накопления усталостных повреждений:

в — при расчете по формуле (13.26); б — по экспериментальным данным [21]

Координата точки перегиба кривой усталости по оси N (N0t v) определяется из условия

0- 1. vN0. v = (1 - v)a” , 0No = const.

(13.28)

Подставив в (13.28) выражение (13.26), получим

 

N0,V = N0(1 - v)1_am(1_A).

(13.29)

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

5-1. v = a - i . о ( l - v ) ,/m.

(13.30)

Графики функции / (v) = 1 — а_1г Ja_lt 0 при

использовании

для определения a_1(V соотношения (13.26) при a = 1 для раз­ личных значений параметра k, а также соответствующие экспери­ ментальные данные приведены на рис. 13.8. Полученные резуль­ таты позволяют сделать вывод о хорошем соответствии расчетных и экспериментальных результатов.

Из соотношения (13.27) получим следующее выражение для определения меры усталостного повреждения за один цикл на­ гружения (или, что то же самое, для скорости накопления уста­

лостных повреждений):

 

 

-^ (a /a _ i,0)M

при

aSso_liV;

V(l) = v =

при

(13.31)

О

п ° - l , v*

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

Из (13.31) следует, что скорость накопления усталостных повреждений зависит от уровня напряжений и величины нако-

142

Рис. 13.9.

Зависимости ресурса

конструкций от истории нагруже­

ния:

 

 

 

 

а — кривые

усталости

исходного

(/)

и поврежденного (2) материала; б, г

процессы

нагружения

с различным

чередованием блоков; а, д — процессы

накопления усталостных повреждений;

rti и П| -

число циклов с напряжения­

ми соответственно о* и а,

 

пленного к данному моменту

времени

усталостного

по-

вреждения. Причем эта ско­ рость не может быть представ­ лена в виде произведения двух функций, одна из которых зависит только от уровня на­ пряжений, а другая — толь­ ко от накопленного повреж­ дения. В этом случае процесс накопления усталостных по­ вреждений не будет обладать автомодельными свойствами, и результат расчета долго­ вечности при переменных амплитудах напряжений бу­ дет зависеть от истории на­ гружения [4].

Рассмотрим для примера простейший случай. Пусть

кривая усталости по мере накопления усталостных повреждений изменяется так, как показано на рис. 13.9, а, а процесс нагруже­ ния состоит из двух этапов, причем на первом этапе высокое напря­ жение 0} сменяется низким напряжением о2 (рис. 13.9, б), а на втором — низкое напряжение аг сменяется более высоким напря­ жением Ох (рис. 13.9, г). Поскольку напряжение а2 выбрано таким, что оно ниже исходного предела выносливости и выше пре­ дела выносливости, устанавливающегося после действия напря­

жений 0Х, то процессы накопления

усталостных повреждений

в рассмотренном случае (рис. 13.9, в,

д) будут существенно раз­

личаться между собой. Это предопределяет и зависимость долго­ вечности от истории нагружения.

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

143

в

V

... г

п

1 —

п

20

40

60 п/Н0

Рис. 13.10. Процессы нагружения с различным чередованием блоков на начальном и конечном этапах нагру­ жения

Рис. 13.11. Процессы накопления уста­ лостных повреждений при экспонен­ циальном законе распределения ам­ плитуд напряжений с учетом посте­ пенного снижения предела выносли­ вости

определения ожидаемого значения скорости накопления устало­ стного повреждения:

 

* =

J

v(i)f(a,s)da,

 

(13.32)

 

 

a- l,v

 

 

 

 

где f (a, s) — плотность

распределения

амплитуд

напряжений

с параметром

s; v<u — усталостное повреждение

за

один цикл

нагружения, определяемое по уравнению (13.31).

через F (v, s),

Обозначив

правую часть

уравнения

(13.32)

получим кинетическое уравнение процесса накопления устало­ стных повреждений:

v = F (V , s),

(13.33)

где черта сверху, обозначающая вычисление среднего, опущена. Решив дифференциальное уравнение (13.33), получим искомую зависимость v = v (п), которая при v = 1 определяет долговеч­

ность конструкции W*.

Подставив соотношения (13.10) при а = s и (1.2) в (13.32)

и (13.33), получим

дифференциальное уравнение, описывающее

процесс накопления

усталостных

повреждений:

 

^

 

= а~тГ +

1, (1 - v) а],

(13.34)

где п = n/N0; а = a_i,0/s; Г (а, г) — неполная гамма-функция. Результаты решения на ЭВМ уравнения (13.34) * при различ­ ных значениях параметра а и т = 4 приведены на рис. 13.11. На рис. 13.12 приведена зависимость NJN^ = f (а), где N1 — долговечность, подсчитанная без учета снижения предела вы­

* Вычисления выполнены Б . И. Петровым.

144

АГ

 

о

 

г

4

в

а

 

 

Рис.

13.12.

Зависимость

отношения

Рис. 13.13. Зависимость параметра х

долговечностей, полученных

без уче­

от а =

о_10/а

та

и

с

учетом

Nt постепенного

 

 

снижения

предела

выносливости при

носливости. Из приведенных на

экспоненциальном законе

распределе­

рис. 13.11 и 13.12 данных сле­

ния амплитуд напряжений от пара­

метра а

 

 

 

 

 

дует,

что неучет снижения пре­

водить к значительным

ошибкам

дела

выносливости может при­

в расчетах усталостной долго­

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

которых его величина

принимается постоянной и равной а_1)к=

= ха_ьо (* — искомый

параметр). Из

соотношения

(13.34) сле­

дует, что параметр х является решением следующего

алгебраи­

ческого уравнеяич:

 

 

 

N, = N*/N0 = ат1Г +

1, х, а).

(13.35)

Результаты решения уравнения

(13.35) приведены на

рис. 13.13.

 

 

 

Полученное выражение для скорости накопления усталостных повреждений (13.34) так же, как и в рассмотренном выше простей­ шем примере детерминированного нагружения (см. рис. 13.10), не может быть представлено в виде произведения двух функций, одна из которых зависит только от накопленного к данному мо­ менту времени повреждения v, а другая — только от параметра интенсивности нагружения s. В этом случае процесс накопления усталостных повреждений также не будет автомодельным и пра­ вило суммирования относительных долговечностей также не будет выполняться. При этом результат расчета усталостной долговеч­ ности будет зависеть от истории нагружения. Пусть, например, эта история состоит в том, что элемент конструкции нагружается

втечение пг циклов с параметром интенсивности воздействий slt

втечение циклов с параметром интенсивности воздействия % и т. д., а принцип сильного перемешивания для режимов нагруже­ ния не выполняется (рис. 13.14, а). Тогда за k блоков нагружения суммарное усталостное повреждение

v = v (tii, si) 4

v (rtai

4 п2>sa) 4

v (пз2 4

п3, S3) -|-

+

• • * 4

V (Пк,к-1 4

nk> sft)>

(13.36)

10 Гусев А. С.

145

 

 

*

 

• •

т Пг ж

л *

т л* 1

 

Рис. 13.14. Изменения параметра s от числа циклов нагружения п (а) и соответствующий процесс накопления усталостных повреждений (б)

Рис. 13.15. Принцип корректировки правила линейного суммирования усталостных повреждений

где v = v

(я, s) — закономерность

накопления усталостных

по­

вреждений

при уровне

воздействия

s; пи г_х (t =

1, 2...... k)

решения уравнений

 

 

 

 

 

 

v (п1. 1-1* sl) — V (Лxi Si) +

v (л2х +

«г* s*) +

 

 

~h • • •

+ v (rtj-x, 2+

S;_2).

(13.37)

Принцип суммирования усталостных повреждений, выражен­ ный соотношениями (13.36) и (13.37), поясняется графиками,

приведенными

на рис. 13.14. Для примера примем, что Ьг = 8,

Лх = 200, Ь2 =

6, л2 =

22. Тогда при прямом порядке приложения

нагрузок v — 0,5, а

при смене порядка приложения нагрузок

v = 0,2.

 

 

Экспериментально установлено, что при нерегулярных процес­ сах нагружения необходима определенная корректировка правила линейного суммирования усталостных повреждений, выраженного соотношением (13.6). Предельное значение накопленного устало­ стного повреждения оказывается меньше единицы, и чем больше процесс нагружения отличается от простого гармонического нагружения, тем больше необходимо снижать это предельное значение усталостного повреждения [21 ]. Корректированное пре­ дельное значение усталостного повреждения определяется по виду интегральной функции распределения амплитуд напряжений (рис. 13.15):

% = = “ i/K - a-i)> (13-38)

146

где щ , о>2— площади фигур,

заштрихованные на рис.

13.15;

а* — максимальное

значение

амплитуд напряжений;

o_j —

расчетное значение

предела выносливости.

 

Вычислив coj и подставив его в формулу (13.38), получим

 

Ук=

,

(13.39)

 

 

а * — a - i

 

где а — среднее значение амплитуд напряжений.

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

ст* = s |/ 2In N,

где N — число циклов нагружения, определяемое ориентиро­ вочно по ожидаемой долговечности; s — среднее квадратическое отклонение действующих напряжений.

В случае релеевского закона распределения амплитуд напря­ жений, описываемого соотношением (14.2), получаем (при а_г = О и JV = 2- 10е) vft « 0,23.

§14. Расчеты при стационарных случайных колебаниях

Как и при потоках случайных импульсных воздействий (см. § 13), при случайных колебаниях будем учитывать возможность разрушения конструкции в момент превышения процессом на­ гружения <т (t) опасного уровня напряжений сг* и в момент, когда накопленное усталостное повреждение достигнет опасного уровня

v* = 1.

Решение первой задачи определяется вероятностью события, что действующие в конструкции напряжения а (/) при заданном времени нагружения t ни разу не превысят опасный уровень а*. Эта вероятность для гауссовских стационарных процессов опре­ деляется по формуле (11.5), которую можно представить в виде

соотношения для определения

меры

надежности

конструкции:

ехр

а2 \

 

 

 

4 )

при

 

(14.1)

 

0

 

при

t > t+y

 

где Я (0 — мера надежности (вероятность неразрушения) для момента времени t; ? — средний период цикла нагружения;

= I exp (<J«/(2so)] —

характерное время,

ю*

147

Расчет усталостной долговечности конструкции при случайных колебаниях рассмотрим раздельно для двух основных случаев: для процессов с простой и со сложной структурой (рис. 14.1, а, б).

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

/(°) = -Згехр(—- ^ - ) ,

(14.2)

где sa — дисперсия процесса нагружения.

Пусть кривая усталости описывается уравнением (13.11). Тогда, подставив (14.2) и (13.11) в формулы (13.8) и (13.9), полу­ чим соответственно выражения для среднего значения и устало­ стного повреждения за один цикл нагружения:

v = 2т/2С~'Г (т/2 -f 1)sm;

(14.3)

4 = 2тС~2[Г (т + 1) - Г2/2 + 1)] s2m.

(14.4)

Коэффициент вариации распределения усталостного поврежде­

ния за один цикл нагружения

 

6V= [Г (т + 1)/Г2(m/2 + 1) - 1],/2.

(14.5)

По сравнению с экспоненциальным распределением распре­ деление Релея приводит к меньшим значениям коэффициента вариации 6V.

Подставив (14.3) и (14.5) в формулы (9.49) и (9.51), получим соответственно следующие выражения для определения среднего значения и коэффициента вариации распределения усталостной долговечности (без учета рассеяния интервала времени между нагружениями. и без учета корреляционной зависимости между повреждениями в циклах нагружения):

Т = 2_т/2(Г (т/2

+ I)]-1 ?Cs_m;

(14.6)

бг = 6v/ / n ,

(14.7)

где 6V— величина, определяемая

по (14.5); п — число

циклов

нагружения.

 

пред­

Графики функций (14.7) при различных значениях т

ставлены на рис. 14.2.

 

 

Рассмотрим теперь случай, когда кривая усталости описы­ вается уравнением (1.2). Подставив (1.2) и (14.2) в формулы (13.8) и (13.9), получим соответственно выражения для среднего значе-

148

<r

v W W v ■,

Рис. 14.1. Случайные процессы с про­

Рис. 14.2.

Зависимость коэффициента

стой (а) и сложной (б) структурой

вариации

распределения долговечно»

 

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

ния и дисперсии усталостного повреждения за один цикл нагру­ жения:

(14.8)

где Г (п,

а) — неполная гамма-функция.

 

 

Коэффициент вариации распределения усталостного поврежде­

ния за один цикл нагружения

 

 

 

6V=

J2L +

1

^

(14.9)

 

2 г

’’

2s*

 

Подставив (14.8) и (14.9) в формулы (9.49) и (9.51), получим соответственно следующие выражения для среднего значения и коэффициента вариации распределения усталостной долговеч­ ности:

Т =

* 0 ° - 1

(14.10)

[- f-+ 1. oil/(2s2)]

2m/V T

 

-- б у / / I 9

где п — число нагружений; 6V— величина, определяемая по (14.9). Аналогично проводятся расчеты долговечности и при других законах распределения амплитуд напряжений и уравнений кри­

вых усталости.

Если в расчете усталостной долговечности учитывается эффект снижения предела выносливости по мере накопления усталостных

149

Рис. 14.3. Процессы накопления усталостных повреждений при релеевеком законе распределения ам­ плитуд напряжений с учетом по­ степенного снижения предела вы­ носливости:

а — при т = 4; б — при т = 6; в — при т = 8

повреждении, можно воспользоваться методикой, изложен­ ной в § 13. Если, например, распределение амплитуд на­ пряжений подчиняется зако­ ну Релея с плотностью (14.2), то тогда кинетическое урав­ нение накопления усталост­

 

 

ного разрушения (13.33) при­

 

 

нимает следующий вид:

 

 

dv _ h- т / Г М + 2

 

 

da ~ 0

1 [

2 ’

 

 

Ь(\ -

V)2] ,

(14.11)

в)

' ’

где п = n/N0; b =<xiI>0/(2s2);

Г (а, г) — неполная

гамма-

 

 

функция.

 

 

Результаты решения на ЭВМ * уравнения (14.11) при различ­ ных значениях параметров b и т приведены на рис. 14.3.

На

рис. 14.4 приведены

графики функций NJN* = f (Ь),

где N,

и AL — долговечности,

подсчитанные соответственно без

учета и с учетом снижения предела выносливости.

Эффект постепенного снижения предела выносливости можно также учесть в обычных методиках расчета путем соответству­ ющего подбора для него постоянного расчетного значения o_lt * = = хо_(>0 (х <; 1), при котором результаты расчета по обеим мето­ дикам будут совпадать. В соответствии с формулой (14.8) для определения параметра х получаем следующее алгебраическое уравнение:

М* = N JN 0 = 6- /2 [ ^ ~ , *?Ь] .

(14.12)

Результаты решения на ЭВМ уравнения (14.12)

приведены

на рис. 14.5.

 

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

* Вычисления выполнены Б. И. Петровым.

150

Рис. 14.4. Зависимости отношения дол­ говечностей, получаемых без учета (NJ и с учетом (N+) постепенного снижения предела выносливости при релеевском законе распределения ам­ плитуд напряжений от параметра b

Рис. 14.5. Зависимость параметра к от b = ai/(2s*)

на долговечность информации о характеристиках сопротивления усталости материалов, полученных в условиях гармонического нагружения в виде кривых и поверхностей усталости (см. § 1), становится затруднительным. Для проведения расчетов необхо­ димо привести процессы со сложной структурой к эквивалентным по повреждающему действию процессам с простой структурой. Полной эквивалентности при этом достичь не удается, и рекомен­ дуемые методы приведения в значительной степени основаны на опыте их практического использо­ вания. В настоящее время приме­ няется более десятка различных ме­ тодов приведения. Некоторые из них приведены ниже.

 

М е т о д м а к с и м у м о в . За

амплитуды

напряжений

принимают­

ся

все положительные

максимумы,

а

за

расчетную

частоту — среднее

число

положительных

 

максимумов

в

единицу

времени.

Схематизация

процесса нагружения

в этом случае

заключается в

условном

переносе

положительных

минимумов и отри­

цательных

максимумов

на нулевую

(среднюю)

линию (точки

I 9 и 2',

рис. 14.6, а, б).

При этом

получаем

наибольшее усталостное повреждение и, следовательно, наименьшее для

Рис. 14.6. Приведение исходного процесса сложной структуры (а) к процессу с про­ стой структурой (б) по методу максимумов

151