книги / Статистические методы в строительной механике
..pdfдля s. Запись большинства предыдущих формул весьма упро щается. Например, формула (6.24) (Принимает вид
|
|
|
|
(уо - 1)а |
|
|
|
Р (х < |
1 |Г ) я г - ^ ехр |
|
п |
mik |
|
(6.28) |
|
|
|
1О |
2 2 |
|
|
||
|
|
[ |
|
J. fe=l |
|
|
|
Рассмотрим иллюстративный числовой |
пример. |
Допустим, |
|||||
что на сооружение действуют два |
вида |
стационарных |
случай |
||||
ных натр узок (для |
определенности |
будем |
говорить о |
ветровой |
|||
и снеговой нагрузках, принимая, |
что опи не коррелировапы). |
||||||
Эффективный период для |
ветровой |
нагрузки qi (0 |
примем рав |
||||
ным 7*1 =0,25 мес., |
а для |
снеговой |
нагрузки qz(t) примем соот |
||||
ветственно 7'2= 12 мес. (рис. 92). |
Далее, пусть тц=0,2, m^z— |
||||||
=0,2. Тогда приведенный период равен |
|
|
|
|
т. е. он существенно больше периода быстро меняющейся ветровой натрузки. Разумеется, з]щсь имеет значение соотноше ние между коэффициентами изменчивости. Имея в виду неко торые опубликованные работы, следует заметить, что утвержде ние о том, что в качестве характерного периода Т0 можно взять минимальный период изменения найрузки, является необосно ванным.
Примем, чго Г=240 мес., х0=4. -Вероятность отказа по формуле (6.28) будет
Р (* < 1 ) |
240 |
ехр Г-------- |
0,0048. |
|
|
0,645 |
I |
2 (0,2 + 0,2) |
J |
66. Применение закона распределения случайных перегрузок для оценки прочности и жесткости трубопроводов
В п. 61 была изложена теория деформаций упругих трубо проводов, прокладываемых в статистически .-неоднородном грунте. Здесь, используя результаты из п. 64, покажем, как можно оценить прочность и жесткость трубопроводов.
Пусть детерминистическое условие прочности сводится к требованию того, чтобы уп-ругий изгибающий момент М(х), яв ляющийся случайной функцией координаты х, ни в одном се чении участка трубопровода длиной L не превысил по модулю предельно допустимого уровня М*. В вероятностной постановке этому условию соответствует требование достаточно малой ве-
оятностисобытия, |
состоящего в том, что |
на |
интервале |
||
5 |
L предельный |
уровень М к будет превышен |
по модулю |
||
хотя |
бы оди-н раз (.рис. 93). Обозначим |
эту вероятность через |
|||
P(\M\>MJ L). Допустим, что процесс |
М{х) |
имеет |
центриро |
ванное нормальное распределение. Тогда для плотностивероят ности р(М) имеем формулу
где в м — среднее квадратичеокое отклонение момента:
Средний 'квадрат кривизны х2 определяется по второй из фор мул (5.125).
Оценка сверху для'малых вероятностей выброса стационар ного гауссовского .процесса на заданном интервале выведена в п. 64. Используя эту оценку, получим
2 L |
( М г \ |
(6.29) |
Р ( |М |> М .|1 ) = — |
e x p l - j ^ j . |
Здесь Lx — параметр, имеющий размерность длины и аналогии
ный эффективному периоду |
(6.8) : |
|
|
Фх (<i>) d о) |
V f |
|
|
|
Lx— 2х |
оо |
(6.30) |
|
|
|
|
Ф% (**) d Cl) |
|
a Ф* (со) — спектральная плотность кривизны оси трубопровода
(5.124). |
Параметр L x характеризует |
изменчивость кривизны |
||||||||||
трубопровода |
(а также изгибающего момента) .по его длине. |
|||||||||||
|
Расчет на |
прочность сводится |
к |
вычислению |
вероятности |
|||||||
Р(|Л4|>Л4*| L) |
по формулам |
(6.29) |
и |
(6.30) |
и « |
их |
сопоставле |
|||||
нию с нормативными значениями. |
|
|
|
|
|
|
||||||
|
Аналогично производится расчет на жесткость по предель |
|||||||||||
но допустимым значениям прогиба |
|
Вместо формулы (6.29) |
||||||||||
получаем |
|
|
|
|
|
|
|
|
|
|
||
|
|
P ( |a 'i ! > o '* U ) ~ |
2L |
|
|
( |
иР |
\ |
(6*31) |
|||
|
|
— |
e x p [ - —v |
) , |
||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Фа,(ш) d<ù |
|
|
|
|
||
|
|
|
Lw= 2 я |
|
|
|
|
|
|
|
(6.32) |
|
|
|
|
|
|
Фа ) (м ) <ü»da» |
|
|
|
|
|||
а |
= |
wï определяется |
по |
первой |
из |
формул |
(5.125). |
|||||
|
Если спектральная |
плотность |
Фг (и) |
Для |
неоднородного |
|||||||
грунта |
задана |
в виде |
(5.133), |
то все |
интегралы, |
входящие в |
формулы (5.125), (6.30) и (6.32), берутся в конечном виде. Со ответствующие формулы -ввиду громоздкости не выписываем.
Используя формулу преобразования (1.33), легко найдем, что плотность вероятности для обобщенной деформации будет
/7(Д|Г) = р[*(Д )|Л |
dà |
(6.34) |
|
|
Математическое ожидание обобщенной деформации Д к концу интервала времени Т определяется по формуле
|
Д (Г )= |
J <p(S)p(S|7)<iS. |
(6.35) |
|
So (Г) |
|
|
Для |
наиболее вероятного |
значения обобщенной |
деформации |
До(Т) |
получим формулу |
|
|
|
MT)æ<p[S(>(7)]. |
(6.36) |
|
Формула применима во всяком случае, если |
|
||
так как при этом условии |
максимум плотности |
вероятности |
|
р(Д| Т) непременно достигается на границе S = S 0(T). |
|||
В |
качестве примера рассмотрим процесс накопления оста |
точных деформаций в жестко-пластической системе с линейным законом упрочнения. Для такой системы
|
-!>(Л)=-Я, + *Д, <p(s)= * ~ * ' . |
|
R |
где |
— предел текучести, k — модуль упрочнения (s !>/?,) |
Для среднего и наиболее вероятного значений деформации лег ко получаем формулы:
00 |
|
-(г>= J -S— Rl |
p(S\T)dS= Jj= Jk . |
St (T) |
(6.37) |
М Л |
S0( T ) - R T |
|
Допустим, что процесс s(t) является стационарным гауссов ским процессом. Вычислим для него среднее 5 (Т)
• S (T) = |
j Sp{S\T)dS. |
|
|
S'(T) |
|
Подстановка сюда формулы (6.13) дает |
|
|
J |
S!exp - (S2~ f ’ 1d S - |
( 6 M ) |
S ' (Г) |
S |
4
1„=е 2 [XQ~2 + (п— 2) лгГ4 + (п - 2) (п - 4) * Г 6 + ... ). (6.43)
Носогласно (6.14) и (6.39)
Рлс. |
96 |
|
|
|
откуда формула (6.41) принимает вид |
|
|
|
|
« (Г) ■* * + 2«. + |
"§+- - |
+ 0 [ - V |
) - |
(6-44) |
|
*о |
\ Х° |
/ |
|
Эта весьма простая формула дает удовлетворительные резуль таты уже начиная сТ= 100 Те.
На рис. 96 представлена зависимость средней деформации
Д(Т) и наиболее вероятной деформации До(Г) от отношения Т/Те. При построении принималось а=0,25Я т, оЛ= 0,25/?т.
68. Накопление усталостных повреждений. Теория суммирования повреждений и ее анализ
Если средний уровень напряжений сравнительно невысок, так что перегрузки с выходом из упругой стадии являются прак тически невероятным событием, а число циклов достаточно ве-
лико, то необходимо считаться с опасностью усталостного раз рушения. Особенно это относится к знакопеременным режимам, для которых, как известно, предел выносливости обычно бывает меньше предела упругости.
Возникает трудный вопрос о том, что следует принять за ме ру усталостного повреждения. Казалось бы, в качестве такой меры можно было бы взять параметр, характеризующий отно сительную глубину усталостной трещины или относительное ос лабление рабочего сечения вследствие развития трещины. Из вестно, однако, что процесс накопления усталостных поврежде ний начинается с первых циклов и состоит из двух стадий [99]. Первая, подготовительная стадия состоит из сочетания двух про цессов: упрочнения наиболее слабых зерен и образования мик роскопических сдвигов, т. е. того «разрыхления», которое пред шествует образованию прогрессирующей микроскопической тре щины. После того как условия для возникновения трещины под готовлены, наступает вторая стадия, в течение которой происхо дят развитие и углубление трещины. Вторая стадия составляет меньшую часть продолжительности жизни образца, между тем как повреждения накапливаются постепенно в течение всей его жизни. Следовательно, мера усталостного повреждения, безус ловно, должна учитывать повреждения, накапливаемые на пер вой стадии.
Как уже указывалось в третьей главе, удовлетворительной статистической теории, описывающей процесс усталостного раз рушения, не существует. Поэтому приходится прибегать к фено менологическому описанию, согласующемуся с более или менее широким кругом опытных данных.
Пока ограничимся случаем, когда прочность конструкции целиком определяется номинальным напряжением s(£) в некото ром наиболее опасном месте, и примем, что напряжение s{t) ме
няется по симметричному циклу со случайной амплитудой. Вве дем априорную меру повреждения D, равную нулю для началь
ного состояния металла и единице при полном разрушении. Ме ра D является, очевидно, неубывающей функцией времени (для
циклов с недогрузкой, не вызывающих необратимых изменений, мера D остается постоянной). Максимальное напряжение цик ла, в отличие от напряжения в любой момент времени s(t), бу
дем обозначать через 5. Приращение |
меры |
повреждения ДD |
||||
при л-м цикле напряжений зависит |
лишь |
от |
состояния кон |
|||
струкции, достигаемого |
к концу |
(п— 1)-го |
цикла |
(т. е. согласно |
||
нашему предположению |
от меры Dn- г), |
и от |
максимального |
|||
напряжения л-го цикла. Следовательно, можно записать: |
||||||
А Ц , * / ( А , - ь |
Sn), |
|
|
|
||
( л = 1 , 2, . |
.£>о=0). |
|
(6.45) |
Здесь f — некоторая функция, подлежащая в дальнейшем кон кретизации.
Поскольку мера повреждения D является весьма медленно меняющейся функцией числа циклов п, то конечную разность Д£>„ в уравнении (6.45) можно с большой точностью заменить производной dD/dn. В результате приходим к «кинетическому уравнению»
~ = H D .S ). |
(6.46) |
Правая часть уравнения (6.46) может быть составлена бес численным множеством способов, один из которых заключается в следующем. Введем вначале меру DS=DS(n/N, S) для одно родного режима (т. е. при S=const). Здесь N(S) — число цик лов, соответствующее разрушению при этом режиме. Взяв про изводную по п и выразив ее через Ds,получим
dDs |
Ds (D . S) |
|
||
~ ~ d n |
N (S ) |
|
|
|
Если принять, что это соотношение |
остается |
справедливым н |
||
для любого неоднородного режима, |
то придем |
к формуле |
||
/ ( А 5 ) = |
Ds (D, S) |
(6.47) |
||
N(S) |
||||
|
|
|||
Допустим, что мера D s зависит только от отношения n/N и |
||||
не зависит явно от амплитуды |
напряжения S |
(рис. 97). Это |
предположение будем называть гипотезой об автомодельности
меры Ds. Действительно, в |
общем случае Ds =DS(л, S) — |
|
=DS (п, N) (N и S являются |
взаимнооднозначными |
функция |
ми). Сокращение числа независимых переменных до одного |
||
Ds = |
Ds ( i ) |
(6.48) |
аналогично сокращению числа независимых безразмерных па раметров при анализе автомодельных решений в механике.
Подставив выражения (6.47) и (6.48) в уравнение (6.46), мы получим полное разделение переменных
|
|
|
|
dD |
(6.49) |
|
|
|
|
(D) |
|
|
|
|
|
|
|
Здесь D — мера повреждений, |
|
накопленных после N циклов. |
|||
Пусть N„ — предельное |
(разрушающее) число циклов. Ему со |
||||
ответствует D —1, откуда |
|
|
|
|
|
N„ |
|
__ |
1 |
|
|
С _ d n _ |
Г |
dD |
|
||
J |
N(S) |
“ |
J |
DS (D) |
|
Правая часть не зависит от программы нагружения и всегда