Problem Of Elasticity Theory Research Articles

First Basic Problem of Elasticity Theory for a Composite Layer with Two Thick-Walled Tubes

The spatial problem of elasticity theory for a fibrous composite in the form of a layer with two thick-walled cylindrical tubes is solved. Stresses are given on the flat surfaces of the layer and on the inner surface of the tubes. The solution to the problem is presented in the form of Lamé equations in different coordinate systems, where the layer is considered in a Cartesian system and the tubes – in local cylindrical ones. To combine the basic solutions in different coordinate systems, the generalized Fourier method is used. Satisfying the boundary conditions and conjugation conditions between the layer and the tubes, an infinite system of integro-algebraic equations is formed, which is reduced to linear algebraic equations of the second kind, and the reduction method is applied. After finding the unknowns, it is possible to obtain the stress-strain state at any point of the elastic combined bodies using the generalized Fourier method to the basic solutions of the problem. According to the results of numerical studies, it can be stated that the problem can be solved with a given accuracy, which depends on the order of the system of equations and has a rapid convergence of solutions to the exact one. Numerical analysis of the stressed state was considered with a variation of the distance between the tubes. The graphs of the distribution of internal stresses in the tubes and the layer are obtained. The results show an inverse relationship between the magnitude of stresses and the distance between the tubes. In addition to the absolute value of stresses, changes in the character of the diagrams and the sign are possible. The proposed method of solution can be applied in the design of a layer with tubes. The obtained stress-strain state makes it possible to preliminarily evaluate the geometric parameters of the structure. Further development of the research topic is necessary for a model where tubes are combined with other types of inhomogeneities.

Tensor linearity of two-dimensional isotropic functions in the plane problem of nonlinear theory of elasticity

It is shown that a nonlinear isotropic tensor function of the second rank in two-dimensional space, which is a power series in its tensor argument, is representable by a finite binomial tensor linear relation. Expressions of two coefficients of this relation are given in terms of an infinite set of coefficients of the original series and two independent invariants of the tensor argument. In relation to continuum mechanics, the reducibility of the constitutive relations in the plane problem of tensor nonlinear elasticity theory to the tensor linear connection of the corresponding second-order minors of stresses and strains is established.

Calculation of resistance force at firm clay digging with a cylindrical bucket

BACKGROUND: The performance of single-bucket hydraulic excavators affects many areas of activity. This parameter of the machine largely depends on the volume of the material being moved. However, the bucket resistance forces do not allow large buckets to be installed on an excavator. AIM: The paper discusses the design of the bucket which resistance to digging is reduced when it embeds into the ground. This statement has to be verified. METHODS: For this aim, the process of embedding the bucket into firm dry clay is considered. This kind of soil is the most difficult to be processed. Since there is no viscosity in this soil, the equations of spatial problems of elasticity theory are used to describe the process. A number of assumptions are made, and the system of differential equations describing stresses in the soil is solved. RESULTS: The result of the solution is the obtained dependence for determining the normal pressure from the soil in the process of its destruction, which helps to determine the overall digging resistance force. The initial parameters of the soil and bucket of the proposed design, developed for the excavator based on the YuMZ tractor, are described. The initial parameters are substituted into the resulting solution. CONCLUSION: The obtained value of the digging resistance force is significantly lower than the force that the hydraulic drive of an ordinary excavator has to overcome.

Analysis of the three‐dimensional elasticity theory problem for a radially inhomogeneous cylinder

AbstractIn this paper an axisymmetric problem of elasticity theory for a radially inhomogeneous cylinder of small thickness is studied. It is assumed that homogeneous mixed boundary conditions are specified on lateral surfaces of the cylinder, and boundary conditions that leave the cylinder in equilibrium are specified on the ends of the cylinder. It is considered that elastic moduli are arbitrary continuous functions of a variable along the radius of the cylinder. The formulated boundary value problem is reduced to a spectral problem containing a small parameter characterizing the thinness of the cylinder walls. To determine its solution, the method of asymptotic integration is used, based on three iteration processes. All solutions of the equilibrium equation are constructed, satisfying the specified homogeneous boundary conditions on lateral surfaces. It is obtained that the solution of the problem consists of a penetrating solution and a solution of the nature of the boundary layer. An analysis of stress‐strain states corresponding to different types of solutions is carried out. Based on the constructed asymptotic solutions, a connection is found between the solutions obtained by the equation of elasticity theory and by the applied theory of shells. It is shown that the penetrating solution determines the internal stress‐strain state of a radially inhomogeneous cylinder. The stress state determined by the penetrating solution is equivalent to the principal vector of stress acting in a cross section perpendicular to the axis of the cylinder.The first terms of asymptotic expansions of the penetrating solution in a small parameter coincide with solutions in the applied theory of shells. New classes of solutions are defined as the character of a boundary layer, which are localized at the ends of the cylinder and decrease exponentially with distance from the ends. These solutions are absent in the applied theories of shells. The first terms of the asymptotic expansion of the solution, having the nature of a boundary layer, are equivalent to the Saint‐Venant edge effect in the theory of heterogeneous plates. Asymptotic formulas for displacement and stresses are constructed, allowing one to calculate the three‐dimensional stress‐strain state of a radially heterogeneous cylinder of small thickness with any predetermined accuracy. Based on the obtained asymptotic expansions, one can estimate the applicability areas of applied theories and construct a refined applied theory for radially heterogeneous cylindrical shells.

