VERIFICATION OF SCAD MODELS OF REINFORCED CONCRETE COFFERED FLOOR BASED ON AN ANALYTICAL CALCULATION METHOD THAT TAKES INTO ACCOUNT SPANS AND STRUCTURAL RIGIDITY
Abstract and keywords
Abstract (English):
In accordance with the modern requirements of urban planning legislation, the design of building structures without the use of BIM technologies is impossible. Strength calculation on an electronic computer is carried out in software complexes implementing the finite element method, in which the forces or stresses calculated in the elements may turn out to be unreliable. There is a number of reasons for this. Analysis of the data of analytical and computer calculations of reinforced concrete coffered structures shows that the forces in the beams may differ significantly depending on the created finite element model and the geometry of the overlap. The work is aimed to find out the most accurate finite element model when calculating reinforced concrete coffered floor. The numerical experiment is based on the work of verification calculations performed in the computing system SCAD of a rectangular overlap with rectangular coffers, modeled by four finite element models. It is concluded that the model consisting of shell finite elements is the most accurate. It addition, the stress obtained is compared with the data of the well-known analytical method for calculating coffered ceilings, which is based on the beam analogy and takes into takes into account only the spans of the structure. In this paper, the analytical calculation of the caisson overlap is carried out both taking into account the spans of the structure and its orthogonal rigidity. The calculation is also performed in the computing system SCAD on a model consisting of rod finite elements of a T-section. The results allow to conclude that the most accurate finite elements models is the bar models. Shell finite element models in the example under consideration show understated results. When performing verification calculations on a computer by the finite element method, in order to confirm the reliability of the stress –strain state obtained, it is necessary to compare with the data of full-scale or model tests of structures. For complex, repeatedly statically indeterminate systems for which there is no an analytical solutions, other methods of studying the convergence of the obtained FEM results have an error.

Keywords:
coffered floors, floor stiffness, verification, the finite element model, SCAD computing complex
Text
Publication text (PDF): Read Download

Введение. Эффективным с конструктивной точки зрения и красивым по архитектуре является часторебристое перекрытие кессонного типа. Пример такой конструкции – павильон № 29 «Цветоводство и озеленение» ВДНХ в г. Москва. Строительство этого объекта осуществлено в 1969–1971 г.г. В настоящее время при строительстве ребристых перекрытий в нашей стране получают распространение опалубочные системы с извлекаемыми пустотообразователями, такие как: SKYDOME, HOLEDECK, ПОБЕДА.

В статью 1 Градостроительного кодекса РФ от 29.12.2004 № 190-ФЗ Федеральным законом от 27.06.2019 № 151-ФЗ введено понятие информационной модели объекта капитального строительства (BIM), принят
СП 333.1325800.2020 Информационное моделирование в строительстве.  Распоряжением Правительства РФ от 20.12.2021 № 3719-р утверждена «Дорожная карта» по использованию технологий информационного моделирования в строительной отрасли. Она касается различных министерств и ведомств, органов государственного строительного надзора, проектных, изыскательских и экспертных организаций. Одним из непростых вопросов
BIM моделирования является процесс создания конструктивной модели [1].

В работах [2, 3] рассматривается важный для проектировщиков вопрос о верификационном моделировании в вычислительном комплексе SCAD при расчете ребристой конструкции прямоугольного железобетонного кессонного перекрытия 9,0×11,55 м с кессонами 1,5×1,65 м (рис. 1). Сравниваются усилия – пролетные изгибающие моменты в балках короткого и длинного направлений, полученные при помощи четырех конечно-элементных моделей. В качестве эталона используется пример 15 работы [4, с 523] аналитического метода расчета. В соответствии с описанием, представленным в [2, 3] по первой расчетной схеме в ВК SCAD моделирование перекрытия заключается в том, что полка плиты задается оболочечными конечными элементами, а ребра стержневыми элементами, отнесенными от полки жесткой вставкой на величину Hc=hр2+hп2  [2, рис. 2.10). По второй схеме моделирование заключается в том, что полка задается оболочечными элементами, а ребра стержневыми элементами в виде тавра, центр тяжести которого совпадает со срединной линией полки. Свесы полок принимаются равными 3hп [2, рис. 2.12]. В третьей модели полка моделируется плитными конечными элементами (№ 11 КЭ), ребро моделируется оболочечными конечными элементами (№ 41 КЭ). Для учета пространственного положения соединение плиты с ребрами осуществляется при помощи жестких вставок [2, рис. 2.14]. В четвертой модели ВК SCAD полка и ребра моделируются оболочечными конечными элементами (№ 41 КЭ), их соединение осуществляется без использования жестких вставок [2, рис. 2.15]. В таблицах 1 работ [2, 3] приведены значения пролетных изгибающих моментов, полученных на ЭВМ, вычислены их отклонения от моментов аналитического метода расчета. По результатам анализа сделан вывод, что наиболее точной является оболочечная конечно-элементная модель. Стержневая модель МКЭ по изгибающим моментам приводит к существенным запасам, которые в рассматриваемом примере расчета доходят до 61,9 %.

