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

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

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

Возможны и другие предположения о виде функции Ds, от­ личные от того, которое дается формулой (6.68). Однако во всех случаях для учета истории нагружения необходимо учиты­ вать явную зависимость меры Ds в правой части уравнения (6.67) от напряжения 5.

Покажем теперь, как применить обобщенную теорию сумми­ рования повреждений к оценке долговечности при напряжениях, представляющих собой узкополосный стационарный случайный ■процесс. Поскольку здесь условие разрушения не может быть представлено в виде, подобном (6.50), будем исходить непосред­ ственно из уравнения (6.67). Разобьем полный интервал време­ ни Т на интервалы Дtk, достаточно малые, чтобы изменение ме­ ры D внутри каждого интервала было мало, но достаточно боль­ шие для того, чтобы реализацию процесса s (t) продолжительно­ стью Дtk можно было считать вполне представительной. Пусть в течение времени Д/Л имеется в среднем Д«1 циклов с напряже­ нием Sj, Дп2 циклов с напряжением S2 и т. д. Тогда приближен­ но можно записать

до*

Ps (Pfe* '-ty)

(6.71)

à П и

/

N ( S j )

 

 

 

С другой стороны, среднее за

интервал Д/Ачисло

циклов, ам

&Sj &Sj

плитуды которых лежат в пределах между о у— — - и оу- г 'у _,

очевидно, будет

Дпу

dV0 (Sj)

à S Л Ь .

(6.72)

dS

 

 

Подставляя выражение (6.72)

в

формулу (6.71) и переходя к

пределу при Д5у- 0, получим

 

 

LDk

Ps (Ри>

S)

dV0 (S)

N (S)

 

dS Д tk.

 

 

dS

Заменяя приращения дифференциалами, получим окончательно

dD

_

?’*

P's (D . S)

dVо (S)

(6.73)

~dt

~

J

N (S)

dS.

dS

 

Предельный срок службы Tnнайдем из условия D(Ta ) = 1:

Г„ =

 

dD

(6.74)

 

-------.

Î

Ds (D . S)

dV0(S)

 

N (S)

dS

 

 

dS

 

Если D s’ (D, S) не зависит явно от S, то мы, очевидно, возвра­ щаемся к формуле (6.54).

71. Теория двух стадий усталостного разрушения

Несмотря на простоту исходных гипотез и основных урав­ нений, практическое применение обобщенной теории накопле­ ния усталостных повреждений затруднительно. Это объясняется необходимостью двухкратного вычисления квадратур в фор­ муле (6.74), причем одна из них должна вычисляться в функ­ ции от параметра D.

Можно выдвинуть другую теорию, приводящую к более удобным расчетным формулам. Хорошо известно, что процесс усталостного разрушения можно представить в виде двух ста­ дий, первая из которых является подготовительной. После того как условия для возникновения усталостной трещины подготов­ лены, наступает вторая стадия, в течение которой происходят развитие и углубление трещины. Вследствие этого имеет смысл описывать процесс накопления повреждений лри помощи двух функций. Одну из них, обозначаемую в дальнейшем через G, можно трактовать как меру разрыхления зерен или как меру подготовки металла к образованию усталостной трещины. Она равна нулю для начального (неупрочненного) состояния метал­ ла и становится равной единице, когда подготовительная стадия закончена. Вторая характеристика, которая по-прежнему обо­ значается через D, является мерой развития усталостной трещины. Равная нулю при G< 4 , она обращается в единицу в момент полного разрушения (рис. 103). Таким образом, вместо уравнения (6.28) получаем два уравнения:

 

 

f

- / К О . *

 

(6.75)

dD

 

 

0,

 

если

G<

1,

dn

 

f (G,

D,

S),

если

G

1,

где fo и f — функции,

которые в дальнейшем подлежат конкре­

тизации1.

 

 

допущения:

1)

предположим, что

Введем два следующих

кинетика развития трещины слабо

зависит от меры разрыхле­

ния G, если G>1;

2)

примем,

что

при

однородном режиме

(S = const) меры Gs

и Ds

не зависят явно

от напряжения S,

а зависят лишь от отношений

соответствующих чисел циклов:

1 Имеются основания подразделить подготовительную стадию в свою очередь на две стадии: стадию упрочнения и стадию разрыхления. Описание процесса усталостного разрушения при помощи трех функций будет, оче­ видно, еще более полным. Вообще теория легко обобщается на любое ко­ нечное число стадий*

Здесь No(S)— число циклов, необходимое для зарождения усталостной трещины, N(.S)— по-прежнему предельное число циклов. В отличие от теории суммирования повреждений здесь предполагается, что автомодельность имеет место лишь для каждой стадии в отдельности (рис. 104).

's," s

