книги / Расчет на надежность композитной балки
..pdfМинистерство образования и науки Российской Федерации
Государственное образовательное учреждение высшего профессионального образования «Пермский государственный технический университет»
Р.Я. Газизов
РАСЧЕТ НА НАДЕЖНОСТЬ КОМПОЗИТНОЙ БАЛКИ
Утверждено Редакционно-издательским советом университета
в качестве учебно-методического пособия
Издательство Пермского государственного технического университета
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 |
|
Ω(S−F >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 = |
∞ z−x |
|
(7) |
|||
∫ |
∫ |
f ( x, y)dy dx. |
||||
|
D |
−∞ |
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 1
Плотность распределения |
f (z) получим, дифференцируя (7) |
||||||
по z, которая входит как параметр в верхний предел интеграла, |
|
||||||
|
∂F (z) |
∞ |
∂ z−x |
|
∞ |
|
|
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 |
|
−∞ |
s−z |
|
|
|
|
|
|
6 |
|
|
|
|
|
Рис. 2
Дифференцируя (11) поz , получим |
|
|
(s, s − z)ds. |
|
||||||||||||||
f (z) |
= ∫ |
∂ ∫ f (s, |
f0 )df ds = ∫ f |
|
||||||||||||||
|
|
∞ |
|
|
|
∞ |
|
|
|
|
|
∞ |
|
|
||||
|
|
−∞ |
∂z |
s−z |
|
|
|
|
|
−∞ |
|
|
||||||
Если случайные величины 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− |
(s−ms )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 |
|
|
|
|
|
∞ |
|
|
− |
(s−ms ) |
2 |
− |
(s−z−m f0 )2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
f (z) = |
|
|
|
|
|
|
∫ e |
|
2 |
σ2s |
|
e |
|
|
2σf0 |
|
ds; |
|
(16) |
||||||||||||
|
2πσsσf0 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
f (z) = |
|
|
|
|
|
|
|
∫ e−As2 ±2B(z)s−C(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 |
|
|
|
|
|
|
|
π |
|
|
− |
AC−B2 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
e |
A |
|
|
, |
|
|
|
|||||||||||||
|
|
|
|
2πσsσf |
0 |
|
|
|
A |
|
|
|
|
|
|
|
|
|||||||||||||||
или после преобразований |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
f (z) = |
|
|
|
|
1 |
|
|
|
|
|
|
− |
(z−mz )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
β= 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