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

книги / Статистические методы в строительной механике

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

Простейшей проверкой условий типа (6.50) является одно­ ступенчатый режим нагружения (рис. 98). Пусть N(Si) и W(S2) — предельное число циклов для однородных режимов с напряжениями Si и S2. Условие (6.50) принимает вид

 

 

П\ , пг

 

(6.52)

 

 

У(50 _г^(5г)

 

 

 

 

 

т. е. оно

выражается

прямой АВ на

рис. 99.

Эксперименты

обычно обнаруживают

отступления от

прямолинейного закона,

характер

которых показан на том же

рисунке.

Предваритель­

ная «тренировка» образца при меньших напряжениях обычно приводит к некоторому упрочнению, т. е. к увеличению суммар­

ной

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

Напротив,

предварительное

нагружение с

напряжениями

S t,

большими,

чем

«контрольное»

напряжение

Ь'г, приводит

к

некоторому раз­

упрочнению.

Если

напряжения

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

ками. Следующий параграф мы посвятим этому использованию. Обобщение «гипотезы суммирования» и анализ погрешностей, связанных с ее применением, будут даны несколько позднее.

69. Применение теории суммирования повреждений к оценке долговечности конструкций

В работах [86, 99] условие (6.50) было применено для оценки долговечности конструкций при стационарном случайном режи­ ме с амплитудой, подчиняющейся нормальному закону распре­ деления. В статье [179] был рассмотрен режим с постоянной ча­ стотой и амплитудой, подчиняющейся распределению Рэлея. Бо­ лее общий случай узкополосного стационарного процесса с нор­ мальным распределением напряжений s(t) был рассмотрен ав­ тором [22], который привел задачу определения долговечности к табулированным функциям и показал, что при режимах с редки­ ми перегрузками формула, приведенная в [179], становится не­ пригодной. Ниже будет дано некоторое развитие результатов, содержащихся в работе [22].

Рассмотрим вначале стационарный случайный процесс s(t). Пусть Vo(S ) — среднее число превышений случайной функцией s ( t ) уровня S в единицу времени. Тогда среднее число циклов, амплитуда которых за время Т оказывается лежащей в преде­ лах между S и S+dS, будет равно

dn(5) = -

TdS.

(6.53)

Для того чтобы выражение (6.53)

можно

было подставить

вместо dn в условие (6.50), необходимо, чтобы процесс накоп­ ления повреждений был достаточно м едленны м процессом . Точ­

 

нее,

необходимо,

чтобы

срок

 

службы

Т можно

было

разбить

 

на интервалы At, достаточно

ма­

 

лые,

чтобы

приращением

ДD

 

внутри этих

интервалов

можно

IgRf

было пренебречь,

но достаточно

большие, для того чтобы реализа­

 

ция

случайного

процесса

s(t)

 

продолжительностью At была до­

l3N,

статочно

представительной. На­

пример, если

предельное число

 

Рис. 100

циклов А/’п=104, то срок

службы

 

можно разбить на 100 интервалов,

в пределах каждого из которых мы имели бы 100 циклов. В пер­ вом приближении реализация, содержащая 100 циклов, может считаться представительной. Из приведенных соображений вы­

текает, что подстановка выражения

(6.53) в условие (6.50) мо­

жет быть оправдана, если Afn ^ 1 0 4.

Для ожидаемого срока

службы (ресурса) мы получаем формулу

т„ =

1

(6.54)

 

 

dV<> (S)

dS

 

dS

N(S)

Допустим, что случайный процесс s ( t ) подчиняется нор­ мальному распределению с математическим ожиданием, равным нулю, и стандартом os. Для кривой усталости N=N(S) возьмем выражение (рис. 100):

N = N1

если S ^ R f ,

(6.55)

N -> со ,

если S < Rf .

 

Здесь R f — предел выносливости

при симметричном цикле, Nt

и т\ — эмпирические константы.

Например, для углеродистой

стали N\ = 10б-т-.107, mi =6-т-12 (99].

Подставляя выражения (6.7) и (6.55) в 'формулу (6.54) и вводя безразмерные переменные

х

S

*0

*/

(6.56)

О •

 

S

 

 

 

получим

 

 

 

 

Тп =

1

NiTetf1

 

 

 

 

 

 

+1 е 2

dx

 

Интеграл, входящий в эту формулу, согласно (6.40) выражается через табулированные функции. А именно:

Гп —

^

(6.57)

 

^(mi + 2) Я(4

mt + 2)

 

Рассмотрим частный случай формулы *(6.57). Допустим, что частота процесса детерминирована и равна tooТогда в формуле

(6.57)

7,е= 2л/(оо. Предположим, что отношение x0=Rj/as доста­

точно мало, чтобы при заданном mj

можно было принять

Р ( 4

rrii+2) « 1. При этих предположениях мы приходим к фор­

муле Майлса [179]

 

 

(6.58)

 

ТП

Нетрудно видеть, что ошибка этой формулы может оказаться весьма большой. Так, если при m, = 8, as =0,5R/ ошибка состав­ ляет около 5%, то уже при а^=0,25 Rf формула (6.58) предска­ зывает примерно в 10 раз меньший срок службы, чем формула (6.57). Причина столь большого расхождения состоит, очевидно, в том, что, полагая нижний предел в интеграле

00 X*

1 * " + ,Г ‘г л«

х.

равным нулю, мы тем самым считаем, что R/-*0. Это равносиль­ но предположению, что даже сколь угодно малые напряжения вызывают усталостные повреждения. Чем ниже средний уровень напряжений, тем меньше вероятность перегрузок и тем, следо­ вательно, значительнее погрешность формулы (6.58).

Рассмотрим числовой пример. Пусть Rf = 20 KHJCM2, N I 10е, nii=8. Тогда при сг^ =10 к‘н/см2ожидаемый срок службы

10в -2а

= 0,703-10е.

384-0,9473

Если среднее

квадратическое

напряжение

уменьшается до

о^=5 KHICM2, то ожидаемый срок

службы увеличивается до

 

 

Т„

ЮМ»

 

1,71 • 10°,

 

 

 

 

Те ~~ 384-0,0996

 

 

 

 

 

 

 

т. е. возрастает примерно в 2500 раз.

Значительный рост долго­

вечности при уменьшении уровня

напряжений

подтверждается

экспериментом. Сошлемся,

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

нию прочности алюминиевых панелей,

помещаемых в сильное

акустическое

поле

работающего

реактивного

двигателя

[168].

При уровне шума

140 дб разрушение панели

наступало

через

Г„ =500 мин. С увеличением уровня до 154 дб этот срок умень­ шался до 14 мин. Наконец, при шуме 166 дб отдельные панели

разрушались спустя 30—40 сек.1

 

существенным

Полученные выше результаты могут быть

образом обобщены, если ввести функцию р0(5)

по формуле

Po(S) = Te

fWpjS)

(6.59)

 

dS

 

Тогда формула (6.54) примет вид

 

Т„ = — — ^

--------.

(6.60)

?' Ро 'S) <is

 