Однако следует заметить, что пример аналитического метода расчета, данные которого приняты за эталонные значения, имеет две существенные ошибки. Первая ошибка заключается в неверном определении погонной нагрузки на балки qx и qy. В соответствии с известной аналитической теорией расчета железобетонных кессонных конструкций, основанной на балочной аналогии, нагрузка на балки qx и qy зависит только от пролетов перекрытия Lx и Ly [4–10]. В работе [11] доказано, что данный аналитический метод определения усилий в балках кессонных перекрытий неверен, он не учитывает ортогональную жесткость конструкции.

Второй ошибкой является неверное значение коэффициента в формулах определения пролетного изгибающего момента. В таблице VII.7 [4, с. 522] для пролетного изгибающего момента балок Мпр шарнирно-опертых по контуру конструкций принимается коэффициент 0,1, должен быть 0,125. Коэффициент 0,1 принят и в примере расчета кессонного перекрытия. Значения коэффициентов для опорных изгибающих моментов Моп в таблице отсутствуют, стоят прочерки.

Таким образом, пример 15 работы [4] не может служить эталоном для верификационных исследований конечно-элементных моделей ВК SCAD работ [2, 3].

Верификационными расчетами железобетонных кессонных и ребристых перекрытий кроме авторов работ [2, 3] занимались и другие исследователи [12–14]. Полученные данные также сравнивались с аналитическим методом, основанном на балочной аналогии. Отклонения компьютерных расчетов от аналитического метода в работе [12] составляют 50 %, в работе [13] 40 %. В работе [14] сравниваются продольные относительные деформации по высоте поперечного сечения балок, значения которых получены при помощи четырех конечно-элементных моделей расчета шарнирно опертой по контуру плиты 6,0×12,0 м. Первая модель представляет собой конструкцию, состоящую из стержневых конечных элементов – балок таврового сечения. Вторая модель состоит из оболочек, моделирующих плиту и стержней, выполняющих функцию ребер. Соединение элементов выполнено абсолютно жесткими телами. В третьей модели как плита, так и ребра представлены оболочками, соединенными абсолютно жесткими телами. Четвертая модель выполнена в виде массивного тела, состоящего из объемных конечных элементов. В соответствии с выводом проведенного численного эксперимента наименьшие деформации (значит и напряжения) получены в четвертой модели, состоящей из объемных КЭ и принятой автором за эталон. Следует акцентировать внимание, что наиболее близкой к четвертой модели оказалась самая простая первая модель, состоящая из стержней. Отклонения деформаций в растянутой зоне составляют +3,9 %, в сжатой
+9,2 %. Наибольшие отклонения от четвертой модели показала вторая модель, в растянутой зоне расхождения составляют +72 %, в сжатой +75,3 %. Отклонения третьей модели в растянутой зоне +44,5 %, в сжатой +42,7 %.

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

Материалы и методы. Численный эксперимент предусматривает сравнение изгибающих моментов, полученных МКЭ в ВК SCAD и аналитическим методом, учитывающим как пролеты конструкции, так и ее ортогональную жесткость.

Для обоснования выбора конечно-элементной модели обратим внимание на ряд работ, посвященных методу конечных элементов и созданию конечно-элементных моделей [15–20].

