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

книги из ГПНТБ / Основы динамики сооружений учеб. пособие

.pdf
Скачиваний:
18
Добавлен:
19.10.2023
Размер:
6.71 Mб
Скачать

Умножив обе части равенства на т(х) и поменяв их местами, получаем

 

=

 

(д)

 

п-1

 

 

Данное

выражение представляет собой разложение импуль­

са s(x),

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

длине балки, в ряд по

функ­

циям Х п (х), характеризующим

главные формы колебаний

и на­

зываемым иногда собственными функциями. Таким образом, по­

стоянные

А„, которые мы должны определить из выражения

(д),

могут быть найдены как

коэффициенты этого ряда. Для этого

поступим

так же, как и

при

определении

коэффициентов

ряда

Фурье. Умножим обе

части равенства (д)

на какую-либо

одну

собственную функцию,

например Х т(х), и проинтегрируем их от

нуля до I:

I

 

і

со

 

 

 

 

 

 

j

s (х) Х т (х) dx =

j‘ Х т (х) 2 A nwnm (х) Х п (х) dx.

 

О

 

 

0

я = 1

 

 

Развернем сумму в правой части этого равенства и проинтегри­

руем каждое слагаемое в отдельности:

 

 

 

I

 

 

I

 

 

 

s (х) Х т (х) dx

j т (х) Х х(.х) Х т (х) dx

 

о

г

 

о

 

 

 

 

 

 

 

•ф Л2Ш2 j* tn (л:) Х%(л) X т(л:) dx -ф ... -ф

0

2

+ -4л $ т (X) х п (х) Х т (х) dx 0

1

A n+\<s>n+\ J tn (jc) Х п+\ (X) Х т(.х) dx -ф ...

0

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

1

 

оJ т

(■ #) х п ( х

) х

т ( х ) dx.

 

 

Однако в силу ортогональности

функций,

описывающих глав­

ные формы колебаний

(см. §

15), все интегралы,

имеющие

в подынтегральном

выражении

произведения

неодинаковых соб­

ственных функций

(если

тфгі),

обращаются

в нуль.

В правой

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

I

I

I s (л) Х п {х) dx =

Ап<лп j tn (x) X n2 (x) dx.

о

0

.170

Отсюда получаем формулу для вычисления постоянной А п:

г

j s(x) Х п(лс) dx

1 о

(е)

^т{х) Х п2 (х) dx

о-

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

 

г

 

*

J s (х) Х п (х) dx

 

