, Россия
УДК 69 Строительство. Строительные материалы. Строительно-монтажные работы
ГРНТИ 67.11 Строительные конструкции
ОКСО 08.03.01 Строительство
ББК 385 Строительные конструкции
ТБК 5414 Строительные конструкции
BISAC TEC009020 Civil / General
Наличие сведений о физико-механических свойствах фибробетона позволяет вести расчет конструкций на его основе с надлежащим приближением к реальности. Этому способствует переход от линейно упругой постановки задачи к учету физической нелинейности. В рамках программных комплексов, основанных на конечно-элементной модели деформируемого тела, приемлем итерационный процесс, обусловленный сходимостью решения. Предложенный метод дополнительных нагрузок в конечно-элементной интерпретации оставляет в силе накопленный опыт автоматизированного проектирования. Итерационная процедура расчета содержит циклы с постепенным приближением физических модулей к их истинному значению в результате сходимости решения задачи. Вычисление модулей ведется на основе зависимостей механики деформируемого твердого тела, что придает постановке задачи фундаментальный характер. Эффект предложенной методики демонстрируется на примере расчета заглубленного цилиндрического резервуара из стеклофибробетона для жидкости
фибробетонные конструкции, физическая нелинейность, конечно-элементная интерпретация метода дополнительных нагрузок
Введение. Успехи в области повышения активности минеральных вяжущих веществ, в первую очередь цемента, позволили значительно повысить предел прочности бетонов на сжатие. В то же время прочность на растяжение оставалась во много раз ниже прочности на сжатие. Это несоответствие в некоторой мере сглаживается использованием волокон в качестве арматуры. Наряду с применением металлических элементов большое распространение получили стеклянные и углеродные волокна [1–11].
Фибробетон, содержащий те или иные волокна, как и бетон, деформируется по нелинейному физическому закону. Механические характеристики (модуль продольной деформации, коэффициент поперечной деформации, модуль сдвига) являются переменными величинами и зависят от деформации. Их значения определяются на основе стандартных испытаний. Автором проведены такого рода исследования для цилиндрических элементов из стеклофибробетона. Серии отличались длиной волокна (20 мм, 35 мм, 50 мм) и процентом армирования (1,5 %, 3 %, 4,5 %) в различных сочетаниях. Результаты испытаний представлены в табличной форме [12].
Для определения прочностных характеристик стеклофибробетона (предела прочности на растяжение, сжатие) также использовалась стандартная методика с применением образов в виде кубиков и призм.
Накопленный материал о физико-механических свойствах фибробетона [13, 14] позволяет вести расчет конструкций на его основе с надлежащим приближением к реальности. Этому способствует переход от линейно упругой постановки задачи к учету физической нелинейности. В рамках программных комплексов, основанных на конечно-элементной модели деформируемого тела, приемлем итерационный процесс, обусловленный сходимостью решения.
Материалы и методы. Известно, что при линейном физическом законе прямая задача содержит в качестве неизвестных 15 компонентов напряжений, деформаций и перемещений. Для их определения используются 15 уравнений по трем законам деформирования твердого тела. При наличии нелинейных деформаций в прямой задаче появляются дополнительные
неизвестные – интенсивности касательных напряжений 
и деформаций сдвига 
, а дополнительными уравнениями являются выражение: и зависимость между и [15].
Возможности современной вычислительной техники привели к разработке общих численных алгоритмов решения нелинейных задач. Одним из наиболее эффективных способов построения таких алгоритмов является линеаризация. Рассматриваемые ниже их варианты предложены российскими учеными А.А. Ильюшиным и И.А. Биргером.
При решении в перемещениях в дифференциальных уравнениях и граничных условиях слагаемые, возникающие из-за отклонения от закона Гука, переносятся в правые части уравнений и считаются известными. Таким образом, получаем систему уравнений линейной теории деформирования относительно компонентов перемещений, но с дополнительными объемными и поверхностными силами. В первом приближении полагаем эти дополнительные нагрузки равными нулю и решаем физически линейную задачу. Найдя перемещения, вычисляем интенсивность , позволяющую определить нелинейные члены в уравнениях Ламе и граничных условиях и решить вторичную физически линейную задачу с дополнительными нагрузками, отличными от нуля и т.д. Такого рода алгоритм получил название метода дополнительных нагрузок.
При решении в напряжениях используется метод дополнительных деформаций. Дифференциальные уравнения равновесия и граничные условия остаются без изменения. Уравнения сплошности вследствие наличия нелинейных членов в выражении физического закона будут содержать дополнительные слагаемые, которые можно интерпретировать как дополнительные деформации и определять их последовательными приближениями.
Имеет место также метод переменных параметров упругости. Систему уравнений представляют в форме уравнений линейной теории с переменными «параметрами упругости» 
(
отражает процесс деформирования), которые в процессе последовательных приближений вычисляются по деформациям.
Природа физической нелинейности фибробетона исходит из механических свойств бетона, проявляющихся в диаграмме «напряжение-деформация», явлении дилатации. Условием создания алгоритма расчета конструкций является наличие экспериментальных данных, касающихся модулей упругости первого и второго рода, коэффициента Пуассона, пределов прочности на растяжение и сжатие при характерных длинах волокон и процентах армирования по массе.
Построение расчета проведено на основе метода конечных элементов (МКЭ), используемого в инженерной практике для решения задач различного типа. Матрица жесткости конечного элемента содержит матрицу упругости D как характеристику связи между напряжениями и деформациями [15]:
где G – модуль сдвига, 
– коэффициент поперечной деформации (коэффициент Пуассона).
Для физической нелинейности характерны переменные величины этих характеристик. Диапазон изменения коэффициента поперечной деформации мал, в связи с чем, его часто принимают неизменным. Фактический модуль сдвига имеет выражение:
(2)
где и 
– интенсивность деформаций сдвига для фактического состояния и ее предельная величина для рассматриваемого вида напряженного состояния, то есть
(3)
где 
– второй инвариант девиатора деформаций, равный
(4)
где 
элементы тензора деформаций.
При трещинообразовании в фибробетоне имеет место увеличение объема тела – дилатация 
где 
– модуль дилатации, зависящий от свойств бетона, типа и содержания армирующих волокон. При чистом сдвиге предельная объемная деформация 
зависит от предельной интенсивности деформации сдвига [16]:
(5)
откуда вытекает значение модуля:
(6)
Линейные деформации дополняются величиной 
. Сохраняется предположение о пропорциональности и соосности девиаторов напряжений и деформаций.
Величина зависит от инвариантов тензора 
и девиаторов 
, 
деформаций, в которых находят отражение прочность фибробетона на сжатие 
и растяжение 
Эти характеристики, как и модуль , зависят от свойств бетона, типа и содержания армирующих волокон. В качестве математической модели можно использовать аналогичную модель для бетона с заменой характеристик 
и 
на и 
[12].
Основная часть. Как следует из изложенных выше положений о физической нелинейности, синтезирующее уравнение метода конечных элементов содержит вместо постоянной матрицы жесткости 
ее переменный аналог 
так, что
(7)
где 
– корректирующая матрица, равная
(8)
где 
– корректирующая матрица отдельного конечного элемента, число которых e.
Эта предпосылка заставляет заменить решение упомянутого уравнения МКЭ итерационным решением рекуррентного уравнения:
(9)
или
(10)
где 
– s-e приближение q;
– фиктивная дополнительная нагрузка, равная
(11)
Сходимость итерационного процесса подтверждена решениями типовых задач. Критерий сходимости имеет энергетическую природу:
, (12)
где 
– потенциальная энергия деформации конечного элемента; 
– средняя величина для e величин ; 
– некоторое малое положительное число.
Алгоритм решения физически нелинейной задачи включает: начальный модуль (модуль упругости) и модуль дилатации; формирование матрицы жесткости для отдельных элементов; формирование матрицы жесткости для тела и вектора узловых сил; решение уравнений МКЭ; определение модулей по формуле (2); формирование новых матриц жесткости конечных элементов; определение ; формирование ; решение уравнения (10); проверка удовлетворения критерию (12) и другим неучтенным ограничениям (по прочности, деформативности и т. д.). При необходимости следующий цикл начинается с определения модулей упругости по формуле (2).
Изложенная теория, касающаяся определения напряженно-деформированного состояния тела при наличии физической нелинейности, использована при расчете цилиндрических резервуаров для нефтепродуктов. При расчетах бетонных резервуаров без учета влияния растяжения бетона допускалось раскрытие трещин до величины 0,2 мм, а при учете влияния растяжения бетона – до 0,1 мм [17, 18]. Включение стекловолокон в бетон увеличивает сопротивляемость трещинообразованию, и допускаемую величину раскрытия трещин можно уменьшить до 0,05 мм.
В качестве примера рассмотрен резервуар, моделью которого является шарнирно опертая оболочка в грунте радиусом 4 м, длиной 6 м, толщиной стенок 0,05 м. Модуль упругости стеклофибробетона 2∙104 МПа, коэффициент Пуассона – 0,17. Толщина прилегающего деформируемого слоя основания 5м, модуль деформации грунта 60 МПа, коэффициент Пуассона 0,3. Давление на боковую поверхность оболочки составляет
500 кН/м.
Исходные данные взяты из решения задачи в линейно-упругой постановке [17]. Учет физической нелинейности потребовал трех приближений расчета. Перемещения срединной поверхности оболочки увеличились на 10−12 %, а изгибающие моменты на 14–16 %. В результате армирование стекловолокнами увеличено на 1,1 %.
Рассмотрим частный случай кручения цилиндрического бруса круглого поперечного сечения радиуса R из материала с билинейной диаграммой сдвига (рис. 1).
Установим зависимость между крутящим моментом Т и кривизной кручения. При 
брус имеет внутреннюю зону деформирования, связанную с участком диаграммы 0–1, и наружную зону, связанную с участком 1–2. Граница между этими зонами определяется радиусом
. (13)
Рис. 1. Диаграмма сдвига: 
деформация сдвига; τ – касательные напряжения;
arctg G (G – модуль сдвига), 
arctg 
( – модуль упрочнения)
Выражение крутящего момента имеет вид
Т
(14)
В данном случае
(15)
или
(16)
В итоге получаем
Т
(17)
Предположим, что пределу несущей способности бруса соответствует появление на контуре сечения напряжений 
.
При этом
(18)
При замене переменной r на 
формула (17) принимает вид
Т
(19)
Формула (19) выгодно отличается от формулы (17) своей непосредственной связью с диаграммой сдвига.
Кроме рассмотренной физической нелинейности существуют также и другие виды нелинейности. Геометрическая нелинейность имеет место, во-первых, когда перемещения конструкции вызывают значительные изменения в ее геометрии, так, что уравнения равновесия приходится составлять с учетом изменения формы и размеров конструкции (по деформированной схеме), во-вторых, когда зависимость между перемещениями и деформациями нелинейная.
Конструктивная нелинейность возникает вследствие конструктивных особенностей системы, вызывающих изменение расчетной схемы в процессе деформирования (например, изменяются условия закрепления). Она имеет место в конструкциях с односторонними связями, встречается также в контактных задачах.
При наличии физической и одного из указанных видов нелинейности в комплексе уравнений для решения такого рода задач появляются дополнительные, специфические уравнения, например, новые зависимости между деформациями и перемещениями, уравнения для условий контакта и другие. При этом остаются в силе предложенные в статье методы решения физически нелинейной задачи.
Выводы. Учет физической нелинейности при расчете конструкций из фибробетона позволяет в определенной мере приблизить физическую модель к реальным обстоятельствам. Предложенный метод дополнительных нагрузок в конечно-элементной интерпретации оставляет в силе накопленный опыт автоматизированного проектирования [19]. Итерационная процедура расчета содержит циклы с постепенным приближением физических модулей к их истинным значениям в результате сходимости решения задачи. Вычисление модулей ведется на основе зависимостей механики деформируемого твердого тела, что придает постановке задачи фундаментальный характер. Эффект предложенной методики демонстрируется на примере расчета заглубленного цилиндрического резервуара из стеклофибробетона для жидкости.
1. Волков И.В. Проблемы применения фибробетона в отечественном строительстве // Строительные материалы. 2004. № 6. С. 12–13.
2. Волков И.В., Газин Э.М., Бабекин В.В. Инженерные методы проектирования фибробетонных конструкций // Бетон и железобетон. 2007. № 4. С. 20–22.
3. Олюнин П.С. Дисперсное армирование цементных композитов полимерными волокнами// Бетон и железобетон. 2009. № 1. С. 21–24.
4. Проектирование и основные положения технологий производства фибробетонных конструкций. (ВСН 56-97). М.: НТУ НИЦ «Строительство», 1997. 91 с.
5. Рабинович Ф.Н. Композиты на основе дисперсно армированных бетонов. Вопросы теории и проектирования, технология, конструкции. М.: Изд-во АСВ, 2004. 560 с.
6. Сталефибробетонные конструкции (СП 52-104-2006). М.: ФГУП НИЦ «Строительство», 2007. 61 с.
7. Grace N.F., Singh S.B. Design approach for carbon fiber-reinforced polymer prestressed concrete bridge beams // ACI Structural Journal. 2003, May-June. P. 365–376.
8. Lorenzis L. de., Miller B., Nanni A. Bond of FRP laminates to concrete // ACI Material Journal. V. 98. 2001. №3. P. 256–264.
9. Reda Taha M.M., Shrive N.G. Enhancing fracture toughness of high-performance carbon fiber cement composites // ACI Materials Journal. 2001. №2. P.168–178.
10. Panchenko L.A. Design of transformable beam systems. Innovations and technologies in construction. BUILDINTECH BIT 2021. Civil Engineering Lecture Notes. 2021. Vol. 151. https://doi.org/10.1007/978-3-030-72910-3_47.
11. Panchenko L.A. Concrete and fiber-reinforced concrete in a fiber-reinforced polymer framework. BUILDINTECH BIT 2020. Civil Engineering Lecture Notes. 2021. Vol. 95. https://doi.org/10.1007/978-3-030-54652-6_20.
12. Панченко Л.А. Строительные конструкции с волокнистыми композитами. Белгород: Изд-во БГТУ им. В.Г. Шухова, 2013, 183 с.
13. Мэттьюз Ф., Ролингс Р. Композитные материалы. Механика и технология / Пер. с англ. С.Л. Баженова. М.: Техносфера, 2004. 408 с.
14. Юрьев А.Г., Панченко Л.А. Исследования деформирования элементов стеклофибробетонных конструкций. Наука и инновации в строительстве SIB-2008. Матер. Междунар. Конгресса. Воронеж. 2008. Том. 2. С. 254–257.
15. Юрьев А.Г. Вариационные принципы строительной механики. Белгород: Изд-во БелГТАСМ, 2002. 90 с.
16. Гениев Г.А., Киссюк В.Н., Тюпин Г.А. Теория пластичности бетона и железобетона. М.: Стройиздат, 1974. 316 с.
17. Панченко Л.А. Железобетонные цилиндрические резервуары с дополнительным армированием стекловолокном // Вестник БГТУ им. В.Г. Шухова. 2009. № 2. С. 68–71.
18. Панченко Л.А., Шевченко А.В. Железобетонный цилиндрический резервуар с дополнительным дисперсным армированием // Вестник БГТУ им. В.Г. Шухова. 2017. № 9. С. 38–42.
19. Матус Е.П. Применение компьютерного моделирования при расчете прочности при растяжении фибробетона // Известия вузов. Строительство. 2005. № 6. С. 103–105.