Используя высказанные гипотезы, придем к уравнениям ти­ па (6.75) :

dG

__

Gs (G)

 

dn

 

N0(S)

 

dD

0,

если G <

1, }

D’.(D)

 

(6.76)

dn

 

 

1.

N ( S ) - N 0(S) * если G>

Проинтегрируем первое из них от нуля доМ,, (Nn„—число циклов, необходимое для зарождения трещины при сложном режиме нагружения), а второе — в пределах от N„0 до Ntt. Тогда по­ лучим:

о

dn

Г

______ dn_____

= 1.

(6.77)

N0 (S)

N ( S ) - N 0(S)

 

 

 

 

 

 

ло

Из первого условия найдем N n„ после чего из второго — предель­ ное число циклов Nп .

Применим условия (6.77) « одноступенчатому режиму. Если Hi<W0(Si), т. е. если в течение первого этапа нагружения ин­ кубационная стадия еще не закончилась, то вместо формул (6.77) получим

 

п\

JVn0 — ni

_ j

ni + л2 — Np0

 

N0(Si) +

N0(Si)

N (Sa) — N0(Sa)

 

Если же ni >

JVo(«Si), то

 

 

 

N*о

, .

ni

N n 0______ j _________ n*

__ t

N0(So

- l >

N {Si)-No (Sx)

^ N{S2) - N o { S t)

 

N(S,)N0(S2)

ni

п2

= 1,

если /ii< W 0(Si),

 

 

N0(SI)N (cS)a

tf(S,)

+ N(S2)

N (S^

 

n2

 

 

 

 

 

■——----- 1- '

N(St) - N 0( S i )

 

=

1,

 

 

 

N (Si)

N (S2) y—1Nо1^2)

N (5i)

 

N(SZ)

 

 

 

 

 

 

если ri! > N0(5x).

 

 

 

 

 

 

 

 

В результате получаем

диаграмму

 

изменения

суммарной

долговечности, приведенную на рис. 105.

Вместо плавных кри­

 

 

 

вых, которые

вытекают

из

 

 

 

обобщенной теории

 

сумми­

 

 

 

рования

повреждений

(см.

 

 

 

рис. 99), мы

получаем

ло­

 

 

 

маные линии; однако харак­

 

 

 

тер

влияния последователь­

 

 

 

ности нагружения

остается

 

 

 

прежним. Это

влияние,

сог­

 

 

 

ласно данной теории, являет­

 

 

 

ся простым следствием

того

 

 

 

обстоятельства,

что

при

бо­

 

 

 

лее

высоких

напряжениях

 

 

 

подготовительная стадия за­

 

 

 

нимает

относительно

мень-

 

 

'

шую часть от общей

долго-

 

 

вечности, чем

при

 

низких

 

 

 

напряжениях (см. рис. 104).

 

 

 

Заметим, что наличие точек

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

Если надлежащим образом выбрать параметры, то обе пред­

лагаемые теории дают примерно

одни и те

же

результаты.

В настоящее время трудно судить,

которая

из

этих теорий

лучше описывает реальный процесс. Правильнее было бы ска­ зать, что обе они являются лишь приближенным феноменоло­ гическим описанием весьма сложного процесса накопления усталостных повреждений, учитывающим в отличие от теории суммирования влияние истории нагружения. Заметим, что в рамках этих теорий можно описать также и явление «аномаль­ ного» упрочнения (сравните штриховые кривые на рис. 99 и 103).

Выведем формулу для ожидаемого срока службы, отправ­ ляясь от условий (6.77). Подставляя в первое из них вместо dn

 

 

dn(.S) = - dv°f>

T,dS,

 

 

 

 

 

dt>

 

 

 

 

a во второе

 

 

 

 

 

 

 

 

dn(S) = -

dV° f - -(T n- T

0)dS

 

 

 

 

 

flo

 

 

 

 

(T0— продолжительность

инкубационной стадии)

и исключая

затем Т0, получим

 

 

 

 

 

 

Тп

<*MS)

 

+

00

1_______

 

-. (6.78)

n

 

dS

Cl dV0 (S)

I

dS

J

dS

N { S ) - N 0{S)

J |

dS

I

N0(S)

0

 

 

 

0

 

 

 

При N0-+0, что соответствует переходу к теории суммирования, мы возвращаемся к формуле (6.54).

72.Применение теории двух стадий усталостного разрушения. Сопоставление с теорией суммирования повреждений

Поскольку, как уже указывалось, обобщенная теория сум­ мирования й теория двух стадий усталостного разрушения дают примерно одни и те же результаты, мы воспользуемся второй из них, как приводящей к более удобной для вычислений фор­ муле (6.78). Рассмотрим стационарный случайный процесс с гауссовским распределением для мгновенных напряжений s(0 - Среднее число превышений уровня 5 для этого процесса опре­ деляется по формуле (6.7). Для кривой усталости N=N(S) примем аппроксимацию (6.55) и, кроме того, положим, что:

 

N0 = N2

)

