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

книги / Механическая усталость в статистическом аспекте

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

О ПРИМЕНЕНИИ УСКОРЕННЫХ МЕТОДОВ ПРО И ЭНОМОТО ДЛЯ ОЦЕНКИ ХАРАКТЕРИСТИК УСТАЛОСТИ ЛЕГКИХ СПЛАВОВ

ИИХ РАССЕЯНИЯ

М.Н. СТЕПНОВ, Ю. А. ШУХМИН'

Всвязи с отсутствием горизонтального участка кривой усталости лег­ ких сплавов часто возникает необходимость определения для них ограни­

ченного предела выносливости на базе 50, 100 миллионов циклов и бо­ лее, что значительно увеличивает длительность испытаний.

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

Большая трудоемкость усталостных испытаний определяет необхо­ димость их форсирования.

Для сравнения результатов ускоренных и обычных усталостных ис­ пытаний для большого числа алюминиемых и магниевых сплавов были получены кривые усталости. Каждую кривую строили на основании ре­ зультатов испытания 100 и более образцов диаметром 8—10 мм. Уста­ лостные испытания проводили при чистом изгибе с вращением со скоро­ стью 6000 об/мин.

В работе [1] для ускоренного определения предела выносливости

предлагается использовать уравнение, предложенное Вейбуллом

[2];

для правой части кривой усталости уравнение имеет вид

 

N , = А > * - с гГ \

(1)

где К и т — параметры, зависящие от материала; Л'* — число циклов до разрушения при напряжении а*; сгг — предел выносливости.

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

в ^ Ь - Ы Г

(2)

