Open Access
BIO Web Conf.
Volume 17, 2020
International Scientific-Practical Conference “Agriculture and Food Security: Technology, Innovation, Markets, Human Resources” (FIES 2019)
Article Number 00133
Number of page(s) 6
Published online 28 February 2020

© The Authors, published by EDP Sciences, 2020

Licence Creative Commons
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

The development of the agro-industrial complex of Russia requires the introduction of modern computing digital technologies in engineering and calculation practice of designing structures of the agro-industrial complex, many of which belong to the class of thin-walled. The examples of such structures include irrigation, drainage, watering pipelines, tanks, bunkers for various purposes and much more.

The search for optimal shapes and sizes of such thin-walled structures leads to the need to simulate the deformation processes of the latter based on the use of finite element shear models for the analysis of their stress-strain state (SSS).

In most common currently, computer complexes, such as Ansys and the like, use simplified finite element models constructed on the basis of the standard FEM interpolation procedure, a separate component of the displacement vector components as scalar fields [125]. The above interpolation procedure allows obtaining correct solutions only when using a Cartesian coordinate system.

At the same time, when calculating thin-walled structures of water management systems, which are shells or their fragments, it is convenient to use curvilinear coordinate systems, which introduce additional problems into computational algorithms, such as the problem of accounting for displacements of a finite element as an absolutely rigid body. The solution of this problem is possible by using the interpolation procedure components of the displacement vector as components of vector fields.

The article describes the main stages of construction of finite element models of shear of thin-walled constructions calculation of AIC based on the use of quadrilateral element discretization, the stiffness matrix which is based on invariant interpolation procedure component of the displacement vector and the component vector of angles of rotation of normals as components of vector fields.

2 Materials and methods

2.1 Geometric relations

When obtaining the dependencies between the components of the strain tensor and the components of the displacement vector in the calculation of thin-walled structures of water systems, the hypothesis of a straight normal [26] is usually used, according to which the displacement vector of the point M, spaced from the median surface by some distance ζ, is determined by the formula


where a and a° are orts of normals to the median surface in the deformed and initial States.

Included in (1) the displacement vector of the point M0 of the median surface can be represented by the components of the local basis of this point


The use of (1) entails the appearance in the Cauchy relations of the second derivatives of the normal component of the displacement vector [26, 27]. In addition, in such relations there are no deformations that take into account the transverse shear, which is necessary in the calculations of short-span structures, as well as structures with rigid support, sliding and other types of seals that prohibit the angle of rotation. Therefore, the most appropriate is to write the displacement vector V of the point M0 ζ, spaced from the median surface at a distance ζ, in the following formulation [28]


where is a vector of angles of rotation normal [3]

After executing the vector product, the formula (3) takes the form


where .

The position of the point of the thin-walled structure before and after deformation is determined by the radius vectors


where R0 is the radius vector of the point of the median surface of the thin-walled structure.

Differentiation (5) on global coordinates a, P vectors of bases of any point of a thin-walled construction are defined


By scalar products (6) covariant components of the metric tensor can be obtained before and after deformation of the thin-walled structure


Deformations at an arbitrary point of a thin-walled structure can be calculated by the continuum mechanics ratio [29]


Here and below, the subscripts a and P consistently take the values 1, 2, 3. And 1 and 2 correspond to the surface curvilinear coordinates, and the Figure 3-the coordinate in the direction perpendicular to the surface.

1.2 Interpolation procedure in a quadrangular finite element

The thin-walled construction is modeled by a set of quadrangular fragments of the middle surface with nodes l, m, n, p. The components of the displacement vector (2) and their partial derivatives with respect to the curvilinear coordinates of the surface α and β, as well as the components of the rotation angle vector of the normal γ (4) were chosen. Thus the column of nodal unknowns in the global curvilinear coordinate system has the following form



Here, w refers to the component of the displacement vector v1, v2 or v.

When determining the SSS of thin-walled structures with curvilinear surface forms, it is most correct to apply the vector form of the interpolation procedure, according to which the following interpolation dependences for the displacement vector and the normal rotation angle vector are used

where – matrices – rows of form functions

whose elements are represented by products of Hermite polynomials of the third degree and bilinear functions of local coordinates -1 ≤ ξ, η ≤ 1, respectively; [M] – transition matrix from nodal unknowns in the local coordinate system to nodal unknowns in the global coordinate system.

