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

книги из ГПНТБ / Попов, Н. Н. Динамический расчет железобетонных конструкций

.pdf
Скачиваний:
81
Добавлен:
20.10.2023
Размер:
8.05 Mб
Скачать

ся при значениях у вплоть до первого максимума у т, если наименьшим по модулю корнем функции ф (у) является

положительное

число, которое

равно

ут. При

проверке

этого условия

сходимости полезно

рассмотреть

функции

из некоторого класса U, которым мы обозначим множество

аналитических

функций и (у),

вещественных на

действи­

тельной оси и обладающих следующими свойствами:

а) функция

и (у) положительна на

некотором отрезке

[О, е] (е > 0) положительной оси,

т. е.

 

 

и (0) >

0, и {у) > 0 при

0 <

у < е;

(2.145)

б) при всех

комплексных у

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

 

и { \ у \ Х \ и ( у ) \ .

 

 

(2.146)

Здесь справа стоит модуль значения функции в некоторой точке у плоскости, слева — значение этой функции в точке действительной оси с координатой, равной модулю числа у.

У функций из класса U среди корней с наименьшими модулями обязательно имеется вещественное положитель­ ное число. Действительно, пусть и (h) = 0. Тогда из (2.146)

следует,-

что и (| h |) ^

0. Если и (| h |) <С 0, то вследствие

условия

(2.145)

и непрерывности и (у) функция и (у) об­

ращается в нуль

при

некотором 1ц, находящемся в интер­

вале (0, I h I). Если и (I h |) = 0, то | Іг| является положитель­ ным корнем и (у). Таким образом, у функций и (у) наимень­ шим по модулю корнем может быть положительное число или несколько комплексных чисел с равными модулями, среди которых обязательно имеется положительное число.

Из общего условия сходимости вытекает следующее: ряды (2.116) и (2.117) сходятся вплоть до первого макси­ мума у т, если ф (у) принадлежит классу U.

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

(2.146):

 

Н \у \)> \ѵ {у)\ .

(2.147)

Отметим для функций из классов U и V ряд свойств, которые позволяют получать конкретные функции из этих классов: 1) произведения функций из U (У) также принадле­ жат классу О (У); 2) сумма любого конечного числа функций

89

из V и сумма сходящегося ряда функций из V принадлежит V; 3) интеграл от ѵ (у) из класса V вида

у

 

w (у) = j V (у) dy

(2.148)

принадлежит У;

4) функции

 

Ui(y) = u(y) — v(y)

(и(у)^ѵ(у) при

0 < г /< е ) ;

 

 

и (У)

 

 

Щ (У) =■ ѵ(у)

 

принадлежат классу

U.

 

 

Используя эти свойства, установим для примера один

вид функций из класса

U. Вначале отметим, что уп (п = 1,

2, ...), и положительные числа а принадлежат одновремен­

но

классам

U и V, так

как \уп \ = \ у\п и а = \а\

при

а >

0. Из

1-го и 2-го

свойств вытекает, что сумма

схо­

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

.ства

получим,

что функция

 

 

 

 

 

СО

 

 

 

'

 

О— S

ßml/'",

(2.149)

где

 

m= 1

 

 

ßo >

0; ßm > 0

 

 

 

 

принадлежит U. Отсюда следует, что при нулевой началь­

ной

скорости,

т. е. при ф0 = 0,

сходимость

рядов (2.116)

и (2.117) вплоть до Ут обеспечивается, если

 

 

 

Ф о>0,

