
книги из ГПНТБ / Бакиров Р.О. Применение современных электронных вычислительных машин при расчете и проектировании конструкций инженерных сооружений учебное пособие
.pdfВеличины я, Я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 |
|
Ък.б |
||
|
||
7м |
т/м 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, + |
• • • + |
85Б Х ъ+ |
о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 |
1Л |
|
|
|
|
|
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