В настоящее время на территории Украины и Молдовы располагается большое количество хранилищ с легковоспламеняющимися либо ядовитыми жидкостями. Если такие хранилища подвергнутся действию наземного взрыва, вызванного террористическим актом, катастрофой авиалайнера, либо сейсмической волной, это приведет к экологической катастрофе. Поэтому проблема оценки состояния таких резервуаров и моделирование аварийных ситуаций является весьма актуальной на сегодняшний день. Полученные в результате численного моделирования данные позволят оценить предел прочности конструкции при ударном либо сейсмическом воздействии. Так как описанные выше хранилища представляют собой тонкостенные конструкции, заполненные жидкостью (частично или полностью) и погруженные в грунт, то становится актуальным построение моделей, учитывающих влияние воздействия жидкости или газа на напряженно-деформированное состояние, частоты и формы их собственных колебаний. Для теоретического исследования динамики взрывного нагружения используются сложные двумерные модели упругопластической среды [4,5]. Будем рассматривать движения сплошной среды в двумерной постановке в лагранжевой системе координат. Математическая модель представляет собой систему нелинейных дифференциальных уравнений в частных производных, для численного решения которой был предложена модифицированная конечно-разностная схема Уилкинса [2]. В данной модели существует возможность проведения численного моделирования детонации заряда взрывчатого вещества, расположенного внутри или вне цилиндрической оболочки. Заряд ВВ в расчетах занимает часть объема (оставшаяся часть заполняется жидкостью) или находится за пределами резервуара. В данном случае для моделирования материала пластины (алюминия, стали) использовались уравнения состояния в форме Ми-Грюнайзена [2]. В качестве теста, для жидкости, было выбрано уравнение состояния воды. Уравнение состояния для ВВ было выбрано в форме закона Тета при упругопластическом нагружении. В нескольких ячейках задавалось условие инициализации детонации. После достижения в ячейках, занятых ВВ некоторого критического значения, осуществлялось переключение уравнения состояния на уравнение состояния продуктов детонации. Таким образом, моделировалось распространение детонационной волны по ВВ. Основные результаты: моделируется немгновенная детонация взрывчатого вещества внутри закрытого с торцов контейнера, заполненного жидкостью и находящегося внутри различных сред (соли, кварца, грунта). Для соли, кварца и воды использовалось уравнение Ми-Грюнайзена (для неметаллов). В качестве модели грунта была выбрана модель Ляхова [6]. Из соображений симметрии производится расчет правой верхней четверти контейнера.