«Решения, получаемые по МКЭ, представляют в общем приближенные или аппроксимативные решения. Как и всегда, когда идет речь о приближенных решениях, ставится вопрос об их точности, устойчивости и конвергенции. Кроме того, с практической точки зрения важно знать точность решений, которая соответствует заданной надежности» [15]. «Какой бы подробной не была компьютерная модель, она всегда будет оставаться моделью, т.е. будет только моделировать различные свойства конструкции, причем, как правило, с различной степенью точности» [16]. Имеется два способа, на основании которых можно судить о точности полученного решения МКЭ, одним из которых является сравнение с задачей, имеющей известное аналитическое решение. Для многих реальных проблем, которые решаются, трудно найти соответствующую параллельную задачу, имеющую аналитическое решение [15]. Например, про расчет пластин в работе [21] сказано: «Задачи изгиба пластин сводятся к решению дифференциальных уравнений в частных производных, часто с переменными коэффициентами и при сложных краевых условиях. Точное решение известно только для отдельных сравнительно простых задач». В этом случае на ЭВМ для поиска сходимости результата применяется второй способ, который основан на многократном решении задачи с различной плотностью сетки конечных элементов [15]. Одним из важнейших моментов применения МКЭ, является выбор элементов модели и интерполяционных функций. Конечные элементы могут быть одномерными или линейными, двумерными или плоскими, трехмерными или пространственными (объемными). В работе [17] отмечается, что для конечных элементов стержневого типа постоянной жесткости по их длине для статической задачи постановка вопроса о сходимости МКЭ лишена смысла. О получении точных решений МКЭ при применении стержней указывается и в работах [15, 16]. В 1988 году отмечалось, что «В расчетах плит в инженерных конструкциях по МКЭ существенное место принадлежит выбору элементов, с помощью которых получается довольно точное и экономичное решение. И хотя этими вопросами занималось много авторов, общий ответ, который бы подходил ко всем отдельным случаям, дать невозможно» [15]. На современном этапе развития компьютерных расчетов в литературе 2022 года для ребристых конструкций можно встретить подобные выводы: «До сих пор не предложена оптимальная расчетная схема, с одной стороны, обладающая достаточной простотой для проведения инженерного анализа, а с другой стороны, позволяющая с большой точностью отражать особенность работы элементов перекрытия» [14]. «Вряд ли можно считать целесообразным использование одной очень подробной компьютерной модели, на основании которой специалист надеется рассмотреть все нюансы работы конструкции. Исследование работы плиты, подпертой ребрами, может быть проверено на основе стратов (моделей). Каждый из этих стратов имеет свои преимущества и недостатки. Исследование всех стратов даст обильную и добротную пищу для размышления» [16, с. 138]. «Следует учитывать, что ряд простых расчетных моделей может дать более точное знание, чем одна сложная (переусложненная) расчетная схема, результаты расчета которой достаточно трудно осмыслить [17]. О трудностях при создании пространственных моделей из объемных конечных элементов отмечается в работах [14, 15]. При проектировании мембранного покрытия центральной базы конькобежного спорта с искусственной конькобежной дорожкой в г. Коломна, представляющей собой овал размером в плане 199,89×110,25 м применялись тестовые компьютерные модели, созданные из различных конечных элементов с варьированием параметров степени дискретизации (крупность сетки, жесткость сталежелезобетонного опорного контура, характеристики упругого основания) [18].  Сопоставление различных конечно-элементных моделей показало, что моделирование опорного контура балочной моделью имеет наименьшую погрешность по сравнению с моделями, состоящими из плоских и объемных элементов. Также отмечается, что в последнее время расчетчики тяготеют к все большей детализации и учету большого числа подробностей. В действительности супердетализация приводит в лучшем случае к усложнению анализа работы конструкции, в худшем – к неизбежному накоплению погрешностей вычислений [18].

Метод конечных элементов позволяет практически полностью автоматизировать расчет стержневых систем. Он позволяет распространить принципы расчета стержневых систем на случай непрерывных тел и сложных конструкций [19]. Определение перемещений (деформаций и напряжений) по МКЭ в случае стержневых систем базируется на использовании технической теории растяжения, изгиба и кручения бруса, позволяющей выразить перемещения и напряжения в любом сечении бруса через узловые перемещения. Для двухмерного (пластина) или трехмерного (массив) сплошного тела эта задача может быть решена только приближенно [19].  Несмотря на свою простоту «…одномерные конечные элементы широко применяются в расчетах таких плоских и пространственных конструкций, как пластины и оболочки, где появляются линейные опоры, как прогоны, ребра и крайние балки» [15, с. 127]. В соответствии с инструкцией к версии 21 ВК SCAD++ [22] расчет армирования железобетонных элементов предусмотрен только для стержневых и пластинчатых элементов.

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

