книги / Статистические методы в строительной механике
..pdfПростейшей проверкой условий типа (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] |
|
|
|
2к |
(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 |
= f—f |
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 |
|
Приводимые в литературе данные дают обычно большой разброс значений рх/?г- • Этот разброс обусловлен, очевидно, ярко выраженным вероятностным характером разрушения от усталости. Здесь нужны систематические массовые эксперимен ты с надлежащей статистической обработкой результатов.