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

книги / Расчет на надежность композитной балки

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
268.16 Кб
Скачать

Министерство образования и науки Российской Федерации

Государственное образовательное учреждение высшего профессионального образования «Пермский государственный технический университет»

Р.Я. Газизов

РАСЧЕТ НА НАДЕЖНОСТЬ КОМПОЗИТНОЙ БАЛКИ

Утверждено Редакционно-издательским советом университета

в качестве учебно-методического пособия

Издательство Пермского государственного технического университета

2011

УДК 539.3 Г12

Рецензент

канд. физ.-мат. наук А.В. Зайцев (Пермский государственный технический университет)

Газизов, Р.Я.

Г12 Расчет на надежность композитной балки: учеб.-метод. пособие / Р.Я. Газизов. – Пермь: Изд-во Перм. гос. техн. ун-та, 2011. – 23 с.

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

Рекомендуется для самостоятельных занятий студентам аэрокосмического факультета.

УДК 539.3

© ГОУ ВПО «Пермский государственный

технический университет», 2011

2

Рассмотрим «традиционный» расчет на прочность в детерминированной постановке.

Например, в методе расчёта по предельным состояниям максимальная действующая нагрузка (точнее, напряженно-деформирован- ное состояние, вызванное этой нагрузкой) сравнивается с нагрузкой, соответствующей предельному состоянию, которое определяет несущую способность конструкции. Несущая способность конструкции – наступление предельного напряженно-деформированного состояния, которое соответствует потере работоспособности конструкции или её разрушению. К таким предельным состояниям относятся: появление пластических деформаций, потеря статической устойчивости, разрушение.

В общем случае можно записать

S F > 0,

(1)

где S – функция, характеризующая предельные возможности конструкции; F – функция, характеризующая реальное состояние конструкции.

Например, пусть S зависит от механических характеристик материала конструкции,

S = S (σТ , σВ, µ, E, G) ,

где σB – предел прочности; µ – коэффициент Пуассона; E, G

модули упругости соответственно первого и второго рода. Тогда функция F зависит от напряжений, возникающих в элементах конструкции при нагружении,

F = F (σx , σy , σz , τxy , τxz , τyz ).

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

Традиционные методы расчета по предельным состояниям и по допускаемым напряжениям в явном виде не учитывают возможные

3

случайные разбросы, т.е. не учитывается вероятностный характер предельных состояний конструкции или вероятностный характер реального состояния конструкции. Поэтому оценивать работоспособность конструкции логичнее не по детерминированным неравенствам типа (1), а по вероятности выполнения этих неравенств, т.е.

P (S F ) >0 ,

(2)

 

 

 

где P – вероятность безотказной работы.

 

Введем понятие надежности системы H ,

которая оценивается

вероятностью выполнения неравенства (1):

 

H = P (S F ) >0 .

(3)

 

 

 

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

ции [1].

Чтобы найти вероятность P, надо знать совместный закон распределения f (z) случайной величины Z = S F при известных за-

конах распределения S и F.

Если

f (z) известно, то

 

P(Z > 0) = f ( z)

dz =

∫∫

f (S, F )dSdF.

(4)

0

 

(SF >0)

 

Для получения численных значений H надо знать (определить) закон распределения случайной величины Z, функционально зави-

сящей от двух непрерывных случайных величин X и Y:

 

Z ( X , Y )

(5)

при известном их совместном законе распределения плотности вероятности f (x, y).

В общем случае закон распределения случайной величины Z имеет следующий вид:

4

 

 

∫∫

f (x, y)dxdy.

 

F ( z) = P(Z < z) = P ( X , Y ) D =

 

(6)

 

D (Z <z)

 

 

Математическая задача определения F ( z) сводится, таким об-

разом, к решению двукратного интеграла

 

 

 

F ( z) = ∫∫

f (x, y)dxdy.

 

 

D (Z ,z)

 

 

 

 

Для случая, когда Z = X +Y

интеграл берётся по области

D

(рис. 1), где x + y < z, поэтому, фиксируя Z

и полагая y = z x,

по-

лучаем конкретные пределы интегрирования, т.е.

F ( z) = ∫∫ f (x, y)dxdy =

zx

 

(7)

f ( x, y)dy dx.

 

D

−∞

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1

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

f (z) получим, дифференцируя (7)

по z, которая входит как параметр в верхний предел интеграла,

 

 

F (z)

zx

 

 

f (z) =

 

=

 

f (x, y)dy dx =

f ( x, z x)dx.

(8)

z

 

 

−∞

z −∞

 

−∞

 

5

Поскольку случайные величины X и Y равноправны, то можно было исключать не y, а x и получить ещё одно выражение для

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

f (z) :

f (z) =

f (z y, y)dy.

 

(9)

 

−∞

 

 

 

Если X и Y независимы, то плотность распределения Z опре-

деляется как

 

 

 

 

f (z) =

( x) f2 ( z x)dx,

 

f1

(10)

или

−∞

 

 

 

 

 

 

 

 

 

 

 

f (z) =

f1 ( z y) f2 ( y)dy.

 

 

−∞

 

 

 

Проинтегрировав правую часть (10) по

x,

получим функцию,

зависящую от z.

 

 

 

 

При определении вероятности безотказной работы требуется найти функцию распределения F (z) случайной величины z, равной разности двух случайных величин s и f:

z = s f .

При известном совместном законе распределения f (s, f0 ),

т.е.

F ( z) =

∫∫

f (s, f0 )ds df0.

 

D (Z <z)

 

 

 

 

 

На рис. 2 показана заштрихованная область, где s f0 < z.

По-

