Creep of Concrete in Shell Structures: Nonlinear Theory

The creep of concrete is one of the main problems threatening concrete structural development and the stability and safety of structures. However, the nonlinear theory is the key to solving the problem of taking into account the physical and mechanical properties of concrete creep in shell structures. To create such a theory, the original shell is replaced by a continuous equivalent elastic shell. To determine the stress–strain state of the structure, the equations of nonlinear creep and crack growth are derived, and a deformation model of the section is created. The behavior of the structure at all stages of the life cycle is investigated by solving the solving systems of differential equations of equilibrium, motion, and perturbation of the equivalent shell. The values of the ratios of dependence of long-term and short-term critical loads on deformations, forces, cracks, etc., are given. The accuracy of the solution of the developed nonlinear theory is compared with the linear theory of concrete creep as well as experimental data. The results show that, according to the linear theory, for the values for the short term and long term, up to 56% and up to 39% of critical loads are overestimated, respectively. The creep process in practical engineering can be effectively controlled by the results of the proposed theory.


Introduction
The phenomenon of metal creep at high temperatures is one of the reasons for the destruction of the Twin Towers in New York (USA) [1]. Due to the creep of materials, the structures of Basmanny Market in Moscow (Russia) and a bridge in Genoa (Italy) collapsed [2].
The development of creep increases the stress loss in the concrete structure and redistributes the internal force of statically indeterminate structures. In turn, excessive deformation of the structure significantly reduces the overall strength and even causes loss of bearing capacity. With the modern, wide range, high-quality, and common use of concrete materials, concrete creep has become one of the main problems threatening structural composition, the quality characteristics of concrete, as well as the safety and long-term stability of concrete structures [3][4][5].
The creep of concrete elements can weaken the greatest tensile stress, which reduces the risk of the early cracking of concrete [6][7][8][9][10][11][12]. Evaluation of the crack resistance of concrete should take into account the creep of concrete [13].
Taking into account the creep of concrete under axial compression in calculations within the framework of the classical theoretical approach was first proposed by A.R.  [14][15][16][17][18][19][20][21][22]. Distefano I.N., Prokopovich I.E., and Linnik V.S. applied this method using the nonlinear creep equation of Harutyunyan N.H. and two significantly simplified assumptions: an idealized section in the form of two thin strips; in the function of the nonlinear creep of concrete, the stress was considered to be the same for both bands (when bending the section) [20,23,24]. Such a simplified model in a narrow range allows us to obtain in the final form qualitative nonlinear estimates of the limit of long-term stability.
Creep deformation has a linear dependence on the applied stress in the case of a low stress level (compressive strength less than 30%). However, creep is an unstable phenomenon, since time and stress increase in the case of a high stress level (from 30% to 80% compressive strength); that is, it is a nonlinear dependence that is found in specific creep theories and experimental research data [25][26][27][28][29].
A review of the literature shows that the development of effective unified computational algorithms for solving physically nonlinear problems of the dynamics of short-term and long-term deformation of shells and plates made of composite materials is an urgent problem today [30,31].
The value of long-term resistance R i is established on the basis of local experiments and depends on the choice of empirical correction coefficients: where R l is long-term resistance (by Rzhanitsin) [14], E is a modulus of the deformation of construction materials, I is the moment of inertia of the section of the structure, L is the length of the column for the core elements, and ϕ ∞ is the standard normalised dimensionless creep characteristic of concrete. The value of long-term resistance contains deviations from the results. With the prolonged action of the load, the empirical method dangerously overestimates the conditional critical force by several times concerning the value of the limiting elasticity of concrete. Formula (2) shows the modulus of deformation of concrete and the numerical values of the ultimate creep characteristic ϕ b,cr : where ϕ b,cr is the concrete creep coefficient. The principles of the general nonlinear theory of the calculation of structures are the basis of European standards [29,32]. Figure 1 shows that according to this concept, on the top line of the diagram, the design is considered to be fully compliant with safety requirements until the load E, increasing continuously from a point 0, will not reach a certain level E 0 , beyond which the design does not meet the safety requirements.  [20,23,24]. Such a simplified model in a narrow range allows us to obtain in the final form qualitative nonlinear estimates of the limit of long-term stability. Creep deformation has a linear dependence on the applied stress in the case of a low stress level (compressive strength less than 30%). However, creep is an unstable phenomenon, since time and stress increase in the case of a high stress level (from 30% to 80% compressive strength); that is, it is a nonlinear dependence that is found in specific creep theories and experimental research data [25][26][27][28][29].
A review of the literature shows that the development of effective unified computational algorithms for solving physically nonlinear problems of the dynamics of short-term and long-term deformation of shells and plates made of composite materials is an urgent problem today [30,31].
The value of long-term resistance is established on the basis of local experiments and depends on the choice of empirical correction coefficients: where is long-term resistance (by Rzhanitsin) [14], is a modulus of the deformation of construction materials, is the moment of inertia of the section of the structure, L is the length of the column for the core elements, and ∞ is the standard normalised dimensionless creep characteristic of concrete.
The value of long-term resistance contains deviations from the results. With the prolonged action of the load, the empirical method dangerously overestimates the conditional critical force by several times concerning the value of the limiting elasticity of concrete. Formula (2) shows the modulus of deformation of concrete and the numerical values of the ultimate creep characteristic , : where , is the concrete creep coefficient. The principles of the general nonlinear theory of the calculation of structures are the basis of European standards [29,32]. Figure 1 shows that according to this concept, on the top line of the diagram, the design is considered to be fully compliant with safety requirements until the load , increasing continuously from a point 0, will not reach a certain level , beyond which the design does not meet the safety requirements. is the elastic behavior of the structure, 0 is the limit state of the structure, ℎ is a plastic hinge in the design, and ℎ is the corrected plastic hinge. Theories for calculating the limit state of the bearing capacity of structures: U e is the elastic behavior of the structure, U 0 is the limit state of the structure, U h is a plastic hinge in the design, and U ch is the corrected plastic hinge. The theory of limit equilibrium excludes the indicator U 0 from the design scheme of the limit state of a compressed reinforced concrete structure [33]. This concept is contrary to the principles of the nonlinear theory of the calculation of structures and is shown on the bottom line of the diagram. The limiting state of the structure is assigned by a volitional method at a point U h > U 0 and corresponds to the formation of a plastic hinge (U = U h ). The value U h may be several times greater than the value U 0 . The process of continuous static loading of the structure is carried out only on the interval (0, U e ) of the elastic operation of the structure. At the point U e , either the first crack in the concrete is formed or plastic deformation begins. From point U e in norms, there is an instantaneous jump to the point U h , bypassing the point U 0 , or to the point U ch , corresponding to the plastic hinge.
Contradictions between the limit state of the structure, established by the nonlinear theory of the calculation of structures, and the limit state are presented in the building codes: 1.
Point U 0 -there is a real point completing the process of static loading of a compressed reinforced concrete structure; point U h -there is an unreal point, but it is used in the building codes and has nothing to do with the process of static loading of any compressed reinforced concrete structure; 2.
The nonlinear theory of the calculation of structures recognizes the unsatisfactory state of the structure, corresponding to the site ( U 0 , U h ); building codes recognize the same deformation site (U 0 , U h ) satisfactory condition for construction; 3.
The nonlinear theory of the calculation of structures establishes the rules for calculating the value U 0 and characterizes the limiting state of the structure; 4.
The nonlinear theory of calculating structures naturally rejects the existence of a deflection at the column, with it having no length; the building codes give deflection to the column, with it having no length. This deflection reaches infinite values in absolute magnitude; 5.
Numerical values U 0 and U h differ from each other up to 100%.
The review shows that the existing theories of the reinforced concrete shell have significant shortcomings.
The nonlinear theory of concrete creep within the framework of instantly elastic models has been developed in scientific works [34][35][36].
The nonlinear creep equation for concrete, which describes the biaxial stress state, has the form [34]: where ϕ(t) is the current value of the creep characteristic, E(0) is the modulus of deformations at the initial moment of time, and f (σ i ) is the nonlinearity function, which is selected based on the experimental data depending on the grade of the concrete and other factors. The equation of the nonlinear creep of concrete has the form [35]: where ε(t, t 0 ) is the total deformation of concrete under the action of constant stress applied at time t 0 , S 0 [dσ (t)] is the nonlinearity function of instantaneous deformation, σ(t) is alternating stress, and E M (t) is the modulus of instantaneous deformation. Deformations at time t, C * 0 (t, t 0 ) и C * 0 (t, t) (Equation (6)) are creep measures at the time of loading and observation in the form of N. Harutyunyan's record [37]. C * 0 (t, τ) is a creep measure in the form proposed by S. Alexandrovsky [38].
The equation of nonlinear creep of concrete [36] for the first time takes into account the nonlinearity of the instantaneous deformation of concrete in the framework of heredity hypotheses: where ε(t) is the relative creep deformation of the material at time t, E(t) is the modulus of elastic-instantaneous deformations, F M is a nonlinear stress function for elasticinstantaneous deformations, σ(τ) is an alternating stress, and F n is a nonlinear stress function for creep deformations.
The nonlinear creep of concrete [36] is written as follows: The equation of the nonlinear creep of concrete takes into account the nonlinearity of the instantaneous deformation of concrete and has the form [39]: where . ε b (t) is the creep deformation rate of concrete at time t, E b (t) is the concrete deformation modulus for the current time, . σ b (t) is the rate of stress change in concrete at time t, ϕ Φ t is the creep characteristic of concrete in the form of a Sanzharovsky [34] notation, α Φ is a constant in the form of a Sanzharovsky notation, and E b0 is the modulus of elasticity.
A general theory of calculation of composite physically and geometrically nonlinear thin-walled systems is proposed, extending to a fairly wide class of smooth and reinforced shell and plate structures, including during their reconstruction, operating in a complex stress state at real loading levels, as well as under conditions of nonlinear creep of the material and the presence of cracks.
A mathematical model of the behavior of the mentioned structures at all stages of loading is constructed within the framework of unified systems of differential equations and an algorithm for their solution at different points in time.
A new method for calculating inelastic shell structures with cracks is proposed, and relations for equivalent elasticity parameters are given, which implement quite strictly the problem of taking into account plasticity, half-strength, and cracks in stretched and compressed cross-section zones. The main idea is to replace an inelastic shell with cracks-a solid equivalent elastic shell structure with equivalent elasticity parameters determined by comparing its deformations with similar deformations of a nonlinear model.