J

N(S)

 

Для узкополосных процессов функция Po(S) имеет смысл плотности вероятности максимальных значений процесса s(/). Приняв, что функция Ро(5) совпадает с распределением Пир­ сона

Po(S) =

S1-1

ехр

 

 

 

 

(6.61)

 

о“ V (а)

 

 

 

 

 

мы сможем описать весьма широкий класс

стационарных

слу­

чайных процессов (см. рис. 8).

 

 

 

 

 

 

Подстановка выражений

(6.61) и (6.55)

в формулу

(6.54)

приводит к формуле

 

 

 

 

 

 

 

7 \,=

N\TeX.™1V (g)

 

 

 

 

(6.62)

1* (Ш! +

а) Р (*£, пц +

а)

 

 

 

 

 

 

1 Уровень акустического излучения L измеряется

здесь в

децибеллах

(дб). Связь между уровнем шума

в децибеллах

и

акустическим

давлением

р в барах дается формулой

 

 

 

 

 

 

 

L = 201g - 2 - ,

Ро

в которой давление ро соответствует порогу слышимости. Отсюда следует, что при увеличении интенсивности от 140 до il66 дб акустическое давление уве­

личивается примерно в 20 раз.

Полагая в формуле (6.62)

а= 1 , получим

формулу для про­

