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

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

.pdf
Скачиваний:
19
Добавлен:
29.10.2023
Размер:
4.93 Mб
Скачать

Величины я, Я0, (а, (Ад jaj,

Е задаются в качестве исходных дан­

ных

усилия в произвольном

сечении стенки

 

 

б)

 

 

 

Му = Ху +

0,5 Syhx - 0,5 е% у2—0,25 (Де« +

Дв«) у2 -

 

 

1

3

1 >>

 

 

 

(51)

 

 

Ny =

°Л-Л дг +

5y — -Gy;

 

(52)

 

Qv =

^ — ерТУ — 0,5 (Дг" +

Д е ") • у — o>y,

(53)

 

•S’, =

i*, [e^-y +

0,5 (betg +

Дe ^ - y +

a»y];

 

 

 

дест =

Двст (Лу—у) : Ay;

 

 

 

 

 

 

 

л-Н

 

 

 

 

 

°у = °0

Лу

 

 

 

 

 

Уц=У

л

 

1

 

(54)

 

 

п

2

 

 

 

 

 

 

 

Gy

 

Л, '.У ' Тжб>

 

 

 

 

_У =У 1, i “ 1, 2, . . . , 5;

 

 

 

 

Уа = 0;

у,+! =

У/ + 0,2 Ау.

 

 

В приведенных выше формулах все величины в т

и м.

4.

Подбор арматуры

 

 

 

 

 

 

 

а) Расчетные формулы:

 

 

 

 

 

 

 

Случай первый х < 0 ,5 5 с ? 0:

 

 

 

 

 

 

tnaFaRa + N <

maFaRax + m6-x-b-R„;

(55)

 

N -e ~-Cm6-x-b-R„

 

 

 

+

maF a'Ra.c (d0 - а ),

(56)

где

 

df)

 

d$

1 j

 

 

 

 

 

 

 

 

 

 

 

a' =

1;

 

 

 

 

 

e =

eQ+

~

— a';

 

 

 

 

е й

 

\ M

\

 

 

 

 

 

 

1*1 ‘

 

 

 

 

 

 

 

 

 

63

Величины т а, т 6, /?а,

/?и, b

задаются как исходные

данные;

F а, F a , х — неизвестные,

определяемые согласно последователь­

ности, приводимой ниже.

 

 

 

Случай второй .x:>■ 0,55^70•

 

 

N-e<C m6RnpS 0 +

F aRa.с (d0—a ’) т а,

(57)

где

 

 

 

 

 

b d 02 .

 

 

/?пр =

0 , 8

 

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

б) При подборе арматуры нужно придерживаться такой после­ довательности действий:

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

т а, т 6, # а = /?а.с, /?,„ й, a', d0, Мь

— согласно (55) определяется х\

— производится анализ

•к < 0,55 d0-,

(д)

при удовлетворении этому условию, величина

Fa—Fa опреде­

ляется согласно формулы (56); в противоположном случае — пло­ щадь арматуры определяется согласно (57);

— далее осуществляется проверка условия

 

F a< 0;

(е)

— если это условие удовлетворяется, то принимается Fa= F a' = =0,002 bd0 и записывается в результат с выдачей специального признака; в противоположном случае производится проверка до­ полнительного условия

F a < 0,002 bd0,

(ж)

— при удовлетворении этому условию принимается

Fa= F a =

=0,002 bdo и записывается в результат с выдачей другого специаль­

ного

признака; в противоположном случае расчетное

значение

Fa= F a записывается в результат без всякого признака.

 

5.

Экономические показатели обделки

 

а)

Стоимость бетона в деле

 

где

Рь — Vb'Ptn

(58)

V = [ 2 hyhx + 0,5 (d3+ dn) L];

(59)

 

 

y7* + 5,33 / 2.

 

64

б)

Стоимость арматуры

 

 

 

 

где

Ри

 

 

(60)

6\, = [1,6 hyFtl +

2,4 hyFti +

(F t, + F a.,) L) T„;