В связи с вышеизложенным, для изучения НДС ребер кессонных плит компьютерная модель представляет собой стержневую конструкцию, состоящую из балок таврового поперечного сечения. Высота балок 460 мм, ширина ребра 200 мм, ширина полки балок, установленных вдоль оси X - 1650 мм, вдоль оси Y-1500 мм, толщина полки 60 мм. Конструкция выполняется из бетона класса В25. Нагрузка на балки прикладывается по законам треугольника и трапеции [8, 9]. Поперечное сечение балок опорного контура принято аналогично пролетным балкам, но с половинным значением ширины полки. По причине отсутствия информации в работах [2, 3] об изменении модуля упругости бетона компьютерный расчет выполнялся как без коэффициентов редуцирования начального модуля упругости бетона, учитывающих ползучесть бетона и его трещинообразование [25], так и с коэффициентом 0,2 для всех балок, как пролетных изгибаемых, так и опорных по причине возникновения в них значительных крутящих моментов [23, 24]. Изменение модуля упругости бетона не привело к изменению усилий в конструкции, но увеличило ее прогиб в 5 раз, с 9,49 мм до 47,43 мм. При сравнении результатов кроме данных вновь рассчитанной модели ВК SCAD используются данные компьютерных расчетов четырех конечно-элементных моделей работ [2,3].

Основная часть. Кессонное перекрытие является плитой, опертой по контуру, подкрепленной снизу ребрами жесткости. В работе [21 с. 419] для пластины, усиленной в двух направлениях взаимно перпендикулярными равноотстоящими друг от друга ребрами, установленными с одной стороны, жесткость конструкции предлагается определять, как сумму, состоящую из цилиндрической жесткости плиты и относительной жесткости ребер:

Dx=E × d312 × (1-n2)+E'× I1d1  ,               (1)

Dу=E × d312 × (1-n2)+E''× I2d2  ,               (2)

где E, E’, Е” – модули упругости материала плиты и ребер; d – толщина плиты; n– коэффициент Пуассона; I1 и I2 – моменты инерции ребер жесткости, установленных вдоль осей X и Y, относительно линии, проходящей через центр тяжести таврового сечения; d1 и d2 – расстояния между ребрами.

Поэтому в данной работе жесткость перекрытия по направлениям X и Y определяется как сумма цилиндрической жесткости плиты и относительной жесткости отдельных ребер.

Первым действием определяем погонную нагрузку qx и qy, приходящуюся на балки и зависящую от пролетов перекрытия Lx, Ly и ортогональных жесткостей конструкции Dx, Dy.

Перепишем формулы расчета железобетонных кессонных перекрытий с учетом жесткостей перекрытия Dx и Dy:

qx=q × Ly4 × DxLx4 × Dy+Ly4 × Dx  ,                    (3)

qy=q × Lx4× DyLx4 × Dy+Ly4 × Dx  .                    (4)

Далее находим максимальные изгибающие моменты:

Mxmax= α1 qxa Lx2nx ,           (5)

Mymax= α2 qybLy2ny ,           (6)

где α1 и α2 – коэффициенты, зависящие от характера распределения нагрузки и вида опорных закреплений. В нашем случае α1=α2=0,125 ; a и b – шаг балок; nx  и ny  – коэффициенты, зависящие от геометрии перекрытия.

Перекрытие рассчитывается на равномерно-распределенную нагрузку q = 0,9 Т/м2, аналогично работам [2, 3, 4].

Произведем расчет по уточненной теории, учитывающей пролеты Lx, Ly, и жесткость перекрытия Dx, Dy.

 

 

Рис. 1. Схема кессонного перекрытия размером в плане 9,0×11,55 м с кессонами 1,5×1,65 м. К1, К2, К3, Д1, Д2, Д3 – рассчитываемые балки. [X, Y, Z] – связи, установленные в узлах балки опорного контура

 

 

Рис. 2. Схема приложения нагрузки на балки q, Т/м

 

 

