Большое количество исследований посвящено линейному анализу напряженно - деформированного состояния (НДС) оболочек классической формы: цилиндрической, сферической, полусферической и конической. Однако НДС тонких оболочек сложной геометрии исследовано недостаточно. Понятие оболочек сложной геометрии возникает тогда, когда коэффициенты первой и второй квадратичных форм их срединных поверхностей представляют собой довольно сложные функции криволинейных координат. В статье рассматривается материальная нелинейная устойчивость железобетонной синусоидальной велароидальной оболочки с внутренним радиусом r0 =1 м, внешним радиусом R = 20 м и числом волн n = 8. Оболочка нагружалась нагрузкой от собственного веса и снеговой равномерно распределенной нагрузкой интенсивностью 0,252 т/м2. Численные расчеты проводились в программных комплексах LIRA-SAPR 2013 и STARK ES 2015. Конечноэлементная модель оболочки состоит из 6400 элементов и 3280 узлов, общее число узловых неизвестных - 18991. Для моделирования поверхности использовались плоские оболочечные элементы, имеющие шесть степеней свободы в узле. Граничные условия соответствовали шарнирному опиранию по наружному и внутреннему контурам. В результате расчетов были получены значения перемещений и формы потери устойчивости.