Methods
Building codes ensure that the general criteria and design methods comply with the requirements of the nonlinear theory of structural analysis.
Practical adaptation requires overcoming existing inconsistencies in the building codes, the calculated position of which is based on an erroneous model of a plastic hinge [40][41][42].
In the nonlinear theory of structures, the calculation of structures is presented in a simple and convenient form; it corresponds to the deformation model of the section. This measure will allow finding alternative solutions in the design of various structures.
The construction of a nonlinear theory for the calculation of reinforced concrete shells and plates based on generally accepted assumptions is shown. The well-known models of concrete deformation under short-term and long-term loading are considered, disobedient to Hooke's law σ = E·ε.
The nonlinear relationship between stresses and deformations was first considered by G. Bulfinger and F. Gerstner. The power law at k = 1 is a nonlinear dependence, which is written in the following form (8): where A is a constant with the dimension of the stresses and k is an exponent (dimensionless quantity).
The following equations are also used to describe the nonlinear relationship between stresses and strains.
F. Gerstner describes the nonlinearity of materials by parabolic dependence: The functional relationship between stress intensity and strain intensity in the material deformation diagram σ ε is approximated by a cubic dependence [43] in the form: The polynomial function is taken as: Within the framework of the hypothesis of the linear creep of concrete, the creep of particular importance is in the form of: The simplified Maxwell-Kachanov concrete creep formula is taken in differential form: Next, the problems of constructing a nonlinear theory for the calculation of reinforced concrete shells and plates are considered.

