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

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

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

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

/ зам (Я) = 0,7/н (3; 1,22) + 0,25/„(6; 1,22) + 0,05/„(9; 1,22);

/зам(Л) = 0,1/„(130; 12,2) + 0,3/„(160; 12,2) + 0,6/„(190; 12,2).

Для максимальных напряжений S с учетом того, что S = Kq, имеем

/з а м (« = 0,7/„(3/Г; 1.22ЛГ) +0,25/„(6/Г; 1.22А) + 0,05/„(9А ; 1,22А ).

Составим разность Z = R -S, для которой можно записать:

/зам<г ) = 0,07/н(1 3 0 -ЗАГ; V 150+ 1 .5A 1 ) + 0,025/„(1 3 0 -

-6 А ;о 2) +0,005/н (130 -9А ;ог) + 0,21/н (160 - ЗА; о г) + 0,075/н (160 -

“ 6А; о г) + 0,0015/„(160 -9А ; ог ) + 0,42/„(190 -ЗА ; о J + 0,15/1,(190 -

-6 А ;о г) + 0 ,0 3 /„ (1 9 0 -9 А ;%) ,г д е о г = V 150+1,5 А1 .

В этом случае

 

 

 

 

 

 

ее

 

 

 

 

 

 

Н = £ /зам С^) dZ,

 

 

 

 

или

 

 

 

 

 

 

Я = 0 ,0 7 Ф (” ^

)

+0, 0 2 5 Ф ( - ^ 1 )

+ 0,005 Ф( 13- ° - ~ ^ ) +

 

О2

 

 

О2

 

О2

* 0 ,2 !^ “

1 “

)

,0 .0 ,ж ‘“

.Г6*

)

.

 

О2

 

 

Oi

 

0 2

+ 0, 4 2 Ф

( ^ И

)

+ 0,15 Ф(—90

)

.О.ОЗФ!.1-90- ^ - ) ,

 

°г

 

°z

 

Чг

где a z = V 150+ 1,5AJ .

Решив приближенными методами это уравнение относительно А, получим А == = 9,75. Для рассматриваемого трубопровода, как это известно из теории безмоментных оболочек, А = d\1h

d

5 • lO '1

= 0,256 ■10"2 м.

Отсюда А = —

= 2

1.8. ДИСКРЕТНЫЙ СПОСОБ ОПРЕДЕЛЕНИЯ РАЗМЕРОВ

ПОПЕРЕЧНОГО СЕЧЕНИЯ ЗАДАННОЙ НАДЕЖНОСТИ

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

51

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

для нагрузки

Ji

Я»» Я1

Я, ; Яг

P i - 1; Pi

Pm - 1; Pm

Pi

я,

Яг

Pi

Pm

Pi

р 1

Р 2

Pi

Pm

для несущей способности

Ji

*0 >^1

/?,; R,

Ri - U Ri

*n - i; R„

*i

 

*i

 

Si

gi

Вг

gi

Sn

Здесь У,- -

обозначение /-го разряда;

яд или (Л,- - ц /?,■) - грани­

цы /-го разряда;

Р^ивал £/ —частота попадания случайной величины в

t-й разряд;

q или /?,• - "представитель”

i-го разряда; т или п - число

разрядов.

Таблицы строят следующим образом. Всю область изменения случай­ ной величины разбивают на разряды в порядке возрастания и заменяют совокупность значений случайной величины внутри разряда "представи­ телем” разряда, с которым производят все дальнейшие операции. В ка­ честве представителя разряда можно брать средневзвешенное значение случайной величины внутри разряда или среднее значение разряда [9]. Для удобства и в запас надежности в качестве представителя разряда бу­ дем брать: для нагрузки - верхнюю границу разряда, а для несущей спо­ собности - нижнюю границу. Учитывая известную зависимость S = Kq, для закона распределения напряжений можно получить следующую таб­ лицу:

Ji

Kpoi KPi

*<?, ; КЯг

 

KPi- i: RPi

 

Kpm- \\KPm

Si

Kq,

Kq,

 

Kqi

 

KPm

PI

Pi

Pi

. . .

Pi

. . .

Pm

Составим разность Z = R - S = R - K q и запишем ее закон распределения в табличной форме

52

R , Kqm

Jos"

ibe

1

 

 

 

!i

 

' Л

g,Pm

 

KiPm- 1

R, ~Kqm

 

R» -Kq,