ф(т>< 0

( m ^ 2).

(2.150)

Эти условия выполняются, например, если функция / (t) задана в виде (2.137) и коэффициенты ак из (2.107) таковы, что

и

С і> 0 ,

ak ^ 0

(А > 2 )

 

(2.151)

 

 

 

 

 

 

 

 

 

ß2=

<

2ai*

 

(2-152)

Действительно, фо =

У>

0.

Из (2.133)

< 0 и из (2.135)

и (2.151) фо <

0. Найдем

из (2.133) и

(2.135)

величину

ß2 +

Ч о =

ß2---- ß2

а1

ß2

а1-

- -

90

Тогда на основании (2.152) получим

 

 

 

 

 

 

 

ß«+ ^ < 0 .

 

 

 

 

Из (2.134) следует,

что

П'о < О,

и

из

(2.135)

и (2.151) —

ф"' <

0.

Также при

всех

т ^

2 " ф

< 0.

 

 

Требуемая сходимость

рядов

нарушается

при

малой

3 начальной скорости.

Это вызвано тем, что при малом зна­

чении ф0 функция ф (у)

имеет отрицательный

корень г/_,

для

которого I г/_ I <

у т.

 

 

изложенного

метода

При

практическом

применении

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

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

Пример 1. Определим прогиб гибкой балки-полоски, вызываемый динамической нагрузкой р (t) = р (і — g). Кон­ цы балки предполагаем шарнирно-опертыми и несмещае-

мыми. Тогда из (2.75) и (2.83)

£ = 1 и а =

2,73. Движение

конструкции описывается

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

(2.81) при

 

 

 

1 ------------ \

=

Ѳ =

с о * Ѳ .

со*Ѳ

 

0

 

Пусть интенсивность р нагрузкитакова, что у = 10. Определим предварительно безразмерный прогиб пластин­ ки при достаточно длительной нагрузке, когда можно при­

нять Ѳ = оо. Интегрируя уравнение (2.81) при / ((*) = 1, находим

9*)=т Ш !=та*~0'5!'-

= 1 0 ^ - 0 ,5 ^ - 0 ,6 8 2 ^ .

Максимальный прогиб пластинки определим из урав-

нения

10 — 0,5(/* — 0,682уі = 0.

Решая это уравнение, находим наименьший положитель­ ный корень г/*т = 2,35. Для решения уравнения (2.81)

91

при / (/„.) = 1 — применим метод степенных рядов (2.116)

и (2.117) при замене у на у* и t на t%. Функцию / (/*) выби­

раем

в

виде

(2.137),

т. е. / (/*) =

 

^l +

cos— ^ 'j , где

0* =

Ѳ.

Пусть

время

действия 0

 

нагрузки

таково, что

0 =

0* =

 

 

тогда

ß = 2.

 

 

 

 

 

 

 

По формулам (2.133)—(2.135) находим коэффициенты

рядов (2.116)

и

(2.117)

при у = 10;

ß = 2; аг — 1;

а3 =

= 2,73;

ah =

0 {k ф

1

и 3) и Ъ =

 

фо =

10; ÜÖ =

—0,4;

фо =

Ю-0,5 (—0,4) — 1 = —3; QJJ =

0,0133; ф" =

ybQ'ö =

= 10 • 0,5 • 0,0133 =

0,0665; Q" =

—0,48;

ф>/ =

ybQ'"

— 6а3 =

 

5 (—0,48) — 6-2,73 = —18,79;

Й‘ѵ =

—0,265.

Ограничившись в

рядах (2.116)

и

(2.117)

членами до

т = 4,

будем иметь

 

 

 

 

 

 

 

 

 

 

й Ы

=

1 - 0 ,4 0 , + 0,0066^ -

0,08^ -0,011 у * ,;

 

 

 

 

Ф Ы

= 1Оу*-

1,50“ + 0,01 \у \- 0,782у\.

 

 

При у* =

2,07 получаем ф = 0 и Q =

—0,72. Так как

Q >

—1,

то

максимальный безразмерный прогиб

дости­

гается при

t* <

0*

и у...т — 2,07.

Время

его достижения

находим

из

равенства

(2.140):

 

 

 

 

 

 

 

 

 

**

arc cos П (г/Фт) ■= -2’37 0* =

1,185.

 

 

 

 

 

 

 

31

 

 

 

31

 

 

 

 

 

Максимальный

прогиб

пластинки у т =

2,076.

 

 

 

Пример 2. Рассмотрим шарнирно-опертую балку из ма­ териала с диаграммой деформирования вида (1.24). Опреде­ лим максимальный прогиб балки, вызываемый динамиче­ ской нагрузкой

р ( 0 = р ( ! - ф ) .

Прогиб найдем из решения дифференциального уравнения (2.96). Представим это уравнение в безразмерном виде

 

~

^

+ y*— ayl=yf(t.t),

 

(2.153)

где

JL -

_

Зи*

В3 .

4рі3

,

У*

/ ’

 

 

4

Bj.1*

У ’

п &Вх

 

 

 

 

/( /,) =

1 - 4

;

0 = ^ 0 .

 

 

92

Получим решение уравнения (2.153) в виде рядов (2.116) и (2.117), причем функцию f (/*) примем в виде

где Ѳ*=Ѳ. Пусть характеристика материала балки и вели­ чина нагрузки таковы, что а — 50, у = 0,04. Примем время действия нагрузки равным 0* = л, тогда из (2.136) ß = 1.

По формулам (2.133)—(2.135) находим:

ф»ѵ = 0,02 (—2860) — 6 (— 50) = 242,8; П'ѵ = — 110 000.

Оставив в рядах (2.116) и (2.117) только члены до т = 4, получим

Q (t/*) = 1 - 25t/.,.- Ь2уІ - Ш у і - 4580^;

Ф (г/*) = 0,04г/*-0,75г/*-0,347г/°-10,12^.

Решая уравнение ф (у*) = 0, найдем у^т = 0,055. При этом Q (0,055) = 1 — 1,373 — 0,157 — 0,079 — 0,042 = = —0,651. Время достижения максимального прогиба

2 28 Ѳ* = 2,28. Максимальный прогиб балки у т =

= 0,055/.

Г л а в а 3

РАСЧЕТ

 

БАЛОЧНЫХ

 

ЖЕЛЕЗОБЕТОННЫХ

 

КОНСТРУКЦИЙ

§ 10. РЕШЕНИЕ УРАВНЕНИЯ ДВИЖЕНИЯ БАЛКИ В УПРУГОЙ СТАДИИ

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

Рассмотрим движение балки в упругой стадии. В за­ висимости от исходных предпосылок в динамике сооружений находят применение различные дифференциальные урав­ нения колебания балки. Приведем вначале уравнение, считающееся наиболее точным в рамках теории сопро­ тивления материалов. При его выводе учитываются дефор-’ мации сдвига и инерция вращения элементов балки [71]. Кроме того, на колебания реальной конструкции сущест­ венно влияют силы сопротивления, приводящие к постепен­ ному затуханию свободных колебаний [4, 47, 57]. Обычно среди этих сопротивлений существенное значение имеют силы внутреннего неупругого сопротивления, вызывающие так называемые гистерезисные потери энергии деформации. Подобные силы сопротивления проявляются главным обра­ зом при наличии циклических процессов нагрузки и раз­ грузки в материале конструкции. Поэтому они слабо влияют на деформирование конструкции в период ее движения от начала приложения динамической нагрузки до достижения

94

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

Обозначим через (х , f) прогиб балки только от изгиба,

а через

у 2 (х , t) — только от сдвига. Тогда полный прогиб

запишем

так:

 

 

У = Уі + Уг-

(3.1)

Для усилий в балке справедливы следующие зависимости:

 

М

=

-

В

(3.2)

 

 

д х 2

 

 

 

Q^-k' GF

.

 

(3.3)

 

 

дх

 

 

'

Здесь М (X, t) — изгибающий момент; Q (х,

f) — попереч­

ная

сила; В — жесткость

балки;

F — площадь

попереч­

ного

сечения; k’ — коэффициент,

учитывающий

неравно­

мерность распределения касательных напряжений по попе­ речному сечению балки. Значения коэффициента k! для ряда поперечных сечений приведены в книге [73].

В исследованиях В. И. Мурашева [39] показано, что при определении прогибов железобетонных балок, когда действуют статические нагрузки, жесткость балки можно считать постоянной по ее длине и равной жесткости в месте максимального изгибающего момента. Будем принимать справедливым это положение и при расчетах на действие динамической нагрузки. Величина жесткости определяется существующими методами с учетом раскрытия трещин в бетоне растянутой зоны.

Уравнения динамического равновесия балки имеют вид

дМ

Q~\~Jp

а*Уі = 0;

(3.4)

дх

 

д Р д х

 

dQ

pF

д2 у

Р =о,

(3.5)

дх

d t2

где J — момент инерции сечения; р — плотность материала балки.

Системы (3.2)—(3.5) в результате преобразований при­ водятся к одному уравнению, носящему название уравнения С. П. Тимошенко:

В

д* у .

т

д2 у

 

В

\

ді у

2 -

dJjL

——+

d t2

J +

 

 

■ p

dt4

о

д х

4

 

k'р G J дх2 d t2

J -k'G

 

 

= P +

Р

d2 p

_

В

d2 p

(3.6)

 

 

 

k ' FG

H F

k ' GF ' d x 2

 

95

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

дМ

(3.7)

Q

дх

Врезультате исключения уг и г/2 получим уравнение

 

d t 2

-----В

 

----- (3.8)

д х 2

'

дх*

k' G

дх°- d t 1

li'G F

Для большинства балочных конструкций, у которых отношение высоты к пролету не превышает Ѵ5, не сущест­ венны также и деформации сдвига. В этом случае у2 = О, перерезывающая сила определяется из уравнения равнове­ сия (3.7) и уравнение колебаний балки имеет вид

<3-э>

Уравнения движения (3.6) и (3.8) являются волновыми: произвольное возмущение распространяется вдоль балки

согласно уравнению (3.6) со скоростями сх = у ~ и с2=

а согласно уравнению (3.8) — со скоростью с2.

Уравнение (3.9) неволновое, произвольное возмущение

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

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

Кратко изложим метод нахождения решения уравнения (3.9), обычно применяемый при исследовании колебатель­ ных процессов. В этом методе, называемом методом Фурье, любой колебательный процесс представляется в виде сум­ мы простейших решений уравнения (3.9). Эти решения имеют вид

У(х, t)= Т it)X (х),

(3.10)

т. е. являются произведением функций, одна из которых

зависит

только

от

времени

t,

а другая — только от х.

Функция

X (X) называется

формой свободных

колебаний

. или собственной

функцией.

 

 

 

Любое решение уравнения (3.9) представляется в виде

бесконечного ряда:

 

 

 

 

 

 

 

 

П— со

 

 

 

У(х >t) =

и

Tn (t)Xn(x).

(3.11)

 

 

 

/1=1

 

 

Собственные функции находятся из решений однород­

ного уравнения

движения

 

балки:

 

 

 

 

B JL L + т д^ - = 0.

(3.12)

 

 

 

дх*

 

 

dt2

 

Подставляя (3.10) в (3.12) и разделяя переменные, полу­

чаем два уравнения для функций Т (t) и X (х):

 

 

 

 

Г +

со2 7 = 0;

(3.13)

 

 

 

ХІѴ — №Х = 0,

(3.14)

где со — частота

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

связанная

с собственным числом X зависимостью

 

 

 

 

со =

Л2 л / ~ J L t

 

 

 

 

 

 

V

m

 

Общие решения

уравнений

(3.13) и (3.14)

имеют вид

 

Т (t) =

а sin (£>t + b cos со/;

(3.15)

X (x) =

А sin Kx +

В cos Хх С sh Хх +

 

 

 

+ DchXx.

(3.16)

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

Собственные функции при обычных типах опорных за­ креплений балок (шарнирные, защемленные и др.) обладают двумя важными свойствами [33]:

1) свойством ортогональности, т. е.

I

 

j Х п (х) Xm (х) dx = 0 при п =ё-т,

(3.17)-

о

 

где I — пролет балки;

4 Н. Н. Попов, ІЗ. С. Расторгуев

97

2)

свойством полноты,

т.

е. произвольная функция

Ф (х),

которая

является

формой прогибов балки при тех

же граничных

условиях, что

и собственные функции, мо­

жет быть разложена в

абсолютно и равномерно сходящий­

ся ряд по собственным

функциям,

при этом

 

 

 

 

СО

ahX h{x),

 

 

 

ф (*)=

 

2

(3.18)

где

 

 

к=1

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

I

Ф (*) Х к (X) dx

 

 

=

о

 

 

 

 

 

]----------------

 

 

[ Х% (x)dx b

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

2

[ ^ ( 0 + “ ^ ( 0 ] * п М =

Р (X, t)

т

п—I

 

Умножив левую и правую части полученного равенства на X h (X), проинтегрировав по пролету балки и приняв во вни­ мание (3.17), найдем уравнение для определения функций

Tk(t):

Tk{ t) + m T h{t) = Ph ü)

где

] p { x , t ) X h (x)dx

Ph (0 = -

f Xk (x) dx

'o

Решение этого уравнения имеет вид

Th(f)=*Akcosa>k t + Bk sin coht + t -

Г j*/?А(тг)sin caÄ(/—x) dx.

(3.19)

mwh

Произвольные постоянные Ak и Bk находятся из на­ чальных условий. В том случае, когда действующая на бал­ ку нагрузка р (х, г1) изменяется во времени по линейному

98

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