CALCULATION OF FLAT SHELLS BY THE BOUNDARY ELEMENT METHOD

This article describes the iterative process of solving the problem of the stress-strain state of a long flexible cylindrical panel. Stress-strain state and the study of the stability of flat shells taking into account the physical, mechanical properties of the material, the process of deformation change. It is based on sequential approximation methods and the boundary element method (MGE). The results of algorithmic weighting are presented in the form of a table showing what value each band represents. The equations are determined by solving twelve unknown values of functions at the ends of the segment. Linear problems of elasticity theory and plate theory fundamental solutions have a simple form, so the method is widely used here. For flat shells, the matrix of fundamental solutions is determined by complex volumetric expressions, and for flat shells - by special functions. Therefore, there is little research on solving problems in the theory of flat shells by the boundary element method. In this regard, an urgent topic of research is the development of methods of boundary integral equations for solving linear and nonlinear problems in the theory of flat shells based on the application of fundamental solutions defined by simple analytical expressions. The scientific novelty of the work consists in the development of a methodology for assessing the reliability of thin-walled spatial structures using the boundary element method.

An Analysis of the Stress–Strain State of a Layer on Two Cylindrical Bearings

A spatial problem of elasticity theory is solved for a layer located on two bearings embedded in it. The bearings are represented as thick-walled pipes embedded in the layer parallel to its boundaries. The pipes are rigidly connected to the layer, and contact-type conditions (normal displacements and tangential stresses) are specified on the insides of the pipes. Stresses are set on the flat surfaces of the layer. The objective of this study is to obtain the stress–strain state of the body of the layer under different geometric characteristics of the model. The solution to the problem is presented in the form of the Lamé equation, whose terms are written in different coordinate systems. The generalized Fourier method is used to transfer the basic solutions between coordinate systems. By satisfying the boundary and conjugation conditions, the problem is reduced to a system of infinite linear algebraic equations of the second kind, to which the reduction method is applied. After finding the unknowns, using the generalized Fourier method, it is possible to find the stress–strain state at any point of the body. The numerical study of the stress state showed high convergence of the approximate solutions to the exact one. The stress–strain state of the composite body was analyzed for different geometric parameters and different pipe materials. The results obtained can be used for the preliminary determination of the geometric parameters of the model and the materials of the joints. The proposed solution method can be used not only to calculate the stress state of bearing joints, but also of bushings (under specified conditions of rigid contact without friction on the internal surfaces).

The asymptotic solution of the elasticity theory problem for a radially inhomogeneous cylinder

The elasticity theory problem for a radially inhomogeneous cylinder of small thickness, whose elastic moduli are arbitrary continuous functions depending on the radius of the cylinder, is considered. It is assumed that the side surface of the cylinder is stress-free, and the boundary conditions that keep the cylinder in equilibrium are given at its seats. Asymptotic solutions are constructed by the asymptotic integration method. It is shown that the asymptotic solution consists of the sum of the penetrating solution, simple boundary effect and boundary layer solutions. The character of the stress-strain state corresponding to the penetrating solution, simple boundary effect and boundary layer solutions is determined. Asymptotic formulas are obtained for displacements and stresses, which allow to calculate the stress-strain state of a cylinder.The problem of torsion of a radially inhomogeneous cylinder is studied, with its lateral surface free from stress and the boundary conditions keeping it in equilibrium at its seats. By applying the asymptotic integration method, it is determined that the asymptotic solution for the torsion problem consists of the sum of the penetrating solution and boundary layer solutions.Numerical analysis is performed and the effect of material inhomogeneity on the stress-strain state of the cylinder is evaluated.