(61}

 

^a2, F a3, Ca, — площади

сечения

арматуры в расчетных

се­

чениях.

 

 

 

 

в)

Стоимость разработки породы

 

 

где

Р г р ~

^rp'Prpi

(62)

Vrp — 2 [а/гу +

/ ' (/0 +

2 sin сря)].

(63)

 

г)

Полная стоимость 1 пог. м обделки

 

 

Р ц =

Р б

+ Р „ +

РГр-

 

Значения рв, Ра, Ргр задаются как исходные данные. В формулах (58) 4- (63) все величины в м, т, руб.

Исходные данные

Величины исходных данных задачи и их размерности приведе­ ны в табл. 11.

Выходные данные

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

При этом:

исходные данные, усилия, площади арматуры и экономиче­ ские показатели обделки печатаются во всех случаях;

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

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

Последовательность вычислений

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

Схема составлена на основе знания физической сущности зада­ чи и основных сведений об ЭЦВМ, наличия всех математических зависимостей и последовательности их использования для решения задачи в полном объеме.

5 Зак. 1018

65

бл о к -схем а реш ения

задачи

 

Вбод исходных данных

б память

 

-----—________________________ 1_______________________________

 

— - Лнализ конца решения бсех

вариантов

 

Да

 

 

Нет

 

Останов

машины

 

J Вызоб исходных данных L+i варианта

 

 

_______t________________________________

 

Лнализ Дур

(9гР +9сВ)п><£

 

Д а

 

 

Нет

 

------------- ---------------

 

 

 

 

о ^ < й д р >

 

 

-------------------------------- Г

 

 

Лнализ

е%%^0

нет

—|

До

 

 

Т

 

 

 

 

0 + <ес” >, 0 < д е § Р >

 

Определение еРр по(35) ', анализДесрыерт<0J

Да

 

 

Нет

 

— -------------1-------------—

 

 

 

0 — <г Л е%р >

 

 

 

 

 

Вычисление е£%

по(32)> анализАеср -еср <£

 

Да

 

 

Нет

- |

 

0 - ^ < А e f >

 

 

 

 

Г ______________________________

 

Вычисление М °

Н? , М0 согласно (38)^39)^90,

 

O JI

<0,11>

 

 

 

 

}

 

 

Определение для свода MPrt Н и , Ох.

L.

Подготовка н обращению в СБ~1 (блок реализации подбора арматуры согласно

формулам (бб), (66) - (6 7 )

Г____________

Обращение к Сб~1 для определения ЛаСБ и FpB б замке и пяте свода

 

 

Днализ оС ^ ОД

Д а

Нет

________________ 1___________

 

Лнализ оС ^ 1,5

Д а

Нет

_________________1_______________

2

<эС >

1 —- <оС>

0 -*■ <ГоС>

1

Вычисление пром еж ут очны х Величин согласно J(96)<-(50)

1

Определение усилий б стенке

Мр, Nv, О,/ согласна (51) -г(S3)

___________

<

 

Подготовка к обращению б СВ~1

_________________ t

Обращение к С6-1 для определения площади арматуры б двух заданных

сечениях стенки

*

Вычисление Vy} GM) Угр ,Ру, Р„ , Ргр и Р0

согласно (5 8) - (63) i

Анализ необходимости Выдачи бсех данных

ттг ________ЧёЛ1_______________

Печать бсех вы ходны х данных

Печать Выходных данных без промежу­ точных результатов

к

Необходимые восстановления

Рис. 12

5 *

67

пп.

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

Т а б л и ц а 11

 

 

Наименование

I

Условн.

Размер-

 

 

обозн.

ность

 

 

 

I

 

 

 

 

 

I. П араметры динамической нагрузки

 

 

 

 

Вертикальная нагрузка ..................................................

 

ёэ

 

 

т ! м 2

Боковая нагрузка

.........................................................

 

и1

 

-