С точки зрения наибольшего приближения к экспериментальным дан­ ным были исследованы также уравнения Бастенейра и Эпремиана [31. Для ряда легких сплавов были найдены параметры по экспериментальным точкам с наибольшей долговечностью 107 циклов и с помощью этих пара­ метров определен предел выносливости на базе 108 циклов. Отклонения полученных значений от экспериментальных, соответствующих той же базе по уравнению (2), не превышали 10%. В пределах такой погрешности уравнение (2) можно применять для экстраполяции в область более значи­

тельных баз на один порядок; при дальнейшем увеличении базы отклоне­ ния от экспериментальных значений заметно возрастают.

В особую группу можно выделить ускоренные методы Локати [4], Про [5], Эномото [6]. Основная особенность этих методов состоит в том, что амплитуда переменных напряжений в процессе испытания каждого образца увеличивается с ростом числа циклов вплоть до разрушения. Некоторой модификацией этих методов применительно к конкретным ус­ ловиям и материалам являются методы Шварева [7] (для сталей в усло­ виях концентрации напряжения) и Нэдэшана [8] (для углеродистых ста­ лей). Первые три метода проверялись многими авторами применительно к сплавам на железной основе. В ряде случаев результаты получались удовлетворительными. Так, в работе [9] показано, что при наличии пара­ метров сходственных кривых усталости ускоренные испытания по методу Локати дают для черных металлов средние значения пределов выносли­ вости и их дисперсию.

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

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

определяли из уравнения

 

бр = 0»г + Аат,

(3)

. где ос — скорость возрастания амплитуды напряжений; ор — разрушающее напряжение; А и т — параметры, зависящие от свойств материала;

о?! — экстраполированное значение предела выносливости, равное раз­ рушающему напряжению при а = 0.

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 1

Химический состав и механические свойства

исследованных сплавов

 

 

 

 

 

 

Содержание, %

 

 

 

 

°0,2'

°ъ,

Сплав

 

 

 

 

 

 

 

 

 

 

Си

ме

Мп

Ре

81

2п

Т1

Сг

А!

»•, % кГ/мм- КГ/Л1Л1*

 

 

 

 

МЛ5

 

Ост.

0,30

0,02

0,07

0,52

 

 

8,10

7,6

12,6

22,1

АВГ

0,05

0,76

0,73

0,11

0,35

Ост.

16,4

29,1

33,5

1

0,011

0.76

0,60

0,25

0,98

0,03

0,04

»

10,5

30,5

33,0

2

0,042

0,71

0,12

0,27

1,00

0,02

0,035

0,20

»

11,2

32,5

34,4

3

0,095

0,73

0,28

0,25

0,93

0,05

0,03

»

13,0

31,1

33,1

АВ

0,46

0,72

0,25

0,34

0,91

0,05

0,03

»

14,2

33,5

36,4

В91

0,80

1,73

0,43

0,28

0,18

3,89

0,02

0,18

»

10,6

37,6

43,8

В95

1,73

2,46

0,29

0,39

0,32

6,4

0,13

»

8,5

53,0

58,1

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

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

и приблизительно соответствует базе 107 циклов.

Ниже сопоставлены значения ограниченных пределов выносливости на базе 107 циклов, найденные обычным способом и вычисленные по фор­ муле (3):

С плав

 

млз

А В Г

1

2

3

А В

В91

Предел выносливости

на базе

6,7

13,6

13,9

12,9

13,4

13,5

17,3

107 циклов, кГ/мм2

 

Экстраполированный

предел

6,3

13,5

14,0

12,5

13,0

14,1

16,3

в ы н о с л и в о с т и , кГ/мм2

Отклонение, %

 

+ 5 ,9

+ 0,7

—0,7

+ 3,1

+ 3,0

- 4 ,5

+ 5 ,8

Значение параметра т .

0,39

0,39

0,34

0,33

0,31

0,32

0,37

При определении предела выносливости по методу Про для произволь­ ных баз предлагается модифицировать уравнение (3) путем введения па­ раметра К , учитывающего базу. В связи с этим уравнение (3) принимает вид

а.1= К(бр^ Ват),

 

 

 

(4)

где ст^! — значение ограниченного

предела

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

для

заданной

базы. Параметр В связан с параметром А уравнения (3)

соотношением

В = АК. Значения коэффициента К для

различных легких сплавов в

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

рис. 1,

на котором видно,

что этот

коэффициент мало зависит от типа сплава. Для баз 106, 107, 5.107и108 циклов средние значения К соответственно составляют 1,51; 1,02; 0,88 и 0,85 при средней квадратичной ошибке бк, равной соответственно 0,065; 0,035; 0,032 и 0,030. Относительная средняя квадратичная ошибка во всех рассмотренных случаях не превышает 4%.

Метод Про предусматривает испытание образцов не менее чем при трех скоростях возрастания амплитуды напряжения, так как уравнение (3) содержит три неизвестных параметра. Результаты обработки показали, что наибольшая точность в определении предела выносливости достига­ ется испытанием образцов при четырех-пяти скоростях возрастания амплитуды напряжений. При выборе диапазона скоростей необ-

Г,5

7.3

7.2

7,7

/,2

0,3

0,8

V

т 7

70'н

Рис. 1. Зависимость коэффициента К =а__1_19 от базы

1 — МЛ5, оь = 22,1 кГ/ммг; 2 — АВ, = 36,4 кГ/лш2; 3 — АВГ, =

= 33,5 кГ/лш2’, Л— В91, аь =* 43,8 кГ/мм2; 5 — АВ16, <ть = 33,0 кГ/мм2;

6 — АВ22, о^ — 33,1 кГ/мм*; 7 — АВ25, аь = 34,4 кГ/лим2; 8 — В95,

а^ = 58,1 кГ/мм2

ходимо учитывать, что минимальная скорость возрастания напряже­ ния определяет длительность усталостных испытаний и оказывает суще­ ственное влияние на точность определения предела выносливости. По­ этому при испытании легких сплавов целесообразно принимать ее равной (2—4) А0~* кГ/мм2на цикл. При максимальной скорости возрастания ам­ плитуды разрушающее напряжение не должно превышать предела теку­ чести. Анализ экспериментальных данных показывает, что для легких сплавов максимальная скорость возрастания амплитуды не должна пре­ вышать (4—6)-10“5 кГ/мм2 на цикл. Приблизительно такой же диапа­ зон изменения скоростей был намечен Шваревым [7] для сталей.

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

персия амплитуды разрушающих напряжений; 8 2 - обобщенная диспер­ сия амплитуды разрушающих напряжений; д—■уровень значимости крите­ рия однородности дисперсий Бартлета [10].

Дисперсия разрушающих напряжений практически не зависит от на­ чального уровня напряжения и скорости возрастания амплитуды. Уров-

Статистический анализ результатов испытаний

Сплав

ос-105, кГ/мм*

п

°р, кГ/мм2

б8

на цикл

АВГ

0,34

10

16,81

1,021

1,70

10

19,39

0,463

 

 

3,34

10

21,13

1,692

 

0,332

25

18,48

1,156

АВ

1,66

20

21,72

0,875

3,32

20

23,85

1,000

 

 

6,64

10

25,86

1,079

МЛ5

0,4

20

8,18

1,330

2,0

20

9,99

1,432

 

 

6,0

15

11,55

1,293

 

0,332

10

17,21

0,366

1

0,87

10

18,11

0,401

1,66

10

19,44

0,423

 

3,32

10

21,01

0,553

 

6,64

10

22,54

0,550

 

0,332

10

17,78

0,500

 

0,87

10

18,08

0,314

2

1,70

10

19,23

0,351

 

3,32

10

20,72

0,318

 

6,64

10

22,50

0,358

3

0,332

10

16,50

0,540

1,83

10

19,44

0,654

 

 

6,64

10

23,18

0,438

Таблица 2

б2 <7, %

1,0357 19

1,079 40

1,356 97

0,454 99

0,369 94

0,543 85

ни значимости критерия однородности дисперсий Бартлета во всех рас­ смотренных случаях не опускаются ниже 19%; это обстоятельство дает возможность производить оценку обобщенной дисперсии.

Анализ кривых распределения разрушающих напряжений, приве­ денных на рис. 2 для сплавов АВ и МЛ5 (для остальных сплавов они име­ ют аналогичный вид), показывает, что для каждого сплава при различ­ ных скоростях возрастания амплитуды характер кривых распределения одинаков, различаются лишь средние значения разрушающего напряжения. Это означает, что при различных скоростях возрастания амплитуды и од­ ном и том же приращении вероятности разрушения приращение разруша­ ющего напряжения будет одинаковым; следовательно, параметры А и т не зависят от уровня вероятности, для которого определяется значение предела выносливости. Подсчеты показали, что параметр т для сплавов 1, 2, 3 и АВ при изменении вероятности в интервале 1—99% колеблется в пределах 0,31—0,35, а для сплавов АВГ и МЛ5 — в пределах 0,38— 0,41. Средние значения параметра А для сплавов АВГ, АВ, МЛ5, 1, 2 и 3 соответственно равны 430, 260, 240, 230, 190 и 230.

Рис. 2. Кривые распределения разрушающих напряжении для сплавов АВ нМЛ5.

Цифры у кривых—значения а, ?.Т/.и.и* пя цикл

Уравнение (4) и независимость параметров А и т от вероятности раз­ рушения позволяют оценивать дисперсию ограниченного предела выносли­ вости по результатам ускоренных испытаний.

Из теоремы о дисперсии [11] применительно к уравнению (4) следует, что дисперсия ограниченного предела выносливости для различных баз может быть выражена через дисперсию разрушающих напряжений фор­ мулой

(5)

где 8 *, —дисперсия ограниченного предела выносливости; «50а — дисперсия

разрушающих напряжений.

Втабл. 3 сопоставлены средние квадратичные отклонения ограничен­ ных пределов выносливости, найденных обычным методом и расчетным путем по уравнению (5) для сплавов АВ и МЛ5. Как следует из этой табцы, соответствие характеристик рассеивания пределов выносливости

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

Средние квадратичные отклонения ограниченных пределов выносливости, к Г /м м 2

 

Характеристика

 

лв

 

МЛ5

 

107

5-107

10я

107

 

 

Среднее

квадратичное отклонение

1,04*

0,92

0,88

1,16

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

1,07

0,94

0,89

1,05

 

 

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

14,1

12,4

12,0

6,3

13,5

11,9

11,3

6,7

Предел выносливости для вероятно­

11,7

10,5

9,8

4,2

сти Р = 1%

11,8

10,8

9,8

4,1

Нижняя граница предела выносли­

11,1

9,8

9,4

4,1

вости

[12]

10,7

9,3

9,2

3,9

 

 

* В числителе —по ускоренному методу, в

знаменателе — по обычному.

 

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

Характер кривых на рис. 2 и 3 свидетельствует о существовании ниж­ ней и верхней границ разрушающих напряжений. Преобразование рас­ пределения путем введения нижней границы разрушающих напряжений (ХРо и верхней границы <тРк приводит распределение случайной величины

х==^ё ' к нормальному (рис. 4).

Рк °р

Значение <тРо зависит от скорости возрастания амплитуды напряже­ ний, однако разность сгр — аРо так же, как и сгРк— бр, для каждого сплава

является устойчивой величиной. Строить функцию распределения преде­ ла выносливости и определять ее параметры по результатам ускоренных испытаний целесообразно в следующем порядке:

1) по средним значениям разрушающих напряжений для образцов, испытанных при трех = пяти уровнях скоростей возрастания амплиту­ ды, определяются параметры т, А и а ? ;

2) строится кривая распределения отклонений. Аар = сгр — ар и вы­ числяется дисперсия величины Дсгр, равная дисперсии разрушающих напряжений и экстраполированного предела выносливости;

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