Конечно-элементное моделирование нелинейных задач упругости в абсолютных узловых координатах на неструктурированных шестигранных сетках

We consider a finite element approach to solving the problem of elasticity theory in terms of absolute nodal coordinate formulation (ANCF), in which large body displacements are described in a global reference frame without using any local coordinate system. The main feature of the method is the absence of gyroscopic effects and, as a result, constancy of the mass matrix and the vector of the generalized force of gravity. In contrast to the traditional ANCF approach, the sets of nodal degrees of freedom of the finite element are formed only on the basis of absolute coordinates of nodes, which allows us to solve the problem using, in particular, unstructured hexahedral meshes. To construct the stiffness matrix, we apply a second-order automatic differentiation algorithm, which ensures its symmetrical form (Hessian matrix) and is analytically accurate in calculating the derivative. This approach also makes it possible to carry out calculations for models of hyperelastic materials without the corresponding Piola-Kirchhoff tensor. It has been shown that within the framework of discretization of the equation of motion, along with the well-known Newmark numerical integration scheme, it is possible to use the HTT-α scheme, which is unconditionally stable, second-order accurate and dissipative for high frequencies. We present several examples of solving static and dynamic elasticity problems for compressible and incompressible models of hyperelastic materials, in which the functions of internal energy density of the body are specified in terms of the deformation gradient.

A Volterra edge dislocation problem in second-order elasticity theory

A displacement function technique has been developed to study plane strain problems in second-order elasticity theory for incompressible elastic materials. The formulation culminates in the development of a single inhomogeneous biharmonic equation of the form [Formula: see text] for the second-order displacement function [Formula: see text]. In the second-order, the isotropic stress associated with constitutive behaviour is governed by an inhomogeneous harmonic equation of the form [Formula: see text]. Both functions [Formula: see text] and [Formula: see text] depend only on the solution to the analogous problem in classical elasticity. While there are similarities to the Airy stress function, the displacement function approach enables the evaluation of the first- and second-order displacement fields purely through its derivatives, thus eliminating the introduction of arbitrary rigid body terms normally associated with formulations where the strains derived from a stress function approach need to be integrated. The displacement function approach is applied to develop a second-order elasticity solution to the problem of a Volterra edge dislocation.

The method of developing the complex stress tensor by basic states for the construction solutions of spatial boundary problems of the linear elasticity theory

The method of developing the complex stress tensor by basic states for the construction solutions of spatial boundary problems of the linear elasticity theory

Triaxial Load Cell for Ergonomic Risk Assessment: A Study Case of Applied Force of Thumb

To assess the ergonomic risk level in work systems involving tasks performed with hands or fingers, it is necessary to know the exerted triaxial forces. To address this need, a prototype of a triaxial load cell based on principles of linear elasticity theory and mechanical problems of torsion, bending and axial load is presented. This work includes an analytical strain model for each instrumented point and its solution regarding the applied force to a triaxial load cell. The proposed load cell was calibrated and validated by performing different static experimental tests. As a case study, the applied force in three directions while the thumb activates a cigarette lighter was measured. Triaxial forces and resultant forces were obtained and compared with the parameter of 10 N established by the ergonomic standards as reference values for pressing down with the thumb, finding that the applied forces in eight tests were 23.73 N, 43.51 N, 12.69 N, 14.50 N 20.35 N, 21.67 N, 39.74 N and 46.02 N, exceeding the reference values and establishing a direct relationship with Quervain syndrome. In conclusion, the developed load cell is a valid and reliable alternative to measure many forces that cannot be obtained with commercial devices, allowing the level of ergonomic risk to be determined with great precision.

Метод чисельного визначення динамічних характеристик лопаті повітряного гвинта з композиційного матеріалу

