https://doi.org/10.21122/2227-1031-2023-22-2-158-167 УДК 624.072
Применение способа Жемочкина в нелинейном расчете железобетонных плит покрытия автомобильных дорог на упругом основании
Магистр Р. В. Кумашов1)
^ОДО «ЭНЭКА» (Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2023 Belarusian National Technical University, 2023
Реферат. Рассмотрена прямоугольная железобетонная плита с учетом ее физической нелинейности на линейно-упругом однородном основании под действием вертикальной внешней нагрузки. Анизотропия и неоднородность плиты обусловлены свойствами железобетона, а также образованием трещин от действия произвольной нагрузки в процессе эксплуатации. Нелинейную задачу решали способом Жемочкина с использованием итерационного алгоритма метода упругих решений Ильюшина. Для определения коэффициентов разрешающих уравнений способа Жемоч-кина применяли метод Ритца (определение прогибов плиты с защемленной нормалью) и решение Буссинеска (определение перемещений точек поверхности упругого полупространства). На первой итерации плиту рассчитывали как линейно-упругую, ортотропную и однородную, на последующих - как линейно-упругую, анизотропную и неоднородную на каждом участке Жемочкина. Прогибы срединной поверхности плиты от единичной силы определяли в виде ряда по первым пяти частным решениям Клебша. Выполнены экспериментальные и численные исследования. Последние - с помощью компьютерной программы MATHEMATICA. Полученные результаты показали, что предлагаемая методика расчета позволяет точно описать распределение осадок и реактивных напряжений под плитой. Верификацию методики статического нелинейного расчета прямоугольной железобетонной плиты с учетом ее физической нелинейности осуществляли путем сравнения результатов расчетов максимальных осадок и средних давлений под плитой, вычисленных с использованием предлагаемой методики, и результатов, полученных с помощью метода послойного суммирования и современных программных комплексов «Лира» и PLAXIS 3D.
Ключевые слова: физическая нелинейность, железобетонная плита, анизотропия плиты, неоднородность, способ Жемочкина, перемещения, упругое полупространство, защемленная нормаль, осадки, итерационный алгоритм, метод упругих решений, метод Ритца
Для цитирования: Кумашов, Р. В. Применение способа Жемочкина в нелинейном расчете железобетонных плит покрытия автомобильных дорог на упругом основании / Р. В. Кумашов // Наука и техника. 2023. Т. 22, № 2. С. 158-167. https://doi.org/10.21122/2227-1031-2023-22-2-158-167
Nonlinear Calculation of the Reinforced Concrete Road Pavement Slabs of Highways on the Elastic Basis by the Zhemochkin Method
R. V. Kumashov1)
1)ALC "ENEKA" (Minsk, Republic of Belarus)
Abstract. A rectangular reinforced concrete slab is considered taking into account its physical nonlinearity on a linearly elastic homogeneous base under the action of a vertical external load. The anisotropy and heterogeneity of the slab are due to the properties of reinforced concrete, as well as the formation of cracks from the action of an arbitrary load during operation. The nonlinear problem was solved by the Zhemochkin method using the iterative algorithm of the Ilyushin elastic solution method. The Ritz method (determining the deflections of a slab with a pinched normal) and the Boussinesq solution (determining the displacements of points on the surface of an elastic half-space) were used to determine the coefficients of the resolving equations of the Zhemochkin method. At the first iteration, the slab was calculated as linearly elastic,
Адрес для переписки
Кумашов Роман Владимирович ОДО «ЭНЭКА» просп. Независимости, 177-1 а, 220125, г. Минск, Республика Беларусь Тел.: +375 29 306-21-18 Gunner09@yandex.ru
Address for correspondence
Kumashov Roman V.
ALC "ENECA"
177-^, Nezavisimosty Ave.,
220125, Minsk, Republic of Belarus
Tel.: +375 29 306-21-18
Gunner09@yandex.ru
Наука
итехника. Т. 22, № 2 (2023)
orthotopic, and homogeneous; at subsequent iterations, it was calculated as linearly elastic, anisotropic, and inhomogeneous at each Zhemochkin site. The deflections of the middle surface of the slab from a unit force were determined as a series according to the first five particular Clebsch solutions. Experimental and numerical researches have been carried out. The latter - with the help of the MATHEMATICA computer program. The results obtained showed that the proposed calculation method allows one to accurately describe the distribution of settlements and reactive stresses under the slab. Verification of the methodology for static nonlinear calculation of a rectangular reinforced concrete slab, taking into account its physical nonlinearity, was carried out by comparing the results of calculations of maximum settlement and average pressures under the slab, obtained using the proposed methodology, and the results obtained using the layer-by-layer summation method and modern software systems Lira and PLAXIS 3D.
Keywords: physical nonlinearity, reinforced concrete slab, anisotropy of slab, heterogeneity, Zhemochkin method, displacements, elastic half-space, pinched normal, settlements, iteration algorithm, method of elastic solutions, Ritz method
For citation: Kumashov R. V. (2023) Nonlinear Calculation of the Reinforced Concrete Road Pavement Slabs of Highways on the Elastic Basis by the Zhemochkin Method. Science and Technique. 22 (2), 158-167. https://doi.org/10.21122/2227-1031-2023-22-2-158-167 (in Russian)
Введение
При эксплуатации дорог передача нагрузки на железобетонную плиту дорожного покрытия от колес автомобиля всегда будет вне оси симметрии конструкции. Также не исключена вероятность образования выбоин, воронок и других дефектов под основанием плиты, поэтому в ней будут проявляться изгибающие и крутящие моменты, поперечная сила. Данные факторы - основные из влияющих на несущую способность, деформативность и долговечность дорожных покрытий [1].
Существует множество методик расчета жесткостных и деформативных характеристик плит и балок на упругом основании. Однако ни одна из них не определяет физическую нелинейность плит на упругом основании в полной мере. В частности, нет методики расчета прямоугольных плит с трещинами на упругом основании. В некоторых работах описаны только предпосылки и допущения, которые могут способствовать созданию такой методики [2-7]. Это объясняется сложностью решения контактной задачи теории упругости для плиты на упругом основании.
Современные программные комплексы по расчету строительных конструкций широко применяют численные методы, в частности метод конечных элементов. Он основан на представлении реальной конструкции ее дискретной моделью и замене дифференциальных уравнений, описывающих напряженно-деформированное состояние сплошных тел, системой алгебраических уравнений. В качестве модели упругого основания в таких расчетных комплексах используются модели с постоянным или переменным коэффициентом постели. Данные методы не лишены недостатков [8-10]. Таким образом, эта отрасль исследований нелинейных
Наука
итехника. Т. 22, № 2 (2023)
задач требует дальнейшей разработки в плане создания общей методологии, формулировки и решения контактных задач для нелинейно-упругих конструкций, взаимодействующих с линейно-упругим однородным основанием при условиях, накладываемых практикой натурных инженерно-геологических изысканий.
Постановка задачи. Граничные условия
В статье рассматривается прямоугольная, физически нелинейная анизотропная неоднородная плита на линейно-упругом однородном основании, которая находится под действием вертикальной внешней нагрузки. Параметры плиты: длина 26, ширина 2а, высота к, жесткости в различных направлениях Оц, 012, 021, 022, Аз. Физическая нелинейность обусловлена нелинейностью материала, анизотропия и неоднородность плиты - образованием в ней трещин от действия произвольной нагрузки.
Линейно-упругое однородное основание моделируется как для линейно-упругого полупространства с постоянным модулем упругости (деформации) этого пространства Е0. Коэффициент Пуассона упругого пространства у0 принимается постоянным. Считается, что на контакте плиты с упругим основанием возникают только нормальные реактивные давления, силами трения пренебрегают. Для плиты справедливы гипотезы теории изгиба, в контактной зоне справедливо равенство осадок основания прогибам плиты.
За основу методики расчета плиты на упругом основании принимали смешанный метод строительной механики, описанный Б. Н. Же-мочкиным [11]. Плита в плане разбивалась на тхп участков Жемочкина (рис. 1). В свою очередь, каждый участок Жемочкина по высоте разбивался на к элементарных слоев (рис. 2).
За неизвестные принимали реактивные усилия Rj(x, y) в связи 7-го участка Жемочкина.
"о
Фо
Г яхт
Рис. 1. Основная система способа Жемочкина: u0, фо - линейное и угловое перемещения введенного защемления на плите; P, q - произвольные сосредоточенная и распределенная нагрузки; Xi,Xi - усилие в связи 1-го и i-го участков Жемочкина
Fig. 1. Main system of the Zhemochkin method: u0, ф0 - linear and angular displacements of the introduced pinching on the plate; P, q - arbitrary concentrated and distributed loads; X1, Xj - an effort in connection with the 1st and ith sections of Zhemochkin
Рис. 2. Расчетное поперечное сечение и распределение деформаций и напряжений изгибаемого железобетонного элемента: Aci, Asi - площадь элементарной площадки для бетона и арматуры; Z0 - координаты положения нейтральной линии; Zi - координата i-й элементарной площадки; gci, esi - относительная деформация бетона
и арматуры в i-й элементарной площадке; aci, asi - нормальное напряжение в бетоне и арматуре i-й элементарной площадки
Fig. 2. Calculated cross-section and distribution of deformations and stresses of a bent reinforced concrete element: Aci, Asi - area of an elementary platform for concrete and reinforcement; Z0 - neutral line position coordinates; Zi - coordinate of the ith elementary area; sci, esi - relative deformation of concrete and reinforcement in the ith elementary area; aci, asi - normal stress in concrete and reinforcement of the ith elementary area
Так как модуль деформации изменяется с высотой сечения Е(г) и аналитическое выражение данной зависимости неизвестно, то вычисление жесткостей выполняли численно [12]. Тогда выражения для жесткостей с учетом армирования плиты принимают вид:
п 2 А
о„ = к Ес,АЛ<г Е,,,—<г;
£ а
О =УЕ В Иг2 + кЕ •
г=1 ]=1 Ь
О21 = ЕЕеАЛ<,- +Е ; (1)
1=1 1=1 а
п 2 А
О =У Е В кг2 +У Е ,у, ] г2 •
22 ^с,^2,2 г су,г ¿^ ь,] , ¿¡¡у,г>
7=1
j=1
D33 = Ё Ec A,AZ<
2
cx ,7 '
7=1
где , , , гьу- положение нейтральной
линии для бетона и арматуры в двух направлениях; Аьх ], Аьу ] - верхнее и нижнее армирование плиты в двух направлениях; Ее г, Е,, - модуль упругости бетона и арматуры; а, Ь - размер плиты; к - толщина г-й площадки; п -количество элементарных площадок бетона и арматуры; В - параметр, рассчитываемый по формуле
1
V2 -1 V2 -1
V2 -1 V2 -1
0
0 1
2 0 + v).
V - коэффициент Пуассона для бетона.
Полученную статически неопределимую систему решали смешанным методом строительной механики [4]. Составляли канонические уравнения смешанного метода, выражающие условия совместности деформаций плиты и упругого основания. К этим уравнениям добавляли уравнения равновесия. Получили систему линейных алгебраических уравнений [13]:
Наука
итехника. Т. 22, № 2 (2023)
(2)
Ё I Ё 5гкКк + Ф0хУ7 + Ф0уХ + и0 + А1р = 0 7=1 V к=1
~Т*кУк + Мрх = 0;
к=1
-Ё^ЯкХк + Му = 0;
к=1
-ТКк + в = 0,
к=1
где и0, ф0х, ф0у - линейное и угловые перемещения введенного защемления на плите; в, Мрх, Мру - равнодействующая внешних
сил, действующих на плиту, и ее моменты относительно координатных осей; Як - неизвестное усилие в связях; 57к - коэффициент при неизвестных усилиях в связях; х, У - координаты центра участка Жемочкина; А1р - вертикальное перемещение 1-й точки основной системы; п - количество участков Жемочкина.
В ходе решения задачи будут получены следующие параметры: осадки плиты, реактивные давления в контактной зоне плиты с основанием, внутренние усилия в сечениях плиты, относительные деформации и напряжения в элементарных слоях участков Жемочкина.
Итерационный алгоритм расчета
В качестве диаграммы бетона, устанавливающей зависимость между напряжениями и деформациями, принимали нелинейную диаграмму состояния бетона с ниспадающей ветвью, рекомендованную в [14] (рис. 3).
На первой итерации при известном усилии от внешней нагрузки задавалось положение нейтральной поверхности сечения плиты в предположении упругой работы бетона и арматуры
Ё
Z0 x = Z0 y =
г=1
Ё EoA
(3)
где Е07 - начальный модуль упругости бетона или арматуры 1-й элементарной площадки; 7 - расстояние от выбранной оси элемента
Наука
итехника. Т. 22, № 2 (2023)
до центра тяжести 7-й элементарной площадки; А - площадь сечения 7-й элементарной площадки; п - количество элементарных площадок бетона и арматуры.
Рис. 3. Зависимость «напряжение - относительная деформация» для бетона: gc - относительная деформация бетона; gc1 - то же, соответствующая прочности бетона
на сжатие; ecu1, ecu1 - то же, соответствующая предельному сжатию и растяжению бетона; gct1 - то же, соответствующая прочности бетона на растяжение; fcm - средняя прочность на осевое сжатие бетона нормального веса; fctm - то же на осевое растяжение бетона нормального веса, установленная для проектирования конструкций; Ecm - модуль упругости бетона; <c - нормальное напряжение в бетоне
Fig. 3. Dependence "stress - relative deformation" for concrete: ec - relative concrete deformation; ec1 - the same, corresponding to compressive strength of concrete; ecu1, ectu1 - the same, corresponding to ultimate compression and tension of concrete; gct1 - the same,
corresponding to the tensile strength of concrete; fcm - average axial compressive strength of concrete of normal weight; fctm - average axial tensile strength of normal weight concrete specified for structural design; Ecm - modulus of concrete elasticity; < c - normal stress in concrete
Составляли канонические уравнения смешанного метода (2). Коэффициенты при неизвестных усилиях в связях в системе зависят от прогибов плиты в основной системе смешанного метода и осадок упругого основания и определяются по формуле [13]
§гк = Wik + Vik,
(4)
где wik - прогиб плиты в точке Х7 в основной системе от действия усилия Як = 1; У7к - перемещение точки приложения Х7 на границе ос-
i =1
нования от действия равномерно распределенного по участку усилия Як = 1.
Прогибы срединной поверхности плиты от единичной силы определяли в виде ряда по первым пяти частным решениям Клебша [15]
Ж (х, у ) = С
2 / \2 Х 1 ( У
-Сп
4 / \4 Х 1 ( У
С
2 ху
02
Ьа
2 / \2 х 1 ( У
(5)
С х С03 Ь
ЬГ - з (У
■с У
04
а
2 / \2 х 1 ( У
31 ь ] 1 а
где С00, С01, С02, С03, С04 - неопределенные коэффициенты.
Для определения коэффициентов использовали метод Ритца. Составляли функционал полной энергии плиты с защемленной нормалью и действующей силой Р(и, 7)
Л а Ь э=211
Г Я2
А,
а 2ж (х, у )
2
ах2
2
а
22
а 2ж (х, у) ау2
2
.а2ж(х,у) а2ж(х,у) ((2 + А: ~а^ + (6)
4а
33
ах
а 2ж (х, у)
дхду
йхйу - РЖ (и, 7),
где Аь а12, а21, а22, а33 - жесткости плиты, определяемые по формуле (1); Ж (х, у) - прогибы плиты в точке с координатами х, у; Р - единичная сила, действующая в точке с координатами и, 7; Ж(и,7) - прогибы плиты от
действия единичной силы Р.
Из условия минимума функционала составляли систему линейных алгебраических уравнений. В результате ее решения получали выражения для коэффициентов и соответственно значения прогибов срединной поверхности плиты. Перемещение точки М(х, у) поверхности упругого полупространства при загружении участка прямоугольной формы на этой поверхности равномерно распределенной нагрузкой
с равнодействующей, равной 1, после решения задачи Буссинеска [16] имеет вид [13]:
Жм ( х, У ) =
1 -У2
лЕ0А х
у
х - Ь + ^(х - Ь )2 + (у - й )2
:- а + ^(х - а )2 +(у - й )2
У - с
1п
х - а + ^(х - а)2 + (у - с)2 Ау х -Ь + ^(х -Ь)2 +(у - с)2
У -й +^(х -Ь)2 + (у - й)2
х-Ь
(7)
х 1п
х - а, +-1п
А у
У - с + д/(х - Ь )2 +(у - с )2 У - с + ^(х - а)2 + (у - с)2
у - й + ^(х - а)2 + (у-<
Ах = Ь - а; Ау = й - с,
где а, Ь, с, й - координаты границ участка Же-
мочкина; х, у - то же точки М; Е0 - модуль
деформации грунта; у0 - коэффициент Пуассона грунта.
После решения системы канонических уравнений (2) по найденным значениям реактивных усилий Як находили реактивное давление под плитой и распределение осадок «к. По величине «к определяли кривизну нейтральной поверхности плиты в двух взаимно перпендикулярных направлениях, а также относительное кручение и главные радиусы кривизны поверхности [2]. Относительные деформации слоя, отстоящего от нейтральной поверхности на расстоянии г7, вычисляли по формуле
8т О - ■
(8)
'1,2
где г12 - главный радиус кривизны.
Далее по диаграмме деформирования бетона с учетом найденных относительных деформаций вычисляли напряжения и секущие модули деформации в каждом элементарном слое, а по формуле (3) - уточненное положение центра изгиба сечения плиты г0х, г0у. Если напряжения растяжения в бетоне элементарной площадки превышали предельные значения, это
Наука
итехника. Т. 22, № 2 (2023)
свидетельствовало об образовании трещины в данной площадке. В последующих расчетах осевые жесткости этой элементарной площадки принимали равными нулю.
Найденные новые секущие модули деформаций вводили в новом расчетном цикле, в котором уточняли положение нейтральной поверхности плиты, жесткости и осадки плиты. Критерием окончания процесса последовательных приближений являлось сравнение осадок и главных радиусов кривизны на смежных этапах:
vn _ on—1
A—, = i--100% <£;
, —П
(9)
( 1 Y
AI =
( Л \
n—1
(1 >
Г
V , /
•100% <£, (10)
где —П, —П - осадка на текущей и предыдущей
( 1 г ( 1 г—1
итерации;
- главные радиусы
V i / V i /
кривизны на текущей и предыдущей итерации; Е, - критерий сходимости; i - номер участка Жемочкина; n - номер итерации.
После окончания итерационного процесса в результате решения системы канонических уравнений (2) по найденным значениям реактивных усилий Rk определяли окончательные величины реактивного давления под плитой и распределения осадок Sk.
Экспериментальные исследования
Изучение реального распределения осадок под железобетонной плитой проводили в натурных условиях при испытании плиты 2ПП30.18-30 серии 3.503.1-1. Методика испытаний и их результаты приведены в [17, 18]. По результатам натурных экспериментов построены изополя перемещений плиты при различных вариантах ее загружения. На рис. 4 представлены изополя основания плиты при центральном нагружении.
Наука
итехника. Т. 22, № 2 (2023)
0,400t-
0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6
4,55 4,80 5,05 5,30 5,55 5,80 6,05 6,30 6,55
Рис. 4. Изополя осадок основания плиты (мм) при центральном загружении силой F = 100 кН
Fig. 4. Settlement isofields for slab basis (mm) with its central loading by force F = 100 kN
Численные исследования и верификация результатов расчета
После проведения экспериментов выполняли численные исследования для центрального загружения плиты силой F = 100 кН [1]. Расчет производили для каждой ступени нагружения в несколько итераций: нулевая (расчет в упругой стадии), первая и последующие (расчет в нелинейной стадии). Для каждой ступени нагружения условия сходимости были выполнены на второй итерации. По результатам численных исследований на рис. 5 приведены изополя осадок основания плиты, вычисленные на второй итерации для последней стадии нагружения, а на рис. 6 - изополя давлений под плитой.
0,8 0,6 0,4 0,2 0 -0,2 -0,4 -0,6
I
-3,3
-3,4 -3,5
-3,6 -3,7
-1,0 -0,5
0,5
1,0
Рис. 5. Изополя осадок основания плиты (мм) при центральном ее загружении силой F = 100 кН
Fig. 5. Settlement isofields for slab basis (mm) with its central loading by force F = 100 kN
—
0
-1,0 -0,5 0 0,5 1,0
Рис. 6. Изополя давлений под плитой (кПа) при центральном загружении силой F = 100 кН
Fig. 6. Pressure isofields under slab (kPa) with its central loading by force F = 100 kN
Верификацию осуществляли путем сравнения результатов расчетов максимальных осадок и средних давлений под плитой, полученных с использованием усовершенствованной методики статического расчета железобетонных плит на упругом основании, и результатов, полученных с помощью метода послойного суммирования [19] и современных программных комплексов «Лира» и PLAXIS 3D.
Расчет методом послойного суммирования. Осадку при использовании расчетной схемы линейно деформируемого полупространства с ограничением глубины сжимаемой толщи по методу послойного суммирования определяли как сумму осадок элементарных слоев грунта в пределах сжимаемой толщи [19]
n с h
S=РЁ-Ет±, (11)
г =1 Ei
где ß - безразмерный коэффициент, равный 0,8; czpi - среднее дополнительное вертикальное нормальное напряжение в г-м слое основания вдоль вертикали, проходящей через центр подошвы фундамента, МПа; hi - толщина г-го слоя, на которые разбивается сжимаемая толща грунта, м; Et - модуль деформации i-го слоя грунта, МПа.
При расчете плиту принимали абсолютно жесткой: нагрузка на плиту составляла 100 кН, размеры плиты 3x1,75 м, давление под ней 19,05 кПа, осадка упругого основания под плитой 3,28 см. Схема распределения вертикальных напряжений под плитой показана на рис. 7.
Напряжения, кПа -50 -45 -40 -35 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25
Рис. 7. Схема распределения вертикальных напряжений под плитой: pzg - вертикальное природное давление; czp - дополнительное вертикальное нормальное напряжение; z - глубина от низа плиты
Fig. 7. Distribution diagram of vertical stresses under the slab: pzg - vertical natural pressure; czp - additional vertical normal stress; z - depth from the slab bottom
Для сравнения результатов расчетов вводили показатель отклонений, рассчитываемый по формуле
Д f = Lok. .100%, (12)
ft
где Af - относительная разница в значениях величины f вычисленной в каждом расчете; fs, ft - значения f вычисленные способом Жемочкина и рассматриваемым методом соответственно.
Максимальные значения осадки плиты, полученные с использованием способа Жемочки-на и методом послойного суммирования, составили 3,722 и 3,280 мм соответственно, показатель отклонения Af = 13 %. Максимальное, минимальное и среднее давления под плитой, полученные с применением способа Жемочкина, соответственно 41,969, 12,821 и 19,05 кПа. Среднее давление под плитой, вычисленное методом послойного суммирования, составило 19,05 кПа, показатель отклонения Af = 0.
Расчет с использованием программного комплекса «Лира». Исследуемую плиту численно рассчитывали в ПК «Лира». При постановке задачи ориентировались на следующие
Наука
итехника. Т. 22, № 2 (2023)
гипотезы и допущения: в зоне контакта плиты с упругим основанием возникают только нормальные напряжения, силы трения пренебрежительно малы. В расчет принимали упругие характеристики:
- для плиты: конечный элемент КЭ-41, модуль упругости Е = 31500 МПа, коэффициент Пуассона V = 0,167, толщина Н = 0,17 м;
- для слоистого основания в каждом упругом слое: КЭ-36, Ео в соответствии с табл. 1 из [17], Vо = 0,3.
На рис. 8 приведена расчетная модель плиты в ПК «Лира», на рис. 9 - результаты расчета.
Максимальные значения осадок плиты, полученные с использованием способа Жемоч-кина и в программном комплексе «Лира»
в линейной постановке, составили 3,722 и 3,576 мм соответственно, показатель отклонения Af = 4 %.
Расчет с использованием программного комплекса PLAXIS 3D (Demo). PLAXIS 3D представляет собой конечно-элементную программу, предназначенную для трехмерных расчетов напряжений, деформаций и прочности (устойчивости) сложных геотехнических объектов. Программа располагает широкими возможностями для решения различных задач, возникающих при проектировании геотехнических сооружений различного назначения. Позволяет моделировать работу системы конструкционных элементов и их взаимодействие с грунтовой средой с учетом нелинейного, реологического и анизотропного поведения грунтов.
Разрез 1-1
Основание КЭ-36
Плита КЭ-41
/
<! < ; i
<
<
a
i х 500 = 9000
i
< CO
i !
I
!-Av: ЙВ w я
я я ■73 ris 'Г Ж 7! Я W - :
Разрез 2-2
Основание КЭ-36
7 х 250 = 1750
Плита КЭ-41
.14 х 500 = 7000
1
\ J
1
/
X
i A
> 0 0 0 9
Й
i jC 0 0 5 х
CO
w
— - У - m M 35
14 х 500 = 7000
Рис. 8. Расчетная модель плиты в программном комплексе «Лира» Fig. 8. Calculation model of the slab in the Lira software package
/ N
/ \
I \ /
\ /
0,875 0,625 0,375
0,125 □ 3,5-4,0
-0,125 □ 3,0-3,5
-0,375 □ 2,5-3,0 -0,625 -0,875
-1,5
-1,1 -0,7 -0,3
0,1
0,5
0,9
1,3
Рис. 9. Распределение осадок основания плиты при расчете в программном комплексе «Лира» (цифры на рисунке - размеры плиты относительно ее центра)
Fig. 9. Settlement distribution of the slab basis in the calculation in Lira software package (numbers in Fig. - dimensions of the slab relative to its center)
Наука
итехника. Т. 22, № 2 (2023)
Расчет вели в линейной постановке. Модуль упругости грунта принимали переменным по высоте. Характеристики:
- плиты: Е = 31500 МПа, V = 0,167, Н = 0,17 м;
- грунта основания: модель МоЬг-Сои1ишЬ, удельный вес грунта выше уровня грунтовых вод уипа = 17,3 кН/м3, удельный вес грунта ниже уровня грунтовых вод уа = 20 кН/м3, модель деформации грунта Е = 4 МПа, приращение модуля деформации грунта на единицу глубины Е7пс = 31,08 МПа, глубина, с которой учитывается упрочнение грунта, - = -0,8 м, коэффициент Пуассона грунта V = 0,3, модуль сдвига грунта О = 1,538 МПа, одометрический модуль грунта ЕоеС1 = 5,385 МПа, удельное сцепление грунта с = 2 кПа, угол внутреннего трения грунта Ф = 33о, угол дилатансии грунта у = 3о.
Сетка конечных элементов показана на рис. 10, результаты расчета перемещений плиты и напряжений под ней - на рис. 11, 12 соответственно.
Рис. 10. Сетка конечных элементов Fig. 10. Finite element mesh
[1C23m]
1 -2,00 -3,1C
i-3,20 -3,30 -3,00 -3,50
Z
То1а1 displacements uz (scaled up 200 times)
Maximum value = -2,250-10-3 m (Element 5 at Node 11) Maximum value = -3,023-10-3 m (Element 9 at Node 13)
Рис. 11. Перемещения плиты Fig. 11. Slab displacements
Максимальные значения осадок плиты, полученные с использованием способа Жемоч-кина и в программном комплексе PLAXIS 3D в линейной постановке, составили 3,722 и 3,023 мм соответственно, показатель отклонения Af = 9 %. Максимальное, минимальное и среднее давления под плитой, вычисленные с помощью способа Жемочкина, составили соответствен-
но 41,969, 12,821 и 19,048 кПа, а в программном комплексе PLAXIS 3D - 27,0, 11,0 и 16,9 кПа, показатель отклонения А/ = 12,7 %.
Сartesian total stress (scaled up 0,0200 times)
Maximum value = 4,654 kN/m2 (Element 131 at Node 5171) Maximum value = -46,54 kN/m2 (Element 540 at Node 1509)
Рис. 12. Напряжения под плитой Fig. 12. Stresses under slab
ВЫВОДЫ
1. Благодаря проведенным исследованиям усовершенствована методика статического расчета прямоугольных железобетонных плит покрытия автомобильных дорог с учетом их физической нелинейности итерационным способом Жемочкина, принимая во внимание тре-щинообразование.
2. Сравнение результатов экспериментальных и численных исследований показало, что способ Жемочкина, в общем, позволяет довольно точно описать распределение осадок и реактивных давлений под плитой.
3. Верификация предлагаемой методики статического расчета по требованиям ТНПА с использованием современных программных комплексов «Лира» и PLAXIS 3D, а также с учетом результатов экспериментальных исследований показала достаточно точное совпадение теоретических и практических значений осадок и напряжений под плитой.
ЛИТЕРАТУРА
1. Семенюк, С. Д. Экспериментальные исследования осадок железобетонных плит покрытия автомобильных дорог и оценка точности расчетной методики / С. Д. Семенюк, Р. В. Кумашов // Проблемы современного бетона и железобетона: сб. науч. тр. / Ин-т БелНИИС; редкол.: О. Н. Лешкевич [и др.]. Минск, 2020. Вып. 12. С. 185-208.
2. Соломин, В. И. Методы расчета и оптимальное проектирование железобетонных фундаментных конструкций / В. И. Соломин, С. Б. Шматков. М.: Стройиздат, 1986. 206 с.
3. Гвоздев, А. А. Расчет несущей способности конструкций по методу предельного равновесия / А. А. Гвоздев. М.: Стройиздат, 1949. 280 с.
4. Ржаницин, А. Р. Строительная механика / А. Р. Ржа-ницин. М.: Высш. шк., 1991. 439 с.
5. Ржаницин, А. Р. Предельное равновесие пластинок и оболочек / А. Р. Ржаницин. М.: Наука. Главная редакция физико-математической литературы, 1983. 288 с.
Наука
итехника. Т. 22, № 2 (2023)
6. Клепиков, С. Н. Взаимодействие конструкций с основанием: дис. ... д-ра техн. наук: 05.23.01 / С. Н. Клепиков; Госстрой СССР, Научно-исследовательский институт строительных конструкций (НИИСК). Киев, 1970. 329 с.
7. Семенюк, С. Д. Железобетонные пространственные фундаменты жилых и гражданских зданий на неравномерно деформированном основании / С. Д. Семе-нюк. Могилев: Белор.-Рос. ун-т, 2003. 269 с.
8. Маковкин, Г. А. Применение МКЭ к решению задач механики деформируемого твердого тела / Г. А. Маковкин, С. Ю. Лихачева. Н. Новгород: Изд-во ННГАСУ, 2012. Ч. 1. 71 с.
9. Лебедев, А. В. Численные методы расчета строительных конструкций / А. В. Лебедев. СПб.: СПбГАСУ, 2012. 55 с.
10. Ильин, В. П. Численные методы решения задач строительной механики / В. П. Ильин, В. В. Карпов, А. М. Масленников. Минск: Вышэйш. шк., 1990. 346 с.
11. Жемочкин, Б. Н. Практические методы расчета фундаментных балок и плит на упругом основании / Б. Н. Же-мочкин, А. П. Синицын. М.: Госстройиздат, 1962. 240 с.
12. Пискунов, Н. С. Дифференциальное и интегральное исчисления для втузов / Н. С. Пискунов. 13-е изд. М.: Наука. Главная ред. физ.-мат. лит., 1985. Т. 1. 432 с.
13. Босаков, С. В. Статические расчеты плит на упругом основании / С. В. Босаков. Минск: БНТУ, 2002. 128 с.
14. Бетонные и железобетонные конструкции: СП 5.03.012020. Введ. 16.11.2020. Минск: Минстройархитектуры, 2020. 178 с.
15. Тимошенко, С. П. Пластинки и оболочки / С. П. Тимошенко, С. Войновский-Кригер. М.: Физ.-мат. изд-во, 1963. 536 с.
16. Boussinesq, J. Applications des Patentiels a L'etude de L'equilibre et du Movement des Solides Elastiques / I. Boussinesq. Paris: Gauthiers-Villars, 1885. 721 p.
17. Семенюк, С. Д. Железобетонные плиты покрытия автомобильных дорог на упругом полупространстве / С. Д. Семенюк, Р. В. Кумашов // International Journal for Computation Civil and Structural Engineering. 2018. T. 14, Ыо 2. С. 149-157. https://doi.org/10.22337/2587-9618-2018-14-2-149-157.
18. Кумашов, Р. В. Методика и результаты полевых испытаний железобетонной плиты покрытия автомобильных дорог 2ПП30.18-30 / Р. В. Кумашов // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. Могилев: Белор.-Рос. ун-т, 2018. С. 298-299.
19. Фундаменты плитные. Правила проектирования: ТКП 45-5.01-67-2007 (02250). Введ. 01.09.2007. Минск: Минстройархитектуры, 2018. 94 с.
Поступила 02.08.2022 Подписана в печать 11.10.2022 Опубликована онлайн 31.03.2023
REFERENCES
1. Semenyuk S. D., Kumashov R. V. (2020) Experimental Research of the Displacements of the Reinforced Concrete Slabs of Highways Pavement and Assessment of the Accuracy of the Calculation Method. Problemy Sovremenno-go Betona i Zhelezobetona: Sb. Nauch. Tr. [Problems of Modern Concrete and Reinforced Concrete: Collection of Scientific Papers]. Minsk, Iss. 12, 185-208 (in Russian).
2. Solomin V. I., Shmatkov S. B. (1986) Calculation Methods and Optimal Design of Reinforced Concrete Foundation Structures. Moscow, Stroyizdat Publ. 206 (in Russian).
3. Gvozdev A. A. (1949) Calculation of the Bearing Capacity of Structures by the Method of Limit Equilibrium. Moscow, Stroyizdat Publ. 280 (in Russian).
4. Rzhanitsyn A. R. (1991) Structural Mechanics. Moscow, Vysshaya Shkola Publ. 439 (in Russian).
5. Rzhanitsyn A. R. (1983) Limit Equilibrium of Plates and Shells. Moscow, Nauka Publ., Main Office of Physical and Mathematical Literature. 288 (in Russian).
6. Klepikov S. N. (1970) Interaction of Structures with the Base. Kiev, State Committee for Construction of the USSR, Scientific Research Institute of Building Constructions (NlISK). 329 (in Russian).
7. Semenyuk S. D. (2003) Reinforced Concrete Spatial Foundations of Residential and Civil Buildings on an Unevenly Deformed Base. Mogilev, Belarusian-Russian University. 269 (in Russian).
8. Makovkin G. A., Likhacheva S. Yu. (2012) Application of FEM to Solving Problems of Solid Mechanics. Part 1. Novgorod, Nizhny Novgorod State University of Architecture and Civil Engineering (NNGASU). 71 (in Russian).
9. Lebedev A. V. (2012) Numerical Methods for Calculating Building Structures. Saint-Petersburg, Saint Petersburg University of Architecture and Civil Engineering (SPbGASU). 55 (in Russian).
10. Ilyin V. P., Karpov V. V., Maslennikov A. M. (1990) Numerical Methods for Solving Structural Mechanics Problems. Minsk, Vysheyshaya Shkola Publ. 346 (in Russian).
11. Zhemochkin B. N., Sinitsyn A. P. (1962) Practical Methods for Calculating Foundation Beams and Slabs on an Elastic Foundation. Moscow, Gosstroyizdat Publ. 240 (in Russian).
12. Piskunov N. S. (1985) Differential and Integral Calculus for Higher Technical Educational Institutions. Vol. 1. 13th ed. Moscow, Nauka Publ., Main Office of Physical and Mathematical Literature. 432 (in Russian).
13. Bosakov, S. V. (2002) Static Calculations of Slabs on an Elastic Foundation. Minsk, Belarusian National Technical University. 128 (in Russian).
14. SP [Construction Regulations] 5.03.01-2020. Concrete and Reinforced Concrete Structures. Minsk, Publishing House of Ministry of Architecture and Construction, 2020. 178 (in Russian).
15. Timoshenko S. P., Voynovsky-Kriger S. (1963) Plates and Shells. Moscow, Fizmat Publ. 536 (in Russian).
16. Boussinesq J. (1885) Applications des Patentiels a L'etude de L'equilibre et du Movement des Solides Elastiques. Paris, Gauthiers-Villars. 721 (in French).
17. Semenyuk S. D., Kumashov R. V. (2018) Reinforced Concrete Coating Plates of Highways on an Elastic HalfSpace. International Journal for Computation Civil and Structural Engineering, 14 (2), 149-157. https://doi.org/ 10.22337/2587-9618-2018-14-2-149-157.
18. Kumashov R. V. (2018) Methodology and Results of Field Tests of a Reinforced Concrete Slab for Road Pavement 2PP30.18-30. Materialy, Oborudovanie i Resursosbere-gayushchie Tekhnologii: Materialy Mezhdunar. Nauch. -Tekhn. Konf. [Materials, Equipment and Resource-Saving Technologies: Proceedings of Scientific and Technical Conference]. Mogilev, Belarusian-Russian University, 298-299 (in Russian).
19. TKP [Technical Code of Common Practice] 45-5.01-672007 (02250). Slab Foundations. Design Rules. Minsk, Publishing House of Ministry of Architecture and Construction, 2018. 94 (in Russian).
Received: 02.08.2022 Accepted: 11.10.2022 Published online: 31.03.2023
Наука
итехника. Т. 22, № 2 (2023)