COMPUTER MODELING OF THE STRESS-STRAIN STATE OF AN EXPLOITED STORAGE TANK FOR PETROLEUM PRODUCTS
Authors: Конопацкий Е.В. , Шевчук О.А. , Крысько А.А. 
Published: June 2022

Целью исследования является развитие методов многомерной интерполяции и аппроксимации как инструментов математического и компьютерного моделирования напряженно-деформированного состояния тонкостенных оболочек инженерных сооружений. Для достижения поставленной цели усовершенствовано дифференциальное уравнение моделирования напряженно-деформированного состояния упругой цилиндрической оболочки при осесимметричном нагружении для численного анализа напряженно-деформированного состояния цилиндрического резервуара с несовершенствами геометрической формы и получено его численное решение с помощью геометрических интерполянтов. Также реализован новый подход к учёту начальных условий дифференциального уравнения, который заключается в параллельном переносе численного решения в нужную точку, координаты которой соответствуют начальным условиям. Предложенный подход численного решения дифференциальных уравнений основан на геометрической теории многомерной интерполяции, реализованной в точечном исчислении. Для сравнения результатов моделирования напряженно-деформированного состояния резервуара для хранения нефтепродуктов с учётом несовершенств геометрической формы было получено эталонное решение путём аппроксимации значений перемещений от действия гидростатической нагрузки с учётом геометрической и конструктивной нелинейности при моделировании напряженно-деформированного состояния в программном пакете конечно-элементного анализа SCAD. Расчёты были проведены в соответствии с прочностной теорией октаэдрических касательных напряжений (энергетическая теория Губера-Хенки-Мизера). Предложенный подход использования геометрических интерполянтов в качестве конечных суперэлементов, включающих в себя не только геометрическую информацию, но и информацию о физических параметрах, может быть использован для численного решения и других дифференциальных уравнений математического моделирования многофакторных процессов и явлений. Преимуществом такого подхода является то, что он позволяет исключить необходимость согласования геометрической информации в процессе взаимодействия между CAD и FEA системами по аналогии с изогеометрическим методом.

The aim of the study is to develop methods of multidimensional interpolation and approximation as tools for mathematical and computer modeling of the stress-strain state of thin-walled shells of engineering structures. To achieve this goal, we improved the differential equation for modeling the stress-strain state of an elastic cylindrical shell under axisymmetric loading for the numerical analysis of the stress-strain state of a cylindrical tank with geometric form imperfections and obtained its numerical solution using geometric interpolators. Also realized a new approach to taking into account the initial conditions of the differential equation, which is a parallel transfer of the numerical solution to the desired point, the coordinates of which correspond to the initial conditions. The proposed approach for the numerical solution of the differential equations is based on the geometrical theory of multidimensional interpolation realized in the point calculus. To compare the results of modeling the stress-strain state of a storage tank for petroleum products with regard for geometric form imperfections, a reference solution was obtained by approximating the values of displacements from the hydrostatic load with regard for geometric and structural nonlinearity in modeling the stress-strain state in the software package of finite element analysis SCAD. Calculations were performed in accordance with the strength theory of octahedral tangential stresses (Huber-Hancke-Miser energy theory). The proposed approach of using geometrical interpolants as finite superelements that include not only geometrical information but also information on physical parameters can be used for numerical solutions and other differential equations of mathematical modeling of multifactor processes and phenomena. The advantage of this approach is that it eliminates the need to coordinate geometric information in the process of interaction between CAD and FEA systems by analogy with the isogeometric method.