The subject matter of this article is the dynamic characteristics of a composite propeller blade. Determination of the modes and frequencies of natural vibrations is necessary to predict the dangerous resonance modes of aircraft engines and to identify the most stressed local zones of the blade surface. The goal of this study is to develop a verified method for determining the dynamic characteristics of composite rotor parts of aircraft engines based on the known properties of the structural components of the composite material. A general mathematical statement of the problem of elasticity theory for the analysis of the dynamic characteristics of composite structures is described. A complete system of equations that describes the mechanical state of the body within the framework of the continuum mechanics approach was developed. Geometric modeling of the propeller blades was performed. Modeling and numerical investigation of the natural frequencies and mode shapes of the propeller blade were performed using the finite element method using the ANSYS software package. Based on the results of numerical research of the stressed and deformed state of the propeller blade, the first five natural vibration frequencies and the distribution of the local stress field were determined. The most stressed local zones on the blade surface were determined for each form of natural vibration. The propeller blade model was verified using an experimental study of the first five forms and frequencies of natural blade vibrations. The eigenfrequencies of the blade vibration were experimentally determined using the method of free (natural) vibrations. The resonance method was used to experimentally determine the resonant frequencies and vibration modes of the blade. The distribution of the blade deformation field was investigated using the strain gauge method. The highest error in the verification of numerical and experimental research is 4.11% for the fourth vibration frequency. Conclusions. The scientific novelty of the results obtained is that the effective elastic properties of the composite material for calculations should be determined using the procedure of numerical homogenisation of composite materials of different reinforcement structures by the properties of the matrix and fibers. The method does not require experimental determination of the effective elastic constants for the layers of blade components of different weave architecture patterns.

Multi-criteria parametric optimization of the displacement and weight of a shell of minimal surface on a circular contour consisting of two inclined ellipses under thermal and power loading with consideration of geometric nonlinearity

Thin shells are well suited to optimal design problems, and the finite element method and gradient descent method make it possible to solve inverse problems in structural and applied mechanics. The calculation process takes into account geometric nonlinearity. In modern numerical studies, a formulation taking into account geometric nonlinearity in the finite element method is used, namely, a stepwise loading procedure. It becomes necessary to take into account the relations between displacement vectors and their derivatives and strain increments. These relations help to determine the stiffness matrix of the finite element at each loading step, which leads to a qualitative study of real displacements in the structure. The geometrically nonlinear calculation of a shell of minimum surface on a circular contour consisting of two inclined ellipses in a multi-criteria parametric optimization is performed by the finite element method. The finite element method is a universal variational method that is focused on solving the most complex problems of elasticity theory and applied and structural mechanics using separate calculation complexes. Stiffness matrix - for the whole body is formed on the basis of the finite element stiffness matrix. The system of solving equations of the finite element method is formed using the Lagrange's variational principle, according to which the total potential energy P of a finite element model of a body in a state of stability and equilibrium has a minimum value. As part of the sensitivity analysis, gradients of the design variables of the structure, displacements in the form of partial derivatives of these characteristics along the design variables, and shell thickness are calculated. The sensitivity information serves as the basis for building an optimal design algorithm using the gradient descent method of the objective function. Standard multi-criteria parametric optimization allows for an average of 10% steel savings, with geometric nonlinearity increasing to 20%. At this study site, 24.4% of sheet steel was saved, which is a significant relative saving. This methodology of the authors shows high results for innovative design and production of steel thin shells in Ukraine and around the world, and also makes it possible to apply several types of optimization to one research object.

A mathematical model of the gravitational potential of the planet taking into account tidal deformations

Objectives. This paper investigates the gravitational potential of a viscoelastic planet moving in the gravitational field of a massive attracting center (star), a satellite and one or more other planets moving in Keplerian elliptical orbits relative to the attracting center. Celestial bodies other than a viscoelastic planet are modeled by material points. Within the framework of the linear model of the theory of viscoelasticity, the problem of finding the vector of elastic displacement has been resolved. Traditionally, a solid body model is used to determine the Earth’s gravitational field, while tidal deformations are taken into account in the form of small corrections to the coefficients of the geopotential model. In this work, the viscoelastic ball model is used to take into account tidal effects. The relevance of the research topic is associated with high-precision forecasting of the movement of artificial satellites of the Earth, high-precision measurement of the Earth’s gravitational field.Methods. In this study the asymptotic and analytical methods developed by V.G. Vilke are used for mechanical systems containing viscoelastic elements of high rigidity, as well as methods of classical mechanics, mathematical analysis. The graphs were plotted using the Octave mathematical package.Results. After resolving the quasi-static problem of elasticity theory by calculating triple integrals over a spherical area, a formula for the gravitational potential of a deformable planet was obtained. In addition, the gravitational potential of the Earth was also calculated taking into account solid-state tidal effects from the Moon, Sun, and Venus at an external point. Graphs were constructed to show the dependence of the Earth’s gravitational potential on time.Conclusions. The theoretical and numerical results established herein show that the main contribution to the gravitational potential of the Earth is made by the Moon and the Sun. The influence of other planets in the solar system is small. The value of the gravitational potential at the outer point of the Earth, taking into account tidal effects, depends both on the position of the point in the moving coordinate system and on the relative position of celestial bodies.