The columns of vector nodal unknowns included in the right parts (10) have the following form


Representing the displacement vectors and the rotation angle vectors of the normal (11) components of the nodal local bases, the relations (10) can be written as


Included in (12) matrices [Av] and [Aγ], containing vectors of local bases of finite element nodes, can be represented as matrix sums(13)

Taking into account (13), (2) and (4), the ratios (12) will take the following form(14)

From (14) it is possible to obtain interpolation expressions for the displacement vector component and the normal rotation angle vector component when implementing the vector form of the interpolation procedure


Further formation of the stiffness matrix and the external load column of the finite element is carried out by minimizing the Lagrange functional in a standard FEM manner [125].

3 Calculation example

An elliptical cylinder loaded in the middle with a concentrated force P=453.6 N and having a hinged support on the diametrically opposite side preventing vertical displacement was calculated (Fig. 1). At the ends of the cylinder has a sliding seal. The initial data had the following values: cylinder length L=26.29-10−2 m; shell thickness h=0.24-10−2 m; elastic modulus E=0.738-105 MPa; Poisson’s ratio ν =0.3125. Due to the presence of symmetry planes, the % part of the shell was calculated.

The calculations were performed according to two variants. In the first variant, the standard interpolation procedure of the components of the displacement vector and the rotation angle vector of the normal as scalar values was used for the arrangement of the stiffness matrix and the nodal forces column of the quadrangular FE [125]. In the second variant, interpolation of the components of the displacement vector and the normal rotation angle vector as components of the vector fields (10) – (15) was implemented.

Was initially designed in a circular cylinder, i.e. the parameters of the ellipse which is the cross section of the cylinder was taken equal to b=c=12.58-10−2 m. the calculation Results for two variants of the interpolation procedure are presented in table № 1, which shows the variation of the deflection under concentrated force (point A) and the value of the normal σ11, σ 22 and σ 13 shear, stresses at points 1 and 2 at the end of the shell to the inner σ int, and σ out of external surfaces of the shell depending on the density of the mesh discretization of shell.

The analysis of the data presented in Table 1 shows that in both variants of calculation there is a satisfactory convergence of the computational process, and the variation values of deflection and stress are almost the same or close enough to each other for the same sampling grids. In addition, it should be noted the practical coincidence of the stress values at points 1 and 2 in both versions of the calculation, which is due to the symmetry of the chosen design scheme of the shell.

Further, the calculation of the thin-walled structure of the AIC in the form of an elliptical cylinder was performed at the ratio of the cross-sectional parameters b/c = 4/1. All other original data retained the same values.

The results of the calculation of the elliptical cylinder are presented in Table 2, the structure of which coincides with Table 1. An analysis of the data in Table 2 shows that there are very significant differences between the two calculation options. Thus, in the variant of scalar interpolation procedure there is a slow rate of convergence of the computational process, which can not be considered satisfactory.

The values of deflection under the concentrated force, as well as the stress values were significantly lower compared to the option of using the vector method of interpolation procedure. It should also be noted a significant difference between the stress values at points 1 and 2 of the cylinder in the first version of the calculation.

In the second variant of calculation it is possible to observe steady convergence of computing process. In addition, it should be noted that the stress values at points 1 and 2 of the cylinder are exactly the same, regardless of the sampling grid used.

Symmetry of the computational scheme of the elliptic shell (Fig. 1) also assumes the equality of normal stresses at the point A of the application of the concentrated force P and at the point B of the hinge support. Figures 2 and 3 show graphs of changes in the normal stress ratios

and at points A and B, depending on the ratio of the semiaxes of the ellipse b/a, which is the cross-section of the cylinder, with a scalar (Fig. 2) and vector (Fig. 3) variants of the interpolation procedure.

The sampling grid 21x21 was used in the construction of these graphs. As can be seen from Figure 2, when using the scalar interpolation procedure, and and are equal to 1 only when the v/s ratio is one and two. With increasing the ratio v/s, with two to four values and dramatically increase to unacceptable values. Using a vector interpolation procedure and remains equal to 1 at any ratio v/s (Fig. 3).

thumbnail Fig. 1.

Design scheme of the pipeline fragment.

