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

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

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

Рис. 14.7. Приведение исходного

процесса сложной структуры (а)

к процессу

с простой структурой (б) по методу размахов

 

всех методов схематизации значение усталостной

долговеч­

ности.

За расчетную амплитуду напряже­

М е т о д р а з м а х о в .

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

по методу

размахов показана

на рис. 14.7.

М е т о

д р а з м а х о в с

у ч е т о м с р е д н е г о . В отли­

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

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

М е т о д п о л н ы х ц и к л о в с у ч е т о м с р е д н е г о . Каждый цикл нагружения характеризуется парой напряжений: амплитудным и средним. Эквивалентная расчетная амплитуда напряжений определяется по диаграммам предёльных амплитуд (см. § 1).

152

Р а с ч е т у с т а л о с т н о й д о л г о в е ч н о с т и по м е т о д у м а к с и м у м о в . В этом методе расчета плотность распределения амплитуд напряжений принимается в следующем виде:

О при

а < 0;

Л

(14.13)

.

,

,

.

cfшах (<г) при

а >

О,

 

где /шах (<т) — плотность{распределения максимумов, определя­

емая по формуле (11.6); с — коэффициент нормировки,

учитыва­

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

В соответствии с соотношением (11.6) имеем

 

\—1

Г — | j / max (®) d o

= 2Л/(*+1),

где k — отношение среднего числа экстремумов к среднему числу нулей.

В области больших значений максимумов (при о -*■ оо), кото­ рые являются определяющими в расчетах усталостной долговеч­

ности, распределение (11.6) принимает вид

 

/.«<<г)=ТЗгехр

04.14)

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

Подставив соотношение (14.14) в выражение (14.13), получим плотность распределения амплитуд напряжений по методу макси­

мумов

 

(

0

при

о <

0;

 

 

 

 

 

/<*) =

2ст

 

 

/

ста \

(14.15)

 

 

 

(k+ l)sa еХр (

)

 

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

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

максимумов)

 

 

где

 

 

«+ =

Prtmax.

 

 

(14.16)

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р =

1/max (<T)d* = ! / *

=

( * + 1 )/( 2 * )

 

 

о

 

 

 

 

 

 

 

— вероятность события, что наугад

выбранный максимум будет

положительным.

 

 

 

 

 

 

 

 

Плотность распределения амплитуд (14.15), приведенная к ча­

стоте процесса по всем его максимумам,

 

 

 

 

 

О при а < 0;

 

 

ж

=

l F

exP ( — - w )

при

а > ° -

(14.17)

 

 

 

153

Поскольку частота процесса по нулям

 

«о = n j k ,

(14.18)

плотность распределения амплитуд по методу максимумов, приведенная к частоте процесса нагружения по его нулям, будет определяться по формуле Релея:

/ (<0 = - j r ехР ( — 15г) •

04-:19)

Пусть кривая усталости описывается уравнением (1.2). Под­ ставив (14.15) и (1.2) в формулы (13.8), (13.9), (9.49) и (9.51), получим, что среднее значение и коэффициент вариации распре­ деления долговечности по методу максимумов могут быть вычис­

лены по формулам (14.11) и (14.12),

где ? — средний период про­

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

д о л г о в е ч н о с т и по

Р а с ч е т у с т а л о с т н о й

м е т о д у р а з м а х о в . В соответствии с (11.13) плотность распределения амплитуд напряжений при схематизации случай­ ного процесса по методу размахов

г ы - - * ? - * р ( — 5 1 ) . <'«■*»

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

2'n/2s>n

vn =

N0km<j^Jr

( - r + 1- ■*& )■•

 

 

1/2

