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

книги из ГПНТБ / Расчет конструкций убежищ

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

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

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

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

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

Вдальнейшем будем применять следующие названия для введенных предельных состояний: предельное состояние по полной несущей способности — состояние 1а, предельное

состояние по отсутствию

больших остаточных деформа­

ций — состояние

16.

расчета конструкций в пласти­

Динамические

методы

ческой стадии, основы которого были заложены в работах А. А. Гвоздева 1151 и И. М. Рабиновича 1501, в настоящее вре­ мя получили большое развитие. Применение этих методов к расчету различных конструкций изложено во многих.ра­ ботах [16, 22, 43, 48, 52, 57]. Экспериментальные данные о влиянии скорости нагружения на прочностные свойства материалов приведены в работах [6, 9, 13, 18, 30, 31, 67].

2. РАСЧЕТ КОНСТРУКЦИЙ В УПРУГОЙ СТАДИИ

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

60

Вывод этих уравнений и методы их решения изложены в различных курсах по динамике сооружений [5, 34, 51, 611.

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

Рассмотрим балочные конструкции. Уравнение колеба­ ний балки постоянной жесткости, как известно [61], имеет

вид

.. .....

 

 

 

 

 

 

В

сНw +

т д2 ш>

Р(х,(),

( 1)

 

 

дх4

 

 

 

 

 

где В — изгибная

жесткость

балки; пг — масса балки на

единицу длины;

р (х, t) — нагрузка на единицу

длины;

w (х, t) — прогиб балки.

поперечная

сила определяются

Изгибающий

момент и

из выражений

 

 

 

 

 

 

 

М =

- В ^

;

Q =

- B

^ .

(2)

 

 

дх

 

 

дх3

v ;

Предположим, что динамическую нагрузку можно пред­

ставить в виде

 

 

 

 

 

 

 

 

 

Р (х,

t) = р/,

(x)f {(),

(3)

где р — некоторое фиксированное (часто наибольшее) значение динамической нагрузки; /\ (х), f (t) — функции, характеризующие изменение нагрузки по пролету и во времени.

Для прогиба балки примем выражение

w {х, t) = pF (х)Т ((),

(4)

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

В Л у (х) = fx (х)

(5)

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

61

которую будем трактовать как коэффициент динамичности в упругой стадии.

Подставив (4) в уравнение (1), получим ошибку

L(x,t)

=

р [Т (/) Д (х) + тТ (t)F (х) —

 

 

I

-

Л Ш (О].

 

(6)

 

L (х, t) F (х) dx — 0

 

 

Из условия

f

находим

уравнение

 

6

 

 

 

 

для функции Т (t):

 

 

 

 

 

т (0 +

со2Г (0 = со2/

(0,

(7)

 

 

 

i

 

 

 

 

 

J /i (х) F (х) dx

 

 

где

 

со2 — — ---------------.

 

(8)

 

 

 

m| F2(х) dx

 

 

 

 

 

о

 

 

Величина со является параметром, подобным круговой частоте колебаний конструкции, соответствующим динами­ ческой нагрузке вида (3). Уравнение (7) удобно использо­ вать в безразмерном виде. Если положить

s = cof,

=

(9)

то будем иметь

(Ю)

Рассмотрим несколько примеров. Распределение нагруз­ ки по пролету балки примем равномерным, т. е. Д (х) = 1 и функция прогиба F (х) находится из уравнения

BF’V (х) = 1.

(11)

Для балки с шарнирно опертыми концами граничные ус­ ловия будут: при х = 0 и х — I F = 0, F" — 0. Из (11) находим

F (х) =

1

 

 

 

( 12)

12В

 

 

 

 

и из (8) будем иметь

 

 

 

 

 

\2]/Tl

f в

9,876

/

В_

(13)

]/ЗТ/2

V

/2

V

m

 

62

Для балки с жестко защемленными концами должно быть: при х — 0 и х — I F = О, F'= 0. В этом случае получим:

F (х) =

1

 

(14)

12В

 

 

 

 

6 У \А Л Г В ^

22,45 . / _В_

(15)

/2 V т

12 у т

 

Усилия в балке будут равны произведению усилий от статической нагрузки (л:) на функцию Т (t). Таким обра­ зом, динамический расчет конструкции сводится к решению уравнения (7), зависимость которого от свойств конструк­ ции проявляется только в величине частоты со. В рассмот­ ренных примерах частоты колебаний балки (13) и (15) почти совпадают с соответствующими по граничным усло­ виям низшими частотами ее собственных колебаний, рав­ ными соответственно [61]:

со —

я2

I

в_ .

со

22,37

1 Г

В_

( 16)

Z2

m J

1°-

V

m

Это объясняется тем, что статическая форма прогибов близка к форме колебания балки с низшей частотой.

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

Значения частот собственных колебаний могут быть взя­ ты из различных курсов теории колебаний или из Инструк­ ции [27].

Следует иметь в виду, что подобный расчет применим для тех конструкций, перемещения которых могут быть пред-

63

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

Выпишем значения функции динамичности Т {t) при трех законах изменения нагрузки во времени, соответствующих различным условиям'Действия волны на конструкцию (гла­ ва II). :

а) Р(0 = р ( 1-----/ ( 0 = 1 —

(см.

рис. 16 пря­

 

 

 

 

мая 2);

(17)

Г (0 = 1-----

cos at + -s ^ - ( О < /< 0 ) .

(18)

 

0

соО

 

 

Коэффициент динамичности

будет

равен: kR =

Т (tm),

где tm — время достижения функцией Т (t) наибольшего зна­

чения.

Время tm находится из уравнения

 

 

Т(0 = ---- 1—+ со sin со/ +

^

^ - = 0,

(19)

 

 

и

 

 

 

о

 

откуда

[51]

соtm= 2 arctg со0

 

(20)

и

 

 

 

 

 

 

 

 

<21»

 

 

 

 

 

 

 

Формула (21) справедлива

при со0 ^

2,33, так как тогда

 

 

Ыт<

©0.

 

 

(22)

При со0 > 200 можно принять kn =

2.

 

 

 

 

Р -Г

О < / < 0 1;

 

 

б)

р(0 =

 

 

 

 

 

(23)

 

Р

1 _ Ч г )

 

 

 

 

 

01< /< 01+ 0*= 0 (см- рис19>б);

 

r ( 0

= Ti = - ----s-i£i2L;

О < / < 0 , ;

(24)

 

 

et

 

coOj

 

*

\ i

64

т(t) = Тг = 1 — Ц А + f_L _ + _ L

b in со ( / - 0 , ) -

 

 

 

0-,

\

 

CO0J

 

Ш0П

 

 

 

 

 

 

 

sinotf ■;

 

0J, <

t <

0.

 

 

 

(2 5 )

 

 

 

 

СО01

 

 

 

 

 

 

 

 

 

 

Коэффициент динамичности

при 02 =

оо равен:

 

 

 

 

 

 

 

 

.

 

0)01

 

 

 

 

 

 

 

 

Ад=1 _1_

2 s i n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(26)

 

 

 

 

 

СО01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О < / <

01 = ротр);

 

в)

р (0 =■■

 

 

 

 

 

 

 

 

 

 

 

(27)

 

*

 

01Л

01< / < 0 1 + 02 = 0;

 

 

 

Робт ( (

 

 

 

 

 

 

 

02

 

(см. рис.

19,

а)

 

 

 

 

 

 

 

 

 

 

 

п ..

.

п'

 

 

 

01

0 j .

 

Т

_

Робт

<

1;

 

 

 

 

 

 

1 — Д

 

 

 

 

Ротр

 

 

 

 

 

T ( i) =

T l = l -

-J

----- cos (£>t

 

0)01

-

О < ^ < 0 ,;

(28)

 

 

 

0,

 

 

 

 

 

 

 

1

'

 

Т ( 0

= Г 2 =

( 1

-

 

 

)

Д

+ [ a ,

s i n со (/■-

0 J

. ■

 

 

 

—a 2cosco(^ —0Х)];

 

0 i ^ ^ ^ 0 ,

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а х =

 

1—cos 0)01 .

 

Д

 

 

 

 

 

sin со0г---------- =------ |-

 

 

 

 

 

 

 

 

 

 

 

(О0j

 

 

 

 

 

 

 

 

 

а 2 -

cos со©!

s i n (001

 

 

 

 

 

 

 

 

 

о) 01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент динамичности будет равен:

 

 

 

при

 

 

 

(o01^2arctgoo01

 

 

 

 

 

 

 

 

 

 

 

 

 

(29)

 

 

К = 2 \ \ -

 

a r c t g

io 0 ]

 

 

 

 

(30)

 

 

 

 

(0

0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 Зак.

344

 

 

 

 

 

 

 

 

 

 

 

 

6 5

при

 

 

 

 

 

 

 

 

0)0 ! < 2

arctgcoOj

 

(31)

&„ =

( 1 — -

. 9l-l А + [ai sin со (t*—0г) —

 

 

— a2 cos со (t *— Bj)],

 

(32)

 

 

где

величина

t* находится из вы­

 

 

ражения

 

 

 

 

 

 

tg

со (/*—01)

 

 

 

 

 

 

 

 

 

 

а . + j /

а |

Д°-

 

 

 

 

+ а j

 

 

 

 

 

 

(со02Г

.(33)

 

 

 

 

 

 

Рис. 26. Коэффициент ди­

 

«1+- со0

 

намичности в

упругой

На рис. 26 приведены графики

стадии (нагрузка при об­

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

динамичности

для

текании)

 

 

 

нагрузки вида (27) при А = 0,5.

При А =

0,5 условия (29) и (31) эквивалентны соответ­

ственно следующим:

 

 

 

 

 

 

со0! >

2, 785;

co0x <

2,785.

(34)

3.ДИНАМИЧЕСКАЯ ПРОЧНОСТЬ АРМАТУРНОЙ СТАЛИ

ИБЕТОНА

В настоящее время в строительных конструкциях наи­ более широко используются малоуглеродистые стали клас­ сов A-I (СтЗ), А-П (Ст5), А-Ш (35ГС). Результаты экспе­ риментов, например [18, 31, 67], свидетельствуют о суще­ ственном влиянии скорости деформирования на прочност­ ные свойства этих сталей. В дальнейшем при расчете кон­ струкций, армированных такими сталями, будем для опреде­ ления динамического предела текучести применять крите­ рий, предложенный Д. Кэмпбеллом [72], в том виде, как он представлен в работе [32]. В ней для случая одноосного на­ пряженного состояния при известном законе изменения на­ пряжения в упругой стадии a (t) динамический предел те­ кучести принят равным:

ад = о (т),

(35)

66

где т — время запаздывания пластических деформаций (вре­ мя конца упругой стадии), определяемое из уравнения

Т

 

^[о(/)]“ d/-= f0a“ при [а (т) > а0]

(36)

о

 

и условии, что в течение некоторого промежутка времени a (t) > сг0. Здесь а0-— величина статического предела теку­ чести стали, t0 и а — некоторые параметры, зависящие от свойств стали и от температуры.

Параметр t0 равен времени запаздывания для случая, когда напряжение прикладывается мгновенно и равно ста­ тическому пределу текучести. На основе обработки опытных данных в работе [32] получены значения параметров для

сталей классов A-I

и

А-П

при комнатной

температуре:

ос = 17, t0 = 0,895

сек.

Для

стали класса

A-III влияние

скорости деформирования несколько меньше, чем для ста­ лей классов A-I и А-П. При практических расчетах это можно учесть путем уменьшения величины расчетного сопро­ тивления для стали класса A-III (на 5—10%).

Выражение для напряжения в арматуре ст (t) в упругой

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

 

о (/)

= асТ (t),

(37)

где а 0 — напряжение от

нагрузки, приложенной

статиче­

ски и равной некоторой величине динамической нагрузки; Т (t) — функция динамичности.

Из (36) и (37) следует

(38)

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

Для случая нагрузки вида (17) примем Т (t) = kt. Коэф­ фициент k и время запаздывания т находятся из (38) и ус­ ловия kx = Т (т).

3 *

67

Исключая k, получаем следующее уравнение для опре­

деления т:

1_

[to(« + !)]“

=

Т(х).

(39)

ч0

Для сталей классов A-I, А-П, А-Ш имеем

1,1776

= т '/ 17 Т (т).

(40)

Из (39) и (40) следует зависимость динамического пре­ дела текучести од = осТ(т) от времени его достижения:

j _

__i_

 

 

од = [^о (а + 1)1 т

 

сто>

 

ад = 1,1776т-*/17 а0

(т в

сек)

(41)

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

при т =

t*,

когда

 

 

 

da

0.

 

 

IF '=r =

При

t >

t* напряжение уменьшается, вследствие чего

в арматуре могут возникнуть лишь ничтожно малые пласти­

ческие деформации и можно считать, что конструкция еще

не вышла из

упругой стадии. Будем называть напряжение,

равное а* =

а ( t*), минимальным динамическим пределом

текучести.

 

Пусть на конструкцию действует постоянная во времени

(т. е. 0 = оо),

мгновенно возрастающая динамическая на­

грузка. Тогда,

как следует из

(18),

Т (t) = I— cos a>t

и Г = —

 

 

и из (38) для этого случая получим точное решение

Л/(0

^ (1 —cos со/)а dt =

о

а я/ш

а + 1

 

 

sin2a xdx.

(42)

 

1

 

68

Последний интеграл точно вычисляется при заданном значении а:

 

 

 

 

si n2a xdx — (2а—1)11

я

 

 

 

 

 

 

 

 

 

 

 

 

2а11

Т '

 

 

 

Из (42) найдем максимальное значение стс (обозначим его

ос), при котором

конструкция

работает

только

в

упругой

стадии:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

(га/0)1/а (Jo

 

 

 

(43)

 

 

 

 

0 ~

 

а± 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

а

/

 

 

 

 

 

Минимальный

динамический

предел

текучести

равен:

 

 

о,, =

2ос =

(21)~1/а (©/„Г ' /а°о,

 

 

(44)

откуда

при а

=

17 имеем

 

 

 

 

 

 

а* =

1,0513 (со/о)1/17

а0;

ас =

0,5256 (сог'0)1/17 сг0

(45)

и при t0 =

0,895 сек

 

 

 

 

 

 

 

 

 

 

 

 

а* =

 

1,0445©'' / ' Ч

 

 

 

(46)

(здесь

под

со следует

понимать

безразмерную

величину

со • 1 сек).

 

 

пользования

формулами

(45)

и (46)

Для

облегчения

на рис. 27 построен график функции со|/|7 .

 

 

 

Если в формуле

(41) принять т = t*

=

то получим

 

 

 

 

о*

=

1,1 со1717 а0.

 

 

 

(47)

Сравнение с (46) показывает, что приближенный способ определения динамического предела текучести дает погреш­ ность ~ 5%.

Рассмотрим теперь действие нагрузки вида (17). Для нее функция Т (t) и ее максимум определяются по формулам

(18) и (21).

Если подставить (18) в уравнение (38), то решить его затруднительно. Поэтому заменим (18) следующим выраже­

нием:

 

Г* (0 = ( 1 — -аГС^ — ) U ~ cos to* Q.

(48)

69

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