Файл: Д’Анжело, Г. Линейные системы с переменными параметрами. Анализ и синтез-1.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 71

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

где Wi—диагональная матрица квадратов собственных частот,

R- 1CR = M~RCR =

Система уравнений сейсмических колебаний будет иметь вид:

43,67 л'і + 53,24 л 2 + 27,28 х ь + 54,53 х ±+ 25,50 л2 +

+ 13,02л'з + ІО4 x L= — 150,15 w0 (t),

53,24 + 107,18 л'., + 58,82 Зс3 + 28,99 jq + 76,63 х 3 +

+ 33,65 ^з+ Ю 1X-2 = - 275,21 w0 (t),

53,24Ä-j+ 114,80л', -f 85,84^3 + 22,33 x, + 53,07 x ,+

-f 79,40 хз+104х 3 = — 171,93 IWQ (0-

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

§7. Системы с распределенны ми параметрам и

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

без

затухания.

 

 

 

 

 

 

 

ня

Уравнение свободных изгибных колебаний однородного стерж­

без затухания имеет вид

 

 

 

 

 

 

 

El

дку -j- т

д^у

=

0,

(11.49)

 

 

 

дх1

dt-

 

 

 

 

где

Е — модуль

упругости

материала,

/ — момент инерции

попе­

речного сечения, т — масса единицы длины.

разде­

Подстановка

в (II.49) у (х, t) = X

(х) Т (t) приводит к

лению переменных

 

 

 

 

 

 

 

 

T(t) + u T { t )

=

0,

(11.50)

 

 

 

- Р Х

( X )

=

0,

(11.51)

где

 

 

 

 

 

 

 

 

m ш2

"ЕЛ '

73.


Уравнение

(11.50) соответствует гармоническим

колебаниям

-с частотой

ш.:

 

 

 

 

 

 

 

 

 

T(t) = a, sin (ш ^ -о ,).

 

 

(11.52)

Решениями уравнения (11.50) являются балочные

функции

Xj (х) = At sin ß, ■* +

Bt cos ^i x + Cl sh ß, л: + Dt ch ß, Л'.

(11.53)

Собственные значения

ß(

и постоянные

At ,B., С., D{ опреде­

ляются по заданной системе граничных условий.

 

Сезава

Уравнение затухающих

колебаний было составлено

(К. Sasawa)

в

1927 г.*

[147]. Вывод уравнений приведен

в более

известной у нас работе

[148]. В основу положена гипотеза вяз­

кого трения, согласно которой нормальная

вязкость

пропорцио­

нальна скорости деформации:

 

 

 

 

 

 

 

dv{t) ^

 

 

(11.54)

где dv (t ) — сила сопротивления, приходящаяся на элементарную

площадку

поперечного сечения; г — продольная де­

формация;

; — коэффициент вязкости.

Изгибающий момент сил сопротивления, действующих на эле­ мент площади поперечного сечения, будет

dp(t) = ^ z d F ,

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

чощий момент:

М г

£ = £7-

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

dp {t) = -ßT~dTzdF.

Момент сил сопротивления в поперечном сечении равен

дМ

дМ.

Р EI dt

Е dt

сотому моменту

соответствует

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

длины стержня:

 

 

 

 

 

 

do =

дз

g

âM

dx.

 

 

d x J

E

dt

 

* В известной книге С. П. Тимошенко «Колебания в инженерном деле», Физматгиз, 1959 г., как первоисточник указывается работа X. Хольцера, опублико­ ванная в 1928 г.

74


Подставляя в последнее равенство М = - £ / ^ и сокращая на

■dx, получаем силу сопротивления в сечении с координатой х, отнесенную к единице длины стержня:

=

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

+ 5/

д:>х

■ Я /|^ = 0.

dt-

 

dxldt

дхі

(11.55)

эту силу после чего

(11.56)

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

 

 

ч

 

 

 

 

V =

С

 

 

 

 

 

 

 

и записать уравнение затухающих колебаний как

 

ді- + ѴЕІ

дъу

+ Ш

д1у

= 0.

(11.57)

dxlât

дх*

Из уравнения (II. 57) ясно,

что принятые

авторами предпосылки

приводят к пропорциональной зависимости сил сопротивления от поперечной жесткости стержня.

Разделяя переменные подстановкой у(х, t)=X(x)T(t), полу­ чаем

E I ^ ( T ( t ) + Z ' f ( t ) ) + m X T ( t ) = 0,

■что равносильно двум уравнениям:

d*X

ЕІ dx*

t{t) + l a 2t(t) + a ? T { t ) = 0 .

Первое уравнение совпадает с уравнением (II. 51), поэтому собст­ венные функции краевой задачи (II. 56) совпадают с таковыми для уравнения незатухающих колебаний. Второе уравнение имеет

DeuieiHHe

Tt (t ) = ate 2

sin ( ш. t — a.j,

где

(11.58)

7o


Декременты колебаний —

9

 

3.1= ігГ -и)І. - .

(11.59)

I

 

Выражения (II. 58) и (II. 59), где со,-—частота

незатухающих

колебаний, совпадают с выражениями (II. 12) и

(II. 13), кото­

рыми определяются частоты и декременты затухающих колеба­ ний дискретной системы, если принять матрицу рассеяния про­ порциональной матрице жесткости. Это совпадение вполне зако­ номерно, так как уравнения (II. 10) и (II. 56) основаны на одной и той же предпосылке о пропорциональности сил сопротивления поперечной жесткости системы. Выше было показано, что реше­ ния, получаемые на основе этой предпосылки, приводят к логи­ ческим противоречиям и не согласуются с результатами опытов. Остановимся еще на выводах, к которым приходит автор ра­ боты [148].

Исходя из выполненных еще в 1908 г. опытов Г. Оморп, в ко­ торых определен логарифмический декремент первой формы ко­

лебаний кирпичной колонны и получено 7^ = 0,26

сек.,

öi = 0,148,

по формулам (11.58), (11.59) можно установить,

что

наивысшая

действительная частота колебаний не превышает

80

в секунду.

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

Для зданий и сооружений более сложной конструкции, чем кирпичная колонна, декременты первых форм колебаний обычно находятся в пределах 0,2ч-0,3, поэтому для них колебания по высшим формам тем более были бы невозможны.

Работы К. Сюэхиро (К. Suyehiro [149]) по инженерной сейсмо­ логии серьезно повлияли на развитие исследований по сейсмо­ стойкости в СССР, и некоторые специалисты до последнего вре­ мени считали, что колебания сооружений по высшим формам не­ возможны. В наши дни это представление, по-видимому, уже пе­ рестало существовать, главным образом под влиянием выводов гистерезисной теории затухания и многочисленных экспериментов, показавших, что формула (11.58) не отражает действительного положения вещей.

Чтобы к стержню с непрерывно распределенной массой приме­ нить методику, разработанную для дискретных систем, следует отказаться от исследования сил сопротивления, действующих на элементарные площадки поперечного сечения и пропорциональ­ ных скорости изменения нормальных деформаций. Вместо этого будем считать, что силы сопротивления, действующие на элемент длины стержня, пропорциональны скорости его перемещений (воз­ можность такой точки зрения отмечена в работе [120], но в ней воп­ рос рассматривается в другой постановке). Кроме того, примем,

76


что эти силы зависят не только от физических свойств материала, но и от механических характеристик сооружения.

Простейший вид

такой

зависимости — пропорциональность

сил трения массе элемента

длины —приводит к уравнению

 

д*у

em _ду_

■EI ^

 

(11.60)

 

de­

öt

Ол■'

 

 

Применяя метод

разделения

переменных, приходим к двум

уравнениям:

 

 

 

 

 

 

 

д4х

 

 

 

 

 

 

d x {

 

 

 

 

 

причем

Т

c f

 

ш2Т = 0,

 

 

 

 

 

 

 

 

 

 

Р4

/ Я Ш 3

 

 

 

 

~ЕГ '

 

 

Первое уравнение совпадает

с

(II. 51),

второе

аналогично

уравнению (II. 7). Следовательно,

формы

колебаний,

соответст­

вующие уравнению

(II. 60), одинаковы с

формами

колебаний

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

затухающих

колебаний

составят

Декременты

колебаний

 

 

 

с~

убывают обратно пропорционально частоте.

Решение уравнения

(II. 60) обладает теми же свойствами, что

и решение системы (II.

3). Оно находится в лучшем соответствии

■с опытом и, по-видимому, ведет к некоторой переоценке влияния

высших форм

колебаний. В

большинстве случаев это

лучше,

чем недооценка, которая связана с решением уравнения

(11.56).

Учет затухания

по уравнению

(II. 60) предпочтительнее

еще и

потому, что решение не содержит очевидных противоречий с фак­

тами и допускает

представление

затухающих колебаний

в

виде

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

(II.

56)

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

существует.

В

выражении

 

 

 

 

оо

 

 

 

 

У(Л'Д)=

6=1

 

 

 

77