Введение. Цилиндрические резервуары на протяжении ряда лет являются объектом глубокого научного исследования [1–7].  При использовании железобетона актуальной проблемой является ограничение трещинообразования. Помимо уменьшения жесткости конструкции, возникает опасность коррозии арматуры. При наличии же сквозных трещин происходит утечка содержимого из резервуара.В тонкостенных конструкциях эффективным вариантом является сочетание стальной арматуры с короткими волокнами или сеткой из стекловолокна. Арматура воспринимает большую часть как растягивающих, так и изгибающих усилий, а массив стекловолокна (или сетка) у поверхности создает дополнительное противостояние изгибным напряжениям и препятствует раскрытию трещин.Основная часть. Рассмотрим проектирование железобетонного резервуара с дополнительным дисперсным армированием стекловолокнами по предельному состоянию второй группы. Круглое днище резервуара опирается на упругий грунтовый массив.Потенциальная энергия системы в этом случае складывается из потенциальной энергии собственно конструкции и потенциальной энергии упругого основания.Функционал Лагранжа в случае винклеровского основания в конечно-элементном представлении имеет вид:   где  – вектор узловой силы;  и  – векторы узловых перемещений конструкции резервуара и основания соответственно;  и  – матрицы жесткости оболочки и днища соответственно; – матрица жесткости основания.В случае круговой цилиндрической оболочки при осесимметричном нагружении возникают четыре внутренних усилия: изгибающий момент  (ему соответствует нормальное напряжение ), поперечная сила (ей соответствует касательное напряжение ), кольцевое усилие и изгибающий момент  в перпендикулярной плоскости (им соответствует нормальное напряжение ). В связи с этим элемент оболочки может быть получен простой комбинацией элементов для плоского напряженного состояния и изгиба пластины.С точки зрения упрощения решения используем для расчета оболочки нулевой кривизны плоские элементы. Учет геометрических особенностей оболочки обеспечивается рассмотрением вписанного многогранника.При назначении расчетной схемы оболочки необходимо, чтобы плоские конечные элементы вписывались в геометрию оболочки. Для развертывающихся на плоскость поверхностей, какой является поверхность цилиндра, можно использовать прямоугольные конечные элементы.Матрица жесткости прямоугольного плоского элемента оболочки получается простым совмещением матриц жесткости прямоугольных элементов плоского напряженного состояния (с множителем , где  – толщина оболочки,  и  – модуль продольной упругости и коэффициент Пуассона соответственно) и пластинки (с множителем , где  и  – размеры конечного элемента) [8].В теории и на практике встречаются процессы и состояния, при описании которых используют величины, имеющие разрывы. Для их описания можно применить единичную функцию Хевисайда:           (1)которую называют импульсивной функцией нулевого порядка.Остановимся на случае, когда модуль продольной упругости конструкции оболочки терпит разрыв (рис.1). Рис. 1. Ступенчато-переменная функция модуля упругости  При использовании функции Хевисайда функция  получает следующий вид:             (2) Когда имеем несколько разрывов, зависимость (2) принимает следующий вид:                                             (3) где  – модуль на i-м участке, то есть на участке , при этом   Статическую симметрию резервуара обеспечивает, например, гидростатическая нагрузка:      ,                      (4)где  – объемный вес жидкости,  – высота резервуара,  – соответствующая ей координата, отсчитываемая от днища.При условии симметрии из шести уравнений равновесия три удовлетворяются тождественно. Одно из оставшихся уравнений показывает, что продольные силы вдоль оси х постоянны. Ввиду их малости по сравнению с другими внутренними усилиями ими обычно пренебрегают.Оставшиеся два уравнения синтезируются в виде:                    (5)или              (6)где  – цилиндрическая жесткость оболочки.Уравнение (6) по форме совпадает с уравнением балки на упругом основании и имеет аналитическое решение [9].Изгибающий момент  вычисляем по формуле [10]:                      (7)Расчет резервуара с симметричной нагрузкой усложняется при оптимизации его конструктивного решения за счет переменного дискретного и дисперсного армирования. Как было сказано, последнее вводится для усиления противостояния трещинообразованию, интенсивность которого нарастает в направлении днища.С точки зрения технологии производства работ целесообразно интенсивность дисперсного армирования наращивать в этом направлении от одного участка оболочки к другому. При этом наблюдается ступенчато-переменное изменение механических характеристик фибробетона. Задача становится нелинейной.Используем шаговый метод в конечно-элементной интерпретации [11]. Шаг наполнения жидкостью составляет n метров по высоте оболочки. На первом шаге принимаем минимальные проценты дискретного и дисперсного армирования, причем последнее осуществляем лишь в нижней трети резервуара. Определяем упомянутые внутренние усилия , , ,  и напряжения:,   и   Находим также главные напряжения  и .Принимая во внимание расчетные сопротивления стержневой или проволочной арматуры сжатию ( и растяжению  [12], а также такого же рода сопротивления фибробетона  и  [13, 14], проверяем условие прочности для арматуры и крайних волокон фибробетона:                                                                                                                                   (8) В случае неудовлетворения неравенств, то есть образования трещин, необходимо увеличивать процент дискретного и (или) дисперсного армирования.На втором шаге в нагрузку включается следующий слой жидкости. Армирование подлежит пересмотру с учетом соотношений левых и правых частей неравенств (8). При значительном превосходстве правых частей оно может быть оставлено неизменным.Число шагов . На последнем шаге производится корректировка армирования в меньшую сторону, если правые части неравенств значительно превосходят левые части.Для резервуаров с идентичными размерами, отличающихся по дискретному и дисперсному армированию, в качестве оптимальной выбирается система с  минимумом потенциальной энергии [15].В качестве примера рассмотрим железобетонный цилиндрический резервуар (рис. 2) из бетона класса В20 с переменным дисперсным армированием, к которому приложена гидростатическая нагрузка  от жидкости с объемным весом  Стенки резервуара толщиной 10 см имеют трехступенчатую структуру из фибробетона с нарастающим к основанию модулем деформации. Коэффициент Пуассона отличается при этом незначительно и принят единым и равным 0,2. Под днищем предполагается скалистое основание.  Рис. 2. Цилиндрический резервуар, подверженный гидростатическому давлению Для расчета шаговым методом использован вычислительный комплекс ЛИРА [11]. Гидростатическая нагрузка разделена на 10 шагов. В табл.1 представлены внутренние усилия на последнем шаге. Использован конечный элемент с размерами 20×70×10 см.Как следует из таблицы, для восприятия растягивающих усилий необходима установка стальной арматуры (класса А400) в продольном направлении в зоне, примыкающей к наружной поверхности (0<x<1м), в зоне, примыкающей к внутренней поверхности (1м<x<2м) и в кольцевом направлении с процентом армирования0,1 %. Кроме того, осуществляется дисперсное армирование стекловолокном d=15 мкм: при   процент армирования 4,5 % (l - длина волокна, d - диаметр волокна); при  , процент армирования 3 %; при , процент армирования3 %. Таблица 1Внутренние усилия в стенке цилиндрического резервуараx, см,кН×см/смQ, кН /смN, кН /см,кН×см/см0-6,470,27801,031001,12-0,032,120,1792000,03-0.0021,650,005300-0,0101,200,00164000,00200,080500000,0406000000  Расчет круглых плит на упругом основании, моделирующих днища цилиндрических резервуаров, подробно рассмотрен в книге [3].Выводы. Представленный автоматизированный алгоритм расчета железобетонного цилиндрического резервуара с дополнительным дисперсным армированием дает возможность оптимизировать проценты дискретного и дисперсного армирования шаговым методом.