Обобщённая математическая модель основных задач теории упругости для анизотропных тел

Introduction. The assessment of the strength of underground structures is a complex task that requires a comprehensive approach. In many cases empirical methods based on experience are not sufficient. For effective application of analytical methods, it is necessary to use mathematical modelling. This allows to simplify the problem considerably, to consider different variants, to identify cases to obtain a solution in a closed form. Materials and methods. The model is based on a system of singular integral equations with a shift and complex conjugate terms. The solution of the system is the complex elastic potential. Systems of singular equations with a Carleman shift are considered on the space of random functions converging in the mean square to expand the possibility of practical use of the work results. Results. The results of the study are: a mathematical model formulation, the development of an algorithm for reducing a system of singular equations to a system of known Fredholm integral equations of the second kind, a method for determining the solvability conditions and the number of solutions. Discussion. As a result of the study: (1) a general algorithm for finding a solution to a mathematical model is obtained; (2) the methodology for determining the number of solutions and the type of solvability conditions is given; (3) The solution algorithm is extended for the stochastic case; (4) The possibility of wider application of the mathematical apparatus of functions of complex variable for research of strength properties of rocks is proved. Conclusion. The proposed comprehensive approach to assessing the strength characteristics of underground structures, combining a set of analytical methods in conjunction with mathematical modelling, expands the variability of calculation and allows to simplify the problem significantly and obtain a solution in closed form. Resume. The paper develops a mathematical model of the stress state of an elastic body, which allows solving the basic problems of elasticity theory for a variety of media. Suggestions for practical applications and directions for future research. The results of the study make it possible to improve the efficiency of analytical methods in the study of the strength of underground structures. The obtained results are applicable for predicting the strength and safety of underground structures.

An Analytical Mechanics Model for the Rotary Sliding Triboelectric Nanogenerator.

In recent years, global attention towards new energy has surged due to increasing energy demand and environmental concerns. Researchers have intensified their focus on new energy, leading to advancements in technologies like triboelectrification, which harnesses energy from the environment. The invention of the triboelectric nanogenerator (TENG) has led to new possibilities, with the rotary sliding TENG standing out for its superior performance. However, understanding its mechanical behavior remains a challenge, potentially leading to structural issues. This paper introduces a novel analytical mechanics model to analyze the mechanical performance of the stator of the rotary sliding TENG, offering a new analytical solution. The solution also presents an innovative approach to solving axisymmetric problems in elasticity theory since it challenges a traditional assumption that the stress function depends solely on the radial coordinate, proposing a new stress function to derive a more general solution, supplementing the classical approach in the theory of elasticity. Through the obtained solutions, the mechanical characteristics of the rotary sliding TENG during operation are analyzed. A clearer relationship between mechanical characteristics and electrical output is expected to provide a theoretical basis for the design of the rotary sliding TENG.

Variation-Difference Method of Calculation of Layered Rubber-Metal Vibration Isolators Used for Protection of Reinforced Concrete Buildings from Anthropogenic Vibration

In the modern construction complex of the city of Moscow for the protection of buildings and structures from man-made vibration arising from the movement of trains of rail transportation (subway trains, railroad lines and trams). To protect buildings and structures from anthropogenic vibration arising from the movement of rail transport trains (subway trains, railroad lines and streetcars), Moscow uses layered rubber-metal vibration isolators [1]. Most often, to determine their static and dynamic characteristics, the finite element method (FEM) is used, which makes it possible to determine all components of the stress-strain state and frequencies of free oscillations in the loaded state practically for any structural forms of isolators. However, for the most popular software packages that implement FEM, the problem of optimizing the structural shape of the vibration isolator still requires significant time expenditures for multiple changes of the finite element mesh, repeated setting of boundary conditions and implementation of a series of calculations. Only some of the software packages implementing the FEM solve optimization problems of the shape of the product being calculated, most often it is related to foreign software products with universal functionality. Variation-difference method (VDM) is the closest to the finite element method (FEM) in terms of its computational capabilities. It is possible to create program modules that repeatedly and automatically solve three-dimensional problems of elasticity theory taking into account the changed geometry of the vibration isolator: the dimensions of the product, the location of perforations within the rubber layers, as well as the thickness of the rubber layer and other parameters important for obtaining an effective technical solution for vibration isolation of buildings. Further, the article describes the method of implementation of the variational-difference method (VDM) as applied to the solution of the problem of determining the components of the stress-strain state inside a three-dimensional layered vibration isolator with perforations of different sizes having different locations relative to the contour of the vibration isolator, i.e., the solution of the problem of optimizing the three-dimensional shape of the vibration isolator is given.