g,Pm

 

 

 

g3P,

Л, -*<7,

R 2 Kqm

Rt -Kq,

g,Pi

КгРт

g*Pi

 

Rn-Rqm

R n - K q ,

 

SnPm

gnP1

Из полученного закона распределения R - S можно найти искомое К, обеспечивающее заданную надежность, следующим образом. Ищем раз­ ряд, до которого сумма вероятностей попаданий в предыдущие разряды равна (1 - И). Тогда, приравняв нулю представителя этого разряда, имеем

R i - K q , = 0.

(1.113)

Это означает, что представители всех остальных разрядов больше нуля, т.е. R - Kq > 0, и вероятность этого равна Н.

Отсюда

K = Rilqh

(1.114)

где Rj и <7/ - представители разрядов, попавших в граничный разряд.

Пример.

Квадратная пластина со стороной 4 = 1м , два края которой шарнирно оперты, а два защемлены, нагружена равномерной нагрузкой q, величина которой случайна и имеет следующий закон распределения:

Jl, Па

0 ...

1

1

...2

2 ...

3

3

...4

4 .

. 5

5

...6

6 ...

7

7 ...

8

8 ...

9

9 ...

10

qu Па

 

1

 

2

 

3

 

4

 

5

 

6

 

7

8

9

10

 

 

Pi

0,04

0,18

0,36

0,22

0,1

0,04

0,02

0,02

0,01

0,01

 

Несущая способность материала пластины также случайна

 

 

 

 

 

■'1-Ю-3, Па 1

...2

2

... 3

3

... 4

4

...5

5

. . 6

6

...7

7

... 8

8

...9

9...

10

10...

11

 

Па

1000

2000

. 3000

4000

5000

6000

7000

8000

9000

10000

_

gi

0,01

0,01

0,02

0,02

0,06

 

0,1

 

0,2

0,36

0,18

0,04

 

Найти толи1ину пластины И, при которой ее надежность по прочности Н = 0,999.

 

Запишем з<1кон распределения напряжени AS’

 

 

 

 

 

 

 

 

JiK". Па 0 ...

1

1

... 2

2

... 3

3

...4

4

...5

5

...6

6

...7

7

...8

8

...9

9 ...

10

5/, Па

 

К

 

2К

 

ък

 

4К

 

5К

 

6К

 

 

8К

 

10*

__

Pi

0,04

0,18

0,36

0,22

 

0,1

0,04

0,02

0,02

0,01

0,01

53

Будем-искать закон распределения разности R - S в виде аналогичной таблицы, располагая значения R - S в порядке возрастания,

Zjt Па 1000-10*

1000-9*

1000-8*

1000 -7*

1000 -6*

 

0,0001

0,0001

0,0002

0,0002

0,0004

. . . .

Для того что&1 значения R - S были расположены в_порядке возрастания, надо

на _каяодом шаге Л,

~_Я^ сравнивать величину Л, - Kqm путем подстановки К =

= RA lqi с величиной Л, “ *<П- 1 - Пройдя все значения Л, - Kqj, переходят к Л, -

