Karpov_Panin_Matematicheskoe_modelirovanie_i_raschet_elementov_stroitelnykh_konstruktsiy2013
.pdfМатематическое моделирование и расчет элементов строительных конструкций
1400 |
|
|
|
|
3 |
|
|||
|
|
|
|
|
3 |
1 |
|
|
|
1200 |
|
|
|
|
|
||||
|
|
|
|
|
|
||||
1000 |
|
|
|
21 |
|
|
|||
|
|
|
|
|
|
||||
|
|
|
800 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
P |
|
|
|
|
|
|
|
|
600 |
|
|
|
|
|
|
|||
400 |
11 |
1 |
|
|
|
|
|||
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
200 |
|
|
|
|
|
|
|||
0 |
0,2 |
|
|
0,4 |
0,6 |
0,8 |
|||
|
|
|
|
|
|
W
Рис. 4.5. График «нагрузка Р – прогиб W » для оболочки варианта III
На приведенных рисунках кривые с номером 1 соответствуют оболочке без ребер, кривые 2 – оболочке, подкрепленной шестью ребрами, 3 – оболочке, подкрепленной 18 ребрами. Кривые без индекса соответствуют прогибу в центре оболочек, с верхним индексом 1 – в четверти оболочки.
Какпоказалирасчеты, наличиеребервоболочкахсущественно понижаетвеличинуих прогиба.Приподкреплении оболочки18ребрамиснижение еепрогибовпосравнениюспрогибамиоболочкибез ребер при одной и той же нагрузке составляет для оболочек разных вариантов: варианта I – на 55 %, вариантаII –на 65 %, вариантаIII –
на 40 %.
Приведенный начальный модуль упругости железобетона Еm принимается в виде (см. параграф 1.9):
E |
m |
= E |
A − As |
+ E |
s |
As |
, |
(4.62) |
|
|
|||||||
|
b |
A |
||||||
|
|
|
A |
|
|
|
где As – площадь сечения арматуры в границах площади единицы длины сечения оболочки А; А – площадь единицы длины сечения
оболочки, A = h s ; s– размерячейки арматурной сетки; h – толщина оболочки.
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
Армирование железобетонныхоболочек толщиной до 5см, как правило, осуществляется одной сеткой; при толщине оболочки 5 см и более – сетками, расположенными в двух уровнях.
Например, при толщине железобетонной оболочки 3 см при одноуровневом расположении арматурной сетки из проволоки клас-
са Вр-I диаметром 3 мм (Es =1,7 105МПа) с размером ячейки сетки 100 мм приведенный модуль упругости железобетона составит (при бетоне класса В30 с Eb = 3,25 104 МПа).
Em = (0,071 10−4 17,0 +(0,03×0,1−0,071 10−4)3,25)104 =
0,03 0,1
= 3,283 104МПа.
При толщине оболочки 9 см и двухуровневом расположении арматурной проволочной сетки (Ж 4 мм, размер ячейки – 100 мм) условный приведенный модуль упругости железобетона составит:
Em = (2 0,13 10−4 17,0 +(0,09 0,1−2 0,13 10−4)3,25)104 =
0,09 0,1
= 3,294 104МПа.
Как следует из анализа полученных результатов, значение приведенного начального модуля железобетона Em незначительно пре-
восходит начальный модуль упругости бетона Eb . Соответственно,
для оценочного анализа НДС железобетонных оболочек можно впрограммном комплексе расчета оболочекзадаваться прочностными характеристиками бетона.
Дляанализапрочностибетонаоболочекприменяетсякритерий (условие прочности) Кулона – Мора
σ1 − |
Rbt σ3 ≤ σдоп . |
(4.63) |
|
Rb |
|
Главные напряжения σ1 и σ3 определяются в бетоне на верхней поверхности оболочки при z = −h /2. Обобщенный коэффициент запаса прочности принимается k = 2. Левая часть критерия Ку-
140 |
141 |
Математическое моделирование и расчет элементов строительных конструкций
лона –Мораобозначаетсятак: σg = σ1 − Rbt σ3 ипредставляетсобой
Rb
максимальное рабочее напряжение в рассматриваемой точке. По результатам расчетов становится целесообразно строить графики
распределения (эпюры) таких напряжений σg по полю оболочек с це-
льюграфическойиллюстрацииизменениянапряженногосостоянияоболочек для различных величин нагрузок (см., например, рис. 4.6–4.10).
Для проведения анализа НДС оболочек значения Rb и Rbt для
бетона принимаются согласно табл. 4.4.
С использованиемформулыперехода к безразмернымпарамет-
рам для напряжения σ = a 2σ, определяются безразмерные значения
E
допускаемых напряжений σдоп для различных вариантов оболочек при разных классах бетона. Выборочные значения допускаемых напряжений σдоп приведены в табл. 4.5.
Таблица 4.5
Безразмерныедопускаемыенапряжения дляразныхвариантов оболочек и разныхклассов бетона
|
|
Безразмерные значения допускаемых |
||||
Класс |
Модуль упругости |
|
напряжений σдоп |
|
||
бетона |
бетона Е, МПа |
|
для оболочек вариантов |
|
||
|
|
I |
|
II |
|
III |
В55 |
4 104 |
7,2 |
|
0,8 |
|
0,2 |
В40 |
3,6 104 |
7,0 |
|
0,78 |
|
0,19 |
В30 |
3,25 104 |
6,65 |
|
0,74 |
|
0,18 |
Для оболочек вариантов I, II, IIIнаходим допускаемуюнагрузку qдоп из условия потери прочности бетона оболочек с использова-
нием формулы перехода qдоп = EPa 4доп . В табл. 4.6 представлены не-
которые результаты расчета допускаемой размерной нагрузки qдоп и допускаемой безразмерной нагрузки Рдоп (цифры в скобках) для
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
оболочек вариантов I, II, III при классах бетона В55, В40 и В30 соответственно [19].
|
|
|
|
|
|
|
Таблица 4.6 |
Допускаемые нагрузки для оболочек вариантов I,II, III |
|||||||
|
|
|
|
|
|
|
|
Номер |
|
σ |
|
103, МПа ( |
Р |
|
доп ) |
варианта |
Число |
доп |
|
|
|
|
|
ребер |
для бетона |
|
для бетона |
|
|
для бетона |
|
оболочки |
класса В55 |
|
класса В40 |
|
|
класса В30 |
|
|
|
|
|
||||
I |
0 |
1,80 (5832) |
|
1,56 (5616) |
|
|
1,33 (5304) |
|
6 |
3,11 (10076) |
|
2,72 (9792) |
|
2,33 (9291) |
|
|
18 |
4,44 (14 386) |
3,89 (14 004) |
|
3,34 (13 319) |
||
II |
0 |
4,13 (1652) |
|
3,63 (1613) |
|
|
3,10 (1526) |
|
6 |
6,40 (2560) |
|
5,62 (2498) |
|
4,80 (2363) |
|
|
18 |
13,33 (5332) |
11,70 (5200) |
|
10,0 (4923) |
||
III |
0 |
6,67 (167) |
|
5,70 (158) |
|
|
4,87 (150) |
|
6 |
10,67 (267) |
|
9,12 (253) |
|
7,80 (240) |
|
|
18 |
21,33 (533) |
|
18,24 (507) |
|
15,60 (480) |
Наоснованиианализарезультатоврасчетадопускаемойнагрузки (см. табл. 4.6) считаем обоснованными выводы о том, что наличие ребер жесткости оказывает существенное влияние на напряжен- но-деформированное состояние оболочек, приводит к увеличению значений допускаемых нагрузок на них, в сравнении с допускаемыми нагрузками на гладкие оболочки. Подкрепление оболочек ребрами жесткости позволяет значительно увеличивать несущую способность оболочек.
Так, увеличение допускаемых нагрузок на оболочки варианта I, подкрепленных шестью ребрами, составляет 75 %, на оболочки варианта II – 55 %, варианта III – 60 %. Аналогично, увеличение допускаемых нагрузок на оболочки при подкреплении их 18 ребрами составляет: для оболочек варианта I – 150 %, для оболочек варианта II – 222 %, для варианта III – 220 %.
Какпоказалипроведенныерасчеты,наличиеребервоболочках существенно понижает величину их прогиба. Например, при подкрепленииоболочки18 ребрамиснижение прогиба в центреоболочки, посравнению с прогибомоболочки без реберпри одной итой же величине нагрузки, составляет: для оболочек варианта I – 55 %,
II – 65 %, III – 40 %.
142 |
143 |
Математическое моделирование и расчет элементов строительных конструкций
4.11.3. Анализ распределения прогибов и напряжений по полю оболочки
На рис. 4.6–4.8 представлены функции прогиба W
и σg = σ1 − Rbt σ3 (названные на рисунках функциями интенсивнос-
Rb
ти) для оболочки варианта I при соответствующих нагрузках (Р =10 000 соответствует q = 2,25 10−3 МПа, Р = 40 000 соответ- ствует q = 9 10−3 МПа, Р = 70 000 соответствует q =1,62 10−2 МПа, Р = 20 000 соответствует q = 4,5 10−3 МПа, Р =120 000 соответ-
ствует q = 2,8 10−2 МПа).
На рис. 4.6 представлены результаты расчета для гладкой оболочки, на рис. 4.7 – для оболочки, подкрепленной шестью ребрами, на рис. 4.8 – то же, подкрепленной 18 ребрами. Как видно из этих рисунков, максимум напряжений находится в угловых точках оболочки.
P = 10 000 |
P = 40 000 |
Рис. 4.6. Функции прогиба W и σg для гладкой оболочки варианта I
144
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
P = 10 000 |
P = 70 000 |
Рис. 4.7. Функции прогиба W и σg для оболочки варианта I с шестью ребрами
P = 20 000 |
P = 12 0000 |
Рис. 4.8. Функции прогиба W и σg для оболочки варианта I с 18 ребрами
145
Математическое моделирование и расчет элементов строительных конструкций
Чтобы снизить напряжения в угловых точках, подкрепим оболочку двумя парами контурных ребер высотой 5h и шириной 12h. На рис. 4.9 представлены функции «W –σg »для оболочки вариантаI.
Из этого рисунка видно, что произошло не только перераспределение напряжений, но и существенное снижение уровня напряжений.
P = 10 000 |
P = 40 000 |
Рис. 4.9. Функции прогиба «W –σg » для оболочки варианта I с контурными ребрами шириной 12h
Какпоказалирасчеты, наличиеребервоболочкахсущественно влияет на ее НДС, повышает величину допускаемой нагрузки относительно оболочек без ребер. Например, при подкреплении оболочек 18 ребрами ее прогибы при одной и той же нагрузке снижаются от40до 65%относительнопрогибовоболочки безребер. Величины допускаемых нагрузок на ребристые оболочки повышаются от 150 до 220 % по сравнению с гладкими оболочками.
Таким образом, используя полученные результаты, можно подобрать соответствующую толщину проектируемой оболочки, размер и число подкрепляющих оболочку ребер, надлежащее армирование по полю оболочки.
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
Для достижения приемлемого времени расчета на ЭВМ одного варианта оболочек, не теряя существенно точности расчета, достаточно брать девять членов разложения перемещений в ряды, т. е.,
принимать N = 9 в формуле (4.48).
4.12. Прочность пологих железобетонных ребристых оболочек при длительном нагружении
В результате развития деформаций ползучести материала во времени начинается бурный рост прогибов оболочек (в 10…15 раз превышающихпрогибыоболочекпри t = 0 ).Время,прикоторомэто
наступает, принимается за критическое время tкр.
На рис. 4.10, 4.11 представлены графики «W −t », полученные при различных значениях нагрузки для оболочки варианта I(о вариантах оболочек см. табл. 4.3).
На рис. 4.10 приведены зависимости «W (0,25;0,25)−t » для не подкрепленной ребрами оболочки. Кривая 1 соответствует нагрузке Р = 20 103, кривая 2 – нагрузке Р = 25 103 , кривая 3 – Р = 30 103 , кривая 4 – Р = 35 103 , кривая 5 – Р = 40 103.
На рис. 4.11 даются зависимости «W (0,5; 0,5)−t » для оболочки, подкрепленной 18 ребрами. Кривая 1 соответствует нагрузке Р = 50 103 , кривая 2 – нагрузке Р = 60 103 , кривая 3 – Р =80 103 , кривая 4 – Р =100 103 , кривая 5 – Р =120 103 .
3 |
|
|
|
|
|
|
|
W |
5 |
4 |
3 |
2 |
|||
2,5 |
|
|
|
|
|
||
|
|
|
|
|
|
||
2 |
|
|
|
|
|
|
1 |
1,5
1
0,5
0
t
Рис. 4.10. Зависимость «W (0,25;0,25) −t » для гладкой оболочки
146 |
147 |
Математическое моделирование и расчет элементов строительных конструкций
3 |
5 |
4 |
3 |
2 |
||
|
|
|
|
|
||
|
|
|
|
|
|
|
2,5 |
|
W |
|
|
|
1 |
2
1,5
1
0,5
0
t
Рис. 4.11. Зависимость «W (0,5;0,5)−t » для оболочки, подкрепленной 18 ребрами
На рис. 4.12 для нагрузки Р = 30 103 представлен характер изме-
нения во времени прогиба W и напряжения σg = σ1 − Rbt σ3 , обозна-
Rb
ченного как функция интенсивности напряжений, для неподкрепленнойребрамиоболочки. Нарис. 4.13представленыаналогичныерезультатыдляоболочки,подкрепленной18ребрами,принагрузкеР =80 103 .
Как видно из этих рисунков, при развитии деформаций ползучести материала оболочки со временем происходит не только рост прогибовинапряжений,ноиперераспределениенапряженийпополю оболочки. Так, например, максимум напряжений из угловых точек оболочки смещается по всему наружному контуру оболочки.
Для рассматриваемых тонких оболочек (h = 600a ) потеря прочности происходит при t > tкр.
|
В табл. 4.7 приведены зависимости « |
|
−tкр»для варианта обо- |
||||||||
Р |
|||||||||||
лочки I. |
|
|
|
|
|
|
Таблица 4.7 |
||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
Оболочка без ребер |
Оболочка, подкрепленная 18 ребрами |
|||||||||
|
|
|
|
tкр, сут |
|
|
|
|
|
tкр, сут |
|
|
|
Р |
Р |
|
|
||||||
|
20 000 |
250 |
50 000 |
|
|
128 |
|
||||
|
|
|
|
|
|
|
|||||
|
25 000 |
135 |
60 000 |
|
|
86 |
|
||||
|
30 000 |
90 |
80 000 |
|
|
48 |
|
||||
|
35 000 |
70 |
100 000 |
|
|
32 |
|
||||
|
40 000 |
55 |
120 000 |
|
|
20 |
|
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
T, сут |
Прогиб |
|
Напряжение σg |
W |
|||
T = 0 |
|
|
|
T = 30
T = 70
T = 100
T = 118
Рис. 4.12. Функции прогиба W и σg для гладкой оболочки, при Р = 30 103
148 |
149 |
Математическое моделирование и расчет элементов строительных конструкций
T, сут |
Прогиб |
|
Напряжение σg |
W |
|||
T = 0 |
|
|
|
T = 20
T = 30
T = 40
T = 48
Рис. 4.13. Функции прогиба W и σg для оболочки, подкрепленной 18 ребрами, при Р =80 103
150
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ... |
|||||||
|
Длябольшей наглядностиэтирезультатырасчетапредставлены |
||||||
графически (рис. 4.14). |
|
|
|
|
|
||
|
140000 |
|
|
|
|
|
|
|
120000 |
|
|
|
|
|
|
|
100000 |
|
2 |
|
|
|
|
P |
80000 |
|
|
|
|
|
|
P |
60000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40000 |
|
|
|
1 |
|
|
|
20000 |
|
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
|
|
|
|
t |
|
|
|
|
Рис. 4.14. Зависимости «Р−tкр » для варианта оболочки I |
|
Кривая 1 на рис. 4.14 соответствует оболочке без ребер, кривая 2 – оболочке, подкрепленной 18 ребрами.
ДлянеподкрепленнойребрамиоболочкивариантаIIнарис.4.15 представленызависимости«W (0,5; 0,5)−t »,полученныеприразлич-
3
2,5 |
4 |
|
|
|
|
5 |
3 |
1 |
2 |
||
2 |
|
|
W 1,5
1
0,5
0
t
Рис. 4.15. Зависимости «W (0,5;0,5)−t » для гладкой оболочки варианта II
151
Математическое моделирование и расчет элементов строительных конструкций
ных значениях нагрузки: кривая 1 соответствует нагрузке Р =1500, кривая 2− Р = 2000, кривая 3 – Р = 2500, кривая 4 – Р = 2700, кри-
вая 5 – Р = 2800. На рис. 4.16 – зависимости «W (0,25;0,25)−t » для оболочки, подкрепленной 18 ребрами.
3
5 4 3
2,5
2
2
W 1,5
1
1
0,5
0
t
Рис. 4.16. Зависимости «W (0,25;0,25)−t » для оболочки варианта II, подкрепленной 18 ребрами
На рис. 4.16 кривая 1 соответствует нагрузке Р = 4000, кривая 2− Р = 5000, кривая 3 – Р = 6000, кривая 4 – Р = 7000, кривая 5 –
Р =8000.
На рис. 4.17 при Р = 2500 представлен характер изменения во
времени прогиба W и напряжения σg = σ1 − Rbt σ3 , обозначенного
Rb
как функция интенсивности, для неподкрепленной ребрами оболочки. На рис. 4.18 представлены аналогичные результаты расчета при
Р = 6000 для оболочки, подкрепленной 18 ребрами.
Из рисунков видно, что и для варианта оболочки II при развитии деформаций ползучести бетона оболочки со временем происходит рост и перераспределение напряжений по полю оболочки. Для
варианта оболочки II потеря прочности произойдет при t < tкр (примерно при 0,75tкр).
Глава 4. Расчетпологих ребристых оболочек при учете физической нелинейности ...
Т, сут |
Прогиб |
|
Напряжение σg |
W |
|||
T = 0 |
|
|
|
T = 25
T = 40
T = 60
T = 66
Рис. 4.17. Функции прогиба W и σg для гладкой оболочки, при Р = 2500
152 |
153 |