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

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

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

симумов ГССП Райса, распределения размахов Гусева [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 Ху по явно содержащемуся

времени

ф:

 

 

 

 

 

 

 

 

 

 

 

 

 

*i<e-¥>—SrJ

0) Ф (a cos ф, — аа>sin ф) sin ф +

 

I

в»008**

1 ^ _

Д(°>

,

°« .

 

+

2со2а

aV

«

+

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)

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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) при случайном нагружении. Процессы нагружения представляли собой ГССП с различными формами спектральной плотности. В каждой из этих работ приведены формы спектральных плот­ ностей, а также кривые усталости при регулярном и случай­ ном нагружении. Для каждой из использованных форм спек­ тральной плотности методом СМ рассчитывались К э по ме-

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