Коэффициент перегрузки для постоянной нагрузки

 

 

 

 

 

 

 

 

Коэффициент перегрузки для временной нагрузки

%

 

 

 

II.

Геометрические

характеристики

 

 

 

 

Пролет в

с в е т у .................................................................

 

 

 

 

 

Расчетная

стрела

подъема с в о д а ..............................

/

 

 

 

Расчетный пролет свод а.................................................

 

1

 

 

 

Толщина свода в зам к е ..................................................

 

d0

 

 

Толщина свода в п я т е ...................................................

 

dn

 

 

Толщина стенки................................................................

 

 

hx

 

 

Наружная

высота с т е н к и ...........................................

 

 

 

Функция соотношения f : 1 ....................................

 

Pi

 

 

 

III.

Характеристики горной породы

 

 

 

 

Объемный вес ................................................................

 

 

Тгр

 

 

Уго л внутреннего

трения ..............................................

 

?

 

 

 

Крепость по П ротодьяконову....................................

 

/ур

<

 

Модуль продольной деформации

F

о

 

с

 

 

Модуль поперечной д еф орм ац и и .............................

Ро

 

 

Коэффициент трения бетона по грунту . .

Pi

 

 

 

IV. Характеристики

материала обделки

 

 

 

19

Расчетное сопротивление арматуры . . . .

20Расчетное сопротивление бетона на сжатие при изгибе ...............................................................................

21Объемный вес ж е л е зо б е т о н а ....................................

22Объемный вес м е т а л л а .............................................

23Модуль упругости б е т о н а ..........................................

24

Коэффициент Пуассона для бетона . . . .

25

Защитный с л о й ...........................................

....... . .

 

 

 

V.

Прочие характеристики

26

Характеристика необходимости учета Дg, Де .

27

Стоимость

1

г м еталла..................................................

28

Стоимость

1

м3 бетона в д е л е ...................................

29

Стоимость

разработки

1 -н3 породы . . . .

30

Признак необходимости учета бокового давления

 

грунта (при учете бокового давления Пе = 0,

 

в противоположном

случае Пе = 1)....................

Да

кг Iсм*.

т / м 3

Ък.б

 

т/м 3

Е

Р

м

а ’

г

ру б ./ т

Рм

Рб

руб./м 3

Ртр

руб.1м3

^е.

 

68

.Наименование выходных данных

пп.

 

1 Исходные данные задачи

2Двойной интервал

3Расчет свода:

результаты промежуточных вычислений

моменты

нормальные силы

поперечные силы

—Да, Да' для сечения в замке

—Да, Да' Для сечения в пяте

4Двойной интервал

5Расчет стенки:

результаты промежуточных вычислений

моменты

нормальные силы

поперечные силы

— Да

Fа'

для первого

сечения

— Да.

Да'

для второго

сечения

•6 Двойной интервал

7Технико-экономические показатели:

объем;

стоимость

Т а б л и ц а 12

Примечания

Отделяются (друг

. от друга) одинарным интервалом

1

)

Отделяются (другот друга) одинар

ным интервалом -

,

Особые требования к программе

а. Программа должна быть пригодной для многократного ис­ пользования.

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

в. Весь материал задачи желательно иметь набитым на два НПЛ: на одном из них набивается программа (нё изменяемая при повторном использовании); на втором — числовой материал, из­ меняющийся при решении различных вариантов задачи.

г. Выдачу на печать промежуточных данных предусмотреть по желанию заказчика (например, по ключам).

69

д. Программа должна быть пригодной к использованию специ листами, не владеющими практическим программированием для ЭЦВМ.

ПРОГРАММИРОВАНИЕ ЗАДАЧИ

Наиболее совершенным методом применения машин является формулирование задачи и непосредственное программирование ее в широком смысле одним специалистом. Изложение этих вопросов выходит за рамки данной работы. Лица, желающие научиться про­ граммированию, в частности, для ЭЦВМ «Урал-2», отсылаются к обширной специальной литературе [16], [19], [30] и др.

