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

Вероятностные методы расчета конструкций

..pdf
Скачиваний:
33
Добавлен:
15.11.2022
Размер:
2.24 Mб
Скачать

Далее получим вероятностные характеристики смещения:

2

u R

;

2

1

uik2 Rk2 S f ( )d

,

myi mf k 1

( p2 2 )

Dyi k 1 2

( p2 2 )2

 

 

ik k

 

 

 

 

 

 

 

k

 

 

 

 

 

k

 

где S f ( ) – спектральная плотность действующихвнешних сил.

Рассмотренная система является узкополосным фильтром, который пропускает в основном составляющие с частотами, близкими к pk , а потому можно принять, что спектральная

плотность в пределах полос пропускания системы постоянна (рис. 3.9). Предполагая, что спектральная плотность внешних сил не имеет острых пиков и разрывов, для дисперсии перемещения масс получим:

 

2

2 2

 

 

d

 

 

 

Dyi

 

S f ( pk )uik Rk

 

 

 

 

.

2

( p

2

2

2

 

k 1

 

)

 

 

 

 

 

 

k

 

 

 

 

Рис. 3.9

2. Рассмотрим теперь случай, когда действуют различные случайные силы, статистически независимые между собой.

Можно использовать два метода решения.

Метод главных координат. Представляя действующие силы в виде fi (t) Fi (0) f0i (t), перейдем, как и прежде, от уравнений (3.31) к уравнениям в главных координатах:

101

 

 

2

 

 

 

u jk f j

 

 

2

j 1

Fk (t) (k = 1,2).

 

2

qk pk qk

 

 

u2jk mj

 

j 1

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

Методинтегральногопредставленияслучайныхфункций[20].

Будем полагать, что действующие на систему силы стационарны и их можно представить в виде интеграла Фурье:

fk (t) Fk ( )ei t d .

Разыскиваемые переменные yk (t) также запишем в виде

yk (t) yk 0 ( )ei t d .

Подставляя оба выражения в систему (3.37), получим:

(1 2m1 11 ) y10 2m2 12 y20 11F1 12 F2 ;

2m1 21 y10 (1 2m2 22 ) y20 21F1 22 F2 .

2

Решение можно представить в виде yk 0 Hki ( )Fi ,

i 1

 

H

( )

 

 

(1 2m

22

)

2i

2m

 

;

где

1i

 

2

 

2 12

 

 

 

 

 

 

 

 

 

 

 

1i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H2i ( )

 

2i

(1 2m )

 

2m

;

 

 

1 11

1i

 

1 12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

102

(1 2m1 11 )(1 2m2 22 ) 4m1m2 12 21.

Находим спектральные плотности решений:

Syi ( ) Hi1 ( ) 2 S f1 Hi2 ( ) 2 S f2Hi1 ( )Hi2 ( ) (S f1 f2 S f2 f1 );

Sy1 y2 ( ) H11 ( )H21 ( )S f1 H11 ( )H22 ( )S f1 f2

H12 ( )H21 ( )S f2 f1 H12 ( )H22 ( )S f2 .

3.2.3. Вынужденные колебания систем с распределенными параметрами

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

Уравнение колебаний системы имеет вид:

L(w(x, y, z,t)) m

2w

q(x, y, z,t),

(3.38)

t2

 

 

 

где w(x, y, z,t) – динамическое перемещение точек системы; m – масса элементарного объема; L(w(x, y, z,t)) – линейный диффе-

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

– дляизгибныхколебаний стержня L(w(x,t))

2

 

2

 

;

 

 

EI

 

 

x

2

 

 

 

 

x2

 

 

103

– для продольных колебаний стержня L(w(x,t))

x GIk x ;

для прямоугольной пластинки L(w(x, y,t)) D w, где D

цилиндрическая жесткость,

D

 

Eh3

;

2

 

2

; E, G,

12(1

2 )

x2

y2

 

 

 

 

 

– константыматериала; I, Ik – моментыинерции; h – толщина.

Решение уравнения (3.38) можно найти в виде разложения по формам свободных (собственных) колебаний [4]:

 

 

y(x,t) uk (x) fk (t),

(3.39)

k 1

где uk (x) – прогиб при k-й форме свободных колебаний, кон-

кретный вид которой зависит от граничных условий;

fk (t) –

функция времени, подлежащая определению.

 

Рассмотрим, например, задачу о вынужденных колеба-

ниях стержня (рис. 3.10) длиной l

с постоянным поперечным се-

 

 

 

 

 

чением при действии на

 

 

 

 

 

него случайной стацио-

 

 

 

 

 

нарной распределенной

 

 

 

 

 

нагрузки q(x,t) с за-

 

 

 

 

 

данным значением кор-

Рис. 3.10

 

 

 

реляционной

функции

 

 

 

Kq (x, t, x1,t1 ).

 

Уравнение движения

 

 

 

 

 

EI

4 y

m

2 y

q(x,t).

(3.40)

x4

t2

 

 

 

 

Умножим (3.40) на uk (x) . Тогда, учитывая условие ортогональности собственных форм колебаний и разложение (3.39), получим уравнения, из которых можно найти функцию fk (t) :

104

f

p2

f

k

F (t).

(3.41)

k

k

 

k

 

Здесь pk k-я собственная частота;

Fk (t) – обобщенная

сила, равная сумме произведений возмущающих сил на перемещения точек их приложения при k-м нормальном колебании, отнесенной к обобщенной массе Mk :

l

 

l

q(x,t)uk (x)dx

Mk m uk (x) 2 dx,

Fk (t)

0

 

 

 

.

l

 

2

 