у ( X , 0 =

---- ------------------------ Х п {х) sin со/.

(9.3)

л- i

" JГ m ( x ) X 2n(x)dx

 

 

о

 

Пусть мгновенный импульс действует не в начальный момент времени, а в момент t = u. Перенесем начало отсчета времени к мо­

менту t = u, т.

е. перейдем к новой переменной t = t и.

В этом слу­

чае в

выражении

(9.3) вместо sin со/ необходимо

подставить

sin со/ =

sin <ол (t — и).

Тогда будет

 

 

 

 

 

I

 

 

 

~

j

s (/) Х п(л;) dx

 

У (х, 0 =

~

1----------------------

Х п (х) Sin CD„ (t и) du. (9.4)

 

 

п=1

[ т (х) Х 2п (х) dx

 

 

 

 

о

 

 

Пример 13. В качестве примера рассмотрим действие на простую балку равномерно распределенного по длине балки импульса s. Балка имеет постоян­ ную жесткость Е І и постоянную по ее длине погонную массу т. Известно

(см. § 16), что главные формы колебаний простой балки имеют вид

*„(•*) = Sin

а частоты колебаний определяются по формуле

я2я2 и ГEl

„п

“ и = -г-

I / — =

пЧ .

I-

I т

 

где га=1, 2, 3, ... ;

 

(при я=1).

<йі — частота первого тона колебаний балки

171

Вычислим интегралы, входящие в знаменатель и числитель выражения (9.3):

I

 

I

 

 

I

ппх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ml

 

 

 

 

 

 

~2

 

 

Г

 

Ппх

si

 

 

 

s (х) Х п (х) dx 5 I sin —j - dx = — (1 —cos me).

 

Так как при нечетных значениях

я

соэпя = —1, а при

четных значениях п

cosnre=l,

то получаем,

что при п = 2,

4,

6 , ...

интеграл в

числителе выраже­

ния (9.3)

обращается в

нуль, а поэтому соответствующие

слагаемые

суммы

в данном выражении выпадают.

 

 

 

 

 

Этого следовало ожидать, так как импульс, распределенный по балке равно­ мерно, является симметричным, а поэтому представляется естественным, что он должен вызывать только симметричные формы колебаний, т. е. формы колеба­ ний с нечетными значениями я. При п 1, 3, 5 ,...

I

I

ппх

si

Г

С

I sX n (х) dx = s

1

sin —j - dx =

2 — .

бо

Подставляя значения вычисленных интегралов, частоты колебаний со„ и соб­ ственной функции Х п (х) в уравнение (9.3), получим выражение для прогиба балки:

00

 

9 si

 

 

У (•*. 0 = V

1

ПП ппх

Я2* 2 Л Г E l

Sin —— Sin Я2<0^ :

 

' "У

 

~ W У m

2

П = 1, 3,

5,

 

4sP

1

ппх

лз у Ш

sin —j—sin я 2м,£.

 

 

/1 = 1 , 3, 5 ,. . .

 

Заметим, что в момент времени, равный четверти периода основного тона свободных колебаний, функция sin n2®d получит одновременно для всех главных форм колебаний максимальное значение, равное единице. В самом деле, пусть

'X = 2«! ’ тогда

sin

= sin «2 = 1 (при п = 1, 3, 5 , ... )

Поэтому уравнение изогнутой линии балки при ее максимальном отклоне­ нии в случае колебаний от действия мгновенного равномерно распределенного импульса имеет вид

 

I1

^

1

ппх

 

4S/2

 

 

 

(ж)

шах у (х) :

шElm,

 

^

sin ~ r -

-3 у

 

 

 

 

/1= 1, 3,

5 , . . .

 

 

Определим максимальный прогиб в середине балки при х = —

4s/2

ѴЧ

1

Пп

max у

Ь

^

sinT -

пз У Ein

/1= 1, 3, 5, . . .

172

Значение вошедшей сюда бесконечной суммы известно из курса математики:

 

1

пк

,

1

1

1

лЗ

2

із sin ~ 2

1 — Зз + 5з

73

; 32 ’

/1 = 1, 3,

5 , . . .

 

 

 

 

 

 

С учетом этого получаем

 

sP

тс3

1

Si2

 

Ушах =

 

 

 

У'Elm

32

8 |

£ /m ’

 

 

* 3

 

 

 

 

 

 

sP

 

 

 

 

) max • : 0,125

 

(з)

 

 

 

 

 

у Elm

 

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

 

М=*

(Ру

 

 

 

 

 

— £/ dx 2 '

 

 

 

 

Продифференцируем дважды уравнение (ж ):

 

 

 

 

[шах у (х)]

 

 

4s

^

1

 

п~х

dx'2-

 

я У Elm

■‘”J

—sin

 

J—

 

 

 

 

 

 

 

 

 

/1= 1, 3,

5, . . .

 

 

 

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

дифференциальное уравнение

изгиба балки:

 

 

 

 

 

 

 

 

max М

тс

г

т

"

*

 

(и)

 

 

 

 

 

 

я —1, 3, 5, .

 

 

 

 

Определим изгибающий момент в середине балки при

х =

2 '

, l £

l / «

2

i

sl„

'

 

 

■к

У

m

 

n

2

 

 

л = 1, 3, 5, . . .

Значение вошедшей сюда бесконечной суммы известно из курса математики:

2

1 Пг.

1 1 1

— sin

= 1 — 3 - + 5 - — т +

••• = Т -

л = 1, 3, 5, . . .

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

ш

 

М - У т

( к )

§ 38. Действие нагрузки, меняющейся во времени

 

по произвольному закону, на систему с бесконечным числом

 

степеней свободы

 

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

173

костью EI, распределенными по длине балки по произвольным за­ конам.

Пусть на балку действует нагрузка, изменяющаяся по длине балки и во времени также по произвольному закону р(х, t). Эту нагрузку можно представить в виде поверхности, описываемой уравнением z = p(x, t) в системе координат г, х, t (рис. 79). Рас-

Рис. 79

смотрим сначала действие на балку элементарного импульса р(х, и) du, приложенного в момент времени и, и определим элемен­ тарное перемещение, которое возникнет от его действия в момент времени і. Считая элементарный импульс s(x) =р{х, и) du мгновен­

ным, применим к рассматриваемому случаю формулу (9.4) преды­

дущего параграфа:

I

“ j J Р (X, и) du Х п (х) dx

dy {х, t ) = 2 j — ■- — i------------------------ X n (x) sin шя (t — u).

n=1 f m (x) XI (x) dx

Представим нагрузку p(x, t) как сумму бесконечного мно­ жества элементарных импульсов р(х, и) du. Тогда перемещение оси балки в момент времени t от действия этой нагрузки можно найти, если предыдущее выражение проинтегрировать в пределах от нуля

I

Х„( х )

*

j'Р( х ,

u ) X n{x)dx

(х, 0 = 2

I

I------------------

------ sin шп (t и) du. (а)

П~1

о

I' т(х) Х 2п(х) dx

 

 

о

 

174

Заметим, что дробное выражение, стоящее под знаком интегра­

ла по

и,

представляет собой коэффициент

разложения

нагруз­

ки р(х,

t)

в ряд по собственным функциям

Х п (х). Чтобы убе­

диться в этом, представим функцию р (х, t) в виде ряда

 

 

 

со

 

 

 

 

р ( х , t) = % p n( t ) m ( x ) X n(x).

(9.5)

 

 

л - 1

 

 

Определим значение коэффициентов разложения pn (t),

приме­

няя прием предыдущего параграфа при нахождении коэффициен­ та Ап разложения в ряд импульса s(x). Умножим левую и пра­

вую части предыдущего

равенства

на Х т (х) и проинтегрируем

обе части по длине балки в пределах от нуля до /:

I

I

со

J Р (■*, 0 х т (х) dx

= J Х т(je) ТіРп(І)т (х) Х п(X) dx.

О

0

/2=1

Записав сумму в правой части почленно и используя свойство

ортогональности собственных функций Х п (х), получим

I

 

I

j р (х, t) Х п(je) dx = рп (t) J т (х) XI (х) dx,

о

 

о

откуда

I

 

 

 

 

J р (х, f) Х п (х) dx

Рп V) =

о

(б)

I

 

j т (х) Х п2(х) dx

 

о

 

Таким образом, выражение для у(х, t) может быть представ­ лено в виде

V х п W

t

 

 

рп{и) Sin (ü„ (t

и) du,

(9.6)

У (х, 0 = 2 j — т г - ^

п=1

где р п (и) определяется по формуле (б).

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

Анализируя формулу (9.6) и сопоставляя ее с формулой (7.17) для определения перемещения в случае действия произвольной нагрузки на систему с одной степенью свободы, приходим к естест­ венному заключению, что каждая главная форма колебаний систе­

175:

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

 

 

рп (х,

t) = pn( t ) m ( x ) X n(x)

формулу

(7.14 в), получим соответствующее слагаемое ряда (9.6):

У п (*>

t) =

1___

\рп (и) т {х) Х п(л:)] sin u)„ (t и) du =

т (д:) шп

 

 

 

 

 

о

t

 

 

 

 

 

Х п ( х )

рп (и) sin шп (t — а) du.

В расчетной практике

распространен случай, когда нагрузка

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

р ( х , t) = p (x )f( t) .

(9.7)

Определим по формуле (б) коэффициент р п (t)\

 

I

 

S Р (х ) х п ( х ) dx

 

P n ( t ) = f { t ) - l---------------------= /(0 Р „ -

( в)

I т (д:) Х 2п (х) dx

 

о

 

Здесь множитель р„ представляет собой коэффициент

разло­

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

р ( х )

=

£ р пт ( х ) Х п(х),

(9.8)

то

 

п =1

 

I

 

 

 

 

j

Р (х ) Х п (X) dx

 

Р п =

Г

(г)

 

I т (х) X 2 (х) dx

 

Таким образом, подставив (в) в формулу (9.6), получим выра­ жение для определения перемещения:

ооt

У (X, 0 = S Р пЛ,'и ^ J / ( “ ) sin шп (t - и) du.

(9.9)

Л«1

 

176

Для ряда частных случаев нагрузки значения интегралов, вхо­ дящих в слагаемые выражений (9.6) и (9.9), были получены в главе 8. Следует иметь в виду, что для каждого слагаемого значения интегралов будут различными, так как подынтегральные выражения не одинаковы для разных п.

Пример

14. Рассмотрим действие на

простую

балку постоянного

сечения

£7 = const с

постоянной погонной массой

от= const

кратковременной

нагрузки,

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

р (х , t) = р f(t),

где р= const.

Главные формы и частоты собственных колебаний простой балки известны:

Хп (х) = sin птех

Т

<1)

п2я2 -|

Г E l = П2шJ .

п

ОТ

 

І2 Г

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

нагрузки в ряд по собственным функциям:

I

Р пт.х

рI sin j - dx

Рпъх

от I sin2 —у—dx

о

Интегралы, входящие в числитель и знаменатель, нами вычислены в преды­ дущем примере 13. Тогда

Ря

Iот/

т

лтсот '

Подставив в формулу (9.9) значения р„,

°°

ХГЛ

УІ*'

 

4 р / 2

 

 

 

7Т/

Elm

1

 

 

 

 

яЗ )

5, . . .

 

 

 

 

 

Я - 1 , 3,

Отсюда

 

 

 

 

 

 

М = — ЕІ

дЪ<

( X ,

t) =IP i f

§l

V

 

 

 

я г

от

/

 

öx2

 

 

 

 

Х п (х),

шп, получим

 

пкх

t

 

sin I

f ( u ) sin шп (t — и) du.

(д)

 

" 3

( У

пт.х

sin

f (u)s'mb3n ( t а) du. (е)

л - i , з, 5, . . .

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

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

12 Основы динамики сооружений

177

нии (9.7) необходимо принять f(t) = 1. Для определения переме­ щения у(х, t) вычислим интеграл, входящий в выражение (9.9):

t і

J f ( u ) sin (o„ (t ti) da — J

1 • sin con (t u)da —

 

о

 

 

 

0

 

 

 

— — COS ti)

(t

ll) f = — (1— COS

Л '

 

(O

Я

n

 

'o

u>„

 

 

 

 

 

Я

 

 

Тогда выражение (9.9) приобретает вид

 

 

У {X, t)

=

J ]

 

^ (1 — cos шя0.

(а)

 

 

/1= 1,

3, 5.

• * •

л

 

 

 

 

 

 

 

В этой формуле общий прогиб получается как сумма прогибов по всем формам колебаний. Каждое слагаемое этой суммы, опре­ деляющее уравнение изгиба оси балки по каждой форме, состоит из двух множителей: дробного и в круглых скобках. Рассмотрим физический смысл каждого из этих множителей. Дробный множи­ тель в (а) представляет собой уравнение изогнутой оси балки при статическом действии соответствующего члена, ряда (9.8) разло­ жения нагрузки. Действительно, дифференциальное уравнение изгиба оси балки при статическом действии какой-то распределен­ ной нагрузки q{x) имеет вид

 

d2 Г Ff d2y(x)

=

Я (X).

 

 

 

 

dx2

dx 2

 

 

 

 

 

 

 

 

 

 

 

Дифференциальное уравнение (4.5) § 13, определяющее формы

колебаний, можно представить в виде

 

 

 

 

 

 

 

p , d2X n (хУ

т (*) ЩгХ„ (х).

 

 

 

dx 2

dx 2

 

 

 

 

 

 

 

 

 

 

Сопоставляя оба уравнения, приходим к выводу, что если на­

грузка q(x) вызывает прогиб у(х), то

точно

так

же

нагрузка

т (х) ш%Хп (х) вызывает

прогиб

Х п(х).

Изменим

в

послед­

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

р„ иы

раз, тогда по закону

сложения действия сил прогиб Х п (х )

изменится в такое лее число

 

\>пХ п (х)

 

 

 

 

 

 

 

раз и будет равен ——"

.

 

 

 

 

 

 

 

Этим

мы показали,

что

статическое

действие

рассуждением

нагрузки

 

 

 

 

 

 

 

 

 

 

-£г И (х) ЩгХп (*)] =

р„т (X) Х п (х),

 

 

 

П

 

 

 

 

 

 

 

 

которая

является одним

из членов

ряда

(9.8),

вызывает

статиче­

ский прогиб балки, описываемый выражением

 

^

, что и

требовалось доказать.

 

 

 

 

 

шл

 

178

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

Пример 15. Рассмотрим действие на

простую

балку постоянного сечения

£ / = const с постоянной массой m = const

внезапно

приложенной равномерно рас­

пределенной нагрузки q. Главные формы, частоты собственных колебаний простой балки, а также коэффициенты разложения нагрузки в ряд по собственным функ­ циям будут такими же, как и в примере 14;

Х п (x) = s i n ^ ± ,

2 n W

шп ~~ ц

Ц_

РП-

пкт '

Подставляем эти значения в формулу (а)

4 ql^ У (х , 0 = —5- ß /

EI

тп

п%х

■(1 — COS iänt).

(б)

л = 1 ,3 , 5, . . .

Получим затем выражение для определения изгибающего момента

 

ш х

 

&2У (х , t)

4qP

sin ■ I

 

М {х, 0 = — ЕІ дх'2

 

rfi (1 — COS U>nt).

(B)

/z = l , 3, 5,

Ряды, содержащиеся в выражениях для у(х, t) и М(х, t), имеют быструю сходимость, и достаточно ограничиться одним членом ряда, чтобы получить точ­ ность, требуемую для практических расчетов. Однако в данном примере при получении максимальных значений прогиба и изгибающего момента во времени и по пролету, как сейчас будет показано, нет необходимости прибегать к усече­ нию рядов, так как в этом случае они суммируются в конечном виде. Заметим,

что выражение в круглых скобках при tM=-X = __ получает максимальное зна­

чение одновременно по всем главным формам. Действительно,

cos и„1 м =

cos rflсо, = cos rfin = — 1 (при п — 1, 3, 5, ...).

 

ші

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

1 —cos &ntu = 2 . С учетом этого обстоятельства эпюры про­

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

 

8 qE

у і

1 . пкх

 

шах у (х) 1ШІ

2j

nbsm ~~Г ’

(r>

 

 

л = і , з ,

5 ____

 

 

шах М (х)

8qll

s

 

nitx

(Д)

тіЗ

n3

sin — .

 

 

 

л = 1, 3, 5, . . .

12*

179

Соседние файлы в папке книги из ГПНТБ