Всоответствии с приведенным выше заданием автором совмест­ но с Самариным В. В. была разработана «Типовая программа расчета и конструирования монолитной железобетонной обделки практического очертания на ЭЦВМ «Урал-2». Эта программа, на­ шедшая применение в учебном процессе и научно-исследователь­ ской работе, приведена в приложении 4.

Взаключение заметим, что в этой программе реализована по­

становка задачи по проектированию обделки подземного сооруже­ ния, рассчитываемой по частному случаю. Когда условие (28) несоблюдается, расчет подземной конструкции производится по об­ щему методу проф. Давыдова С. С. [20], с помощью механических или электрических клавишных счетных приборов. При ручных спо­ собах расчета проектировщик не гарантирован от возможных оши­ бок вычисления. Попытка механического переноса ручных методов расчета военно-инженерных сооружений на ЭЦВМ He1 приводит к желаемому результату из-за малой эффективности получающихся программ.

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

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

Автоматизации расчета подземных конструкций должна пред­ шествовать работа по развитию, уточнению и видоизменению ранее существовавших методов расчета, а также по разработке новых ме­ тодов, ориентированных на использование ЭЦВМ. Для стати-

70

ческого и динамического расчета подземных конструкций следует шире использовать матричные формулировки задачи, которые срав­ нительно просто программируются. Матричные формулировки за ­ дач строительной механики в настоящее время интенсивно разраба­ тываются. Здесь следует отметить прежде всего труды профессоров Смирнова А. Ф., Филина А. П., доктора технических наук Лужи­ на О. В., Аргирос Дж. и др.

Рассмотрим матричную формулировку расчета одноэтажной об­ делки подземного сооружения как упругой системы в упругой сре­ де методами теории упругости и строительной механики. Схема рамы приведена на рис. 10. Расчетные размеры: /р— пролет; /р — стрела подъема; /гр— высота стенки; d3, hx, йл — толщина криволи­ нейного свода, стенки и лотка. Расчетные симметричные нагрузки: g p — вертикальная, ер — боковая. Прочие характеристики: Е 0, ро. Е, р — модуль упругости и коэффициент Пуассона соответственно для породы и материала обделки; уб — объемный вес бетона.

Исследованиями проф. Жемочкина Б. Н. [22], Давыдова С. С. [20] и Синицына А. П. [32] доказано, что касательные напряжения по подготовке и боковой поверхности заглубленного сооружения в большинстве случаев оказывают незначительное влияние на ре­ зультаты, поэтому при расчете подземных сооружений ими можно пренебречь. При этом может быть применена модель упругого слоя по методу Давыдова С. С. или модель двухслойного основания Си­ ницына А. П., которая позволяет учитывать наличие деформируе­ мой прокладки или разрыхленного грунта по наружному контуру обделки.

Для получения -расчетной схемы между обделкой и породой введем упругие стерженьки: 10 стерженьков в лотковой части и по 5 стерженьков по боковой поверхности стенок. Основную систему получим, введя заделку в середине лотка, сделав разрез в замке свода и заменив стерженьки сосредоточенными силами.

В результате, при симметричных конструкции и нагрузке, полу­ чим 13 основных неизвестных: Ai-f-As — в лотке; Хб-г-Аю — по бо­ ковой поверхности; Хц, Х ,2 — в замке, уо — в месте условной за ­ делки.

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

нонических уравнений смешанного метода:

 

 

8,, X, +

• • • +

815Х 5 +

816 Хв -)-••• +

Sn2 Х 12 —Уо +

Д)Р =

0;

З51X, +

• • • +

8Х ъ+

о66Х е +

• ■ • +

8512Х*,2— Уо+

Д5Р =

0;

3<л X, +

• ■ • +

866 Х ъ+

806 Х й+

• • • +

8612 Х 12 + 0 +

Д6р =

0;

..................... ........................................................................................(64)

З1 2 1

X, +