0

 

m uk (x)

dx

 

 

 

 

 

 

 

0

 

 

 

 

Таким образом, задача сведена к задаче, рассмотренной ранее (колебания с n степенями свободы). Некоторые особенности ее решения состоят лишь в определении обобщенных сил.

Имеем граничныеусловия:

 

y(0, t) y(l, t ) 0;

2 y

0 .

 

x2

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

x l

Для рассматриваемого случая закрепления собственные

формы и частоты имеют вид: u

 

(x) sin kx ,

p2 EI

4k4 .

 

 

 

 

k

 

l

k

ml4

Это означает, что решение уравнения (3.40) можно пред-

ставить в виде

y(x,t) fk (t)sin k x .

 

 

(3.42)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

l

 

 

 

 

Обобщенные силы из уравнения (3.41)

 

 

 

F (t)

2

l q(x,t)sin k x dx.

 

 

(3.43)

 

 

 

k

 

ml

0

 

 

l

 

 

 

 

 

 

 

 

 

 

Взаимную корреляционную функцию обобщенных сил определим по формуле

 

 

 

2 2 l l

KFi Fj

(x,t, x1

,t1 )

 

 

K f

 

 

 

ml

0 0

(x,t, x1,t1 )sin i lx sin jl x dxdx1.

105

Для стационарных случайных функций взаимная корреляционная функция зависит от разности моментов времени

( t t1 ), т.е. KFi Fj (t,t1 ) KFi Fj ( ) .

Корреляционная функция прогиба

Ky (x, x1

, ) K fi f j

( )sin i x sin

j x1 .

(3.44)

 

 

 

 

 

 

 

 

 

 

i 1

j 1

 

l

 

l

 

 

Для определения взаимных корреляционных функций K fi f j предварительно найдем взаимные спектральные функции S fi f j .

Для этого обратимся к уравнению (3.41).

Распределенная нагрузка q(x,t) является стационарной

функцией, поэтому диагональная матрица передаточных коэффициентов

 

 

 

Hk (i )

 

 

 

1

 

.

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

pk

 

 

 

Следовательно, спектральная функция

 

 

 

S f

( )

 

Hk (i )

 

2

SF ( ),

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

k

 

 

а взаимная спектральная функция

 

 

 

 

 

 

 

 

 

S f

f

(i )

 

H j (i )

 

Hi (i )

 

S F F ( ).

 

 

 

i

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i j

Зная взаимную спектральную плотность S fi f j , определим

взаимную корреляционную функцию:

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

i

 

 

K fi f j ( )

 

S fi f j (i )e

 

d .

 

2

 

Подставляя полученное соотношение в (3.44), получим окончательное выражение корреляционной функции прогиба.

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

106

Например, если ввести спектральные плотности для функций q(x,t) и y(x,t) следующим образом:

Sq (x, x1

 

1

 

Kq (x, x1, )e i d ;

 

, )

 

 

2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

i

 

 

Sy (x, x1

, )

 

 

Ky (x, x1, )e

 

d ,

 

2

 

 

 

 

 

 

 

 

то, рассматривая уравнение (3.40) и применяя к нему дифференциальное преобразование Фурье по времени, получим уравнение связи для спектральных плотностей:

EI x44 m 2 EI x44 m 2 Sy (x, x1, ) Sq (x, x1,t). (3.45)

Имеем граничные условия: Sy 2 S2y 0 при x = 0, x = l.

x

Решение уравнения (3.45) можно представить в виде:

 

 

 

 

Sij ( )

 

 

 

 

 

i x

 

j x

Sy (x, x1, )

 

 

 

 

 

sin

 

sin

1

.

2

2 2

 

2

 

l

 

 

i 1 j 1

( pi

)( pj

)

 

 

 

l

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2 l l

 

 

 

i x

 

 

 

j x

 

 

 

Sij ( )

 

 

Sq (x, x1, )sin

 

sin

 

1

dx dx1.

 

l

 

l

ml

0 0

 

 

 

 

 

 

 

 

 

Корреляционная функция прогиба

Ky (x, x1, ) Sy (x, x1, )ei d .

Задачи для самостоятельного решения

1. Упругий однородный шарнирно закрепленный стержень длиной l, показанный на рисунке, совершает случайные изгибные колебания под действием случайной сосредоточенной силы q1(t), приложенной в точке х = х1.

107

Определить вероятностные характеристики поперечных смещений балки при ее колебании.

2. Упругий однородный шарнирно закрепленный стержень совершает случайные изгибные колебания (как показано на рисунке) под действием поперечной случайной нагрузки q(x,t) и слу-

чайной сосредоточенной силы q1(t), приложенной в точкех = х1.

Определить вероятностные характеристики поперечных смещений балки при ее колебании. Принять, что нагрузка q(x,t) дельта-коррелирована как по координате, так и по вре-

мени:

Kq (x,t, x1,t1 ) s (x x1 ) (t t1 ).

3. Струна, изображенная на рисунке, возбуждается сосредоточенной силой q(t) в точке А [7].

108

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

иметь дополнительные прогибы y и дополнительные скоростиy , которые внесут искажения в снимаемый сигнал.

Определить дисперсию скорости y в точке В при стацио-

нарных колебаниях струны, если спектральную плотность случайной составляющей силы q можно представить в виде

S

q

( )

a

.

b2 2

 

 

 

109

4. ОПРЕДЕЛЕНИЕ НАДЕЖНОСТИ КОНСТРУКЦИЙ НА СТАДИИ ПРОЕКТИРОВАНИЯ

4.1.Вероятностные основы задач надежности конструкций

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

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

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

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

110