книги / Многоцикловая усталость при переменных амплитудах нагружения
..pdfсимумов ГССП Райса, распределения размахов Гусева [59], Ковалевского [234], распределения полных циклов, предло женного А. С. Гусевым (2.62). Для узкополосного режима нагружения Р = 1 все распределения превращаются в рас пределение Рэлея и Кэ определяется одинаково для всех
расчетных случаев:
К„ = К 2 |
(4.33) |
( - S - + |
* ) Г |
Согласно теории выбросов распределение выбросов для ГССП является рэлеевским, а их число совпадает с числом пересечений среднего уровня процесса. Поэтому по распре делению выбросов зависимость К ъ от Р следует рассчитывать
по формуле
Кг = Я„/р1/\ |
(4.34) |
В методах Шефера — Ежова и |
Трофимова не содержит |
ся непосредственно способов для вычисления Кд. Рассмотрим
в качестве примера соотношение, предложенное О. Ф. Трофи мовым (4.25), причем положим р# = 1:
h
Теперь получим выражение для Кэ, при этом отметим, что согласно формуле (4.21) i V p — (Ка^/К31)~ь<, откуда с учетом формул (4.33), (4.35) следует К3р = Каi/p2/b и
Я э Р = / 2 [Г (4 - + 1)' |
1/Ь |
2/Ъ |
(4.36) |
|
/р: |
Аналогичные выражения получены для формул Ежова — Шефера и Рамеша (4.26) и (4.27). Результаты расчетов Ка
по распределению Райса представлены в табл. 2. Расчеты величины Кэ по теории выбросов и по формуле Трофимова
представлены в табл. 3, а по распределению полных циклов Гусева — в табл. 4. Выражение для К а по гипотезе спект
рального суммирования легко получить, сопоставив равен ства (4.21) и (4.28):
ЛГ„ = [L (Ь) со* v M„/Mt |
(4.37) |
Результаты расчетов Кв методами статистического мо
делирования представлены на рис. 79 для 2>, равного 3; 6 и 10. Следует отметить, что расчеты Ка по методу СМ дают
лишь приближенную оценку, что и учитывается при представ лении результатов расчетов — на рисунках приведены оцен ки Ка и 90 %-ный доверительный интервал.
U
ь
|
1,0 |
1.2 |
1,4 |
1.0 |
1,0 |
2.0 |
2.5 |
3.0 |
3.5 |
4.0 |
1 |
1,253 |
1,070 |
0,9534 0,8713 0,8103 0,7632 |
0,0819 |
0,6302 0,5943 |
0,5681 |
||||
2 |
1,414 |
1,296 |
1,212 |
1,149 |
1,100 |
1,061 |
0,9899 |
0,9428 0,9091 |
0,8839 |
|
3 |
1,555 |
1,465 |
1,396 |
1,343 |
1,301 |
1,266 |
1,202 |
1,159 |
1,127 |
0,1103 |
4 |
1,682 |
1,607 |
1,549 |
1,502 |
1,464 |
1,433 |
1,374 |
1,333 |
1,303 |
1,280 |
5 |
1,798 |
1,734 |
1,682 1,640 |
1,606 |
1,577 |
1,522 |
1,483 |
1,455 |
1,432 |
|
6 |
1,906 |
1,849 |
1,803 |
1,765 |
1,733 |
1,706 |
1,654 |
1,617 |
1,589 |
1,568 |
7 |
2,007 |
1,955 |
1,913 1,878 |
1,948 |
1,823 |
1,774 |
1,739 |
1,712 |
1,691 |
|
8 |
2,102 |
2,055 |
2,616 |
1,983 |
1,956 |
1,932 |
1,885 |
1,851 |
1,825 |
1,805 |
9 |
2,193 |
2,149 |
2,113 |
2,082 |
2,056 |
2,033 |
1,989 |
1,956 |
1,931 |
1,911 |
10 |
2,279 |
2,238 |
2,204 2,175 |
2,150 |
2,150 |
2,128 |
2,086 |
2,030 |
2,010 |
|
11 |
2,361 |
2,322 |
2,290 |
2,262 |
2,239 |
2,218 |
2,177 |
2,147 |
2,123 |
2,104 |
12 |
2,439 |
2,402 |
2,372 |
2,346 |
2,323 |
2,304 |
2,264 |
2,235 |
2,211 |
2,193 |
Т а б л и ц а 3. Расчеты К э по теории выбросов (над чертой) и форму
ле Трофимова (под чертой)
Р
ь |
1.2 |
1.4 |
1.6 |
1,8 |
2.0 |
?.5 |
3,0 |
3,!> |
|
|
|
||||||||
3 |
1,463 |
1,390 |
1,329 |
1,278 |
1,234 |
1,146 |
1,078 |
1,024 |
0,980 |
|
1,377 |
1,243 |
1,137 |
1,051 |
0,980 |
0,844 |
0,748 |
0,675 |
0,617 |
4 |
1,607 |
1,546 |
1,495 |
1,452 |
1,414 |
1,337 |
1,278 |
1,230 |
1,189 |
|
1,535 |
1,421 |
1,330 |
1,254 |
1,189 |
1,064 |
0,971 |
0,890 0,841 |
|
5 |
1,734 |
1,681 |
1,637 |
1,599 |
1,565 |
1,497 |
1,443 |
1,400 |
1,363 |
|
1,672 |
1,572 |
1,490 |
1,421 |
1,363 |
1,246 |
1,159 |
1,089 |
1,033 |
6 |
1,849 |
1,802 |
1,763 |
1,728 |
1,698 |
1,636 |
1,587 |
1,547 |
1,513 |
|
1,794 |
1,704 |
1,630 |
1,561 |
1,513 |
1,405 |
1,322 |
1,256 |
1,201 |
7 |
1,956 |
1,914 |
1,878 |
1,846 |
1,819 |
1,762 |
1,716 |
1,679 |
1,647 |
|
1,906 |
1,824 |
1,756 |
1,698 |
1,647 |
1,545 |
1,467 |
1,404 |
1,351 |
8 |
2,057 |
2,017 |
1,984 |
1,955 |
1,929 |
1,876 |
1,834 |
1,799 |
1,769 |
|
2,010 |
1,934 |
1,871 |
1,816 |
1,769 |
1,673 |
1,599 |
1,538 |
1,488 |
9 |
2,151 |
2,115 |
2,084 |
2,057 |
2,033 |
1,983 |
1,934 |
1,910 |
1,882 |
|
2,108 |
2,037 |
1,978 |
1,927 |
1,882 |
1,791 |
1,720 |
1,662 |
1,613 |
10 |
2,241 |
2,207 |
2,178 |
2,152 |
2,130 |
2,083 |
2,045 |
2,014 |
1,987 |
|
2,201 |
2,134 |
2,078 |
2,029 |
1,987 |
1,900 |
1,832 |
1,777 |
1,730 |
11 |
2,327 |
2,295 |
2,267 |
2,243 |
2,222 |
2,177 |
2,141 |
2,112 |
2,086 |
|
2,289 |
2,226 |
2,173 |
2,127 |
2,086 |
2,003 |
1,938 |
1,884 |
1,839 |
12 |
2,410 |
2,379 |
2,353 |
2,330 |
2,310 |
2,267 |
2,233 |
2,204 |
2,180 |
|
2,374 |
2,314 |
2,263 |
2,219 |
2,180 |
2,100 |
2,038 |
1,986 |
1,942 |
Р
ь
|
1.2 |
м |
1.6 |
1,8 |
2,0 |
2,й |
3,0 |
3.5 |
4.0 |
1 |
1,223 |
1,177 |
1,127 |
1,077 |
1,029 |
0,921 |
0,831 |
0,755 |
0,691 |
2 |
1,380 |
1,328 |
1,271 |
1,215 |
1,161 |
1,040 |
0,938 |
0,852 |
0,780 |
3 |
1,518 |
1,460 |
1,398 |
1,336 |
1,277 |
1,143 |
1,031 |
0,937 |
0,858 |
4 |
1,642 |
1,579 |
1,512 |
1,445 |
1,381 |
1,236 |
1,115 |
1,013 |
0,928 |
5 |
1,755 |
1,688 |
1,617 |
1,545 |
1,376 |
1,322 |
1,192 |
1,083 |
0,922 |
6 |
1,961 |
1,790 |
1,714 |
1,658 |
1,565 |
1,402 |
1,264 |
1,149 |
1,051 |
7 |
1,960 |
1,885 |
1,805 |
1,725 |
1,649 |
1,476 |
1,331 |
1,210 |
1,107 |
8 |
2,054 |
1,976 |
1,892 |
1,808 |
1,728 |
1,547 |
1,395 |
1,268 |
1,160 |
9 |
2,143 |
2,061 |
1,974 |
1,886 |
1,803 |
1,614 |
1,455 |
1,323 |
1,211 |
10 |
2,228 |
2,143 |
2,052 |
1,961 |
1,874 |
1,678 |
1,513 |
1,375 |
1,259 |
11 |
2,318 |
2,222 |
2,127 |
2,033 |
1,943 |
1,740 |
1,569 |
1,426 |
1,305 |
12 |
2,388 |
2,298 |
2,200 |
2,103 |
2,009 |
1,799 |
1,622 |
1,474 |
1,350 |
|
Результаты расчетов |
показывают |
хорошее соответствие |
между расчетами при схематизации методом максимумов по распределению Райса и по методу СМ также при схемати зации методом максимумов: только для одного вида спектра (при р = 4 , 1 ) и Ь > - 8 значения ЛГЭ, рассчитанные по форму ле Райса, лежат внс90%-ных доверительных интервалов, построенных методом СМ.
Расчеты по формуле Шефера — Ежова дают результаты, близкие к формуле Трофимова, и поэтому на рисунках не показаны. Также не приведены результаты расчетов по фор мулам работ Рамеша и Ковалевского, так как соответству ющие значения К9 оказываются заниженными. Анализ ре
зультатов, представленных ва рис. 79, свидетельствует о том, что метод выбросов дает оценку сверху по сравнению с методом СМ. Формула Трофимова дает результаты, хоро шо согласующиеся с расчетами по методу СМ в диапазоне Р = 1 -г- 2. Паилучшео соответствие с результатами расчета по методу СМ дает формула (4.37) спектрального суммиро вания и расчет с использованием распределения полных циклов Гусева. В табл. 5 приведены результаты количествен ного сопоставления расчетов по формуле (4.37) и по распре делению полных циклов Гусева. В качестве меры отклонения Д принимали усредненную сумму квадратов разностей ме жду значениями К&, рассчитанными по методу СМ и по каж
дому из сравниваемых соотношений для всех вариантов спектров.
Как видно из таблицы, соответствие расчетов по гипотезе спектрального суммирования и по распределению полных циклов примерно одинаковое.
Рис. 79. Зависимость К э от |
Р при d = 3 |
(о), d = 6 |
(б) и d — 10 (в): |
1 — расчет по формуле Р ай са для |
распределения |
максимумов |
случайного процесса; 2 — расчет по теории выбросов; з — расчет по фор |
м уле Трофимова; |
4 — расчет по ф орм уле Г усева д л я распределения разыахов; 5 — расчет по формуле Гусева для распределения полны х |
|
циклов; 1 — расчет по |
гипотезе спектрального суммирования; 11 — расчеты методом СМ и схем атизации по методу полны х циклов; 111 — |
|
расчеты методом |
СМ и |
схем атизации по м етоду максимумов. |
Т а б л и ц а 5. |
Данные |
расче |
|
тов значения |
Д |
|
|
Гипотез» |
Распределе |
||
сп е ктр а л ь |
ние |
полны х |
ьного сумми ц и кл ов Гу
|
рования |
|
сева |
|
|
|
|
|||
3 |
|
0 ,1 3 5 |
|
|
0 ,1 2 6 |
|
|
|
|
|
6 |
|
0 ,1 5 0 |
|
|
0 ,1 1 7 |
|
|
|
|
|
8 |
|
0 ,1 4 3 |
|
|
0 ,1 3 0 |
|
|
|
|
|
10 |
|
0 ,1 5 0 |
|
|
0 ,1 8 4 |
|
|
|
|
|
Для |
проверки |
предпо |
|
|
|
|||||
ложения |
о том, что Р пол |
|
|
|
||||||
ностью определяет К0, бы |
|
|
|
|||||||
ли |
рассчитаны |
значения |
|
|
|
|||||
Ка |
для |
трех |
вариантов |
|
|
|
||||
спектральных |
плотностей |
|
|
|
||||||
при р « |
1,4 и для трех |
ва |
|
|
|
|||||
риантов |
W (/) |
при р « |
2. |
|
|
|
||||
На |
рис. |
80 |
представлены |
|
|
|
||||
зависимости |
К ь от b для |
|
|
|
||||||
разных |
|
спектральных |
|
|
|
|||||
плотностей. |
При |
р й* 1,4 |
Рис. 80. Зависимость |
К 0 от |
при |
|||||
значения К а для |
разных |
|||||||||
различных значениях 0. |
|
|||||||||
спектров практически сов- |
|
|||||||||
|
|
|
пали, при Р й? 2 они различаются значимо, так как довери тельные интервалы не пересекаются, что качественно согла суется с расчетами по формуле спектрального суммирования. Расчеты по этой формуле несколько проще, чем с помощью распределения полных циклов Гусева, что позволяет реко мендовать формулу (4.37) для расчетов Кэ с учетом формы
спектральной плотности. В тех случаях, когда форма спек тральной плотности неизвестна, но известен параметр Р, ре комендуется применять распределение полных циклов Гусе ва. При нагрузках с параметром широкополосности Р < 2 целесообразно применять формулу Трофимова (4.36) как наиболее простую.
5. РАСЧЕТ ДОЛГОВЕЧНОСТИ ПРИ СЛУЧАЙНЫХ ОДНОЧАСТОТНЫХ КОЛЕБАНИЯХ С УЧЕТОМ АМПЛИТУДНО-ЗАВИСИМОГО ТРЕНИЯ
В предыдущих параграфах рассмотрена зависимость пара метров усталостного разрушения при случайном нагружении от статистических характеристик нагрувок с использо ванием предположения о нормальности процессов пагруже-
ния. Здесь же рассмотрен вопрос о том, насколько существен но отклонение от гауссовского закона распределения мгновен ных значений процесса нагружения влияет на прогнозируе мую долговечность в практически важном случае, что по зволяет с некоторым приближением установить рамки приме нимости общепринятой модели случайной нагрузки.
Предположим, что прочность случайно колеблющейся конструкции определяется напряжением а (£) в наиболее
опасном месте, и допустим, что в конструкции возбуждает ся одна резонансная форма колебаний. В этом случае а {£)
является узкополосным случайным процессом. Рассматри вая а (t) как обобщенную координату, уравнение движения
по одной форме колебаний можно записать в виде
о -f- со2а + еФ (сг, а) = g (t), |
(4.38) |
где со — собственная частота системы в отсутствие демпфи рования; еФ (а, а) — функция, описывающая рассеяние
энергии в конструкции при колебаниях на данной резонирую щей форме; 8 — малый параметр, отражающий относитель ную величину рассеиваемой за цикл энергии; g (t) — обоб
щенное возбуждение типа нормального белого шума, по Р. Л. Стратоновичу [17, 1631 (производная от винеровского
процесса), с нулевым математическим ожиданием и интенсив-
о
ностыо as (имеющей порядок в).
Как показано в работах Г. С. Писаренко, рассеяние энер гии в материале упругой системы можно описать с помощью гипотезы Н. Н. Давиденкова [123, 1241. Возможность приме нения этой гипотезы для изучения случайных колебаний упругих систем с нелинейным рассеянием энергии обоснова на В. М. Гончаренко в работе [481.
Для нахождения плотности вероятности р (а) амплитуд
процесса применим теорию марковских процессов. Посколь ку методика нахождения функции р (а) изложена в работах
[15, 17, 48], ниже приводятся только основные этапы вычис лений.
В уравнении (4.38) сделаем замену переменных: |
|
||||
|
a (t) = А (t) cos ф (f); |
и (t) — — A (f) to sin ф (t); |
|
||
|
|
ф (0 = |
ctf + |
9 (0 , |
(4.39) |
где A |
(t) |
и 0 (/) — огибающая |
и фаза процесса a (t). |
Для |
|
функций А и 0 получим следующие уравнения: |
|
||||
А = |
1 |
[еФ (A cos ф, — A tosin ф) sin ф — g (t) sin ф]; |
|
||
— |
|
||||
|
" |
|
|
|
(4.40) |
0 = [еФ (A cos ф, — Аы sin ф) cos ф — g (t) cos ф].
Плотность распределения вероятностей р (а, ф) удовлет
воряет уравнению Фоккера |
Планка — Колмогорова |
||||
|
|
dsip |
( а д + 4 - - £ - ( хиР) + |
|
|
|
1 |
д‘ |
а2 |
|
(4.41) |
+ |
Т |
<х22?) + |
дадф (Xl2Р). |
|
|
где |
|
|
|
|
|
Xi («» Ф, *) = |
— |
Ф (a cos ф, — а© sin ф) sin -ф + |
О^ COS2 ф |
||
*ы ч |
; |
||||
|
|
|
|
|
(4.42) |
|
Р |
^ I -----1- |
— |
^.f.V « >I- |
Хг (а, ф, t) = |
*■* |
|||
саа |
(a cos ф, |
— а© sin ф) cos ф — |
(Лsin ф cos ф
а _______ •
2М а-(0
|
|
|
|
|
|
|
|
(4.43) |
|
|
о? sin 3 ф |
|
|
COS2 ф |
|
||
|
|
— _ i |
; |
(а , ф , |
*) = |
со2а2 |
; |
(4.44) |
Хп (а* Ф> 0 — |
||||||||
|
|
X,, (а, ф, |
() = |
о? sin ф cos ф |
|
|
(4.45) |
|
|
|
- « -------------- |
|
|
||||
|
|
|
|
дат |
|
|
|
|
Отметим, |
что уравнение |
Фоккера — Планка — Колмо |
||||||
горова |
при |
интерпретации g |
(I) как белого шума |
в смысле |
||||
Ито не |
отличается от |
соотношений |
(4.41) — (4.45), |
если |
учесть при выводе уравнений (4,40) формулы Ито для стоха стических дифференциалов [112].
Применяя принцип усреднения Р. 3. Хасьминского к дифференциальному оператору параболического типа (4.41), усредняем коэффициенты Xi Ху по явно содержащемуся
времени |
ф: |
|
|
|
|
|
|
|
2я |
|
|
|
|
|
|
*i<e-¥>—SrJ |
0) Ф (a cos ф, — аа>sin ф) sin ф + |
|
|||||
I |
в»008** |
1 ^ _ |
Д(°> |
, |
°« . |
|
|
+ |
2со2а |
aV |
« |
+ |
2оа>2 » |
|
|
|
2л |
|
|
|
|
|
(4.46) |
Х2 (а,ф) |
|
соа |
Ф (асоэф, |
— а© 8тф )созф — |
|||
|
|
||||||
|
|
|
|
|
|
|
|
|
2а2а>2 |
|
а со |
|
|
||
|
|
|
|
|
•> |
|
|
Хц(а, ф) |
|
Х22(«. Ф) = b)'Jaa |
f Ф12 (а, ф) = О, |
|
R (a) = |
— L j |
ф (a cos ф, |
— aco sin ф) sin |
|
(4.47) |
|||
|
|
о |
|
|
|
|
|
|
|
|
2л |
|
|
|
|
|
|
Q (a) = |
~2^- j |
Ф (a cos ф, |
— a<s>sin ф) cos ф d\p. |
(4.48) |
||||
|
|
о |
|
|
|
|
|
|
Рассматривая |
стационарные |
колебания |
= 0), |
полу- |
||||
чаем, что фаза распределена равномерно о р (а, |
ф) = |
X |
||||||
X р (а), откуда |
следует уравнение для р (а): |
|
|
|||||
— j r 1д <“> р/«>] - |
i L |
& |
\ |
, 1 |
J |
= 0. |
(4.49) |
|
2oj2a |
} ' |
2 |
d a 1 |
|||||
Функция aR (а) — величина |
рассеиваемой |
энергии за |
цикл колебаний (на единицу массы), и поэтому для большин ства конструкционных материалов допустима степенная ап
проксимация этой |
функции: |
|
|
R (а) = — га”, |
(4.50) |
где п — параметр |
петли гистерезиса материала |
(для случая |
амплитудно-независимого трения п — 1); г — величина, за
висящая от параметров петли гистерезиса и от формы коле баний конструкции.
Из дифференциального уравнения (4.49) получим следу
ющее выражение для р (а) |
с учетом граничных условий на |
|||
бесконечности р (а) |
0 и |
0 при а |
сю: |
|
р (а) = са ехр (— pan+J), р. = |
2шг |
(4.51) |
||
|
||||
|
|
(я + 1)o2g |
|
|
Из формулы (4.51) |
находим |
среднеквадратическое значение |
||
процесса a (t): |
|
|
|
|
о?к = Е С<42 cos2 ф) = |
~ J а2р (a) da = -у- ^ |
а3 ехр (— рап+1)с?а. |
||
|
о |
г о |
|
|
|
|
|
|
(4.52) |
|
|
|
|
оо |
Используя выражения (4.51), (4.52) и условие J р (a)da =
о
= 1, выразим с и |х через п и ис„:
(4.53)
В |
случае, если и — 1 (ам- |
Рис. 81. Распределенве амплитуд |
плитудно-независимое тре- |
узкополосного процесса при раз- |
|
иие), |
распределение (4.53) |
личных значениях л. |
превращается в распределе ние Рэлея, справедливое для амплитуд узкополосного гаус
совского процесса [14]. Для большинства конструкционных материалов п изменяется от 1 до 3. На рис. 81 представлены
графики р (—— ) для различных значений п.
Рассмотрим вопрос о влиянии на прогнозируемую устало стную долговечность отклонения распределения амплитуд согласпо выражению (4.53) от распределения Рэлея [31, 32]. Для этого в формулу линейного суммирования подста вим плотность распределения (4.53) и получим следующее выражение для числа максимумов до разрушения:
|
|
N п — |
(п - f 1) Яц}?+2)/(Т1+1) |
(4.54) |
|
|
6 + 2 |
||
|
|
|
|
|
|
|
|
п + 1 |
|
где Г (и, |
v) — неполная гамма-функция. При п = 1 форму |
|||
ла (4.54) приобретает вид (4.30). |
|
|||
На рис. |
82 показана |
зависимость N J N j от а при п = 2 |
||
и п — 3 |
и |
различных |
значений Ъ в полулогарифмической |
системе координат. Анализ графиков показывает, что форму ла (4.54) при п = 1 предсказывает меньший срок службы конструкции, чем при п > 1, причем это расхождение мо жет быть весьма значительным. Величина N J N j возрастает при увеличении параметра петли гистерезиса п и при уве
личении се. Это объясняется тем, что существенный вклад в повреждаемость материала дают выбросы случайного процес са нагружения с малой вероятностью появления, которая, как видно из рис, 81, уменьшается с увеличением п.
Полученные результаты показывают, что для более точ ного прогнозирования усталостной долговечности конструк ции необходимо привлекать экспериментальные данные о распределении плотности вероятностей амплитуд напряжений
в материале конструк ции, а при проведении усталостных испытаний при случайном нагруже нии на установках резо нансного типа необхо димо учитывать нели нейные эффекты, связан ные с конструкционным гистерезисом и рассея нием энергии в матери але резонансных узлов испытательной устано вки.
Следует отметить, что полученные резуль таты, точность которых соответствует точности
Рис. 82. Зависимость Л/п/Л/, от а при
первого приближения
различных значениях b u n .
метода Крылова — Бо
голюбова, не зависит от формы петли гистерезиса и основываются только на степен ной зависимости рассеиваемой энергии от амплитуды.
6.СОПОСТАВЛЕНИЕ РЕЗУЛЬТАТОВ РАСЧЕТНЫХ МЕТОДОВ
ИЭКСПЕРИМЕНТА ПРИ ОПРЕДЕЛЕНИИ ДОЛГОВЕЧНОСТИ
Впредыдущем параграфе рассмотрен вопрос о вычислении
величины К д, используемой в расчетах по линейной гипоте
зе долговечности при случайном нагружении с учетом формы спектральной плотности нагрузки. Однако применимость гипотезы линейного суммирования не всегда обоснованна, поэтому в данном параграфе представлены результаты со поставления экспериментов и расчетов по различным гипоте зам накопления повреждений для выбора расчетных формул, удовлетворительно прогнозирующих долговечность при слу чайном нагружении.
В работах [202, 221, 236] приведены результаты испыта ний алюминиевого сплава 2024 (аналог Д16) при случайном нагружении. Процессы нагружения представляли собой ГССП с различными формами спектральной плотности. В каждой из этих работ приведены формы спектральных плот ностей, а также кривые усталости при регулярном и случай ном нагружении. Для каждой из использованных форм спек тральной плотности методом СМ рассчитывались К э по ме-