Файл: Морозов, В. А. Регулярные методы решения некорректно поставленных задач-1.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

в колебательных системах, — необратимое поглощение (рассеяние) энергии колебаний. Детальное рассмотрение закономерностей это­ го сложного явления и методов его описания не входят в наши за­ дачи [101, 137]. Некоторые данные о рассеянии энергии в мостах приведены в гл. III, здесь же мы только изложим способ учета это­ го фактора в уравнениях колебаний.

В зарубежных работах по теории сейсмостойкости рассеяние энергии обычно учитывается по теории вязкого сопротивления (тео­ рия Фохта). С этой целью в дифференциальные уравнения сейсми­ ческих колебаний вводятся дополнительные диссипативные силы, пропорциональные скоростям колебаний и направленные в обрат­ ную сторону [65]; количественной характеристикой меры рассеяния энергии служат коэффициенты пропорциональности, именуемые ко­ эффициентами вязкого сопротивления. Теория Фохта имеет ряд недостатков [137]- В СССР за последние годы наибольшее распро­ странение получила теория Е. С. Сорокина, основанная на рас­ смотрении неупругой работы материала. Мы используем первона­ чальный, более простой вариант этой теории, где в качестве меры рассеяния энергии принят коэффициент неупругого сопротивления у, смысл которого будет пояснен ниже. Для учета рассеяния энер­ гии по теории Сорокина в дифференциальных уравнениях колеба­ ний все смещения и силы представляются комплексными функция­ ми и к силам вводится дополнительный комплексный множитель

(1 + Y/)> гДе / — мнимая единица [137].

Учет рассеяния энергии по обоим вышеуказанным методам в задаче сейсмических колебаний приводит к тому, что в подынте­ гральных выражениях решений (II.4), (И.10) возникают дополни­ тельные экспоненциальные множители. Мы будем придерживаться

в дальнейшем теории Е. С. Сорокина,

которая

физически-

более

обоснована. На основе этой теории решения (11.4), (НЛО)

задачи

сейсмических колебаний, дополненные

за

счет

учета

рассеяния

энергии, могут быть представлены в таком виде1.

 

 

Для дискретной расчетной схемы

 

 

 

 

 