Establishing a Functional Relationship between Stresses and Strains
In the first stage, a deformed thin-walled element and its cross-section of unit length are considered, receiving displacements with these U, V, W rotation angles, this ω 1 , ω 2 , ω 3 , elongation strain, and these ε 1 , ε 2 , longitudinal γ and transverse γ 1 , γ 2 shifts [44].
The hypotheses of the theory of plasticity are accepted, and specific equations of fiber deformation in the differential form are obtained with increasing load for the case of short-term loading in the form: .
where E * j σ ij (t) are variable modules of deformations, E c σ ij (t) are secant modules, and E k σ ij (t) are tangent modules of deformations from the concrete deformation diagram σ ij ∼ ε ij . For research on the processes of the creep of materials in structures, any of the equations of creep theories are used; exactly, equations of the theory of ageing in differential form, which have the form: . where ] is the nonlinearity function; and ϕ(t) is a creep characteristic.
The combination of these two laws of short-term and long-term deformation is applicable at all stages of the structure's existence.
Geometric relations are obtained based on a Timoshenko-type shift model: Transverse shear stresses are distributed according to the parabola law in the following form: In this case, the following conditions are met: where k sh is the shape coefficient at shift and z oj (t), z mj (t) are the upper and lower limits of integration within the solid part of the section, depending on the presence or absence of cracks.
The shape coefficient at shift is derived to obtain shear stresses at the center of gravity of the shell cross-section with cracks, depending on the shape of the section.

