книги / Расчет элементов конструкций заданной надежности при случайных воздействиях
..pdfчто совпадает с результатом (1.34), полученным ранее. В безразмерной форме уравнение (1.38) имеет вид
и
3. Нагрузка распределена по закону распределения наибольших зна чений (двойное зкспоненциальное распределение), несущая способ ность - по нормальному закону
/з (Я) = |
. 1..._J |
Ч —ОСЗ |
|
|
|
|
|
~ 7 Г |
|
|
|
|
1 |
r |
1 |
R ~м р |
а |
/ Н * ) = |
,___ |
ехР [- |
- |
( --------— ) |
]• |
|
oR s/ г п |
|
2 |
oR |
|
С учетом того, что S= Kq, для закона распределения S имеем
U (s ) = |
exp |
S - O i |
, |
s - |
|
S T |
' e x p [' |
~ |
|||
|
|
где Pi — Kp3 ; a j — K a 3 .
Используя (1.4), для надежности можно записать
J Л (Л)[ J / , (S)dS]dR = J Л (/?)F, (Л)Л?.
—ОО |
—оо |
_ оо |
Для рассматриваемого закона распределения напряжений
F i (R) = e x p [ - e x p ( - ^ - ^ i )].
Тогда |
|
|
|
|
|
|
|
|
оо |
|
1 |
|
1 |
Л - и » 2 |
|
, |
Л |
Н = J |
------------ |
ехр[ - — ( -------- — ) |
]ехр[ -е х р (---------- )]</Л. |
|||||
— оо |
aR \ГгИ |
|
2 |
aR |
|
|
PJ |
|
Обозначим: |
R |
- m R |
_ |
|
|
|
+ oRy. |
|
|
°R |
= у, тогда dR = aRdy,R = |
||||||
|
|
|
|
|
|
|
|
|
Отсюда для Я имеем |
|
|
|
|
|
|
||
H = |
1 |
~ |
( v J |
|
Op |
m. - a |
|
)] \d y = |
|
|
|
-exp[ - ( — |
— |
- |
|||
>/27 |
- |
L 2 |
. |
Pi |
Pi |
|
|
|
sflU |
J exp |
|
-exp[-(cy + i4)]|rfy, |
|
|
|||
|
|
|
|
|
|
|
||
гдес = |
|
|
^ = |
|
mR - K a 3 |
|
|
|
|
кРэ ’ |
mR ~ tt| |
|
|
|
|||
Pi |
|
|
Pi |
|
|
|
|
21
В безразмерной форме с учетом того, что/я. = а3 + 0,57703, оа = 1,28303,
с = |
1,2834^ л |
1,283 (я - 1 ) |
А = |
+ 0,577. |
Проведя аналогичные выкладки для различных сочетаний законов распределения нагрузки и несущей способности, когда не удается анали тическими методами взять интеграл в выражении для надежности, мож но получить подобные же выражения для определения К (эти результаты приведены в табл- 1-2).
По этим выражениям для различных сочетаний изменчивости нагруз ки Aq и несущей способности AR , определив значение интеграла для ряда значений п методами численного интегрирования, можно построить гра фики зависимости п = /(Я ), которыми удобно пользоваться при рас
четах.
Пример.
Цилиндрическая оболочка радиусом г = 1 м нагружена внутренним давлением q, величина которого случайна, с нормальным законом распределения с параметра ми rriq = 1,8 МПа, Oq = 0,036 МПа. Несущая способность материала оболочки слу чайна и распределена по закону Вейбулла с параметрами 0 = 3, Л0 = 670 МПа, о = 2263 МПа3.
Подобрать толщину оболочки И таким образом, чтобы надежность ее Язад == = 0,9985. Учет случайного разброса толщины оболочки провести с доверительной вероятностью Я/, = 0,9986, т.е. Язаа1Ц, = 0,9999.
Воспользуемся данными табл. 1.2. Для рассматриваемых параметров законов распределения имеем Г = 0,8934, Г * = 0,324,
Од |
AR |
OR |
Г* |
Aq = —L. = о,02, |
= — - = |
------------------ = 0,084. |
|
mq |
|
mR |
R„ + %Ja Г |
Определив значение интеграла для ряда значений п методами численного интегриро вания и построив график зависимости л = /( Я ) , найдем для надежности Я = 0,9999 значение п = 1,3. Тогда для К имеем
К = m R |
= R 0 + $fa Г |
— 670 + 226 • 0,8934 . = 374 |
mqn |
1,8 ■1,3 |
1,8 ■1,3 |
Для рассматриваемой оболочки К = r/h, отсюда/) = r/К = 2,67 ■10"3 м. В предпо ложении нормального закона распределения значения толщины оболочки с коэф фициентом вариации 4^ = 0,033 и доверительной вероятности Щ =< 0,9986 (для которой 7 = 3 ) по формуле (1.12) для номинальной толщины можно получить
*ном = “j----- ~т —2,97 ■10-3 м. 1 - у Ah
Таким образом, для искомой толщины имеем (2,97 ± 0,3) • 10"3 м.
4. Нагрузка и несущая способность распределены по логарифмичес ки нормальному закону
В этом случае |
tlru? - 1п<7о) |
|
/з(<7)= ------ 1 __ ■ехр[- |
]; |
|
o Z3q ТТтг |
2а 2з |
|
22
/И * ) = |
ехр1- . < |
« - « . > |
‘ |. |
у/гп |
2а? |
|
|
о 22R |
|
||
Для этого случая [31] гауссовский уровень надежности равен |
|||
’J L |
J I 1 + л ; |
|
|
К т , |
1 + >4 I |
|
(1.39) |
у = • |
|
|
|
s/inta +^ й Г 7 л | ) 7 |
|
|
|
Разрешив это уравнение относительно К, получим |
|||
К = |
mR J T T A ? |
(1.40) |
|
|
|
||
y / l + A j i ехр^7 v/ln[(l + / ф |
(1 + А \ T jJ |
5. Нагрузка распределена по нормальному закону, а несущая способ ность - по закону Релея
/з (?) = |
г (<? " т |
<?)2 |
1 |
•ехр[---------- |
— |
]; |
|
>/27Г ос |
20л |
|
|
/ а ( Л) = т - е х р [ - ^ - ]• |
|
|
|
Д2 |
2^2 |
|
|
Для этого случая надежность определяют по следующему выражению [31]:
Н = |
1 |
г |
1 |
]■ |
(1-41) |
1ТАд ШдК2 |
ехр[- |
( 2> 4 |
|||
|
|
+ Aq) • 2 |
|
||
|
J T . |
|
|
||
|
2m2R |
|
пт 1к2 |
|
|
После несложных преобразований это выражение можно представить в виде, удобном для графического решения,
H2 (2mR + itntqAq К2) = 2тR exp( - 1 ■)• (1.42)
2"R
+ A i
m iqK 2
Графическое решение данного уравнения позволяет получить искомое значение К. При малых изменчивостях нагрузки Ад для определения К можно использовать приближенную формулу
2тR |
In// |
К = |
(1.43) |
Зная К, легко найти размеры поперечного сечения.
6. Нагрузка распределена по закону Релея, несущая способность - по нормальному закону
23
h ( a ) = - V |
ехР[_ - г т |
1» |
|
|
|
|
a\ |
2al |
(Л - m R )2 |
|
|
fz(R) = |
|
exp[- |
|
|
|
|
]• |
|
|
||
|
■jraoR |
2 ok |
|
|
|
В этом случае [31] |
|
|
|
||
# = 1 - |
|
t Г |
exp[- |
]. |
(1 .4 4 ) |
|
|
Wn%A |
* K 2n?q |
|
|
|
|
R A R |
+ |
2A i |
|
|
|
2K 2m2„ |
|
||
|
|
m 2R |
|
|
После несложных преобразований это выражение можно представить в виде, удобном для графического решения,
( 1 - Я ) 2 (1 + ™ R A R ) = ехр[ • |
2K 2n?q |
■ |
(1 .45) |
2К 2т2„ |
|
+ A i
1Гт%
в результате которого получают искомое значение К. При малой измен чивости несущей способности AR для определения К можно применять приближенную зависимость
2т, |
1п(1 - Н ) |
(1 .4 6 ) |
|
Пример.
Прямоугольная пластина длиной а = 2 м и шириной b = 1 м, шарнирно опер тая по трем сторонам и защемленная по четвертой стороне, нагружена распреде ленной нагрузкой, меняющейся по треугольному закону, величина qa которой слу чайна с релеевским законом распределения с параметром а3 = 0,06 МПа (рнс. 4 ).
При в, = 0,6 МПа имеем mq = 1,253 • 0,06 = 0,075 МПа. Несущая способность материала пластины R случайна и распределена по нормальному закону, имеющему
параметры mR - |
400 МПа; oR = 40 МПа. |
|
||
Необходимо найти толщину пластины А, при которой ее надежность по проч |
||||
ности Н = 0,99. |
|
|
|
|
По (1.45) имеем |
|
|
|
|
|
я - 16 ■Ю4 • 10"1 |
= ехр " |
|
|
(0,99 - 1 ) 1 (1 + |
2 • К* (0,075)* |
|
||
|
|
Д*1 (0.075)" |
+10_, ’ |
|
|
|
|
v 16104 |
- |
64,68 • 104 |
|
|
||
или 10-4 (1+ |
|
) = ехр- |
|
|
^--т-тТТТЛл |
|
+ 0,01 |
|
|
V. |
64,68 • 10* |
|
||
|
|
|
Рис. 4. Схема нагружения прямоугольной пластины рас пределенной нагрузкой, меняющейся по треугольному
азакону
24
Графическое решение полученного уравнения дает К -- 2S43. Для рассматриваемой пластины имеем [28 ]
К = 6 ■0,084а2 / А1, отсюда А = У —-2’^ Z E = 2,8 ■10'* м.
7. Нагрузка распределена по экспоненциальному закону, а несущая способность подчиняется гамма-распределению
/ з ( < 7 ) = Х з е - ^ ;
f*(R) = |
R ‘ |
|
а + 1 |
||
|
Так как S = Kq, то для напряжения S по (1.29) в качестве закона распре деления имеем
Л Ф ^ е - ^ . г д е Х , = \ 3/К.
Тогда используя (1.4) и проведя преобразования, аналогичные сделан ным ранее, можно получить [21 ]
|
|
1 |
|
н —\ - |
(_ЛГ_г_ ) “ + 1 |
(1.47) |
|
Тогда |
|
|
|
|
1 |
|
|
\3 |
0 2 |
= |
|
|
|
||
+ |
1 |
|
|
К |
|
0г |
|
Отсюда для К имеем |
|
||
г _ |
Хз02 а *-УГПГ |
(1.48) |
|
|
|
|
1 - а+У~Гчн ■
Пример.
Стержень растянут силой Р, величина которой случайна и распределена по экспо ненциальному закону, имеющему параметр распределения Л, == 10~4 1/Н. Несущая способность материала стержня также случайна, но подчиняется гнмма-распределе-
нию с параметрами а = |
1 и 0, = |
100 МПа. |
|
|
Найти площадь поперечного сечения стержня, обеспечивающую надежность |
||||
Я = 0,99 . |
|
|
|
|
По (1.48) имеем |
|
|
|
|
М»а+ |
= ю~4 • ю1 • у/Щ • ю6 |
ю4 |
1 |
|
1 - а * \ т г |
1 - ч / а о г |
9 |
м1 " |
|
Так как К = IjF, то F = |
1 .9 |
= 9 • 10-4 м*. |
|
|
- г - |
|
|
25
8. Нагрузка подчиняется гамма-распределению, а несущая способ ность - экспоненциальному закону
h |
{Я) = |
— ?-■ |
, |
; Л (R) = Х2 e~x*R . |
||
|
|
|
e!U?+ l |
|
||
Проведя преобразования, аналогичные предыдущим, получим [21] |
||||||
н = |
( — |
i— _ |
а + 1 |
(1.49) |
||
|
|
|
1 |
+ Pi х2 )и |
||
где |
= 03К. |
|
|
|
||
Тогш |
- V |
7 T , т ^ |
т, |
|||
отсюда имеем |
|
|
||||
r _ |
1 - а +У1Г |
(1.50) |
||||
|
|
0з |
|
Х2 “ +1/н~ |
||
|
|
|
|
Пример.
Прямоугольная пластина, у которой b <а, имеет две шарнирно опертые сторо
ны, одну защемленную и одну свободную (рис. 5 ). Посредине свободной стороны приложена сосредоточенная сила Р, величина которой случайна и распределена по гамма-распределению с параметрами а = 3; 05 = 5000 Н. Несущая способность ма
териала пластинки также случайна с экспоненциальным законом распределения,
параметр которого К, |
1 |
1 |
= ----- |
--------- . |
г^ 400 МПа
Необходимо определить толщину пластины А, при которой ее надежность по прочности Н = 0,99.
В соответствии с (1.50) |
|
|
|
|
|
|||
Л — |
1 + “ +{/ЗГ |
|
( 1 - ^ 9 9 * ) |
- 4 0 0 - 10‘ |
1 |
• |
||
|
|
— > ■ 1 |
■ |
■' —— |
—1о0 ■ |
|||
|
|
|
|
5000 |
|
|
м |
|
Для рассматриваемой пластины имеем [29] |
|
|
|
|||||
К = 3,06/Л1. |
|
|
|
|
|
|
|
|
Отсюда |
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А = |
|
|
3,06 |
= 4,36 • 10- 1 |
м. |
|
|
|
|
|
160 |
|
|
||||
|
|
|
|
|
|
|
|
|
9. Нагрузка |
подчиняется нормальному, |
а несущая |
способность - |
|||||
|
Нагрузка |
усеченному |
нормальному закону распреде |
|||||
|
|
|
|
ления |
|
|
|
|
...f ... ...Т " |
Рис. S. Схема нагружения прямоугольной пластины |
|||||||
сосредоточенной силой |
|
26
В этом случае/з (<7) и / 2 (Ю описываются выражениями вида (П.29)
и (П.34) соответственно. Тогда, |
воспользовавшись |
уравнением |
(1.4) |
|||
и проведя соответствующие преобразования, можно записать |
|
|||||
Н = С<Ф( |
mR ~Ктд |
) - |
- |
/ Л (*)[ { U (S)dS]dR |
(151) |
|
|
ч/аД + Ко; |
|
с |
_ 00 |
^ |
|
гдеС=[Ф(Г2) |
-Ф(г,)]- 1 ; |
г, = (Л, - m R)/aR -, |
|
|
||
t2 = (Л2 - m R)laR \ S = Kq; |
|
|
|
|
||
|
|
(S-Kmq)2 |
|
|
||
Л Р ) = |
■ехр[ ■ |
|
]; |
|
|
|
spTn КО, |
|
|
|
|
|
Л 1 и Л2 —точки усечения слева и справа соответственно. Значение интег рала, входящего в выражение надежности, можно рассчитать лишь чис ленными методами с помощью ЭВМ, пользоваться этим выражением для нахождения искомого К крайне неудобно. Поэтому для высоконадеж ных систем, когда точка усечения слева R\ достаточно близка к mR (рис. 6), в качестве нижней оценки для надежности можно записать
R « |
|
R I -т о |
|
н = J |
/,(3)< ю = ф ( — - — - ) . |
|
|
-о» |
|
as |
|
R t -K n iq |
= 7 |
(1.52) |
|
Отсюда |
К0„ |
||
|
|
|
|
и К = |
У°д |
|
(1-53) |
'? + |
|
|
' Пример.
Найти толщину стенки А трубопровода диаметром d = 5 см, обеспечивающую
надежность Н = 0,999. Трубопровод выполнен из стали, несущая способность кото рой случайна, и нагружен внутренним избыточным давлением q, величина которого случайна с нормальным законом распределения с параметрами та = 1 0 МПа; о„ =
= 1 МПа. |
|
4 |
|
Закон распределения несущей способности — усеченный нормальный с пара |
|||
метрами mR = 200 МПа; |
aR — 20 МПа, параметр усечения слева R i = 180 МПа. |
||
Согласно выражению (1.S3) можно записать |
|||
К |
R. |
|
180 •= 13,7, |
= ------- |
|
||
|
mq + yoq |
10+3,1 • 1 |
|
где 7 = |
3,1 для Н = 0,999. |
|
Как известно из теории безмоментных оболочек, К = djTh. Отсюда А = —~г =
5 io-J
= 0,182 • 10-1 м.
27,4
Для сравнения найдем толщину стенки трубопровода без учета усечения за кона распределения несущей способности.
27
fWifrsj
.f,(5) |
_ f 2W |
Коэффициент К определим по уравнению (1.18), подставив в него следующие значевния:
“ = 90,39 МПа1; д = 3,62 • 104 МПа1; f = 4 • 103 МПа1,
к = 4 ' J° 3 “ / 1 6 • 10* - 4 • 90,39 • 3,62 • 104
2• 90,39
адля толщины стенки получим
Таким образом, учет усечения распределения несущей способности слева (а это означает отбраковку трубопроводов, у которых несущая способность меньше Л ,) приводит к уменьшению толщины стенки на 9 % при той же надежности.
10. Нагрузка и несущая способность подчиняются гамма-распреде лению
Проведя соответствующие преобразования,из (1.4) можно получить[31]
1 - # = / 5[(а2 + 1) ( а 3 + 1)] , |
(1-54) |
где /§ [(а 2 + 1) (аз + 1) ] —функция, определяемая по специальным таб лицам.
Она представляет собой отношение функций Эйлера первого и вто рого рода
«/fi[(a2 + 1) (аз + О ] - |
В&[(а2 + 1) (а3 + Q ] |
|
(1.55) |
|
В[(а2 + 1) (а3 + 1) ] |
' |
|||
|
|
|||
где В[(а2 + 1) (а3 + 1)] - |
г (a2 + 1) г (а3 + 1) |
|
(1.56) |
|
г (а2 + а3 + 2) ’ |
|
|||
|
|
|
8
(1.57)
(1-58)
28
Рис. 7. Схема нагружения круглой пласти ны равномерно распределенной нагрузкой
Схема решения выглядит следующим образом. Из условия /{ [(а 2 + 1) X X (аа + 1) ] = 1 - Н по специальным таблицам [6], зная а 2 и а 3, ищут па
раметр 8, |
обеспечивающий заданную величину 1 -Н . Тогда из (1.S8) |
по найденному 8 легко найти искомое значение К : |
|
К = |
fos |
& а - в ) '
Пример.
Круглая пластина диаметром 1 м нагружена равномерно распределенной наг рузкой q, величина которой случайна и подчиняется гамма-распределению с пара метрами а3 = 10 и 0з = 0,2 МПа (рис. 7 ). Несущая способность материала пласти
ны также случайна и имеет законом распределения гамма-распределение с пара
метрами а , = 9; ft |
= 20 МПа. |
|
|
||
Найти толщину пластины Л, обеспечивающую надежность по прочности 0,997. |
|||||
Из условия У6 (10,11) = |
1 - Н = 0,003 |
по таблице [6, с. 291 находим в, при |
|||
котором вероятность отказа равна 0,003: 6 = 0,2. Тогда |
|||||
К = |
g»6 |
= |
20 ° ’2 |
= 25 |
|
|
ft (1- 6) |
|
0,2 -0,8 |
|
|
Для рассматриваемой пластины |2 ) |
|
||||
|
0,1263ft1 |
|
|
|
|
отсюда |
|
|
|
|
|
А = |
/0.1263Л 1 |
/ 0,1263 0.51 |
„ „ |
||
V ----- - ------ |
= > /--------—---------- = |
3,55 ■10 -1 м. |
Очень часто бывает неудобно использовать в расчетной практике специальные таблицы табулированных функций. Во-первых, они не очень широко распространены и не всегда доступны. Во-вторых, в тех случаях, когда расчет надежности является частью расчета конструкции на ЭВМ, использование табулированных функций затруднено. Поэтому можно ис пользовать приемлемо точную аппроксимацию параметра 8 функции Уб[(а2+ 1) (а3 + 1) ] и гауссовского уровня надежности у*
0 = |
|
а2 + 1 |
|
(1.59) |
|
а2 + l + (а3 + 1) е 2w ’ |
|||
|
|
2а 2VrWf ЗА } ' |
||
|
у |
А |
" ( 2а3 + 1 |
|
где w = |
|
-— - |
|
- L ) - |
|
!( |
1 |
1 |
X = У2 - 3 |
|
2а2 + 1 |
2а3 1 ) |
* Справочник по специальным функциям. М.: Наука, 1979.832 с.
29
Рис. 8. Зависимость относительных размеров поперечного сечения от надежности по проч ности для различных комбинаций законов распределения нагрузки и несущей способ ности:
1 - экспоненциальный - экспоненциальный;
2 — распределение Релея - распределение Релея; 3 - логарифмически нормальный - лога рифмический нормальный; 4 — экспоненци альный - гамма-распределение; S - распре деление Вейбулла - распределение Вейбулла; 6 - гамма-распределение - экспоненциальный
оо + a ,t |
Г |
|
?=г~ттмгг,?-;г= ^1п[ (1 ~Ну ■]; |
||
а0 = 2,30753; |
А, = 0,99229; |
|
в, =0,2706; |
Ь2 = 0,04481. |
|
В безразмерной форме с учетом того, что п = mRIKmq,
mR ~ 02 (<*2 + 1) , mq = 0з (а3 + 1), для п можно получить: п = e2w.
Расчет по этим зависимостям будет выглядеть следующим образом. По заданному Н находят t, затем у, после этого, найдя X и А, подсчиты вают w. Зная vv, определяет п и К = mR/nmq . Имея значение К, легко найти размеры поперечного сечения. В частности для вышеизложенного примера для Н = 0,997 t = 3,40856, у т 2,74979, X = 7,602, А = 19,95, w = 0,6454, п = 3,636, тогда
20 |
10 |
К = |
25,0025. |
0.2 • 11 ■3,636
Для удобства использования результаты, полученные в данном разделе и полезные в практических приложениях, сведены в табл- 1.2.
На рис. 8 показаны графики зависимости относительных размеров поперечного сечения А/А* от надежности по прочности для различных комбинаций законов распределения нагрузки и несущей способности. Здесь А* — размеры поперечного сечения, подсчитанные при значениях нагрузки и несущей способности, равных их математическим ожида ниям. Для наглядности по оси абсцисс откладывается величина -lg (1-И).
30