A new approach to numerical simulation of boundary value problems of the theory of elasticity in stresses and strains

The main parameters characterizing the process of deformation of solids are displacements, strain and stress tensors. From the point of view of the strength and reliability of the structure and its elements, researchers and engineers are mainly interested in the distribution of stresses in the objects under study. Unfortunately, all boundary value problems are formulated and solved in solid mechanics mainly with respect to displacements, or an additional stress functions. And the required stresses are calculated from known displacements or stress functions. In this case, the accuracy of stress calculation is strongly affected by the error of numerical differentiation, as well as the approximation order of the boundary conditions. The formulation of boundary value problems directly with respect to stresses or strains allows to increase the accuracy of stress calculation by bypassing the process of numerical differentiation. Therefore, the present work is devoted to the formulation and numerical solution of boundary value problems directly with respect to stresses and strains. Using the well-known Beltrami-Miеchell equation, and considering the equilibrium equation as ah additional boundary condition, a boundary value problem(BVP) is formulated directly with respect to stresses. In a similar way, using the strain compatibility condition, the Beltrami-Mitchell type equations for strains are written. The finite difference equations for two-dimensional BVP are constructed and written in convenient a form for the use of iterative method. A number of problems on the equilibrium of a rectangular plate under the action of various loads applied on opposite sides are numerically solved. The reliability of the results is ensured by comparing the numerical results of the 2D elasticity problems in stresses and strains, and with the exact solution, as well as with the known solutions of the plate tension problem under parabolic and uniformly distributed load