Введение. Успехи в области повышения активности минеральных вяжущих веществ, в первую очередь цемента, позволили значительно повысить предел прочности бетонов на сжатие. В то же время прочность на растяжение оставалась во много раз ниже прочности на сжатие. Это несоответствие в некоторой мере сглаживается использованием волокон в качестве арматуры. Наряду с применением металлических элементов большое распространение получили стеклянные и углеродные волокна [1–11].Фибробетон, содержащий те или иные волокна, как и бетон, деформируется по нелинейному физическому закону. Механические характеристики (модуль продольной деформации, коэффициент поперечной деформации, модуль сдвига) являются переменными величинами и зависят от деформации. Их значения определяются на основе стандартных испытаний. Автором проведены такого рода исследования для цилиндрических элементов из стеклофибробетона. Серии отличались длиной волокна (20 мм, 35 мм, 50 мм) и процентом армирования (1,5 %, 3 %, 4,5 %) в различных сочетаниях. Результаты испытаний представлены в табличной форме [12].Для определения прочностных характеристик стеклофибробетона (предела прочности на растяжение, сжатие) также использовалась стандартная методика с применением образов в виде кубиков и призм.Накопленный материал о физико-механических свойствах фибробетона [13, 14] позволяет вести расчет конструкций на его основе с надлежащим приближением к реальности. Этому способствует переход от линейно упругой постановки задачи к учету физической нелинейности. В рамках программных комплексов, основанных на конечно-элементной модели деформируемого тела, приемлем итерационный процесс, обусловленный сходимостью решения.Материалы и методы. Известно, что при линейном физическом законе прямая задача содержит в качестве неизвестных 15 компонентов напряжений, деформаций и перемещений. Для их определения используются 15 уравнений по трем законам деформирования твердого тела. При наличии нелинейных деформаций в прямой задаче появляются дополнительныенеизвестные – интенсивности касательных напряжений τi  и деформаций сдвига γi , а дополнительными уравнениями являются выражение: γi  и зависимость между τi  и γi  [15]. Возможности современной вычислительной техники привели к разработке общих численных алгоритмов решения нелинейных задач. Одним из наиболее эффективных способов построения таких алгоритмов является линеаризация. Рассматриваемые ниже их варианты предложены российскими учеными А.А. Ильюшиным и И.А. Биргером.При решении в перемещениях в дифференциальных уравнениях и граничных условиях слагаемые, возникающие из-за отклонения от закона Гука, переносятся в правые части уравнений и считаются известными. Таким образом, получаем систему уравнений линейной теории деформирования относительно компонентов перемещений, но с дополнительными объемными и поверхностными силами. В первом приближении полагаем эти дополнительные нагрузки равными нулю и решаем физически линейную задачу. Найдя перемещения, вычисляем интенсивность γi , позволяющую определить нелинейные члены в уравнениях Ламе и граничных условиях и решить вторичную физически линейную задачу с дополнительными нагрузками, отличными от нуля и т.д. Такого рода алгоритм получил название метода дополнительных нагрузок.При решении в напряжениях используется метод дополнительных деформаций. Дифференциальные уравнения равновесия и граничные условия остаются без изменения. Уравнения сплошности вследствие наличия нелинейных членов в выражении физического закона будут содержать дополнительные слагаемые, которые можно интерпретировать как дополнительные деформации и определять их последовательными приближениями. Имеет место также метод переменных параметров упругости. Систему уравнений представляют в форме уравнений линейной теории с переменными «параметрами упругости» Eλ и νλ  (λ – отражает процесс деформирования), которые в процессе последовательных приближений вычисляются по деформациям. Природа физической нелинейности фибробетона исходит из механических свойств бетона, проявляющихся в диаграмме «напряжение-деформация», явлении дилатации. Условием создания алгоритма расчета конструкций является наличие экспериментальных данных, касающихся модулей упругости первого и второго рода, коэффициента Пуассона, пределов прочности на растяжение и сжатие при характерных длинах волокон и процентах армирования по массе. Построение расчета проведено на основе метода конечных элементов (МКЭ), используемого в инженерной практике для решения задач различного типа. Матрица жесткости конечного элемента содержит матрицу упругости D как характеристику связи между напряжениями и деформациями [15]:  D=2G1-2ν1-ννν000ν1-νν000νν1-ν000000121-2ν000000121-2ν000000121-2ν,                            (1) где G – модуль сдвига, ν  – коэффициент поперечной деформации (коэффициент Пуассона).Для физической нелинейности характерны переменные величины этих характеристик. Диапазон изменения коэффициента поперечной деформации мал, в связи с чем, его часто принимают неизменным. Фактический модуль сдвига имеет выражение:Gγi=G1-γi2γis,                  (2)где γi  и γis  – интенсивность деформаций сдвига для фактического состояния и ее предельная величина для рассматриваемого вида напряженного состояния, то есть γi=2 Iε2* ,                          (3)где Iε2*  – второй инвариант девиатора деформаций, равный  Iε2*= 16εx-εy2+εy-εz2+εz-εx2+14γxy2+γyz2+γzx2 ,                  (4)   где εx , εy , εz , γxy , γyz , γzx –  элементы тензора деформаций.При трещинообразовании в фибробетоне имеет место увеличение объема тела – дилатация g0γi2 ,  где g0  – модуль дилатации, зависящий от свойств бетона, типа и содержания армирующих волокон. При чистом сдвиге предельная объемная деформация ϑc  зависит от предельной интенсивности деформации сдвига γis  [16]:ϑc=g0γiс2 ,                          (5)откуда вытекает значение модуля:g0= ϑcγiс2 .                        (6)Линейные деформации дополняются величиной g0γi23 . Сохраняется предположение о пропорциональности и соосности девиаторов напряжений и деформаций.Величина γis  зависит от инвариантов тензора Iε1  и девиаторов Iε2* , Iε3*  деформаций, в которых находят отражение прочность фибробетона на сжатие Rfbc  и растяжение Rfbt .  Эти характеристики, как и модуль g0 , зависят от свойств бетона, типа и содержания армирующих волокон. В качестве математической модели γis  можно использовать аналогичную модель для бетона с заменой характеристик Rbc  и Rbt  на Rfbc  и Rfbt  [12].Основная часть. Как следует из изложенных выше положений о физической нелинейности, синтезирующее уравнение метода конечных элементов содержит вместо постоянной матрицы жесткости K ее переменный аналог K  так, что K= K+ ∆K ,                          (7)где ∆K  – корректирующая матрица, равная ∆K=e∆ke ,                          (8)где ∆ke  – корректирующая матрица отдельного конечного элемента, число которых e.Эта предпосылка заставляет заменить решение упомянутого уравнения МКЭ итерационным решением рекуррентного уравнения:Kq(s)=F-∆Kqs-1qs-1,             (9)илиKq(s)= F-∆ Fs-1 ,                  (10)где qs  – s-e приближение q; ∆F – фиктивная дополнительная нагрузка, равная ∆ Fs-1=∆Kqs-1qs-1.           (11)Сходимость итерационного процесса подтверждена решениями типовых задач. Критерий сходимости имеет энергетическую природу:We-WmWm≤ε0 ,                          (12)где We  – потенциальная энергия деформации конечного элемента; Wm  – средняя величина для e величин We ; ε0  – некоторое малое положительное число.Алгоритм решения физически нелинейной задачи включает: начальный модуль (модуль упругости) и модуль дилатации; формирование матрицы жесткости для отдельных элементов; формирование матрицы жесткости для тела и вектора узловых сил; решение уравнений МКЭ; определение модулей по формуле (2); формирование новых матриц жесткости конечных элементов; определение ∆ke ; формирование ∆K ; решение уравнения (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). Установим зависимость между крутящим моментом Т и кривизной кручения. При γ1R<kk≤γ2R  брус имеет внутреннюю зону деформирования, связанную с участком диаграммы 0–1, и наружную зону, связанную с участком 1–2. Граница между этими зонами определяется радиусомr1=τ1Gkk .                         (13)  Рис. 1. Диаграмма сдвига: γ –  деформация сдвига; τ – касательные напряжения; α=  arctg G (G – модуль сдвига), α1=  arctg G1  (G1  – модуль упрочнения)  Выражение крутящего момента имеет видТ =2πkk0r1Gr3dr+r1RG(λ)r3dr.   (14)В данном случаеGλ=1γγ1G+(γ-γ1)G1,        (15)илиGλ=r1rG+1-r1rG1.       (16)В итоге получаем Т =πkk12Gr14+23r1R3-r13G-G1+12R4-r14G1.                           (17) Предположим, что пределу несущей способности бруса соответствует появление на контуре сечения напряжений τ2 .  При этом r1=γ1γ2R.                            (18)При замене переменной r на γkk формула (17) принимает вид Т =1kk3G0γ1γ2dγ+γ1γ2γ1(G-G1)+γG1γdγ.                                       (19) Формула (19) выгодно отличается от формулы (17) своей непосредственной связью с диаграммой сдвига.Кроме рассмотренной физической нелинейности существуют также и другие виды нелинейности. Геометрическая нелинейность имеет место, во-первых, когда перемещения конструкции вызывают значительные изменения в ее геометрии, так, что уравнения равновесия приходится составлять с учетом изменения формы и размеров конструкции (по деформированной схеме), во-вторых, когда зависимость между перемещениями и деформациями нелинейная.Конструктивная нелинейность возникает вследствие конструктивных особенностей системы, вызывающих изменение расчетной схемы в процессе деформирования (например, изменяются условия закрепления). Она имеет место в конструкциях с односторонними связями, встречается также в контактных задачах. При наличии физической и одного из указанных видов нелинейности в комплексе уравнений для решения такого рода задач появляются дополнительные, специфические уравнения, например, новые зависимости между деформациями и перемещениями, уравнения для условий контакта и другие. При этом остаются в силе предложенные в статье методы решения физически нелинейной задачи.Выводы. Учет физической нелинейности при расчете конструкций из фибробетона позволяет в определенной мере приблизить физическую модель к реальным обстоятельствам. Предложенный метод дополнительных нагрузок в конечно-элементной интерпретации оставляет в силе накопленный опыт автоматизированного проектирования [19]. Итерационная процедура расчета содержит циклы с постепенным приближением физических модулей к их истинным значениям в результате сходимости решения задачи. Вычисление модулей ведется на основе зависимостей механики деформируемого твердого тела, что придает постановке задачи фундаментальный характер. Эффект предложенной методики демонстрируется на примере расчета заглубленного цилиндрического резервуара из стеклофибробетона для жидкости.