Calculations of building structures, the development and improvement of calculation methods, the use of numerical methods in solving these problems are one of the promising areas. The paper presents on the basis of the theory of plastic flow with hardening the equations of state in increments, the equations of virtual work, the geometric equations in increments for small elongations, shifts, and rotation angles. The calculations take into account, depending on the intensity of the current load, the elastoplastic properties of the material, both the structure and the soil environment interacting with it, with a complex process of loading [1-3]. For numerical implementation, the equations arising from the conditions of Genius – applicable to concrete, Coulomb - Mohr – applicable to soil environments are applicable. The information presented is implemented in the form of software systems that make it possible to obtain numerical results of solutions to problems in the calculation of underground structures. An example of calculating a structure based on these methods using a software package is also given. The elastoplastic properties of the material, both the structure and the soil medium interacting with it, are taken into account during the complex process of their loading.