цесса, амплитуды которого имеют распределение,

совпадающее

с положительной ветвью

нормального

распределения. Такой

процесс рассматривался в работах {86, 99],

где

интегралы вы­

числялись при помощи квадратур. Здесь

задача впервые приве­

дена к табулированным функциям. При а = 2 мы возвращаемся к формуле (6.57) и т. д.

График для отношения Тп fN{Teb логарифмическом мас­

штабе представлен на рис. 101. Если *0< 2 ,

т. е. если <^>0,5./^

зависимость довольно хорошо выражается

прямыми

линиями.

Их

аналитические выражения легко

получить,

полагая

Р(х%,

mi + a ) « l :

 

 

NxTex™' ¥ (а)

(6.63)

W(mx + a)

àia формула пригодна при режимах с высоким уровнем напря­ женности.

Для больших х0 можно получить другое упрощение. Исполь­ зуя формулу (6.43), найдем, что

_4 _

Irrti+ a = 6

2

[xîr,+e 2 “Ь { Щ

+

Я— 2)

4 “Ь

 

+

а — 2) (/их -J- а — 4) х%'+ 6 +

..•!•

Формула (6.62) принимает вид

 

 

 

 

Тп =

 

 

NiTLe 2

х.2—а

 

 

(6.64)

Ш1 + Ч-

2

, (mt + а -

2) (mt +

а — 4)

1 I

о

 

*i

 

4

 

 

Ha рис.

102 показаны значения

Та !N{Te, вычисленные по

точной формуле (6.62), по приближенной формуле (6.63) и по

асимптотической формуле (6.64); в последнем

случае удержи­

валось по одному, по два и по три члена ряда. Для вычислений

было принято, что mi = 10, а=2 . Как

видно из

графика, фор­

мула (6.63) пригодна при малых х0, а

формула

(6.64) — при до­

статочно больших Хо. Формула (6.63) дает в данном случае для

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

приближение снизу,

а формула (6.64) — прибли­

жение сверху.

 

 

который

входит в формулу

Необходимо отметить, что ряд,

(6.64),

вообще

говоря,

расходится

(исключение

составляет

случай

целого положительного m j+o,

при

котором

ряд содер­

жит конечное число членов). Поэтому,

используя его

для вы­

числений, следует строго

придерживаться

известных

правил,

относящихся к асимптотическим рядам. При отсутствии соот­ ветствующих навыков могут быть недоразумения. Однако, если *£<30, то нет необходимости применять формулу (6.64).

В этих пределах изменения х0 имеются доступные таблицы для функции Р(х£, п), а точная формула не сложнее приближен­

ной формулы (6.64). Если же х£ >30, то для грубой оценки может быть использована приближенная формула

4

е 2 ЧГ (а).

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

1 В машиностроении коэффициент запаса для неоднородных режимов переменных напряжений определяется обычно как

До сих пор рассматривалась задача об определении ожидае­ мого срока службы Т„. Можно поставить задачу об определе­ нии меры повреждений D при Т<Тп. Возвратимся к соотноше­ нию (6.49).

Обозначим его левую часть через v:

N

dn

-Jï (S ) ’

Для реализации стационарного случайного процесса тельностью Г ожидаемое значение v будет равно

00 .

dV0 (S)

dS

T

=ГЛ

dS

N (S)

T„

 

 

 

продолжи­

(6.65)

где Тп — ожидаемое значение срока службы. С другой стороны, из соотношения

D

следует, что D = D s (v). Учитывая (Ь.65), получим окончательно простую формулу для нахождения меры повреждений

D = Ds { ~ y

(6.66)

Вообще говоря, мера D является случайной функцией вре­ мени Т. Однако в силу сделанного предположения о крайне медленном течении процесса накопления повреждений ее дис­ персия весьма мала, и поэтому можно считать, что при одном и том жё Т значения D будут весьма плотно группироваться вокруг среднего значения (6.66). В самом деле, параметр v вы­ ражается через сумму N квазинезависимых случайных величин с конечной дисперсией

N

. = 2 — • /Я "<S/>

При больших N стандарт <*„ ведет себя как

и, следовательно, приближенно можно считать v детерминиро­ ванной величиной. Отсюда вытекает аналогичный вывод для D и Г„.

То обстоятельство, что в действительности срок службы Т„, а также зависящая от него мера D обнаруживают большой раз­ брос, объясняется другими причинами. Параметры кривой усталости Rf, N\ и mx являются случайными величинами и по­ этому срок службы Та (Rf , N1, mx) — случайная величина. Зная

совместную плотность вероятности р(Р/,

mx), можно найти

распределение случайных величин Т„ и D.

Некоторые

сообра­

жения, относящиеся к этому вопросу, будут

приведены

в п. 75.

70. Обобщение теории суммирования повреждений

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

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

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

Один из путей построения обобщенной теории состоит в от­

казе от гипотезы

об автомодельности

меры

/)_,(п. 68). Если

принять,

что функция D's

(D, S)

в

правой

части

уравне­

ния (6.46)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dD

_

D'

(D, S)

 

 

(6.67)

 

 

 

~dn

 

 

N (S)

 

 

 

 

 

 

 

 

 

 

 

 

зависит не только от D, но и

от

максимального

напряжения

цикла 5, то мы получим уравнение,

которое,

вообще

говоря,

не интегрируется

путем

разделения

 

переменных.

Отсутствие

интегрируемости

такого

рода

указывает на зависимость пре­

дельного

числа

циклов

Мп,

определяемого

из

условия

D(Nn) = \, от вида функции S = S (я), т. е. от истории нагруже­ ния. Следовательно, задача заключается в выборе надлежащей функции D ' (D, 5). Некоторые частные случаи этой функции

рассматривались в работах [137, 178, 184].

Допустим, что при S = const имеет место зависимость

 

 

 

г »

i *(S)

 

 

 

 

 

{ 6 Щ

гд ер (5 )> -1 — некоторая

функция

напряжения 5. Уравнение

(6.67) принимает вид

 

 

 

 

4P

_

р (S)

P ( S ) - l

 

п ? (S)

(6.69)

dn

~~

N(S)

 

 

 

Проинтегрируем уравнение (6.50) для одноступенчатого ре­ жима (см. рис. 85). На первом этапе

гдер1=р(5!). Интегрируя уравнение на втором этапе (при

N„ , S = S2 const) *и определяя постоянную инте* грирования, найдем

D

= ff

f t

^' N(Ss)O

J

 

U N(SOi)

J

Условие разрушения принимает вид

 

 

р»

 

 

 

При рх = р8 = const мы вновь

получаем

зависимость (6.52).

Если же pi^pa»

то мы получаем

семейство кривых, представ­

ленных на рис. 99. Учитывая обнаруживаемое на опыте неко­ торое увеличение долговечности при предварительном нагру­ жении с напряжениями, меньшими контрольного напряжения, и уменьшение долговечности — в противоположном случае, при­

ходим к выводу, что показатель

р (S)

должен быть

невозра­

стающей функцией напряжения S. Эта

функция

может быть

найдена из опыта на одноступенчатое нагружение

 

 

 

 

 

PL =

lg l l

N (SQ

]

 

 

 

 

 

 

Pa

 

ni

 

 

 

 

 

Очевидно, что показатель p(S)

определяется с

точностью

до

множителя, который может выбираться произвольно.

 

 

В качестве примера в табл. 6

приведены

результаты обра­

ботки одной

из

кривых

Коммерса [56].

Для

этой

кривой

Si —1,29Rj,

S2= \,\R J .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

6

Пг

hr

п±

Пг

1

Г 1

Пг

1

Pi

 

N { S ,)

8 N{St)

N (5,)

e

1 1

mst) J

р*

 

0,1

—1,000

0,78

 

—0,658

 

0,66

 

0 .2

—0,699

0,61

 

—0,409

 

0,62

 

0,3

—0,523

0,52

 

—0,319

 

0,61

 

0,5

-0 ,3 0 1

0,36

 

—0,194

 

0,65

 

 

 

 

 

 

Среднее

 

0,635

 

Приводимые в литературе данные дают обычно большой разброс значений рх/?г- • Этот разброс обусловлен, очевидно, ярко выраженным вероятностным характером разрушения от усталости. Здесь нужны систематические массовые эксперимен­ ты с надлежащей статистической обработкой результатов.

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