если 5 >

Rf ,

(6.79)

 

N0-+ оо, если S < Rf .

 

 

Здесь No и пг2— эмпирические

константы

{N2^CN\,

условие

Тогда,

подставляя формулы

 

(6.11), (6.55)

и (6.79) в

(6.78)

и вводя, как и ранее,

безразмерные

переменные (6.56),

получим

 

 

 

 

 

Тп = TeNlXp

 

1

1

(6.80)

 

 

Я*

Z (тг , x 0)

 

Y

 

 

 

 

- ,

 

 

При этом Z(m2, x0) = 'F(2) P(xl, m2-f 2),

Y ^ти m2,

,

x0

 

xm*+

1 e

dx

(6.81)

 

 

 

 

 

 

 

 

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

Интеграл

(6.81)

должен

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

численными методами.

Лишь в случае, когда т2—т\,

\ _ 2Z (m ь JC0)

 

 

 

У(т „

т„, ^

 

 

 

° ) ~ ~ У г

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— -

 

 

Однако тогда формула

(6.80) переходит

в формулу

(6.57).

Этот результат

является вполне

естественным, так

как при

пц = т2 имеет место

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

меры Ds.

Согласно до­

казанному

в п. 68,

возвращаемся

к

теории

суммирования

повреждений. Следовательно, интерес

представляет лишь слу­

чай

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.

106

 

 

 

 

 

Рассмотрим

численные

примеры.

Пусть R f = 20 кн/см2,

/У1 = 10®, тi= 8. Оставляя эти данные

неизменными,

произведем

вычисления

для

различных as,

N2 и т 2.

Предположим,

что

os =10 mjcM2, W2=0,8 106,

т2= 8

(что соответствует теории

суммирования),

9,

10 и 11. Кривые

усталости

N=M(S) и

N0=No(S) в

логарифмическом

масштабе

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

на

рис. 106. Интегралы

(6.81)

вычислены по

формуле

трапеций,

а интегралы Z(m2, х0) — по таблицам [43, 101].

 

 

Вычисленные значения

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

Тп представлены

на

рис. 107. Существенно, что

учет влияния

истории нагружения

всегда приводит

к некоторому снижению

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

по

сравнению с той, которая предсказывается

теорией суммирова­

ния повреждений.

Так, при

ш2 —11 вычисления по

первой из

формул (6.80) дают Т„ = 3,87

104 сек, в то время как по теории

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

Тп =4,42 - 104 сек.

Разница

со­

ставляет, таким образом, 14,2%. Аналогичные результаты по­

лучаются и при других

значениях

параметров.

Так,

при os=

= 10 кн/см2, N2= 0,9 • 106 и m2= l l

находим, чтоТ„ =3,58 • 104 сек,

в то время как

значение Т п по теории суммирования

остается

неизменным. На этот раз разница составляет около 23,4%.

Были выполнены вычисления и при других

значениях пара­

метров. Для

случая

Rj= 20

кн/см2,

о5 = 5

кн/см2,

N: = 106,

mi = 8, N2 = 0,8- 106 результаты

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

на рис.

108. И

здесь имеет место снижение долговечности.

 

 

 

 

 

 

г„ сек

1

1

 

,

 

 

 

Ч-Ю'|---------

 

 

 

 

 

 

Ыг *0.8-10 6

 

 

 

 

 

 

 

нг =0 84 *

 

 

 

 

 

 

Nt •'00,9'S - t 0 s

"

— .

 

 

 

 

8-Ю7----------------------------

 

 

 

 

 

 

 

5

Ю1----------------------------

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

и-Ю’----------------------------

 

 

 

 

 

 

 

г-ю1----------------------------------

 

 

 

 

 

 

 

о I

_|___ ^

н

mt

 

 

 

 

â

з

ю

Рис. 107

 

 

 

Рис. 408

 

 

Полученные

результаты легко

допускают физическое ис­

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

Формула (6.80) неудобна для вычислений, поскольку ее использование связано с вычислением трудоемких квадратур.

Было бы желательным получить простую оценку для долговеч­ ности снизу. Используя преобразование

В Д 4

и замечая, что при

х0и m2 > т .

получим оценку снизу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.82)

(знак

равенства

имеет место

при m i=m 2). В

отличие

от фор­

мулы

(6.80) эта

формула содержит

лишь

табулированные

функции. К сожалению, формула (6.82)

обычно дает чрезмер­

ную заниженную оценку. Так, при R f = 20 KH/ CM2, Qs —10

кн/см2,

N1=106, m i=8, N2=0,8’ 10е,

т 2= 9 оценка