Creation of a Deformation Model of a Section with Cracks and the Inelastic Properties of Materials
The main vector and the main moment of the stress diagram in a linear section are compiled in the form: where ε oj , ε mj (j = 1, 2), γ 0 , γ m are fiber deformations of the cross-section; σ jj , τ 12 , τ j3 , are the linear stresses in the concrete matrix; σ aj , . . . , σ aj , . . . uF aj , F aj are, stresses and total areas j upper and lower reinforcing bars, respectively; h aj , h aj are coordinates of the reinforcing bars of the location relative to the center of reduction of the internal forces of the section' and a aj , a aj is the thickness of the upper and lower protective layers of concrete. To simplify the epure of the distribution of stresses and strains k, there is an amount of upper and lower layers of reinforcement. To represent the distribution of stresses and strains, Figure 2 shows the location and number of top and bottom reinforcing bars in the structure section.
where , ( = 1, 2), , are fiber deformations of the cross-section; , , , are the linear stresses in the concrete matrix; ′ , … , , … и ′ , are, stresses and total areas upper and lower reinforcing bars, respectively; ℎ′ , ℎ are coordinates of the reinforcing bars of the location relative to the center of reduction of the internal forces of the section' and ′ , is the thickness of the upper and lower protective layers of concrete. To simplify the epure of the distribution of stresses and strains , there is an amount of upper and lower layers of reinforcement. To represent the distribution of stresses and strains, Figure 2 shows the location and number of top and bottom reinforcing bars in the structure section.
where * ( = 1, 2 1 ↔ 2) are fixed ultimate deformations of the composite matrix in tension and compression and п are conditions for the strength of the Equation (21). With numerical integration, the boundary of the occurrence of plastic deformations in the reinforcement corresponds to the achievement of yield strains т on the reinforcement deformation diagram ~ .
When creating a deformation model of the cross-section, the condition for the appearance of a crack during compression and tension was adopted in the form of a condition (criterion) of strength by G. Geniev (21) [45]:  The fixed ultimate deformations in the section are taken subject to the strength condition in the form: where ε * jj (j = 1, 2 1 ↔ 2) are fixed ultimate deformations of the composite matrix in tension and compression and n cond are conditions for the strength of the Equation (21).
With numerical integration, the boundary of the occurrence of plastic deformations in the reinforcement corresponds to the achievement of yield strains ε m on the reinforcement deformation diagram σ a ∼ ε a .
When creating a deformation model of the cross-section, the condition for the appearance of a crack during compression and tension was adopted in the form of a condition (criterion) of strength by G. Geniev (21) [45]: where R cs and R ts are compressive and tensile strength, respectively. To relate the stresses and strains in the fibers, there are equations of the forms (14) or (15). The specific equilibrium equations in the cross-section with cracks are compiled in differential form. A linear system of differential equations with time-varying coefficients is obtained in the form: where [A(x, t)] is a square matrix of coefficients with an unknown system, the order of which depends on the accepted model, and .
x is a vector column of unknown velocities; A p (t, x) ; there is a matrix column of free terms.
The solution of Equation (22) makes it possible to find edge deformations with the help of which all parameters of the stress-strain state of the section are determined.

Derivation of the Integral Relationship between Deformations and Stresses, and Finding Equivalent Elasticity Parameters
Equivalent elastic parameters are obtained based on te calculated edge deformations of the section.
To achieve this, m-the moment relationship between the total deformation of an arbitrary layer-is expressed with a coordinate z + z 0 . This is how the equivalent elasticity parameters are found for different models of a section with cracks and inelastic properties of materials: where L = sinα; n = cosα; . Expressions for elongation, shift, and curvature from the elastic calculation are substituted in the right part of Formula (23), and integration is performed. As a result, equivalent elasticity parameters are found for different models of a section with cracks and the inelastic properties of materials, in particular, for a model of the Timoshenko type: where N rj, M rj, Q r,r are the internal forces in the principal axes; F nj , I nj are the reduced characteristics of the reinforced concrete section; F · S r , S · I r are the integral geometric characteristics of the section; and m is an indicator of the weight of the influence of the extreme fibers of the section. As a result of the integral representation of material properties (24), physical equations of a nonlinear model with load-dependent parameters E, v can be replaced by physical relations with equivalent elasticity parameters.