Procedure of Continuation of Boundary Conditions in the Problems of Elasticity Theory

Procedure of Continuation of Boundary Conditions in the Problems of Elasticity Theory

The Problem of the Reliability of Bending Models for Composite Plates of Medium Thickness

Most refined bending models of medium-thick plates, which consider transverse shear and partial compression deformations, differ little. However, despite a significant increase in the order of the governing differential equations, the results obtained from their equations give mainly a small increase in accuracy compared to the existing theories. On the other hand, such an increase in the order of the constructed systems of differential equations requires a significant increase in the effort required to solve them, complicates their physical interpretation, and narrows the range of people who can use them, primarily engineers and designers. Therefore, developing a plate-bending model that incorporates all the above factors and is on par with previously applied theories regarding the complexity of the calculation equations remains relevant. For example, most of the applied theories that do not consider transverse compression cannot be used to solve problems of contact interaction with rigid and elastic dies and bases because it is impossible to satisfy the conditions at the contact boundary of the outer surface of the plate, as well as the boundary conditions at the edges of the plate. Therefore, to provide guaranteed accuracy of the results, some researchers of these problems have introduced such a concept as “energy consistency” between the functions of representation of the displacement vector components, their number, the order of equations, and the number of boundary conditions. The authors, based on the developed version of the model of orthotropic plates of medium thickness, investigate the problem of taking into account the so-called “energy consistency” effect of the bending model, depending on the order of the design equations and the number of boundary conditions, as well as its usefulness and disadvantages in practical calculations. The equations of equilibrium in displacements and expressions for stresses in terms of force and moment forces are recorded. For rectangular and circular plates of medium thickness, test problems are solved, and the numerical data are compared with those obtained using spatial problems of elasticity theory, as well as the refined Timoshenko and Reissner theories. An analysis of the obtained results is provided.

Application of the Double Approximation Method for Constructing Stiffness Matrices of Volumetric Finite Elements

Introduction. When numerically solving problems of elasticity theory in a three-dimensional formulation by the finite element method, finite elements (FE) in the form of parallelepipeds, prisms and tetrahedra are used. Regularly, the construction of stiffness matrices of volumetric FE is based on the principle of isoparametricity, which involves the Lagrange polynomials to approximate the geometry and displacements. In computational practice, the most widespread FE are the so-called multilinear isoparametric FE with a linear law of approximation of displacements. The main disadvantage of these elements lies in the “locking” effect when modulating bending deformations. Moreover, the error of the numerical solution increases drastically in the case when the structure, in comparison to conventional deformations, undergoes significant displacements as a rigid whole. Long-term experience in solving problems of deformable solid mechanics by the finite element method has shown that existing volumetric FE have slow convergence, specifically, when modeling bending deformations of plates and shells. This study aims at constructing stiffness matrices of multilinear volumetric FE of increased accuracy allowing for rigid displacements based on the double approximation method.Materials and Methods. The mathematical apparatus of the double approximation method based on the principle of a separate representation of the distribution functions of displacements and deformations inside the element, was used to construct the stiffness matrices of volumetric FE. The storage and processing of the resulting system of equations was implemented in algorithmic terms of sparse matrices. Software development and computational experiments were carried out using the Microsoft Visual Studio 2013 64-bit computing platform and the Intel ® Parallel Studio XE 2019 compiler with the integrated Intel ® Visual Fortran Composer XE 2019 text editor. Visualization of the calculation results was performed using the descriptor graphics of the MATLAB computer mathematics package. A large eight-node SOLID185 CE of the ANSYS Mechanical software complex was used as a test sample.Results. Mathematical tool and software were developed to study the stress-strain state of massive structures under various types of external actions. The authorized application software package was verified on test examples with known analytical solutions. It has been shown that the constructed FE accurately satisfy the basic requirements for finite element modeling of spatial problems of elasticity theory.Discussion and Conclusion. The performed testing of the developed mathematical and program toolkit has shown that the finite elements constructed on the basis of the double approximation method can successfully compete with similar SOLID185 volumetric elements of the ANSYS Mechanical software complex. The proposed elements can be integrated into domestic import-substituting software systems that implement the finite element method in the form of the displacement method.