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

книги из ГПНТБ / Фоломеев, А. А. Снижение материалоемкости железобетонных конструкций-1

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

может быть реализована на трехадресных машинах типа М-20,

М-220, М-222, БЭСМ-4.

Программа «корень функции» (КФ) является составной частью библиотеки компилирующих программ (БКП), которая занимает ячейки 5000—5162. В качестве рабочих ячеек использованы

5273—5373, 0160, в качестве ячеек для констант — 5350—5360, 3001—3006.

Ввод информации осуществляется следующим образом.

1.Вводится перфокарта «Вызов ИС-2», которая вызывает

ИС-2 из барабана

МОЗУ

и автоматически

начинает

вводить

БКП. После ввода

БКП

проверяется

контрольная

сумма

Изгибающие моменты

и\ (г, t)

W 2 (х, t)

u ’i (у. 0 v"2 (X, t)

(*. t)

■U[ (г, 0 W'2 (x, t)

(y. t)

v l (x, t)

<f'i (*. 0

 

 

 

 

 

Т а б л и ц а

1

 

 

 

Узлы рамы

 

1 в,

 

A

B 1 в,

■^1

A

В

A t

 

t =

1

 

 

1

 

 

50

 

 

1 — 30 *0

 

+ 2 , 0 0

+ 2,05 - 2 ,0 4

- 1 ,9 0

+ 2 , 2 2

+ 2,18

- 2 ,1 3

- 2 , 0 0

+ 0,92 + 0,75

- 0 , 6 8

- 0 ,7 0

+ 1,04 +0,97

- 0 ,7 9

- 0 ,8 0

- 1 ,4 6

- 1 ,4 0

± 1 .36

+ 1.30 - 1 ,5 6

- 1 ,5 2

+ 1.46 + 1.40

—1,03 - 1 , 0 0

+0,98 + 1 , 0 0

- 1 ,1 5

- 1 , 1 2

+ 1,06 + 1,08

- 0 ,4 3

- 0 ,4 0

+0,38

+ 0,30 - 0 ,5 5

-0,51

+ 0 ,4 0 + 0,32

 

t

1

 

 

 

 

 

 

К) ' 0

 

 

 

 

 

 

 

 

 

 

 

 

+2,46 T 2,50 - 2 ,4 3

—2,40

f 3,13

+3,10

- 3 ,0 2

-3 ,1 1

+ 1.34 + 1,26

-1 ,3 7

- 1 ,2 4

+ 1,67

+ 1,54 - 1 ,5 0

- 1 ,2 3

— 1,90!— 1,88 + 1,80 + 1.82 - 2 ,4 0

- 2 ,3 2

+ 2,26

+ 2,17

1. 20

- 1, 12

f-1 , 1 0

+ 1,14 - 1 ,4 6

- 1 ,4 0

+ 1,29 + 1,16

- 0 ,7 0

-0 ,7 6

+ 0,70

+ 0 , 6 8

-0 ,9 4

-0 ,9 2

+0,97

+ 1,01

К2 172 5242 6264 1741, при совпадении которой происходит ввод ■основной программы.

2.Далее, выполняется передача управления (ПУ) на начало программы и машина начинает счет.

3.Для обращения к программе КФ необходимо, чтобы функция была оформлена как блок со стандартным входом и выходом.

90

4. Перед передачей управления на начало программы КФ (ячейка 5065 — начало программы) нужно записать в ячейку5121 команду: 056 5355 Бл. (дг) х рдля / (д);

здесь

Бл. f ( x ) — начальная ячейка блока счета функции, х р для / (л) — рабочая ячейка аргумента функции,

5355 — рабочая ячейка аргумента функции при нахож­ дении корня.

Изгибающие моменты

U \ (+ 0

W 2 (х, t)

Ul (у, t)

V \ (x, t)

<p! ( ?. 0

U[ (г. 0

W"2 (x , t)

Щ (y, t)

V\ (x, t)

?! (z. *)

 

 

 

 

 

Т а б л и ц а

2

 

 

 

Узлы

рамы

 

 

 

A

В

Si

+ I

A

В

в.

At

 

t =

1

 

 

1

 

 

 

4 to

 

 

i - 3*0

 

+2,70^+2,75 - 2 ,7 2

- 2 ,6 3

+ 2 . 1 0

+ 2,16

- 2 ,1 7

- 2 ,0 3

+ 1,94 + 1,84 —1 , 86

- 2 ,0 5

f 1 , 1 2

+0,78

- 0 ,6 7

- 0 ,7 4

- 2 ,3 2

- 2 ,3 0

+ 2 ,2 9

+ 2,24

-1,61

- 1 .4 7

+ 1,42 + 1.39

—1,40 - 1 ,4 7

+ 1,42 + 1,37 - 1 ,0 3

- 1 ,0 7

+ 1 , 00

+ 1 , 0 1

- 0 ,9 2

- 0 ,8 3

+ 0,87

+0.87

- 0 ,5 8 —0,40 + 0,42

+ 0,38

 

t =

1

 

t = tQ=

0 , 1 сек.

 

2,5 to

 

+ 2 , 0 0

+ 1 , 8 6

—1,79 - 1 ,7 6

+ 1,77 + 1,70 - 1 ,6 7

- 1 ,6 0

+ 1 , 20 +0,98

- 0 ,9 7

- 0 , 8 8

+ 1.07 + 0,90 —0,85 - 0 ,8 0

—1,60 - 1 ,4 2

+ 1.38 ±1.37

- 1 ,4 2

- 1 ,3 0

+ 1,26 + 1 , 2 0

1 , 00

—0,96 + 0,95 + 0,89 - 1 ,0 3 - 0 ,9 8

+ 0 ,9 3 + 0 ,9 0

-0 ,6 0 - 0 ,4 6

+ 0,43 + 0,38

- 0 ,3 9

- 0 ,3 2

+ 0,33 + 0,30

5.Производится обращение к программе «корень функции»

(КФ): Я: 016 Я+1 5065 0007.

6.При выходе из программы КФ окажутся вычисленными пер­

вые пять корней функции, они засылаются в ячейки 3001—3005 со­

ответственно.

ранее обозначения

^2п,

ц,)п и т. д. можно преоб­

Принятые

разовать следующим образом:

 

 

 

^гл =

^2 »

 

^зл =

1*1*

;

91

здесь

Перерезывающие силы

и \ " Ц. 0

w 2" (X, t)

С'з" (у. 0

V's" { х ,t)

?I t)

u \ " (г, t)

W ’2 ' {x, t )

C 3” (У, t)

-

£

—\ /

ii (*. 0

 

 

 

 

 

Т а б л и ц а

3

 

 

 

Узл ы

рамы

 

 

 

 

 

B

в .

■ 4,

 

-

 

в ,

л,

 

 

1

 

 

 

 

 

 

 

t — 50

 

 

t =

35*°

 

+ 1, 21

-1 ,1 7

-1 ,1 6 + 1 , 2 0

т 1,28 —1,24 —1 , 2 0

+ 1,26

—0,87 +0,84 +0,84 -0 ,8 7 -0 ,9 4 +0,91

+0,91

-0 ,9 4

+ 1 , 2 2

- 1 , 2 0

-1 ,1 9

+ 1 , 2 2

+ 1,30 -1 ,2 6

-1 ,2 7

+ 1,31

-0,91

+0,90

f 0.90 -0 ,9 1

-0 ,9 9

+0,97

г0,97 -0 ,9 9

-0 ,6 5

+0,63 +0,61

-0 ,6 4

-0 ,6 5

+ 0,63

+0,59 -0 ,6 4

 

t =

1

 

 

 

1

 

 

 

10

 

 

i =

e

 

+ 1,29 -1,31

-1 ,3 6

+ 1,28 + 1,41 -1 ,4 0

-1 ,3 7

-1 ,4 0

- 1. 01

+0,98 +0,98 - 1, 01

-1 ,1 6

+ 1.Ю + 1 , 1 0

-1 ,1 6

+ 1,30 -1 ,3 3

-1 ,3 2

+ 1,34 + 1,45 -1 ,4 3

-1 ,4 3

+ 1,46

- 1 , 1 2

+ 1 , 1 0

+ 1 , 1 0

- 1 , 1 2

-1,27

+ 1 , 2 2

+ 1 , 2 2

-1 ,2 7

-0 ,4 6

+0,56 +0,54

-0 ,4 9

-0 ,4 3

+ 0,51

+0,48 -0 ,4 3

Для соотношений длин элементов пространственной рамной

системы -j- = С2, -j- = чз принимаем следующие значения: 1,00.

1,25; 1,50; 1,75; 2,00.

Рассмотрим порядок выполнения машинных операций на ЭЦВМ БЭСМ-4 для определения функций смещений четырехстоечной рамной системы без учета кручения стоек.

1. После определения корней частотного уравнения находятся частоты собственных колебаний.

2. Определяются численные значения выражений

■\и’ ^ 2л' \ п ’ ••• ’ ^ 1л’ Л ’л’ ^Лн

92

срезывающие силы

' (*, о w 2 ' (м 0

(у. О

К(•*. О

<h (*. 0

'(*, О

'(х , 0

+' (у<0

У'г ' (х, t)

ь (г, ()

 

 

 

 

 

 

Т а б л и ц а

4

 

 

 

 

Узлы рамы

 

 

 

 

 

 

А

в

В.

А,

А

 

В

в,

 

л.

 

 

 

1

 

 

 

1

 

 

 

 

 

t =

 

 

 

 

^ — 2

 

 

 

 

4

 

 

 

 

" *и

 

 

 

+ 1,27 - 1 ,2 8

— 1,26 + 1,28 + 1 , 2 0

- 1 ,1 7

- 1,11

+ 1 , is

- 1.11

+ 1, 10

+ 1, 10

1,11

- 1 ,0 4

+ 1 , 0 2

+ 1 , 02

- 1 ,0 4

+ 1,33 —1,36 - 1 ,3 6

+ 1,35 + 1,27

- 1 ,3 0

- 1 ,3 3

+ 1,30

- 1 .1 7

+ 1, 12

+ 1,12

- 1 ,1 7

1 , 1 0

+ 1,06

+ 1,06

-

1 , 1 0

- 0 ,3 2

+ 0,42

+0,40 - 0 ,3 5

- 0

33 + 0,39

+0,36 - 0 ,3 4

 

t =

1

 

t

/ , = 0, 1

сек.

 

 

 

 

25*о

 

 

 

 

+ 0,98 - 0 ,9 6

- 0 .9 6

+ 0,98

+ 0,70 - 0 ,7 2

- 0 ,7 2

 

О

О

- 0 ,8 0

+0,81

+0,81

- 0 ,8 0

-0 ,6 1

+ 0,59

—0,59 - 0 ,6 4

+ 1,03 - 1 , 0 0

- 1 ,0 3

+ 1,06 + 0,73 - 0 ,7 0

1О+1о+ 0 .7 2

- 0 ,7 2

+ 0,70

+ 0 ,7 0 - 0 ,7 2

- 0 ,5 7

+ 0 ,5 6

+ 0,56 - 0 ,5 7

- 0 ,4 9

+ 0,45

+ 0,48 - 0 ,5 2 —0,25

+ 0,27

+0,27 - 0 ,2 4

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

4. После нахождения численного значения (

+ ]

для

\

11/ !Ч—Т-1л

каждого отдельного \>in определяют численные значения

для вы­

ражения т.д/Ц

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

Для вычисления оригиналов перемещений, изгибающих момен­ тов и перерезывающих сил на ЭЦВМ БЭСМ-4 сделаем допуще­ ние, что корни частотного уравнения собственных колебаний или известны, или значения +/л необходимо определить.

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

1.Засылка начальной цифровой информации, которая принята

ввиде констант I, Е, J и т. д.

93

2.Вычисление частот собственных колебаний ы.л-

3.Нахождение численных значений выражений/1л ( [^л), / 2п(У‘ч^)

ит. д. и производных от этих выражений при определенном со­ отношении параметров.

4.Представление сейсмической нагрузки в виде аналитичес­ ких выражений и нахождение численных значений ^in(t) при раз­

личных значениях времени t и жесткостных и геометрических параметрах. Использованы численные значения, полученные по формулам (IV.53)— (IV.57).

5.Вычисление членов рядов (IV.53) —(IV.57), а также выраже­ ний и их численных значений для изгибающих моментов и перере­ зывающих сил с последующим суммированием для разных п.

6.Вывод в печать окончательных данных, которые далее пред­ ставлены в виде графиков и табл. 1—4.

Ниже приводятся блок-схемы выполнения машинных операций

на ЭЦВМ БЭСМ-4 (рис. 26, 27). I

§ 6. АНАЛИЗ ДЕФОРМАТИВНЫХ ФАКТОРОВ ЭЛЕМЕНТОВ ПРОСТРАНСТВЕННОЙ РАМЫ

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

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

чения ярко выражено для высших частот

(четвертой и пятой)

(рис. 28).

 

 

На рис. 29—31 показан характер изменения смещений в объем­

ном

кресте по длине стойки (t),

продольного (/) и поперечного

(k)

ригелей.

 

объемного креста

Как видно из иллюстраций, форма изгиба

близка к форме изгиба консоли с массой на конце.

В поперечном ригеле k изменение смещения по длине проявля­

ется

больше, чем в продольном i.

Эта характерность выявляется

94

Рис. 26. Блок-схема машинных операций на БЭСМ-4 (расчет без учета кручения).

Очистка регистра, перевод 10-^2

Засылка начальных значений a, {i, у,

l1! и' s. k'i-

Изменение а, 8, у. Восстановление и,

Изменение р., , восстановление

-----Изменен ie 6, восстановление t

Изменение t,

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

Изменение

k\\ k\ = £ /( : pS

Счет правой части, т. е. вычисление

U, V', W,

с?

перевод 2 - ^ 1 0

и печать

Проверка окончания по kf

Да.

Проверка окончания по t

Даj

Проверка окончания по t

Да!

Нет

Нет

Нет

Нет

Проверка окончания по р.,л

Да!

Проверка окончания по а, ,8 у

Да

Нет

Стоп

Рис. 27. Блок-схема машинных операций на БЭСМ-4 (расчет с учетом кручения).

при меньших размерах поперечного сечения и больших значениях длин элементов.

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

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

На рис. 32—34 представлены эпюры смещений, изгибающих мо­ ментов и перерезывающих сил дл четырехстоечнои рамы, расчлененной на плоские рамы. По ха­

рактеру

изменений

смещения,

 

 

 

изгибающих моментов и перере­

 

 

 

зывающих

 

сил

в портальных и

 

 

 

замкнутых плоских

рамах, полу­

 

 

 

ченных при

расчленении

прост­

 

 

 

ранственной рамы, можно сде­

 

 

 

лать следующие выводы:

 

 

 

 

 

1.

С увеличением времени дей­

 

 

 

ствия

импульса

возрастает

вели­

 

 

 

чина смещения, причем величина

 

 

 

смещений

продольного

ригеля

в

 

 

 

плоскости xoz больше, чем в

 

 

 

плоскости

хоу.

Это

свидетельст­

 

 

 

вует о преобладающей

 

доле де­

 

 

 

формаций

продольного

ригеля

в

 

 

 

плоскости xoz. Точки перегиба в

 

 

 

плоскостях xoz

и

хоу

близки,

 

 

 

но не совпадают.

 

и попереч­

 

 

 

2.

Если

для

стоек

 

 

 

ных

ригелей характер

изменения

 

 

 

смещений

 

аналогичен

случаю,

 

 

 

приведенному ранее (т.

е. в сис­

 

 

 

теме «объемный крест»),

то для

 

 

 

продольных

ригелей

i

изменение

 

 

 

Vn Wn

насколько

 

позволил

Рис. 28. Изменение корне'] частот­

уловить масштаб, в

замкнутых

ного уравнения с учетом и без учета

рамах иное,

чем

в

объёмном

 

кручения.

 

кресте.

 

 

 

 

 

 

 

 

 

 

 

3.

Эпюры изгибающих

моментов

Vt для продольных риге­

лей имеют криволинейное очертание.

 

4.

Эпюры изменения

перерезывающих сил для V’t и W t“ кри­

волинейные,

причем

величина

 

Wt

меняет знак, a I/”

имеет вы­

пуклую пологость

с

уменьшением

значений ближе

к узлам и

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

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

сивной нагрузки продолжительностью t =

t0.

7-207

97

( а

1,0

 

 

 

|

■ ■ ■.

1,5

2,5

3,5

4.5

0,5

1,5

Uj(Z,

Рис. 29. Изменение смещений в объемном крес­ те по длине элемента U^ при действии импуль­

сивной нагрузки, I = 6 ,0 (а) и 3,0 м (б)\

/ - p

F =

0

, 0 0

9

S ,

2

- р

/ ? = 0 , 0 3 8 4 .

Рис. 30. Изменение смещений в объемном крес­ те по длине элементов Vt и Uk при действии

импульсивной нагрузки, I = 6,0 (а) и 3,0 м (б):

/ - p

F =

0 , 0 0 9 6 ,

2

- р

/ ? = 0 , 0 3 8 4 .

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

Характер изменения Vt и Wt был с достаточной степенью

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

моментов и перерезывающих сил. Наибольшая доля

при построе­

нии эпюр смещений с учетом кручения

в стойках приходится

на

стойки и

поперечные

ригели.

На

 

 

 

 

 

рис.

35 показаны эпюры смещений,

 

 

 

 

 

которые с

большой

наглядностью

 

 

 

 

 

характеризуют изменения

конфи­

 

 

 

 

 

гураций рамы во времени и в про­

 

 

 

 

 

странстве.

 

 

было

реализова­

 

 

 

 

 

В дальнейшем

 

 

 

 

 

но предложение

о раздельном

 

по­

 

 

 

 

 

строении эпюр изгибающих

момен­

 

 

 

 

 

тов и перерезывающих сил с учетом

 

 

 

 

 

изгибных колебаний стоек и колеба­

 

 

 

 

 

ний продольного ригеля в плоскос­

 

 

 

 

 

ти xoz, с учетом крутильных колеба­

 

 

 

 

 

ний стоек

совместно

с изгибными

 

 

 

 

 

колебаниями

ригеля

в

плоскости

 

 

 

 

 

хоу. Это схематическое разделение

 

 

 

 

 

оказалось весьма удачным и позво­

 

 

 

 

 

лило выделить две фазы

в колеба­

 

 

 

 

 

ниях рамной системы при действии

Рис.

31. Изменение смещений i

импульса:

фазу изгибных

и

фазу

крутильных

колебаний

стойки.

В

объемном кресте по длине эле

действительности

стойки

испыты­

ментов Vt при действии импуль­

сивной нагрузки,

I = 6 , 0

(a)

i

вают

изгибно-крутильные

колеба­

 

3,0 м

(б) :

 

 

ния, поперечный

ригель — изгибные

 

/ —р/7—0,0096, г-рд^.озв-}.

 

 

колебания в плоскости хоу,

а про­

 

в плоскостях хоу

и xoz.

дольный ригель — изгибные

колебания

Было принято предположение, что изгибные колебания ригеля в плоскости хоу более причастны к появлению кручения в стойках.

Рассмотрим первое явление. Судя по эпюрам изгибающих мо­ ментов (рис. Щ), доля кручения в стойках по сравнению с изгибом была, естественно, меньше, но в пределах 20% от изгиба. Нам ка­ жется, что на данные цифры следует обратить внимание, если

учесть, что бетон плохо работает на срез.

а гг

Характер изменения эпюр V{ и Wt в обоих случаях криво­ линейный, с несовпадающими точками перегиба. Даже простое

99

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