Description of the Shape of the Deformed Scheme and Nonlinear Analysis of Structures
The shell calculation model is constructed according to a deformed scheme based on differential equilibrium equations.
The features of deformation and the existence of the shell are accounted for. Functional dependencies of changes in external influences and/or the rheological properties of materials corresponding to the real process of loading and deformation are set: where ∂M 1 ∂x + ∂M 12 ∂y − Q 1 = 0; ∂M 12 ∂x + ∂M 2 ∂y − Q 2 = 0. To trace the entire loading process, the law of increasing load is presented in the form: q(t, x, y) = q 0 (x, y) + ∆q(x, y) · t. (26) A resolving system of quasi-static equations of shell motion due to changes in external load or geological processes is constructed. All relations used are differentiated once in time, including physical relationships with equivalent elasticity parameters, taking into account variability. As a result, a linear system of differential equations is obtained concerning the generalized velocities of the deformation parameters: . y J (t) · a jj (t) = a JP (t), (27) where . y J (t) is the generalized velocities of the deformation components; a J J (t) is the variable coefficients, depending on the condition and age of the structure; and a JP (t) is variable functions of changing the external load and/or creep parameters of materials.

Evaluation of the Stability of the Equilibrium under Study and the Establishment of the Criterion for the Loss of Bearing Capacity
To check the stability of the equilibrium of the structure, its varied state is considered. A system of linear equations in variations of the desired quantities is derived. As a result, a homogeneous system of linear equations in variations is obtained with coefficients recalculated at each calculation step: The equality of the determinant of this system to zero determines the condition for the loss of stability of the equilibrium of the structure.
where δy J (t)is the perturbation components and a J J (t)is the virtual equivalent stiffnesses when the structure is perturbed. The proposed mathematical model and algorithm for solving nonlinear problems are implemented in the form of a compiled set of computer programs. The developed software systems have been tested on solving test problems, with satisfactory comparison of shell calculation results with known solutions [46,47] and experimental data [48,49]. They passed the state registration of copyright in the authorized bodies of Kazakhstan [50].
Based on the developed software systems, the operation of structures at all stages of the life cycle is studied. Various effects of the influence of nonlinearities on the behavior of structures are revealed. A study established the dependence of the equivalent parameters of elasticity on the level of loading, the percentage of reinforcement, the cracks formation scheme, and other factors [44].
Solutions for the development of practical recommendations have been obtained on the effective calculation of composite shells and plates in a nonlinear formulation in the concept of European standards. Figure 3 shows a graph of the build-up deflection at the center of the plate. The results of the calculation of the plate according to the linear elastic scheme are presented without taking into account crack formation, according to the proposed theory and also according to experimental research. An insignificant discrepancy between the calculation results according to the proposed theory and experimental data (7-12%) confirms the reliability of the developed nonlinear calculation theory. Taking into account the physical nonlinearity gives a significant correction in bearing capacity values for plates with high content of reinforcement. of the calculation of the plate according to the linear elastic scheme are presented without taking into account crack formation, according to the proposed theory and also according to experimental research. An insignificant discrepancy between the calculation results according to the proposed theory and experimental data (7-12%) confirms the reliability of the developed nonlinear calculation theory. Taking into account the physical nonlinearity gives a significant correction in bearing capacity values for plates with high content of reinforcement.    Short-term ultimate loads for P-1 with the following parameters: a = b = 55 cm, 3 Mpa, R cs = 14.5 Mpa, R ts = 1.45 Mpa, and R a = 327.6 MPa. Figure 4 shows the results of the calculation of reinforced concrete shells according to a linear elastic scheme without taking into account cracks and also taking into account crack formation and the inelastic properties of materials.
taking into account crack formation, according to the proposed theory and also accord to experimental research. An insignificant discrepancy between the calculation resu according to the proposed theory and experimental data (7-12%) confirms the reliabi of the developed nonlinear calculation theory. Taking into account the phys nonlinearity gives a significant correction in bearing capacity values for plates with h content of reinforcement.    In Figure 4a, the triangles show the curves of a numerical and experimental study of the behavior of shells O-1 and O-2 according to a nonlinear elastic scheme. A comparison of them with the rectangle curve obtained according to the proposed theory shows their insignificant discrepancy for particularly flat shells and more significant discrepancy for shells with a large rise. Circles represent curves calculated without taking into account the nonlinear components of deformations and crack formation. The higher the lift and the load level, the greater the discrepancy in the values (and other parameters) obtained from the linear and nonlinear theory of elasticity, taking into account plasticity and crack resistance.
A comparison of the values of short-term limit loads for the shell and plate under consideration, calculated with and without taking into account the nonlinear properties of concrete, is shown in Table 1.
At the same time, short-term limit loads P-1 have the following parameters: a = b = 55 sm, 3 MPa, R cs = 14.5 MPa, R ts = 1.45 MPa, and R a = 327.6 MPa.
Numerical studies have shown that according to the linear elastic scheme, the calculation of the instantaneous critical load q M gives inflated values up to 56%; according to the linear theory of concrete creep, without taking into account cracking, the calculation of the long-term critical load q overestimates its value up to 39%. With an increase in the lifting of the shell, the correction made taking into account the nonlinearity of concrete deformation is most significantly manifested. The dependence of the ratio of the long-term critical load to the corresponding short-term α 0 = q/q M on: parameters of plasticity, creep, cracking, geometric characteristics, initial deflection, support conditions, etc., is revealed.
As a result of numerical experiments, the dependences of equivalent elasticity parameters on the loading level of E эк (η), υ эк (η), the percentage of reinforcement of E эк (µ), υ эк (µ), schemes of crack formation and development, the redistribution of forces, changes in stiffness, etc., are established.

