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

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

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

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

Добавлен: 23.10.2024

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

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

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

 

п—т

 

 

хг = ш.

п оператор F может

быть выражен в виде

отрицатель­

ной, а

оператор С — положительной рациональной

степени опе­

ратора

D.

колебаний однородного стержня

Для

сдвиговых поперечных

можно написать дифференциальные уравнения с учетом затуха­ ния, аналогичные (II. 57) и (II. 60):

д”-у

с Y n t â F

 

д 2у

- G'F

д:у

 

0;

 

(11.76')

т d/2

дх -ді

дх-

=

 

+ с V m G ’F ^ - G ' F

дх?

=

0.

 

 

(11.76")

дг-

 

 

 

 

 

 

 

 

 

 

 

Через G' обозначен модуль сдвига с учетом

коэффициента фор­

мы поперечного сечения, площадь которого равна F. Для

(11.76')

метод разделения (переменных приводит к уравнениям

 

 

 

 

г" +

 

ßr =

0;

 

 

 

 

 

 

 

 

ср -ф CCüßcp +

(іГ® = 0.

 

 

 

Решением первого из них являются собственные функции

уравнения незатухающих колебаний, а решение

второго —

 

 

С«ц

ßf

 

 

 

 

 

 

 

 

®l {t) = al e

2

 

sin

t — azj,

 

 

 

где

 

 

 

 

 

саза

 

 

 

 

 

 

 

°'i = 0,i

1'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

G'F

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как величины

ß.

пропорциональны

частотам

,

решение

 

/

 

 

уравнений (11.76)

качественно

 

аналогично

решению

уравне-

ния (11.57).

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (11.76") эквивалентно следующим двум:

 

 

 

 

г" +

3V = 0;

 

О ,

 

 

 

 

 

 

Ш

 

 

о

-

 

 

 

 

 

 

ѵО-f- С Н ‘ со - } - o r со =

 

 

 

 

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

-

І І І - 1

с?l ( t ) = a i e

2fi sin (ш' t — аг j,

82


где

шI

 

1 -

С3

 

 

2 *

 

 

 

I

Декремен'г колебаний ог =

с у - уменьшается с увеличением

номера гармоники. Из двух

уравнений

(11.76) затухающих коле­

1/

 

 

 

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

Чтобы улучшить соответствие решений уравнения с опытными

данными, в данном случае, так же как и в уравнении

(11.61), можно

использовать линейные комбинации

диссипативных

членов урав­

нений (11.76) с двумя различными

коэффициентами затухания

С\ и с2.

 

 

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

Нх) сходный с оператором затухания, предложенным В. Колоуше-

ком [44], с той разницей, что в данном случае оператор действует на функцию переменной х — координаты стержня, а не на функ­ цию времени t, как в работе [44].

Оператор Гильберта Н определяется следующим соотношением:

= 4-J іт_х d\.

(11.77)

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

Н sin ßx = — cos ßx

(11.78)

H cos ßx — sin ßx

Уравнение

_ G 'F —

(11.79)

u

dx?

 

удовлетворяет поставленным требованиям. В этом можно убедить­ ся, используя подстановку

со

 

 

 

У (х , t) = 2

(Л, sin ß, X +

Bt cos ß, x) <?l (t),

(11.80)

1

 

 

 

что приводит к уравнению для <р((£)

 

 

?i +

7ші <?i + w\ t

= °*

 

83


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

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

переписав уравнение (11.79) в виде

 

 

 

<Э- у -р GTF

н

д-у

-

G ' F ^ = 0.

(11.81)

т

 

дхді

 

д х 1

 

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

Г л а в а III

ИССЛЕДОВАНИЕ СЕЙСМИЧЕСКИХ ВОЗДЕЙСТВИЙ МЕТОДАМИ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ И НАТУРНЫХ ИСПЫТАНИЙ

§ 1. Определение реакций сооружений на сейсмические воздействия с помощью весовых функций

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

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

Lx = ш0(t ) Ж,

(III.1)

где L—дифференциально-матричный оператор;

X—вектор обобщенных координат, которыми в большинстве случаев являются горизонтальные перемещения точек, несущих сосредоточенные массы в расчетной схеме;' Ж—матрица сосредоточенных масс;

ш0(^)—вектор-функция ускорения основания.

Число строк матрицы L равно числу степеней свободы системы. Решение системы (III. 1) получается в -главных координатах, при этом разделение переменных достигается только при опреде­ ленных дополнительных предположениях о характеристиках рас­ сеяния энергии ів системе. В результате каждая координата х к

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

рых при нулевых начальных условиях зависят от wo{t)

и закона

внутреннего трения:

 

П

 

x k = ^ l x lkTl (t,B),

(III.2)

г-і

85

 


где е — параметр, характеризующий внутреннее трение.

Для расчета сооружений необходимо определить максималь­

ные значения

координат и связанных с ними силовых факторов

за время землетрясения, которое принимается равным

длине

обрабатываемого интервала записи на акселерограмме.

этапа"

Решение

задачи обычно выполняется на ЭЦВМ в два

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

Второй этап решения заключается в определении максималь­ ных значений суммы (ІП.2). Акселерограммы обычно задаются в табличной форме и для сокращения вычислений при подготовке материала и определении значений х к max интервал обработки

акселерограммы сокращают до 4э-б сек. Максимальные значения суммы вида (III.2) приходится определять отдельно для всех необходимых кинематических компонент и силовых факторов, при­ чем структура формулы (ІП.2) каждый раз меняется.

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

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

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

Переходной функцией g(t)

мы будем называть реакцию си­

стемы на воздействие в виде

единичной

ступенчатой функции

о(7); определяемой равенствами

 

 

o(t) =

1, t > и,

 

о ( * ) = 0 , * < 0 ,

86


а весовой функцией h(t) — реакцию на воздействие единичной импульсной функции 8(t), обладающей свойствами

8(0 = 0,

О,

j 8 (t) dt = 1.

Если L — оператор системы, то

gj(t) = L - l o.(t)

(ІИ-3)

J j ( t ) = L - 18.(^)

где gj {t) — столбец переходных функций, а lij (t) — столбец весо­

вых функций, соответствующих воздействию на /-м входе, при ра­ венстве нулю воздействия на остальных входах [43, 76, 86]:

 

1, (0 = [e*j (0);

 

 

 

 

 

hj(t) = [hk j (t)],

Ä = l,?, ... ,

n.

 

 

Линейная

система с т входами и ѣ выходами

может

быть

полностью описана матрицей

переходных

функций

нли

весовых функций || hkj |J

, причем

 

 

 

ë kj (0

= hkt (0

(k = 1,2, ... , ѣ\ j =

1,2, ... ,

т).

 

Как известно, реакция на k-u выходе при воздействии w0j (х) на у'-м входе определяется выражением

x kj (0 - J woj (') ё кі (t

і) dr = j w0j (x) hkj (t — x) dz.

(III.4)

U

U

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

няется вектором с

составляющими gk{t)

или соответственно

hk(t), где

т

 

 

 

ё к (*) = У>£ы

 

 

j - 1

(III.5)

 

 

м * )

= 2 А*у{t)' k = h 2

п

 

/=•1

 

Здесь движение всех точек основания сооружения рассматрива­ ется как синхронное. Весовыми функциями можно пользоваться

87