S - [ r Gr+'- 4 H / I*(T + 1-

В этом случае

получаем

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

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

Т Р —

? э / ^ р 5

бгр =

бур/у^/Г,

где п = Т/1Э — ожидаемое число нагружений за время Т; ?э =

= 1/k — интервал времени между экстремумами процесса нагру­ жения; ? — средний период процесса нагружения по его нулям.

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

гиба, получаем для плотности распределения средних

[в соответ­

ствии с формулой (11.14)] следующее выражение:

 

/ (°т) —

ехр (

<

(14.21)

У2п&(1 —А-*)

\

2s2 (1 — А-*) /

 

154

Совместную плотность распределения амплитудных и средних напряжений циклов можно представить в следующем виде:

/ К . от) = / Ю / КО,

где f (оа) и / (а т ) — плотности, определяемые соответственно по формулам (14.20) и (14.21).

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

(14.22)

о о

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

лости

(1.17).

 

 

После несложных преобразований получаем

- = — Ц— j J I К») Г " f

(°а + ФхСГт)т f (Оа) <*Оа d(Jm +

 

« л

 

[ ’- А

}

 

°А

Г ат

|ат. сж1

(Ра - № т)т f (оа) doa dam 4-

+

J

/ (°m) I

J

 

—os

|_

 

}

 

 

 

 

ОГ °т. р-1

+ р ( о . )

j

(оа - ь°т)т f (оа) doa I dam }, (14.23)

<ТД

|_

 

где пределы интегрирования приняты в соответствии с диаграммой предельных амплитуд, показанной на рис. 1.4.

Внутренние интегралы в соотношении (14.21) приводятся

к виду

в

j x = 2b J (х+ а)т х ехр (—Ьх2) dx.

Так как

(* + а)т =

храт- р,

-т

где

о-

m(m— 1) ... (от — р+ 1)

 

 

С )-' (14.24)

то

1-2 ... р

л - 2 ( Г ) “"“р,’^/![г (-г + |’ ь■4' ) - г ( т + 1’ »**)]■

155

Для

первого интеграла а = \|зхсгт ,

b = /s2/(2s2), А = а_г

В =

аТ. р — ат ;

для

второго

интеграла

а =

—ф2сгт ,

Ъ = k2/(2s2),

А

= о_! — t|?acrm,

В = сгт

+ | <гТ. сж |;

Для третьего

интеграла а =

—ф2ат ,

b — k2/(2s2),

Л = а_г — ф2ат ,

В =

= <*т. р

От.

 

 

 

 

 

 

Подставив (14.24) и (14.23) в формулу (14.21), приходим к сле­

дующему интегралу:

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

Л = j

exp (—сх2) хпГ [-у- + 1.

Ь (о_х — х)2] dx,

(14.25)

 

 

с

 

 

 

 

 

 

где С,

D, с,

b — константы, определяемые по заданным параме­

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

териала.

д о л г о в е ч н о с т и по

Р а с ч е т у с т а л о с т н о й

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

Для расчета долговечности

по этому методу можно также воспользоваться формулами (13.8) и (13.9). При этом плотность распределения амплитуд напря­ жений вычисляется либо по точной формуле (11.24), либо по приближенной (11.33). Если использовать (11.24), то результаты могут быть получены только численными методами с использо­ ванием ЭВМ. В случае, когда используется соотношение (11.33), ожидаемая долговечность

lN0arH_l2~mf2s - m

- § г )

 

j r ( i +

I.

k2a2а!г \

 

2sa

J

- r ( » + , ) i - L ]

 

 

 

Т —

 

 

 

 

 

 

при

<

a*;

 

 

W0a«i2- m/2s~m

при o_x >

<r*,

 

c r f - J + l ,

aij/(2s2)]

 

 

 

 

 

где a* = sx*; ? — среднее значение периода процесса нагружения по его нулям.

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

начального значения k0 до единицы. В

результате средние зна­

чения

циклов

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

уменьшаются

до

нуля,

и их

плотность

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

1 переходит

в

дельта­

функцию [см. (14.21)). Таким образом, параметр сложности структуры процесса нагружения k в методе полных циклов будет функцией значений амплитуды напряжений. С учетом этого

156