• • • + 8 125 ХЬ+ 8 Ш Х 6+

• • ■ + 8 1 2 1 2 Х , 2 + 0 +

Д,2р =

0;

- х

, ----------

х 6

 

 

+

£ + = о .

71

Обозначим матрицы:

 

 

 

 

 

 

 

8 п

• S . 5

°16

'

°112

-

1

 

 

 

 

 

 

 

<\

 

 

2,5 ®5б ’ ‘' '

 

-

1

 

 

 

 

 

 

 

°51

8 512

 

 

 

 

 

 

 

 

 

 

°65

°бв

■'

°612

 

0

 

 

 

 

 

 

«S

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

°121

'

'

C’l25

° 1 2 б ‘

 

''

* 0 1212

 

 

 

 

 

 

 

 

 

 

 

 

• -

1

0

 

 

 

0

 

0

8 ц

'

8 15

°16

'

'

8

 

 

■-

1

 

 

 

 

 

 

 

«51

'"

 

<>53

°50

'

'’ '

°512

— 1

 

 

 

 

 

 

\ л

■ '' •

 

\ ь

8ш>

°612

 

0

 

 

 

 

 

 

 

12!

'

° 1 2 6 ° 12«

'

'

°1 2 1 2 6

 

 

 

 

 

 

 

- 1 • • • - 1

0 „ 2

 

 

 

0

 

0

 

 

 

 

 

 

 

а также

, (65)

Х

г

-

 

-

.

,

 

 

 

 

 

 

,

х

=

 

 

-

 

-

У

»

 

 

X ,

X »

^1 р

-

^12 р

.

L Ъ Р \

Д 1Р

(66)

L- "^12 р -

где 6 — полная матрица с учетом влияния упругого основания и прокладки; б — полная матрица коэффициентов, но без учета упру­

гого основания и прокладки; б — неполная матрица коэффициен­ тов и без учета упругого основания и прокладки; 2Р — проекция внешней нагрузки на вертикальную ось.

С учетом (65) и (66) система (64) примет вид:

г * + Д р = 0.

(64')

Каждый элемент основной системы разделим на 4 равные меж­ ду собой части и их длину обозначим: li—l2= . . . ~ l 4 — в лотке; k —k = ... = k — в стенке; /9—/ю = • • • =/12 — в своде.

72

Введем трехдиагональную матрицу податливости [26]:

- 2 U

 

 

 

 

 

6 Е1Х 6 £7,

 

 

 

 

 

 

2 /,

_ 2 /з_

/2

 

 

 

6 E li 6 Е1Х 6 £У2 6 £У2

 

 

 

= [^/] 1з,1з =

 

 

 

 

 

 

 

 

*-11

~ -11 I

~ *12

*42

 

 

6 £УИ 6Е1и

6Е112 6 Е /13

 

 

 

 

^12

2

/12

 

 

 

 

6 £/„ 6 £/r, J

 

 

 

 

 

 

(67)

Кроме того, имеем матрицу

моментов в основной системе от

действия единичных сил

 

... = X i2= l :

 

 

 

Ж ' — [- М // ] 12,13

^^12 >••• » Л4пз

( 68)

- - ^ 1 2 1

122» ••• 5 ^ 1 2 1 3

Теперь матрица 6 будет равна

F = ( l _ ,х*)М ' В М ,

где

М = ( М ' У .

(69)

В соответствии с (65) дополним матрицу 6 последней строкой и

столбцом и получим матрицу 6. Учтем влияние упругого основания и прокладки и сформируем матрицу коэффициентов 6.

Матрицу грузовых моментов в основной системе от действия на­

грузок g p и ер обозначим

 

 

 

М р ~[У И 1рЛ*1р,...,Л * 13р].

(70)

Тогда неполная матрица свободных членов будет равна

 

Лр =

0

- Г ) М ' В М Р,

(71)

где

__

•_

 

м , = ( л у у .

Отсюда легко получить матрицу Др.

73

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