х = щ1 - __СРо нормальному;

°Рк ср

4) для ряда вероятностей разрушения (например, при Р , равных 1,10, 30, 50, 70, 90 и 99%) с кривой распределения снимаются соответствующие значения разностей величин разрушающих напряжений (Дсгр)р;

- / о г г

Рис. 3. Обобщенные кривые распределения отклонений разру­ шающих напряжений от среднего значения

1 — сплав АВ* п = 75;

2 — сплав АВ16, п = 90;

3 — сплав АВ22,

?г = 50; 1 — сплав АВ25,

п = 00; 5 — сплав МЛ5, п =

55

■ол о ол о,в ц 16 -!,г -о»

о,о 0,0 0,8

Рис. 4. Преобразование распределения (см. рис. 3) путем введе­ ния верхней (ар и нижней (ар о ) границ разрушающих на­

пряжений

1 — сплав АВ, а р — <Тро = 3 кГ/мм*,

сгр.к

— йр в 2 кГ/мм8; 2 — сплав

АВ22, Ир—сГр о =

1,5 кГ/иш5, Яр.к._—^р =

 

*»5 кГ/ммг; 3 — сплав МЛ5,

<тр — Ср_о = 2,2

кГ/мм'-, ар#к — ар =

2,3

кГ/мм*

5)вычисляются значения ограниченного предела выносливости для