Рис. 14.8. Интегральные функции распределений амплитуд при схематизации случайного процесса сложной структуры различными методами размахов (/), полных циклов (2) и максимумов (3):

а — при kt = 2; б — при к = 2,5; в — при *<, =3; г — при к0 — 3,5

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

/ [®а> tXml = / О^а* ^о) / 1®т> ^ (*7а)Ь

(14.26)

где функция k — k (оа) — решение уравнения (11.23). Подставив (14.26) в формулу (14.22), получим возможность

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

С о п о с т а в л е н и е м е т о д о в р а с ч е т а д о л г о ­ в е ч н о с т и . Для сопоставления полученных результатов на рис. 14.8 приведены графики интегральных функций распре­ деления амплитуд напряжений при схематизации случайных процессов нагружения по методам максимумов, размахов и пол­ ных циклов при использовании формулы (11.33). Метод полных циклов дает для больших квантилей промежуточные значения вероятностей по сравнению с методами максимумов и размахов. Сопоставление значений долговечностей, получаемых этими мето­ дами, производится с помощью следующих соотношений:

Гм (1 + Фо) = = 1 Л »

где Ти, Тр и Т0 — значения долговечностей,полученных по методу максимумов, размахов и полных циклов; <р0 и <р — параметры.

Зависимости <р и <р0 от параметра k сложности структуры про­ цесса нагружения приведены на рис. 14.9 и 14.10.

Расчеты на сопротивление усталости элементов конструкций, находящихся в условиях интенсивной коррозии *. Коррозия ме-

* Разработка выполнена совместно с Н. А. Черновой.

157

Vo

Vo

Vo

 

1

1

 

 

/

 

 

w

1

1

1 1

»

t

»____ L _ I ____

3 5

7 9 k

t

3

5

7 9 к

'

3

5

7 9 k

 

a)

 

 

 

6)

 

 

0}

 

Рис. 14.9. Зависимость <р0 от

параметра к:

 

 

 

 

а — при

т = 3;

б — при т — 4;

в — при

т = 5; 1 — a_1/s = 1 ;

2

= 2

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

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

Снижение прочностных характеристик металлов в результате их коррозии может быть описано с различной степенью детали­ зации. Так, если кривая усталости описывается степенным урав­ нением (1.2), то в общем случае необходимо выявить в экспери­ менте три функции, описывающие эволюцию параметров кривой

Рис. 14.10. Зависимость <р от параметра к:

а — при т = 3; б — при т = 4; в — при т = 5; / -■ o ^ /s = 1; 2 —a ^/s = 2

158

усталости: а_г = a_t (t), m = m (t) и N0 = N0 (t). Рассмотрим только случай, когда предел выносливости постепенно пони­ жается, а параметры т и N0 остаются постоянными.

Описание эволюции предела выносливости, обусловленной на­ растанием коррозии металла, основывается на накопленном экс­ периментальном материале [21 ]. В достаточно общем случае можно считать, что скорость изменения предела выносливости зависит от качества металла (высокопрочные металлы обычно более чувствительны к коррозии, чем менее прочные) и от уровня действующих напряжений. В этом случае кинетическое уравнение, описывающее изменение предела выносливости со временем, можно записать в следующем виде:

<j-i = / (о, a j),

где / (a, a_i) — функция, зависящая от уровня действующих напряжений в данный момент времени и от пределов выносливости материала для данного момента времени.

Полагая, что функцию / (о, a_t) можно представить в виде произведения двух функций, одна из которых а (о) зависит только от уровня напряжений, а. вторая Р (о.!) — только от предела выносливости, получаем

0_i = a (а) р (о.!).

(14.27)

Разделив переменные в дифференциальном уравнении (14.27) и проинтегрировав его, получим искомую зависимость a_j = = сг_х (/). Для случая, когда р (ст^) = —a_lt имеем экспоненциаль­ ный закон изменения предела выносливости со временем