Момент инерции отдельных балок вдоль осей X и Y относительно центра тяжести тавров:

Ix=Iy=bb × hb312+bb × hb× a2 ,                                                         (7)

где bb – ширина балки (ребра); hb – высота ребра; a – расстояние от центра тяжести ребра до центра тяжести тавра.

 

Ix=0,2 × 0,4312+0,2 × 0,4 × 0,1272=0,002357 м4 ,                                       (8)

Iy=0,2 × 0,4312+0,2 × 0,4 × 0,1222=0,002257 м4 .                                       (9)

Жесткость перекрытия вдоль оси X:

Dx=E × d312 × (1-n2)+E× Ixa=3058104 × 0,06312 × 1-0,22+0,0023571,65=4426 Тм  .                     (10)

Жесткость перекрытия вдоль оси Y:

Dy=E × d312 × (1-n2)+E × Iyb=3058104 × 0,06312 × 1-0,22+0,0022571,5=4659 Тм  .                     (11)

Составляющие общей нагрузки, приходящиеся на балки:

qx=q × Ly4 × DxLx4 × Dy+Ly4 × Dx=0,9 × 11,554 × 442694 × 4659+11,554 × 4426=0,648 Тм2  ,                                        (12)

qy=q × Lx4 × DyLx4 × Dy+Ly4 × Dx=0,9 × 94 × 465994 × 4659+11,554 × 4426=0,252 Тм2  .                                       (13)

Рассматриваем балки короткого направления, установленные вдоль оси X.

Балка К1.

Коэффициент пропорциональности, учитывающий расположение балки от опорного контура вдоль оси Y:

hy= yLy=1,6511,55=0,143 ,                                                                    (14)

nx= 165 × hy- 2 × hy3 + hy4= 165 × 0,143-2 × 0,1433+0,1434=0,44 .              (15)

Максимальный пролетный изгибающий момент:

MК1max= α × qx × a × Lx2 × nx=0,125 × 0,648 × 1,65 × 92 × 0,44=4,76 Тм .             (16)

Балка К2.

Коэффициент пропорциональности, учитывающий расположение балки от опорного контура вдоль оси Y:

hy= yLy=3,311,55=0,286 ,                                                                    (17)

nx= 165 × hy- 2 × hy3 + hy4= 165 × 0,286-2 × 0,2863+0,2864=0,787 .            (18)

Максимальный пролетный изгибающий момент:

MК2max= α × qx × a × Lx2 × nx=0,125 × 0,648 × 1,65 × 92 × 0,787=8,52 Тм .            (19)

Балка К3.

Коэффициент пропорциональности, учитывающий расположение балки от опорного контура вдоль оси Y:

hy= yLy=4,9511,55=0,429 ,                                                                    (20)

nx= 165 × hy- 2 × hy3 + hy4= 165 × 0,429-2 × 0,4293+0,4294=0,976 .            (21)

Максимальный пролетный изгибающий момент:

MК3max= α × qx × a × Lx2 × nx=0,125 × 0,648 × 1,65 × 92 × 0,976=10,57 Тм .           (22)

Рассматриваем балки длинного направления, установленные вдоль оси Y.

Балка Д1.

Коэффициент пропорциональности, учитывающий расположение балки от опорного контура вдоль оси X:

hx= xLx=1,59=0,167 ,                                                                    (23)

ny= 165 × hx- 2 × hx3 + hx4= 165 × 0,167-2 × 0,1673+0,1674=0,507 .             (24)

Максимальный пролетный изгибающий момент:

MД1max= α × qy × b × Ly2 × ny=0,125 × 0,252 × 1,5 × 11,552 × 0,507=3,2 Тм .       (25)

 

Балка Д2.

Коэффициент пропорциональности, учитывающий расположение балки от опорного контура вдоль оси X:

hx= xLx=39=0,333 ,                                                                    (26)

ny= 165 × hx- 2 × hx3 + hx4= 165 × 0,333-2 × 0,3333+0,3334=0,869 .             (27)

Максимальный пролетный изгибающий момент:

MД2max= α × qy × b × Ly2 × ny=0,125 × 0,252 × 1,5 × 11,552 × 0,869=5,48 Тм .       (28)

Балка Д3.

Максимальный пролетный изгибающий момент:

MД3max=α × qy × b × Ly2 × ny=0,125 × 0,252 × 1,5 × 11,552 × 1=6,3 Тм .                  (29)

 

 

Рис. 3. Эпюры изгибающих моментов М, Тм в балках модели ВК SCAD, состоящей из стержневых конечных элементов таврового поперечного сечения с шириной полки, равной расстоянию шага балок, перекрытия
размером в плане (Lx x Ly) 9,0
×11,55 м с кессонами 1,5×1,65 м, работающего на равномерно – распределенную нагрузку q = 0,9 Т/м2

 

 

Данные аналитического расчета и компьютерных моделей кессонного перекрытия представлены в таблице 1.

Выводы

  1. В формулах аналитического метода расчета прямых шарнирно – опертых по контуру ребристых кессонных железобетонных перекрытий, основанного на балочной аналогии, при определении составляющих общей нагрузки, необходимо использовать величины пролетов и ортогональную жесткость перекрытия, состоящую из суммы цилиндрической жесткости плиты и относительной жесткости ребер. Условием применения данной жесткости является одинаковое расстояние между параллельными балками.
  2. Значения изгибающих моментов в балках прямоугольного в плане шарнирно-опертого по контуру ребристого перекрытия с прямоугольными кессонами, определенные уточненным аналитическим способом, основанном на балочной аналогии, а также методом конечных элементов, реализованным в ВК SCAD с применением стержневой конечно-элементной модели таврового поперечного сечения с шириной полки, равной шагу балок, имеют близкие значения. Отклонения компьютерного расчета от аналитического метода составляют от +5,0 до +9,7 %, что можно объяснить погрешностями расчетов, как аналитического, так и компьютерного.

3.  Значения изгибающих моментов в балках прямоугольного в плане шарнирно-опертого по контуру ребристого перекрытия с прямоугольными кессонами, определенные уточненным аналитическим способом, основанном на балочной аналогии, а также методом конечных элементов, реализованным в ВК SCAD с применением модели, состоящей из полки, смоделированной оболочечными конечными элементами и ребер, представляющих собой стержни таврового поперечного сечения со свесами полок, равными тройной толщине плиты, имеют близкие значения. Максимальные отклонения компьютерного расчета от аналитического метода составляют -5,7 %, что можно объяснить погрешностями расчетов, как аналитического, так и компьютерного.

 

Таблица 1

Сравнение значений изгибающих моментов в балках шарнирно опертого по контуру
кессонного перекрытия размером в плане
(Lx x Ly) 9,0×11,55 м, полученные аналитическим
методом с учетом величин пролетов и жесткостей перекрытия, а также при помощи компьютерных моделей вычислительного комплекса
SCAD

 

Расчетная модель

Изгибающий момент в балках М, Тм

К1

К2

К3

Д1

Д2

Д3

Аналитическая балочная

с учетом величин пролетов и

жесткости перекрытия

 

4,76

100%

 

8,52

100%

 

10,57

100%

 

3,2

100%

 

5,48

100%

 

6,3

100%

Балочная модель

Тавровые балки - пространственный стержень тип КЭ 5

 

 

 

5,18

108,8%

 

 

 

9,33

109,5%

 

 

 

11,6

109,7%

 

 

 

3,36

105%

 

 

 

5,76

105,1%

 

 

 

6,64

105,4%

Модель 1 [1,2]

Полка – оболочечные КЭ, ребро – пространственный стержень тип КЭ 5, соединение элементов жесткими вставками

 

 

 

5,137

107,9%

 

 

 

9,155

107,5%

 

 

 

11,153

105,5%

 

 

 

3,709

115,9%

 

 

 

6,674

121,8%

 

 

 

7,691

122,1%

Модель 2 [1, 2]

Полка перекрытия – оболочечные КЭ, тавровое ребро – пространственный стержень тип КЭ 5, соединение элементов в центре тяжести

 

 

4,49

94,3%

 

 

8,1

95,1%

 

 

10,07

95,3%

 

 

3,03

94,7%

 

 

5,17

94,3%

 

 

5,95

94,4%

Модель 3 [1, 2]

Полка – плитные КЭ тип 11, ребро – оболочечные КЭ тип 41, соединение абсолютно жестким телом

 

 

3,634

76,3%

 

 

