книги / Сопротивление усталости и живучесть конструкций при случайных нагрузках
..pdfРис. 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Г Г М + 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