The work considers rods with a constant cross-section. The deformation law of each layer of the rod is adopted as an approximation by a polynomial of the second order. The method of determining the coefficients of the indicated polynomial and the limit deformations under compression and tension of the material of each layer is described with the presence of three traditional characteristics: modulus of elasticity, limit stresses at compression and tension. On the basis of deformation diagrams of the concrete grades B10, B30, B50 under tension and compression, these coefficients are determined by the method of least squares. The deformation diagrams of these concrete grades are compared on the basis of the approximations obtained by the limit values and the method of least squares, and it is found that these diagrams approximate quite well the real deformation diagrams at deformations close to the limit. The main problem in this work is to determine if the rod is able withstand the applied loads, before intensive cracking processes in concrete. So as a criterion of the conditional limit state this work adopts the maximum permissible deformation value under tension or compression corresponding to the points of transition to a falling branch on the deformation diagram level in one or more layers of the rod. The Kirchhoff-Lyav classical kinematic hypotheses are assumed to be valid for the rod deformation. The cases of statically determinable and statically indeterminable problems of bend of the rod are considered. It is shown that in the case of statically determinable loadings, the general solution of the problem comes to solving a system of three nonlinear algebraic equations which roots can be obtained with the necessary accuracy using the well-developed methods of computational mathematics. The general solution of the problem for statically indeterminable problems is reduced to obtaining a solution to a system of three nonlinear differential equations for three functions - deformation and curvatures. The Bubnov-Galerkin method is used to approximate the solution of this equation on the segment along the length of the rod, and specific examples of its application to the Maple system of symbolic calculations are considered.