6,57

77,1%

 

 

8,177

77,4%

 

 

2,433

76%

 

 

4,151

75,8%

 

 

4,781

75,9%

Модель 4 [1, 2]

Полка и ребра – оболочечные КЭ

тип 41, соединение в узлах

 

 

3,815

80,2%

 

 

6,893

80,9%

 

 

8,578

81,2%

 

 

2,554

79,8%

 

 

4,36

79,6%

 

 

5,021

79,7%

 

4.  Значения изгибающих моментов в балках прямоугольного в плане шарнирно-опертого по контуру ребристого перекрытия с прямоугольными кессонами, определенные уточненным аналитическим способом, а также методом конечных элементов, реализованным в ВК SCAD с применением плитных и оболочечных конечных элементов, существенно отличаются. Максимальные отклонения компьютерного расчета от аналитического метода в рассматриваемом примере составили -24,2 %.

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

6.  Для всестороннего изучения работы сложной многократно статически неопределимой конструкции необходимо создавать несколько конечно-элементных моделей.

7.  При выполнении верификационных расчетов на ЭВМ методом конечных элементов для подтверждения достоверности полученного напряженно – деформированного состояния необходимо проводить сравнение с данными натурных или модельных испытаний конструкций. Для сложных многократно статически неопределимых систем, для которых отсутствуют известные аналитические решения, другие способы изучения сходимости получаемых результатов МКЭ имеют погрешность.

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

References

1. Perelmuter A., Kriksunov E. BIM in construction design - white spots. CAD-MASTER [BIM v stroitel'nom proektirovanii - belye pyatna. CAD-MASTER]. 2021. No. 1(95). Pp. 96-99. (rus)

2. Loskutov I.S. Monolithic reinforced concrete coffered floors [Monolitnye zhelezobetonnye kessonnye perekrytiya]. URL: https://dwg.ru/lib/2046 (rus)

3. Efimtseva E.E. Methods of modeling coffered ceilings // Actual problems of humanities and natural sciences. 2015. No. 11-5. URL: https://cyberleninka.ru/article/n/sposoby-modelirovaniya-kessonnyh-perekrytiy (accessed: 26.09.2022) (rus)

4. Linovich L.E. Calculation and construction of parts of civil buildings [Raschet i konstruirovanie chastey grazhdanskikh zdaniy]. Kyiv, Budivel'nik Publ., 1972. 644 p. (rus)

5. Zaliger R. Reinforced concrete, its calculation and design. Translation from German by professor P. Kamentsev [Zhelezobeton: ego raschet i proektirovanie. Perevod s nemetskogo pod red. P.Ya. Kamentseva]. M., Leningrad, GNTI Publ., 1931. 671 p. (rus)

6. Vakhnenko P.F., Khilobok V.G., Andreyko N.T., Yarovoy M.L. Calculation and construction of parts of residential and public buildings: designer's handbook. [Raschet i konstruirovanie chastey zhilykh i obshchestvennykh zdaniy: spravochnik proektirovshchika]. Kyiv, Budivel'nik Publ., 1987. 424 p. (rus)

7. Ulitskii I.I., Rivkin S.A., Samoletov M.V., Dykhovichnyi A.A., Frenkel' M.M., Kretov V.I. Reinforced concrete structures [Zhelezobetonnye konstruktsii]. Kyiv, Budivel'nik Publ., 1972. 992 p. (rus)

8. Ivanov-Dyatlov I.G. Reinforced concrete structures [Zhelezobetonnye konstruktsii]. Moscow, Leningrad, Ministerstvo kommunal'nogo khozyaystva RSFSR Publ., 1950. 296 p. (rus)

9. Karpukhin N.S. Reinforced concrete structures [Zhelezobetonnye konstruktsii]. Moscow, Gosudarstvennoe izdatel'stvo literatury po stroitel'stvu i arkhitekture Publ., 1957. 442 p. (rus)

10. Murashev V.I., Sigalov E.E., Baykov V.N. Reinforced concrete structures. General course [Zhelezobetonnye konstruktsii. Obshchiy kurs]. Moscow, Gosudarstvennoe izdatel'stvo literatury po stroitel'stvu i arkhitekture i stroitel'nym materialam Publ., 1962. 662 p. (rus)

11. Mozgolov M.V., Kozlova E.V. Use of analytical method for calculating stiffnesses of straight waffle slabs. Bulletin of Science and Research Center of Construction. 2022. 33(2). Pp. 122-138. doi:https://doi.org/10.37538/2224-9494-2022-2(33)-122-138 (rus)

12. Malakhova A.N. Monolithic caisson floors of buildings. Vestnik MGSU. 2013. (1). Pp. 79-86. (rus)

13. Shibaeva V.D. Investigation of the stress-strain state of monolithic coffered ceilings. Molodoi uchenyi [Young Scientist]. 2021. (16). Pp. 119-123. (rus)

14. Nikitin K.E., Kirsanov O.A. Comparative study of finite element methods of calculation of ribbed reinforced concrete floors. Construction mechanics of engineering structures and structures. 2022. 18(3). Pp. 242-254. doihttps://doi.org/10.22363/1815-5235-2022-18-3-242-254 (rus)

15. Sekulovich M. Finite element method [Metod konechnykh elementov]. Translated from Serbian by Yu.N. Zuev. Edited by V.Sh. Barbakadze. Moscow: Stroyizdat Publ.; 1993. 664 p. (rus)

16. Gorodetskiy A.S., Evzerov I.D. Computer models of structures [Komp'yuternye modeli konstruktsiy]. Kiev: Fact Publ.; 2005. 344 p. (rus)

17. Perel'muter A.V., Slivker V.I. Design models of structures and the possibility of their analysis [Raschetnye modeli sooruzheniy i vozmozhnost' ikh analiza]. Moscow, DMK Press Publ., 2007. 600 p. (rus)

18. Kancheli N.V., Batov P.A., Drobot D.Y. Realized membrane shells. Calculation, design and construction [Realizovannye membrannye obolochki. Raschet, proektirovanie i vozvedenie]. Moscow, Association of Construction Universities Publ. 2009. 120 p. (rus)

19. Perel'muter A.V. Conversations about construction mechanics [Besedy o stroitel'noy mekhanike]. Moscow, SCAD Soft Publ., ASV Publ., 2016. 304 p. (rus)