снизу

оказывается

равной Т п =0,495* 10е Те вместо точного значения

0,671

10е Те.

Поскольку Rf, Nu Ns, гп\

и m2 являются

случайными вели­

чинами, а долговечность оказывается чувствительной к их ма­ лым изменениям (это особенно относится к пределу выносли­ вости Rf), то во многих случаях расхождения в 10—20% сле­ дует признать несущественными. Тогда может быть использо­ вана оценка для долговечности, которая дается теорией сумми­ рования повреждений.

73. Оценка долговечности при квазистационарных случайных режимах напряжений

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

Например, средняя скорость ветра и среднее давление ветра на поверхность меняются медленно по сравнению с пульсационными составляющими ветра и давления. Можно привести

и другие примеры, из которых видно, что режимы напряжений,

близкие

к стационарным, являются

типичными режимами для

условий

эксплуатации. Поэтому их

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

практический интерес.

 

Мерой близости процесса s(t) к

стационарному случайному

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

теристик за время т, достаточное для

того, чтобы корреляция

между значениями функции

s(t)

стала

достаточно малой,

а также для того, чтобы реализацию

процесса продолжитель­

ностью х можно было считать

достаточно

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

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

Если исключить режимы с очень высоким уровнем напря­ жений, когда разрушение наступает после небольшого числа циклов, то можно считать, что изменение меры повреждения за время At также является малым. Таким образом, при квазистационарных режимах случайных напряжений имеют место три процесса: быстрое изменение напряжений, медленное изменение вероятностных характеристик .режима напряжений и медлен­ ный рост усталостных повреждений. Два последних процесса должны быть рассмотрены совместно.

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

тервала некоторыми средними значениями G и D. Пусть в те­ чение времени At имеется Дпх циклов с напряжением S|, Дп2 циклов с напряжением S2 и т. д. На основании (6.75) можно приближенно записать:

A G ^ f 0(G,D, Sj)Atij,

J

Л£>~ 2 /(G, D, SJ)

Anj.

/

 

Подставляя сюда вместо Дп} выражение

(6.72), зависящее от

t как от параметра, переходя к пределу

при Д5/-Ю и заменяя

приближенно конечные разности соответствующими дифферен­ циалами, получим:

 

Д

S)

dV0 (S; t)

dS,

 

dS

 

о

 

 

 

 

 

(6.83)

 

 

 

 

dD

= ] н о , D,

S)

dV0 (S; t) dS.

dt

 

 

dS

 

Простейшее предположение о виде функций / и /о приводит к выражениям:

 

1

О,

если

G<

1,

 

/о —

1

 

_

,

(б-8*)

. / = {

если

 

No (S)

 

G>

1.

 

N ( S ) - N 0{S)

Рассмотрим гауссовский случайный процесс, близкий к ста­ ционарному. Пусть Те (t) и os (t) — эффективный период и стан­ дарт напряжения, которые согласно предположению являются функциями времени, медленно меняющимися по сравнению с s(t). Подстановка выражений (6.55), (6.79) и (6.84) в урав­ нения (6.83) дает

 

 

 

dG

__ г \щ , *0 (/)]

 

 

 

 

 

 

 

dt

N%Te (t) XQ1* (t)

 

 

 

 

 

 

 

 

О,

 

если

G <

1,

 

 

 

dD

 

щ ,

, x0(t)

 

 

 

 

(6.85)

 

 

dt

 

 

если

G>

1.

 

 

N,Te(t)xp(/)

 

 

 

 

 

 

 

 

 

 

Задача

определения

ожидаемого

срока службы

сводйтся

к интегрированию

уравнёний (6.85)

при

начальных

условиях

0(0) =D(0) =0. Ожидаемый

срок службы

Тп найдется из ус­

ловия

£>\ТП) = 1.

Уравнение,соответствующее

теории

суммиро­

вания,

получится

из первого

уравнения,

если заменить в нем

G на D,

на Ni и т2на т\\

 

 

 

 

 

 

 

 

 

dD

 

Z\tnx, *о(0]

 

 

 

(6.86)

 

 

 

dt

 

NxTe

(О*

 

 

 

 

 

 

 

 

 

 

 

Для стационарного случайного процесса Те и х0 постоянны; от­ сюда легко получается формула (6.57).

Рассмотрим иллюстративный пример. Пусть R f = 20 кн/см.2, ЛЛ1= IО6, mi =8, #2=0,8 • 106, /«2=9. Предположим, что эффек­ тивный период постоянен и равен 0,01 сек, а процесс s(t) со­ стоит из двух реализаций узкополосного стационарного случай­ ного процесса со средними квадратическими напряжениями = = 10 кн/см2 и « ,= 5 кн/см2 (рис. 109). Строго говоря, в окрест­

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