выбранного ряда вероятностей разрушения и базы по формуле

(3-1)р = К + (А<5р)р] . (6 )

На рис. 5 сопоставлены кривые распределения ограниченных преде­ лов выносливости сплавов АВ и МЛ5, построенные по результатам обыч­ ным и ускоренных испытаний. ;

Приведенные данные показывают на удовлетворительное соответст­ вие характеристик сопротивления усталости и их дисперсии, найден­ ных указанными методами. Расхождение в значениях пределов выносли­ вости не превышает 0,5 — 0,8 кГ/мм2, что соответствует средней квадра­ тичной ошибке определения предела выносливости обычным методом по результатам испытаний серии из 50 образцов.

Рис. 5. Сопоставление кривых распределения ограниченных пределов выносливости ускоренным (сплошные линии) и обыч­ ным (пунктирные) методами

1 — сплав МЛ5, N — 107; 2 — сплав АВ, N = 10*; 3 — сплав АВ, N = 107

При определении медианного значения предела выносливости уско­ ренным методом средняя квадратичная ошибка б может быть найдена по формуле, вытекающей из теоремы о распределении выборочного средне­ го [И]

где Оо.х — среднее квадратичное

отклонение предела

выносливости;

п — число испытанных образцов;

с — минимальное число

образцов,

не­

обходимое для оценки предела выносливости ускоренным методом;

сле­

дует принимать с = 3, что определяется числом параметров уравнения

(3).

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

нюю квадратичную ошибку б/о_г

следует принимать равной 4—6%,

что соответствует обычной ошибке

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

Анализ показал, что для сплавов со стабильными усталостными свой­ ствами (коэффициент вариации предела выносливости у 10%) число образцов п должно составлять 8—12, а для сплавов с повышенным рассея­ нием предела выносливости = 15—20%) — 12—20 образцов.

Указанный объем испытаний при использовании ускоренного метода обеспечивает такой же уровень случайной ошибки, как и обычный метод испытания того же объема. Однако надо иметь в виду и систематическую ошибку ускоренного метода, равную 4—6%, как это следует из резуль­ татов обычных и ускоренных испытаний на усталость большого количе­ ства образцов из ряда легких конструкционных сплавов.

Приведенные данные показывают, что модифицированный метод Про не дает экономии в числе образцов при опредёлении медианного значе­ ния предела выносливости.

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

С

 

 

7,9

 

7,8

 

о

 

 

7,7

 