Table 1.

Stresses in characteristic sections of a circular cylinder

Table 2.

Stresses in characteristic sections of an elliptical cylinder

thumbnail Fig. 2.

Comparative diagram in the conventional interpolation procedure

thumbnail Fig. 3.

Comparative chart when the vector forms the interpolation

3 Conclusion

The following conclusions can be drawn from the analysis of tabular data.

  1. The use of the standard scalar interpolation procedure [125] allowed obtaining the correct VAT parameters only when calculating the circular cylinder,i.e. the simplest version of the thin shell.

  2. When calculating thin-walled constructions agriculture, non-circular cylinders (for example, elliptical cylinders), to obtain the correct values of VAT needed to apply the vector version of the interpolation procedure because the interpolation component of the displacement vector and the component vector of the rotation angle of the normal components of vector fields allows considering not only changes in the aforementioned component, but also a change of the local basis vectors point in the middle surface of a thin shell in the process of its deformation when using curvilinear coordinate systems.


The study was supported by the Russian Foundation for Basic Research and the Volgograd Region Administration as part of research project No. 19-41343003 r_mol_a.


  • K.-J. Bathe, F. Brezzi, L.D. Marini, The MITC9 shell element in plate bending: mathematical analysis of a simplified case Comput. Mech. 47(6), 617–626 (2011) [Google Scholar]
  • A.I. Golovanov, O.N. Tyuleneva, A.F. Shigabutdinov, Finite element Method in the statics and dynamics of thin-walled structures (Fizmatlit, Moscow, 2006) [Google Scholar]
  • R.A. Kayumov, I.Z. Muhamedova, F. R. Shakirzjanov, B.F. Tazyukov, Parameter determination of hereditary models of deformation of composite materials based on identification method J. of Phys.: Conf. Ser. Current Probl. 012006 (2018) [CrossRef] [Google Scholar]
  • O.C. Zienkiewicz, R.L. Taylor, The Finite Element Method for Solid and Structural Mechanics (Elsevier, Oxford, 2005) [Google Scholar]
  • V.K. Belov, L.P. Zheleznov, T.S. Ognyanova, Nonlinear deformation and stability of a noncircular cylindrical shell under combined loading with bending and twisting moments J. of Applied Mech. and Techn. Phys. 59(2), 352–358 (2018) [CrossRef] [Google Scholar]
  • Nguyen Nhung, Waas Anthonym, Nonlinear, finite deformation, finite element analysise ZAMP. Z. Angew. math. and Phys. 67(9), 35/1–35/24 (2016) [CrossRef] [Google Scholar]
  • S.L. Paznanova, G.P. Vasilev, P.S. Dineva, G. D. Manolis, Dynamic analysis of nanoheterogeneities in a finite-sized solid by boundary and finite element methods Int. J. Solids and Struct. 80, 1–18 (2016) [CrossRef] [Google Scholar]
  • P. Wriggers, B. Reddy, W. Rust, B. Hudobivnik, Efficient virtual element formulations for compressible and incompressible finite deformations Comput. Mech. 60, 253–268 (2017) [Google Scholar]
  • A.S. Solodovnikov, S.V. Sheshenin, Numerical study of strength properties for a composite material with short reinforcing fibers Moscow Univer. Mech. Bull. 72(4), 94–100 (2017) [CrossRef] [Google Scholar]
  • V.N. Paimushin, S.A. Kholmogorov, M.V. Makarov, D.V. Tarlakovskii, A. Lukaszewicz, Mechanics of fiber composites: forms of loss of stability and fracture of test specimens resulting from three-point bending tests ZAMM. Zeitschrift fur Angewandte Mathem. und Mech. 99(1), e201800063 (2019) [CrossRef] [Google Scholar]
  • L.U. Sultanov, Analysis of finite elasto-plastic strains. Medium kinematics and constitutive equations Lobachevskii J. of Mathem. 37(6), 787–793 (2016) [CrossRef] [Google Scholar]
  • A.I. Abdrakhmanova, L.U. Sultanov, Numerical modelling of deformation of hyperelastic incompressible solids Mat. Phys. and Mech. 26(1), 30–32 (2016) [Google Scholar]
  • L.P. Zheleznov, V.V. Kabanov, D.V. Boiko, Nonlinear deformation and stability of discrete reinforced elliptical cylindrical composite shells under torsion and internal pressure Russian Aeronautics 61(2), 175–182 (2018) [CrossRef] [Google Scholar]
  • I.V. Stankevich, S.S. Volkov, Algorithms for solving MDTT boundary value problems taking into account creep deformation Mathematics and mathematical modeling 1, 1–14 (2018) [CrossRef] [Google Scholar]
  • A.M. Belostotsky, P.A. Akimov, I.N. Afanasyeva, T.B. Kaytukov, Contemporary problems of numerical modelling of unique structures and buildings Int. J. for Computational Civil and Structural Engin. 13(2), 9 (2018) [CrossRef] [Google Scholar]
  • A.M. Belostotsky, A.L. Potapenko, P.A. Akimov, Universal software system “Stadyo” for the numerical solution of linear and nonlinear problems of the field theory, statics, stability and dynamics of spatial combined systems: general parameters and superelemental features Int. J. for Computational Civil and Structural Engin. 14(3), 26–41 (2018) [CrossRef] [Google Scholar]
  • V.N. Paimushin, I.B. Badriev, S.A. Kholmogorov, Consistent equations of nonlinear multilayer shells theory in the quadratic approximation Lobachevskii J. of Mathem. 40(3), 349–363 (2019) [CrossRef] [Google Scholar]
  • V.V. Lalin, E.V. Zdanchuk, Conditions on the surface of discontinuity for the reduced cosserat continuum Mat. Phys. and Mech. 31 (1-2), 28–31 (2017) [Google Scholar]
  • A.V. Ignatiev, V.A. Ignatiev, E.A. Gamzatov, Calculation of thin plates by the finite element method in the form of the classical mixed method with the exclusion of displacements of finite elements as a rigid whole in Proc. of higher educational institutions. Construction 3(711), 5–13 (2018) [Google Scholar]
  • V.V. Galishnikova, P.Ja. Pahl, Analysis of frame buckling without sidesway classification Structural mech. of engin. Struct. and struct. 14(4), 299–312 (2018) [Google Scholar]
  • K.A. Skoptsov, S.V. Sheshenin, V.V. Galatenko, A.P. Malakho, O.N. Shornikova, V.V. Avdeev, V.A. Sadovnichy, Particle simulation for predicting effective properties of short fiber reinforced composites Int. J. of Applied Mech. 8(2), 1650016 (2016) [CrossRef] [Google Scholar]
  • L. Stupishin, A. Kolesnikov, K. Nikitin, Variable form forming investigation for flexible shallow shells on circular base Asian J. of Civil Engin. 18(2), 163–171 (2017) [Google Scholar]
  • Yu.Ya. Tyukalov, Finite element models in stresses for plane elasticity problems Civil Engin. J. 1(77), 23–37 (2018) [Google Scholar]
  • G.A. Genius, K.P. Pyatikrestovsky, Issues of long durability of precious reinforced concrete shells Struct. Mech. and calculation of struct. 1(270), 14–19 (2017) [Google Scholar]
  • V.G. Bazhenov, M.N. Zhestkov, Applicability of a structurally orthotropic model in problems on tension, bending, and stability of densely perforated plates and shells Mech. of Compos. Mat. 53(2), 159–164 (2017) [CrossRef] [Google Scholar]
  • V.V. Novozhilov, Theory of thin shells 380 (St. Petersburg university, St. Petersburg, 2010) [Google Scholar]
  • R.B. Rickards, Finite element Method in shell and plate theory 248 (Zinatne, Riga, 1988) [Google Scholar]
  • K.F. Chernykh, Linear theory of shells, 1, 374 (LSU publ. House, Leningrad, 1962) [Google Scholar]
  • L.I. Sedov, Continuum Mechanics, 1, 536 (Nauka, Moscow, 1976) [Google Scholar]

All Tables

Table 1.

Stresses in characteristic sections of a circular cylinder

Table 2.

Stresses in characteristic sections of an elliptical cylinder

All Figures

thumbnail Fig. 1.

Design scheme of the pipeline fragment.

In the text
thumbnail Fig. 2.

Comparative diagram in the conventional interpolation procedure

In the text
thumbnail Fig. 3.

Comparative chart when the vector forms the interpolation

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.