\ ? ^ х)е Т‘ (‘ ' ) s i n ^ - ( t - T ) d x

( 11. 12)

о

 

 

 

 

( £ = 1 , 2 ,.. .„ п).

 

 

 

 

Для схемы с распределенными параметрами

 

 

 

^ L - T iX i (x)^V0(x)e

HLТ;

(/-■О

 

—x)dt.

 

sin

о

7'/

(11.13)

 

1 Порядок составления и решения дифференциальных уравнений сейсмических колебаний по теории Е. С. Сорокина мы здесь не приводим, так как этот вопрос детально изложен в работе [46].

44


Нужно отметить, что учет рассеяния энергии по теории Соро­ кина (так же, как и по теории Фохта) не отражается на формах собственных колебаний, т. е. на коэффициентах X,-;t и функциях Х{(х). При обычных значениях коэффициента, у, характерных для строительных конструкций, периоды (частоты) колебаний также практически не меняются, поэтому все эти величины, а также ко­ эффициенты Di в формулах (11.12), (11.13) имеют тот же смысл, что и раньше.

Выражения (11.12), (11.13) полностью определяют закон сейс­ мических колебаний сооружения при заданном законе движения основания. Как видим, структура этих выражений в обоМх случаях одинакова: смещения точек сооружения представлены в виде сово­ купности однотипных членов; при схеме с распределенными пара­ метрами число этих членов бесконечно и они образуют сходящий­ ся ряд. При дискретной схеме имеем конечную сумму членов, число которых равно числу степеней свободы системы. Каждый член в формулах (11.12), (IL13) описывает так называемое нормальное колебание (нормальную составляющую колебаний), п.ри котором все точки сооружения смещаются синхронно и форма колебания не меняется во времени. Для г'-го нормального колебания эта фор­

ма описывается собственной функцией Xi(x)

в случае континуаль­

ной схемы

и

совокупностью амплитудных

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

(6=1, 2, ..., п)

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

ординаты

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

прикрепления

масс (см.

рис.

11.4). Закон изменения1 смещений во

времени при

i-м колебании

(общий для всех точек) описывается выражением,

носящим название интеграла Дюамеля:,

 

 

(11.14)

О

Относительная роль отдельных нормальных составляющих в сумме (11.12) или ряде (11.13) определяется соотношением произ­ ведений DiTi и интегралов h ( t ) . Как правило, основное значение имеют первые несколько членов.

Таким образом, формулы (II.12), (II.13) представляют смеще­ ния точек сооружения в процессе сейсмических колебаний в виде разложения (ряда или суммы) по нормальным составляющим. В та­ ком же виде представляются и другие факторы сейсмических коле­ баний. Приводим без вывода выражения для сил инерции [46]. При дискретной схеме сосредоточенные силы инерции, приложенные в точках 6, выражаются формулой

ип

2 n D ^ /i(t)

(II. 15)

 

(6 = 1 , 2,. .. , п).

45


 

 

При континуальной схеме для интен­

 

 

сивности

 

распределенных

сил инерции

 

 

имеем:

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

s(x,

t ) =

S, {х,

0 =

 

 

 

 

 

 

 

/Si

 

 

 

 

= - V

 

т (Л-) 2jlDi* ‘ {х)

/, (0-

(Н.16)

 

 

/Si

 

 

 

 

 

 

Здесь

 

 

Sj — сосредоточенные и рас­

ма осциллятора

пределенные

сейсмические

силы

по от­

 

 

дельным нормальным колебаниям.

 

Представляет интерес частный

случай

дискретной системы —

линейный

осциллятор,

имеющий

одну

сосредоточенную

массу

(система

с одной степенью свободы,

рис.

II.3). Обозначим его мас­

су

и период колебаний через

т и Т. В этом случае согласно форму­

ле

(II.5)

можно принять Ь =

1 и Z u = l. Тогда на основе формулы

(11.15)

сейсмическая сила

осциллятора запишется

таким об­

разом:

 

 

 

 

 

S(t) =

- m ^ y I{t).

(11.17)

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

/

т-

,

W ( t , T , y ) = ^ \ Y 0(x)e

т

X

о

 

 

Xsin - у - (t — x)dx.

 

(11.18)

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

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

Укажем принципиальный путь определения и некоторые свой­ ства этих факторов.

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

46


п г 1о 11

m 2&i2.

■<Пп\ п

 

 

^ ’1^21

/7?2^2^ •

•m А л

[m/A,]*" =

 

(11.19)

- т А г

ш 28л 2 .

• ' « А л -

Амплитудные коэффициенты Хщ, описывающие формы собст­ венных колебаний, являются* координатами собственных векторов этой матрицы; круговые частоты колебании (без учета рассеяния энергии) определяются ее собственными значениями Xf.

<Р? =

- Т = .

(‘ = 1. 2 , . . . ,

п)

(11.20)

 

V h

 

 

 

Круговые частоты с

учетом

рассеяния

энергии

выражаются

формулой

 

 

 

 

Для строительных конструкций значения у не превосходят 0,20. Поэтому учет .рассеяния энергии практически не изменяет частот

и можно принять <Р/= <Р/; сообразно с этим для периодов колебаний будем иметь

Т 1 = — = 2 л У Т 1.

(11.21)

<Р/

При схемах с распределенными параметрами собственные функции Xi(x), описывающие формы собственных колебаний, полу­ чаются решением (при соответствующих граничных условиях) сле­ дующего обыкновенного дифференциального уравнения:1

1 [А (

х )]-ср 0; т ( л Щ х ) = 0.

(11.22)

Собственные частоты

ср*° здесь определяются из характеристи­

ческого уравнения описанной краевой задачи. Частоты ср*

(с уче­

том рассеяния энергии) и периоды Ti вычисляются по приведенным выше формулам.

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

Отметим, что амплитудные коэффициенты Xik и собственные функции А (х ), как решения однородных задач собственных зна.- чений, определяются с точностью до постоянного множителя. Од­ нако, как нетрудно усмотреть из формул (II.15), (II.5) или (П.16), (11.11), это не отражается на однозначности полученных решений.

47


Рис. 11.4. Формы собственных колебаний:

я — дискретная система; и — система с распределенными параметрами

Представляет интерес следующая зависимость периодов собст­ венных колебаний от массы и жесткости системы: при увеличении массы системы собственные периоды, как правило, возрастают и лишь некоторые из них могут не изменяться. При увеличении-жест­ кости системы собственные периоды, как правило, уменьшаются и лишь некоторые могут не изменяться [6]. В частности, 'период

основного тона

(i= l)

всегда уменьшается, с возрастанием жест­

кости системы.

В этом

смысле период основного тона Т\ часто

принимают в качестве мерила динамической жесткости соору­ жения.

Формы собственных колебаний отличаются друг от друга чис­ лом полуволн (экстремумов, нулевых точек). Наиболее простая форма с минимальным числом полуволн соответствует первой (ос­ новной) нормальной составляющей (г = 1). С ростом номера i чис­ ло полуволн возрастает. На рис. 11.4 для примера показаны первые три формы собственных колебаний дискретной и континуальной систем, изображенных на рис. II.1, 11.2.

Для форм собственных колебаний справедливы условия орто­ гональности, выражающие свойство независимости нормальных ко­ лебаний [44, 46]. Для дискретных и континуальных систем эти усло­ вия соответственно имеют вид:

2

mbXibX )* = ®’ j m { х ) Х [( х ) Х } ( x ) d x = 0 (i Ф j). (11.23)

fe=l

о

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

48

Рис. II.5. Диаграмма дефор­ мирования не вполне упру­ гого тела

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

_ 1 ДП

(11.24)

2л П

Как показывают опыты, величина у не зависит от частоты колебаний; в диапа­ зоне эксплуатационных напряжений она не зависит и от амплитуды последних. Все предыдущие формулы получены- в предположении y=const.

Численные значения у, обусловленные внутренним трением в материале, уста­

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

§ 11.2. ПРАКТИЧЕСКИЕ МЕТОДЫ РАСЧЕТА НА СЕЙСМОСТОЙКОСТЬ. НОРМИРОВАННАЯ МЕТОДИКА ОПРЕДЕЛЕНИЯ СЕЙСМИЧЕСКИХ СИЛ И УСИЛИЙ ПО СПЕКТРАЛЬНОЙ КРИВОЙ

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

та. Функция WQ{t) = Y^{t), представляющая закон изменения ускорений грунта в расчетных формулах, носит, как отмечалось в § 1.3, нестабильный, случайный характер и не поддается удобному аналитическому описанию. Имеющиеся в настоящее время акселе­ рограммы далеко не достаточны для полной характеристики коле­ баний грунта в различных региональных и местных условиях.

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

по максимальным ускорениям грунта:

 

.Sft = max |S k(t)\ = mk max \F0= ^ ^ ? Q * = / C cQ*.

(11.25)

g

 

49