- Kq( и т.д. Для рассматриваемого случая после 1-го шага К = 1000/10 = 100; 1000 -9* = 100, а 2000—10* = 1000, поэтому переходим ко второму шагу. После второго шага К = 1000/9 = 111,1 опять 1000—8* < 2000 -10* и т,д.

Как только сумма частот попадания в разряды достигнет величины 1 - Н, процесс можно прекратить и для последнего разряда приравнять/? - 5 = 0 [для раз­

ряда (1000-6/Q

сумма частот равна 0,0001 + 0,0001 н 0,0002 + 0,0002 + 0,0004 =

= 0,001 = 1 -Н\. Поэтому 1000 - 6 * = 0

 

К = 1000/6 = 166,7.

 

 

Для рассматриваемой пластины [29]

 

0,0697л2 -6

 

 

Отсюда следует

 

 

 

/ 0 , 0697' 1 -6

/ 0 ,0697 1 -6

,

" v —

*-----------v

 

10-1 м.

— 1ббл--------- 5

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

Схема решения здесь выглядит следующим образом. Пусть задан за­ кон распределения нагрузки

Ji

 

qi \ ч*

Qi- i; <?i

 

qm - 1; Ят

4i

4\

Яг

4i

. . .

Ят

Pi

Pi

Рг

Pi

 

Pm

Тогда закон распределения деформаций w = K*q будет иметь вид

wi

K*qt

К*Яг

**<?,'

. . .

К*Ят

Pi

P,

Рг

Pi

. . .

Pm

Ищем к -й разряд, для которого выполняется условие

 

P i

+ Рг + — + РЛ:-1 ~ н -

 

 

(1 .115)

54

Рис. 18. Схема нагружения прямоугольной пластины распределен-

Ч

ной нагрузкой

0

Тогда можно записать

 

К*Як = "’зад!

10

что означает, что P(w < w ^ ) = Я.

______

Отсюда для К* запишем

|^ а .

К * = . 1^5_.

(1.116)

Чк

 

Найдя К*, по табл. 1.1 легко определить искомые размеры поперечного сечения.

Пример.

Шарнирно опертая прямоугольная пластина (а = 1 м; Ъ= 2 м; д = 0,3) нагру­ жена гидростатическим давлением (рис. 18); q0 - случайная величина со следую­

щим законом распределения:

/,-ЮМПа

0 ... 1

1 ...

2 2 ... 3

3 ...4

4 ... 5

5 ...6

6 ...7

7 ...8

T IP 9 ... 10

<7/10МПа

1

2

3

4

5

6

7

8

--------- 1.......

"

9

10

Pi

0,05

0,1

0,35

0,25

0,14

0,07

0,03

0;005

0,004

0,001

Найти толщину пластины А, при которой ее надежность по жесткости Н = 0,99; и'зад == 0,3 - 10-* м. Закон распределения величины максимального прогиба с уче­ том того, что wmax = K*q, будет иметь следующий вид:

* 1

wmax

Pi

0 ...

1

1

...2

2 ...

3

3...

4

4 ...

5

5 ...

6

6 ...

7

7 ...

8

8 ...

9

9 ...

10

0,1**

0,2К*

0,3**

0.4**

0,5** 0,6** 0,7** 0,8** 0,9**

1,0**

0,05

 

0,1

0,35

0,25

0,14

6,07

0,03

0,005

0,004

0,001

А- 1

Найдем к-А разряд, для которого

р,- = 0,99. Внашем случае это 8-й разряд,

так как для него

 

 

, = 1

 

0,05 + 0,1 + 0,35 + 0,25 + 0,14 + 0,07 + 0,03 =

0,99.

Следовательно,

 

 

 

 

** =

‘'зад

0,3 • ю - 1

 

 

 

= 3,75 • 10-3 1/МПа.

 

0,8

 

0,8

 

 

Для рассматриваемой пластины [29]

 

 

К* =

0,00206л*

 

0,00206а* • 12(1 - ц 1)

 

D

 

Eh*

 

 

 

 

 

,/0 ,0 0 2 0 6 -

1 • 12 - 0,91

 

io-j м.

Отсюда А =

У 2 ■Ш5 - 3,75 ■10~

" 3>И

55

Рис. 19. Схема нагружения треугольной пластины равно­ мерной распределенной нагрузкой

Также легко применим дискретный метод для расчета иЗ условий надежности по устойчивости.

В этом случае из закона распределения нагрузки

Ji

4»; <7.

Ч\ ; Чг

 

4 /- i; 4i

 

Чт- i ; Чт

—»

_

 

 

_

ж

Ч1

ЧJ

. . .

4i

 

Чт

Pi

р.

Рг

. . .

PI

. . .

Pm

следует определить к-й разряд, для которого выполняется условие

Pi

+ Рг + -

+ P k - l =Н.

 

 

 

(1.117)

Приравняв <7кр = qk, что означает,что вероятность непревышения нагруз­ кой q величины <7кр равна Я, легко найти необходимые размеры попереч­ ного сечения.

Пример.

Треугольная пластина, у которой размер а 2 м, сжата со всех сторон распре­ деленной нагрузкой q, величина которой случайна и имеет следующий закон распре­ деления (рис. 19):

Jt-10-s НД*

0 ... 2 2 —4

4...6

6 - 8

8 ... 10 10... 12

1 2 -1 4

14...16

1 6 -1 8

1 8 -2 0

A,"10-S H/M

2

4

6

8

10

12

14

16

18

20

Pi

0,05

0,1

0,35

0,25

0,14

0,07

0,03

0,005

0,004

0,001

Найти толщину пластины, при которой ее надежность по устойчивости Н =

= 0,99.

 

 

 

 

 

 

 

 

 

 

Ищем А-й разряд, дЛя которого выполняется условие

 

 

 

Р, + Рг + ~ + Р к - 1 = #•

 

 

 

 

 

 

 

В нашем случае это 8-й разряд, так как для него

 

 

 

 

0,05 + 0,1 + 0,35 + 0,25 + 0,14 + 0,07 + 0,03 = 0,99.

 

 

 

Тогда flKp =

16 • 10s Н/м.

 

 

 

 

 

 

 

Для рассматриваемой пластины [29]

 

 

 

 

 

 

 

An1 Eh*

 

 

 

 

 

 

 

 

ffKp_

12(1-д г)аг '

 

 

 

 

 

 

 

О т с ю д а __________________________________________

 

 

 

 

12(1 -м а)о*Акр

 

/1 2 -0 ,9 1 - 4 -1 6 -1 0 *

 

 

 

 

 

4 п*Е

 

~

4 • 3,142 • 2 • 10"

=

2’®7

м‘

56

ГЛАВА 2

ДИНАМИЧЕСКИЕ ЗАДАЧИ

Сначала поясним, что мы будем понимать под термином ’’динами­ ческие задачи”, так как обычно этим термином обозначают задачипроектирования и расчета конструкций с учетом сил инерции. Но как мы ви­ дели, ряд задач, в которых учитываются силы инерции, с успехом могут решаться квазистатическими методами и могут быть отнесены к квазистатическим. Поэтому в данной работе под термином "динамические задачи” мы будем понимать задачи, для решения которых необходим аппарат теории случайных функций.

2.1. РАСЧЕТЫ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ ЗАДАННОЙ НАДЕЖНОСТИ ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ, ОПИСЫВАЕМЫХ СЛУЧАЙНОЙ ФУНКЦИЕЙ

Найдем

вероятность того, что в течение данного времени будет

не более

заданного числа выбросов. Особый интерес представляет-

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

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

Под мерой надежности Н будем понимать вероятность того, что ни разу за срок службы Т максимальное напряжение S не превысит несущей способности/?.

Пользуясь некоторыми результатами теории выбросов случайной функции за случайный уровень [34, 41], для среднего числа выбросов случайной функции 5 (0 за случайный уровень R в течение срока службы

Т можно записать

 

V(R/Т ) = ТЦ Sf(R, Б/ t)dSdt,

(2.1)

где V(RjT) — среднее число превышений уровня/?; S(t) —производная функции S (г) ; /(/?, S/t) совместный двухмерный закон распределения S и ее производной в момент времени t при S = R.

Найдем вероятность того, что в течение данного времени будет не более заданного числа выбросов. Особый интерес представляет частный случай, когда появление последовательных выбросов можно считать независимыми ’’редкими” событиями. При этом принимаем, что число выбросов в течение времени подчиняется закону Пуассона.

57

Так как вероятность Рт появления события т раз в этом случае за­ висит только от V(RjT), то вычисление Рт может быть доведено до конца.

Например, для вероятности Р0 того, что за время Т не произойдет ни одного выброса (а это и есть принятая надежность), получим [34]

Н = Р0 = ехр[ - Ц Sf(R, S/t)dSdt].

(2.2)

Подставим в формулу (2.2) известные вероятностные характеристики нагрузки <7(0 и несущей способности R. Принимая во внимание, что должно выполняться равенство Н -- Н^ц, для определения размера по­ перечного сечения А получим выражение

й = |/г(д1( Лг,..., ли, Г,/^ззд),

(2.3)

■где 0\ , «2.......о„ — известные заранее

вероятностные характеристики

нагрузки <7 и несущей способности R.

 

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

* = Л, = е х р [ - ^ wf(w3m, w /t)dwdt].

(2.4)

Подставив в уравнение (2.4) известные вероятностные характеристики нагрузки, получим выражение для определения искомых размеров по­ перечного сечения h, обеспечивающих заданную надежность:

А = Ф*(Р1, ®2, •••» ®и> w3au> ^ з ад)»

(2-5)

где в |,

..., ап —известные заранее вероятностные характеристики наг­

рузки <7 (0 • Аналогично решается задача проектирования элементов конст­

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

Таким образом, надежность по устойчивости будет

 

Н = Р0 = ехр [ - ^ < 7/(<7кр. <7/ 0<*7<*]-

(2.6)

Решая это уравнение с учетом того, что должно выполняться условие Н ~ Я 3ад, получаем необходимый уровень <7кр, зная который легко найти размеры поперечного сечения, которые и обеспечат заданную надежность.

S8

2.2.РАСЧЕТ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ

ЗАДАННОЙ НАДЕЖНОСТИ ПРИ СЛУЧАЙНЫХ ВОЗДЕЙСТВИЯХ БЕЗ УЧЕТА СИЛ ИНЕРЦИИ

Рассмотрим, как и в гл. 1, элементы конструкций,

максимальные нап­

ряжения S в которых линейно зависят от нагрузки q

 

S = Kq,

(2.7)

где коэффициент К связан простыми зависимостями с размерами попе­ речных сечений элемента конструкции (см. табл. 1.1).

Предположим, что характер действия нагрузки q (t), которая пред­ ставляет собой случайную функцию, таков, что силами инерции при опре­ делении напряжений можно пренебречь, т.е. напряжения определяются по (2.7). Подставив в (2.2) уравнение (2.7), получим выражение для определения К. Зная К, легко найти размеры поперечного сечения.

Если S (Г)

- нормальный стационарный процесс, для

V(R/T) имеем

 

 

То*

 

г, -tn of

 

(2.8)

V(R/T)= ------ ехр[ -

 

-----£ _ ] .

 

 

 

2т г ч / о | + о £

2 (о |+ а£)

 

 

Отсюда для надежности получим [34]

 

 

 

Н = exp [-

------- T°s

—■ ехр[ -

im R ~ms)

])

(2.9)

 

[

2л >/ с$ + а \

2 (ct2s +o2R)

J

 

где о | -

дисперсия напряжений; о% - дисперсия несущей способности

материала

конструкции;

ms — математическое ожидание напряжений;

тд - математическое ожидание несущей способности материала конст­ рукции; Т - срок службы.

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

ция нагрузки может быть аппроксимирована формулой [27]

 

Кя (т) = Оде “ 1т1(со$0т+ - j - sin(31т | ).

(2.10)

Константы а и 0 подбираются так, чтобы экспериментальная кривая Kq (г) совпадала с теоретической кривой, построенной по формуле (2.10).

Для этого случая имеем

= *S (0) = - *S ( 0 1 г = 0 = °s («* + Р У= к2 о2ч (ОС2 + р2). (2.1D

(В случае других вариантов корреляционных функций нагрузки расчет Можно вести по табл- П.4.)

Запишем выражение для надежности

Н = ехр

f

Т ч/ а 2 +

рг Koq

( m p ~ K m q )2 1

( 2. 12)

------

,

■ехр[ - -------=------------]>,

 

2 n - J K 2 Oq + Од

2 ( K 2 Oq + 0 2д ) )

 

из которого можно определить искомое К.

59

В общем случае уравнение (2-12) удобно решать графически, для чего перепишем его в виде

------------------------- - I n

--------------

—---------- .

(2.13)

2 ( K 2 0q + Од)

Т у / a 2 +

К 0 q

 

Для случая aR = 0 (т.е. уровень, выбросы за который запрещены, детер­ минирован) уравнение (2.12) удается разрешить относительно К:

К =

т о

.

(2.14)

2.-------

mq + oq у/Тл

 

 

где А = - In

2тг(-1пЯ)

 

(2.15)

Ту /а 1 + р2

Втом случае, когда уровень, выбросы за который запрещены, является случайным с законом распределения Релея

/( /? ) = 4 - р ( - ~ ) ,

 

 

а

 

 

 

для V(R/T) имеем [16]

 

 

V(R/T)=

o l + а2

(2.16)

 

 

2п

 

Отсюда по (2.2) определим надежность

 

Н —ехр

Т

os а

 

(2.17)

2 п

о |+

а2

 

 

 

Подставим вместо

 

 

 

os = Koq; Og = Koq / а 2 + 02 ,

(2.18)

тогда

 

 

 

 

 

I

TKOq у / a 2 + f a

(2.19)

H = e x p < - ------ 2_ _

----- ---

 

I

2 l T( K2 0 q + a 2 )

 

Решая относительно К, получим

 

К2 -

в

 

 

 

 

 

2П0Й ( - Ь Н )

 

 

Обозначив

' у/ а2 + Г

_

 

 

 

= В, получим

 

 

2я( - In//)

 

 

*«-£<4 *^4 -о-

( 2.20)

 

60

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