этому интегрирование по области

D,

как и в предыдущем случае,

можно заменить интегрированием

по

f0

при фиксированном

z ( f0 = s z), а затем по s:

 

 

 

 

 

 

 

∫∫ f (s, f0 )ds df0 =

f (s,

 

 

 

f

0 )df0 ds.

(11)

D

 

−∞

sz

 

 

 

 

 

 

6

 

 

 

 

 

Рис. 2

Дифференцируя (11) поz , получим

 

 

(s, s z)ds.

 

f (z)

=

f (s,

f0 )df ds = f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

z

sz

 

 

 

 

 

−∞

 

 

Если случайные величины s и

f0 независимы, то имеем

 

 

 

 

 

 

 

f1 (s) f2 (s z)ds,

 

 

 

 

 

 

f (z) =

(12)

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

или, изменяя порядок интегрирования,

 

 

 

 

 

 

 

z+ f0

 

 

 

 

 

 

f (z) =

 

 

 

 

f (s, f0 )ds df0 =

 

f

(z + f0 , f0 )df

=

 

 

 

 

 

 

−∞ z

−∞

 

 

 

 

 

 

 

−∞

z

 

 

 

 

 

 

 

 

 

+ f0 ) f2 ( f0 )df0.

 

 

 

 

 

 

 

 

= f1 (z

 

(13)

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

Рассмотрим случай, когда s и

f0

имеют нормальные распре-

деления:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

e

(sms )2

 

 

 

 

 

 

 

 

f (s) =

 

 

 

2σ2s ;

 

(14)

 

 

 

 

 

 

 

2πσs

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

( f0 m f0 )2

 

 

 

 

 

 

 

f ( f0 ) =

 

 

 

2σ2f0

 

 

 

 

 

 

 

 

 

 

e

 

.

(15)

 

 

 

 

 

 

2πσf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

Закон распределения случайной величины z (12) в этом случае имеет вид

 

 

 

 

 

 

1

 

 

 

 

 

 

 

(sms )

2

(szm f0 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (z) =

 

 

 

 

 

 

e

 

2

σ2s

 

e

 

 

2σf0

 

ds;

 

(16)

 

2πσsσf0

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (z) =

 

 

 

 

 

 

 

eAs2 ±2B(z)sC(z)ds,

 

 

 

 

 

2πσsσf

 

 

 

 

 

 

 

 

 

 

0

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z +mf0

)2

A =

σ2s 2f

 

 

 

m

 

 

 

z +mf

0

; C =

 

m2

+

 

 

0

; B =

 

s

 

+

 

 

 

 

 

 

 

 

 

 

s

 

.

2σ2s

 

2σ2s

 

 

 

2σ2f

 

 

 

 

 

2σ2f

 

σ2f

 

 

 

 

 

 

 

0

 

 

 

 

 

 

2σ2s

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

Воспользовавшись табличными значениями для определённых

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (z) =

1

 

 

 

 

 

 

 

π

 

 

ACB2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

A

 

 

,

 

 

 

 

 

 

 

2πσsσf

0

 

 

 

A

 

 

 

 

 

 

 

 

или после преобразований

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (z) =

 

 

 

 

1

 

 

 

 

 

 

(zmz )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

2σz

,

 

 

 

 

 

(17)

 

 

 

 

 

 

2πσz

 

 

 

 

 

 

 

 

 

 

где mz = ms mf

0

; σz =

σ2s

2f

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определив

 

f (z),

находим

 

 

вероятность

безотказной работы

(надёжность):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H = P(Z > 0) = f ( z)dz.

 

 

 

 

 

(18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Введя новое обозначение

8

Рис. 3

β= z (ms mf0 ) ,

σz

получим

 

1

e

β2

(β0

 

z=0 ).

 

H =

2 dβ

 

(19)

 

 

2π−β0

 

 

 

 

 

 

График подынтегральной функции показан на рис. 3. Функция f (β) симметрична относительно вертикальной оси, поэтому

1

 

 

 

 

 

1

 

 

β0

 

 

1

0

1

β0

(...)dβ =

 

 

 

 

(...)dβ =

 

 

(...)dβ+

(...)dβ =

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

−β0

 

 

 

 

 

2π

−∞

 

 

2π −∞

2π

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0,5 +

 

 

0

(...)dβ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

Окончательно

 

 

числовое

 

 

 

 

 

 

 

 

значение надёжности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ms m f0

 

 

β2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

σz

 

 

 

 

 

 

 

 

 

 

H = 0,5 +

 

 

0

e

2

dβ. (20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл, входящий в правую часть (20), является табличным интегралом («интеграл вероятности»).

Пусть распределение несущей способности (прочности) подчиняется нормальному закону с плотностью вероятности f1 (s), математическим ожиданием прочности ms и среднеквадратическим (стандартным) отклонением σs , и распределение сил (нагрузки) подчиняется нормальному закону с плотностью вероятности f2 (s0 ), математическим ожиданием усилия mf0 и среднеквадратическим (стандартным) отклонением σf0 (рис. 4).

9

Рис. 4. Распределение плотностей вероятностей усилий и прочности в узлах металлоконструкций

Функцию надёжности такой системы определяют, используя зависимость [1]:

 

m m

 

 

 

P = Ф

1

2

 

,

(21)

σ2

2

 

 

 

 

 

1

2

 

 

 

где Ф(up ) – нормированная нормальная функция распределения

(см. таблицу приложения).

Оценим надёжность двухопорной балки сложного сечения (рис. 5). На рисунке представлена схема нагружения с эпюрами поперечных сил Q и изгибающих моментов M .

10

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