Conclusions
A deformation model of a section with cracks and inelastic properties of materials has been created. A new nonlinear theory for calculating a reinforced concrete shell is proposed based on the replacement of the original shell with a continuous equivalent elastic shell. The behavior of the shell over the entire range of loading and existence is investigated by solving resolving systems of differential equations of equilibrium, motion, and perturbation of an equivalent structure. New scientific results on the work of inelastic composite shell structures have been obtained, and they are in good agreement with the known solutions [46,47] and experienced data [48,49]:

•
To research the stability of an equilibrium state, it is not enough to use geometrically and physically nonlinear dependencies, based on the nonlinear theory of elasticity. Calculation of short-term and long-term critical load according to the linear elastic and nonlinear elastic scheme in comparison with the developed theory gives an overestimation of the values q s and q l ; • The linear elastic calculation scheme overestimates the values q s up to 56%; calculation according to the linear theory of creep without taking into account crack formation overestimates the values q l up to 39%; • The installed dependency α = q l /q s from the load level, the parameters of plasticity, creep, reinforcement percentage, the crack formation scheme, geometric characteristics, initial deflection support conditions of support and other factors; • Taking into account that physical nonlinearity gives a significant correction in the values of the ultimate load q u for shells with high lift and reinforcement content, for very shallow shells, it is essential to take geometric nonlinearity into account.
Numerical studies are given on the basis of the proposed theory, algorithms, and software package. The results obtained in the course of numerical experimental studies have an insignificant error of 5-10%, which is in satisfactory agreement with the data of field experiments on the models of other authors. Consequently, the application of the theory gives quite acceptable results in solving problems of the theory of concrete creep.
The developed nonlinear theory of calculation of composite physically and geometrically nonlinear thin-walled systems can be applied to the calculation of structures made of concrete, metal, and other structural materials operating in a complex stress state at real loading levels, as well as in conditions of nonlinear creep of the material and the presence of cracks.
Other possible areas of coverage of the proposed theory can be the calculation of reinforced shells and plates on a nonlinearly deformable base, the calculation of reinforced concrete prismatic systems, the calculation of shells and plates of through-sections made of structural materials under conditions of nonlinear deformation and creep, and the calculation of reinforced shells and plates reinforced with stiffeners.
The circumstances listed above show that the codes of reinforced concrete structures must be brought into line with the principle general nonlinear theory of the calculation of structures (European standard), which will require significant material costs and organizational efforts [51]. However, this is justified, as it facilitates the flow of foreign investment into the country, improving the quality of manufactured materials and structures and objects under construction, which means improving the quality of life of people [52][53][54].

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.