20. Mozgolov M.V., Kozlova E.V. Creation of a SCAD verification model for the design calculations of a reinforced-concrete waffle slab floor system. Bulletin of Science and Research Center of Construction. 2022. 32(1). Pp. 128-140. doi:https://doi.org/10.37538/2224-9494-2022-1(32)-128-140 (rus)

21. Vaynberg D.V., Vaynberg E.D. Calculation of plates [Raschet plastin]. Kyiv, Budivel'nik Publ., 1970. 436 p. (rus)

22. Karpilovskiy V.S., Kriksunov E.Z., Malyarenko A.A., Fialko S.Yu., Perel’muter A.V., Perel’muter M.A. SCAD Office. Version 21. The SCAD ++computing complex [SCAD Office. Versiya 21. Vychislitel'nyy kompleks SCAD ++]. Moscow; SKAD SOFT Publ.; 2015. 844 p. (rus)

23. Mozgolov M.V., Kozlova E.V. On the Unloading Action of Torque Moments in the Beams of Reinforced Concrete Casson Floors. Gradostroitel’stvo i arhitektura [Urban Construction and Architecture]. 2022. Vol. 12, No. 3. Pp. 11-20. DOI:https://doi.org/10.17673/Vestnik.2022.03.2 (rus)

24. Mozgolov M.V., Bryl S.V., Kozlova E.V. On the influence of the support contour beam on the stress -strain state of beams of straight coffered reinforced concrete floors. System technologies. 2022. No. 43. Pp. 31-40. doi:https://doi.org/10.55287/22275398_2022_2_31 (rus)

25. Ministerstvo stroitel'stva i zhilishchno-kommunal'nogo khozyaystva Rossiyskoy Federa-tsii, Federal'noe avtonomnoe uchrezhdenie

26. «Federal'nyy tsentr normirovaniya, standarti-zatsii i otsenki sootvetstviya v stroitel'stve» flat girderless reinforced concrete floors [Ploskie bezbalochnye zhelezobetonnye perekrytiya: Moscow, 2017. URL: https://www.faufcc.ru/upload/methodical_materials/mp60_2017.pdf (rus)


Login or Create
* Forgot password?