о_1 (t) = 0^-1 exp [—a (0) t],

(14.28)

где cr^i — начальное значение предела выносливости.

Если темп снижения предела выносливости со временем не зависит от уровня напряжений, то а (о) = a = const.

В результате эволюции предела выносливости мера устало­ стного повреждения v за один цикл даже при стационарном про­ цессе нагружения будет изменяться от цикла к циклу и для не­ которого бесконечно малого интервала времени dt

dv = v dt/t0,

где t0 — средний период процесса нагружения; v = v (/) — уста­ лостное повреждение за один цикл нагружения, соответствующий моменту времени t.

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

стное повреждение за время t

 

t

(14.29)

** = -% -$ * (О Я.

О

 

159

 

 

0,01

0,03

0,05

0,07

0,09

0,11

0,13 2act

Рис. 14.11.

Графики функции JA=

Рис.

14.12.

Графики функции / в =

= У4 (2а/)

при различных значени­

=

(2а/)

 

при

различных значениях

ях x_i

 

 

 

 

 

 

 

 

 

Используя в соотношении (14.29) формулы (1.2), (13.8) и

(14.19), получаем

 

 

=

1

0-1 (0/(2s2) | dt. (14.30)

о

'

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

vt = *^т/(2оо^о^о)<

где

Ут =

J 2-m/2- i г [(т -)_ 2)/2, 2] dz,

(14.31)

/_,е-2а<

*-i = a!Li/(2s2).

Интеграл (14.31) может быть выражен через табулированные функции. После некоторых преобразований получим

J . -

ехр ( - А ) - < А ) - ,! Г

А ) -

 

— ехр £—JC^_I exp (—2a/)] -f (*-i)_m/2 x

X

exp (atm) Г ( m ~ ^2- , x!Li exp (—2a/)) .

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

деляется из уравнения vt = 1.

протабулированы.

Функции

Jm для различных значений m

Соответствующие таблицы приведены в прил. III. На рис. 14.11— 14.14 представлены графики этих функций.

160

Гх.,==*5

г О у/

о M

s = e = j _ j ------------

0,09

1., .. ------

2а t

0,01

0,03

0,05

0,07

0,11

0,13

Рис.

14.13.

Графики

функции У8 =

= JB(2аt)

при

различных

значени­

ях х_х

 

 

 

 

 

 

 

X . J

о

>

0,05

---- 1---- 1.

0,13 2a t

0,0/

0,03

0,07

0,09

0J1

Рис.

14.14.

Графики функции У10 =

= J10 (2аt)

при

различных

значени­

ях

 

 

 

 

 

 

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

Из всех возможных математических моделей нестационарных случайных колебаний а (I) рассмотрим только такие, которые задаются детерминированными функциями изменения со време­

нем их

основных параметров: п = n(t) — средней частоты; s2 =

= s2 (/)

— дисперсии; д = д (() — среднего значения процесса

нагружения и т. п. В некоторых моделях таких процессов изме­ няется только один из этих параметров, в других — могут изме­ няться одновременно несколько из этих параметров.

Усталостное повреждение за один цикл нагружения, соответ­ ствующий моменту времени t:

*■<'> = 1 У

<|51)

О —оо

 

 

где / (<х0, ат, t) — совместная плотность распределения амплитуд­ ных и средних напряжений в циклах нагружения, соответству­ ющих моменту времени t; N (аа, от, t) — число циклов до раз­ рушения при амплитуде <т0 и среднем значении цикла нагруже­ ния ат , соответствующего моменту времени / (имеется в виду, что в результате, например, старения материала его прочностные характеристики могут изменяться со временем).

Усталостное повреждение, накопленное в конструкции за

время

t,

 

 

t

 

 

v (t) = j V(1) (t) n (i)dt.

(15.2)

 

о

 

Долговечность T конструкции определяется из решения урав­

нения

v (Г) = 1. Рассмотрим несколько частных

случаев.

11 Гусев А. С.

161