А

 

 

/,8

о

А

 

 

 

!

д

о

/

 

7,5

 

/ х

 

 

 

 

 

2

 

7,Ь

А

 

 

Д 3

 

А

 

А

4

 

 

У

 

X

/

 

7,3

V«К А

 

 

 

 

 

5

 

 

&

 

 

7

 

 

 

 

 

 

1

 

 

 

 

 

 

5

7

ос, кГ/ммг/цикл

Рпс. 6. | Зависимость коэф­

фициента С от скорости возрастания амплитуды на­ пряжений для сплавов

1 — АВГ; 12 — АВ; 3—16; 4— АВ22; 5—25; 0 — В91; 7 — МЛ5

напряжений от среднего значения, как показали результаты статисти­ ческого моделирования, для сплавов со стабильными усталостными свойствами число образцов должно составлять 30—40, а для сплавов с высокой дисперсией усталостных свойств — до 40—60. Случайная ошибка определения предела выносливости для малых вероятностей разрушения будет такой же, как и при обычном методе испытаний 100—150 образцов. Таким образом, при построении кривой распределения пределов вынос­ ливости наряду с многократной экономией машинного времени достига­ ется двух-трехкратная экономия образцов.

Ускоренный метод усталостных испытаний, предложенный Эномото [6], предусматривает испытание одной серии образцов при постоянной скорости возрастания амплитуды напряжений. Определение предела вы­ носливости по этому методу основано на предположении, что отношение разрушающих напряжений к величине предела выносливости при посто­ янной скорости возрастания амплитуды напряжения является постоян­ ной величиной для всех материалов, т. е.

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

Проверка этого метода для различных марок стали и чугуна [6] при значениях а в пределах 10_б—105 кГ/мм2 за цикл показывает, что сред­ няя квадратичная ошибка составляет 7—12%.

Для оценки возможности применения метода Эномото к легким спла­ вам определялся коэффициент С при различных скоростях возрастания амплитуды напряжения для ряда алюминиевых и магниевых сплавов. За исходную базу испытаний было принято 107 циклов. Зависимость зна­ чений коэффициента С от скорости возрастания амплитуды показана на рис. 6.

Из приведенных результатов видно, что при низких значениях а — в пределах (0,3—1)-10~5 кГ/мм2 за цикл — относительная средняя квад­ ратичная ошибка составляет для исследованной группы сплавов около 5%, а при скоростях выше 10 ~5 кГ/мм2 за цикл указанная ошибка дости­ гает 10%. Эти ошибки следует считать систематическими, так как на одном уровне скорости возрастания амплитуды напряжений испытыва­ лось достаточно большое число образцов, обеспечивающее определение медианы разрушающих напряжений с относительной средней квадратич­ ной ошибкой не менее 1—2%.

Таким образом, для ориентировочной оценки медианы предела вынос­ ливости для базы 107 циклов можно использовать метод Эномото. При одной скорости возрастания амплитуды достаточно испытать 4—5 образ­ цов. Для определения предела выносливости на большей базе можно вос­ пользоваться коэффициентом /*Г, применяемым в модифицированном ме­ тоде Про. Систематическая ошибка в этом случае возрастает до 8—12%. Предел выносливости для произвольной базы определится по формуле

Отмеченные выше закономерности рассеивания разрушающих напря­ жений при ускоренных испытаниях дают возможность по результатам испытаний 20—30 образцов при одной скорости возрастания амплитуды на основании формул (5) и (9) производить оценку дисперсии и функции распределения ограниченного предела выносливости.

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

Выводы

1. Использование уравнения Вейбулла для ускоренного определе­ ния предела выносливости показало, что для легких сплавов с ростом базы предел выносливости приближается к нулю. Уравнение кривой усталости для этих сплавов приобретает обычную степенную форму. Та­ кое уравнение соответствует результатам усталостных испытаний лег­ ких сплавов. Если найти его параметры по экспериментальным точкам с наибольшей долговечностью 107 циклов, то с точностью до 10% путем экс­ траполяции степенного уравнения можно определить ограниченный предел выносливости на базе 108 циклов. Экстраполяция при помощи ука­ занного уравнения на 1,5—2 порядка приводит к большим погрешностям.

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

3. Дисперсия разрушающих напряжений не зависит от скорости воз­ растания амплитуды напряжения и по дисперсии разрушающих напря-