Issue 
BIO Web Conf.
Volume 17, 2020
International ScientificPractical Conference “Agriculture and Food Security: Technology, Innovation, Markets, Human Resources” (FIES 2019)



Article Number  00200  
Number of page(s)  5  
DOI  https://doi.org/10.1051/bioconf/20201700200  
Published online  28 February 2020 
Application of finite elements of various dimensions in strength calculations of thinwall constructions of agroindustrial complex
^{1}
Volgograd State Agrarian University, 400002 Volgograd, Russia
^{2}
Lomonosov Moscow State University, 119991 Moscow, Russia
^{*} Corresponding author: Klotchkov@bk.ru
The article presents a comparative analysis of the effectiveness of the use of finite elements of various dimensions in the study of the stressstrain state (SSS) of objects of the agroindustrial complex (AIC). To determine the strength parameters of the AIC objects, which can be attributed to the class of thinwalled, it is proposed to use a twodimensional finite element in the form of a fragment of the middle surface of a triangular shape with nodes at its vertices. To improve the compatibility of a twodimensional finite element at the boundaries of adjacent elements, it is proposed to use the Lagrange multipliers introduced in additional nodes located in the middle of the sides of the triangular fragment as additional unknowns. It is proposed to use a threedimensional finite element in the form of a prism with triangular bases to study the SSS of agricultural objects of medium thickness and thickwalled. To improve the compatibility of the prismatic element, Lagrange multipliers in the middle of the sides of the upper and lower bases are also used. On the example of calculating a fragment of a cylindrical pipeline rigidly clamped at the ends loaded with internal pressure, the effectiveness of the developed twodimensional and threedimensional finite elements with Lagrange multipliers was proved. The validity of the use of a twodimensional element for researching the SSS of agricultural objects belonging to the class of thinwalled was proved.
© The Authors, published by EDP Sciences, 2020
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
Currently, the agroindustrial complex (AIC) is one of the highest priority areas for the development of the economy of our country. In line with the policy of import substitution, a course has been taken to intensify agricultural production, which provides for the construction of new and reconstruction of existing engineering systems and agricultural facilities. Such systems and objects include irrigation, watering, drainage and drainage systems, tanks, bunkers, storage, silos and other structures, the need for which will only grow At the same time, considerations of material and resource conservation necessitate the development of modern computing technologies to determine the strength parameters of the aforementioned systems and AIC facilities. Such technologies include numerical methods for determining the stressstrain state (SSS) of systems and objects [1–9], in particular, the finite element method (FEM) [10–16].
The construction and troublefree operation remains quite relevant despite a wide selection of finite element computing systems, mainly foreign ones, the task of developing and improving domestic computing algorithms for determining the strength parameters of AIC objects for their rational design.
The article presents a comparative analysis of the effectiveness of the use of finite elements of various dimensions in determining the SSS fragments of agribusiness systems.
2 Materials and methods
2.1 Twodimensional finite elements
When modeling the processes of deformation of thinwalled objects of the agroindustrial complex, such as pipelines for various purposes, tanks, reservoirs, bunkers and others, twodimensional discretization elements, for example, of a triangular shape (Fig. 1), can be used [17, 18].
Consider a triangular fragment of the middle surface, for example, a cylindrical pipeline, with nodes i, j, k at its vertices. To organize the procedure of numerical integration over the area of an element, each such fragment will be mapped onto a right triangle with a local coordinate system O < ξ,n < 1. As nodal variable parameters, we take the components of the displacement vector and their derivatives. Thus, the column of nodal unknowns in the local coordinate system ξ, n will have the form(1)
where .
Here, by q we mean the tangential u, v or normal w component of the displacement vector.
Taking into account (4), the conditional Lagrange functional is formed
In order to solve the compatibility problem of the triangular discretization element at the boundaries of adjacent elements, it is proposed to introduce Lagrange multipliers at nodes 1, 2, 3 located in the middle of the sides as additional nodal unknowns. Using the mentioned Lagrange multipliers, the equality of the derivatives of the normal components of the displacement vectors calculated along the normals to the sides of the triangular element at nodes 1, 2, 3 is ensured(2)
where the index l takes the value 1, 2, 3, and the prime indicates the values calculated from the side of the adjacent finite element (Fig. 2).
Based on (2), for the separately considered triangular discretization element, the following condition can be written(3)
which can be represented in matrix form(4)
where matrix contains form functions consisting of complete polynomials of the third degree.
Taking into account (4), the conditional Lagrange functional is formed(5)
where – matrix rows containing deformations and stresses at an arbitrary point of the thinwalled structure of the agroindustrial complex; – row matrices containing components of the displacement vector and the external surface load vector; – column of unknown nodal unknowns in the selected global curvilinear α, β coordinate system; – transition matrix from column to column .
Performing the minimization operation over and over (5), we obtain the stiffness matrix and the column of nodal forces of the triangular discretization element.(6)
where – stiffness matrix of a triangular element; – column of unknown nodal unknowns; – column of nodal forces of external surface load.
Fig. 1. Design scheme of the fragment of the pipeline using twodimensional finite elements 
Fig. 2. Twodimensional finite element 
2.2 Threedimensional discretization element
In order to solve the problem of determining the strength parameters of AIC constructions in a threedimensional formulation (Fig. 3), it is proposed to use a finite element in the form of a prism with a triangular base with nodes i, j, k, m, n, p located at its vertices. In order to perform numerical integration over the volume, each such element will be mapped onto a prism, the bases of which are rectangular triangles with local coordinates 0 ≤ ξ, n ≤ 1. Vertical coordinate ζ varies within −1 ≤ ζ ≤ 1.
The column of the desired nodal unknowns of the prismatic element in the local coordinate system includes the components of the displacement vector, as well as their partial derivatives of the first order with respect to ξ,n, ξ (7)
where .
The interpolation procedure for the components of the displacement vector of a point belonging to the inner region of the prismatic finite element is determined by the expression(8)
where G_{1}(ξ, η)… G_{9}(ξ, η) are the complete polynomials of the third degree;H_{1}(ϛ)…H_{4}(ϛ) – Hermite polynomials of the third degree.
In order to solve the compatibility problem of a volumetric finite element in the form of a prism, we use the technique described in paragraph 1 above. To do this, we introduce the Lagrange multipliers at nodes 1, 2, 3 and 4, 5, 6 located in the middle of the sides of the lower side as additional nodal unknowns and upper bases respectively (Fig. 4). Thus, for the bases of the prismatic element, the following two additional conditions can be formed(9)
which can be represented in matrix form(10)
where {λ_{S}}^{T} = {λ_{4} λ_{5} λ_{6}};[Φ_{1}] и [Φ_{2}] – matrices containing form functions representing the products of the full polynomial of the third degree and Hermite polynomial of the third degree.
The conditional Lagrange functional for a prismatic finite element, taking into account (10), can be formed as follows(11)
where – a column of nodal unknowns of a prismatic element in global α, β, γ coordinate system; – transition matrix from column to column .
Performing the minimization procedure over (11) with respect to , , and , we can compose the stiffness matrix and the column of nodal forces of the prismatic finite element, standard for FEM [19–30].
Fig. 3. Design scheme of a fragment of a pipeline using threedimensional finite elements 
Fig. 4. Threedimensional finite element 
3 Calculation example
As an example, the problem was solved to determine the strength parameters of the design of the AIC in the form of a fragment of a cylindrical pipeline rigidly clamped at the ends and loaded with internal pressure. It should be noted that the calculated fragment of the pipeline, by its geometric characteristics, belongs to the class of thinwalled structures, to which the theory of thin shells with KirchkoffLove hypotheses is applicable.
The calculations were performed in two versions: in the first embodiment, a fragment of the pipeline was modeled by an ensemble of twodimensional finite elements described in paragraph 1; in the second embodiment, a prismatic finite element was used, oriented so that the triangular bases were located on the external and internal surfaces of the pipeline, and the ribs were oriented in the direction normal to the surface of the pipeline.
The analysis of the calculated values of normal stresses on the inner and outer surfaces of the pipeline fragment showed the following. In the span section, normal stresses in the meridional and annular directions have positive values (i.e., tensile strain is observed). The numerical values of stresses in magnitude in this section turned out to be close in both versions of the calculation. In the reference section, where a moment stress state is observed, the results of the variant calculation differ significantly from each other.
Despite the stable convergence of the computational process, in both versions of the calculation, the numerical values of the meridional stresses in the second version of the calculation were underestimated by about 20% compared to the first version.
Moreover, further thickening of the grid of nodes in the second embodiment did not affect the values of the calculated meridional stresses. The values of ring stresses, which are several times less than the meridional stresses, in the second embodiment were approximately 5 and 10% higher than in the first embodiment. However, as noted above, the meridional stresses on the internal and external surfaces of the calculated structure play the most important role in the VAT picture.
Therefore, in the second version of the calculation, a sampling grid was used, which provides for dividing the cylinder in thickness not into one, but into two or even three layers of volume elements. This approach made it possible to achieve the values of meridional stresses in the reference section similar to the values of the first calculation option.
4 Conclusion
The total number of unknowns when using the second version of the calculation turned out to be about five times larger than in the first version, which makes the first option the most preferable when analyzing the SSS of agricultural structures, which can be attributed to the class of thinwalled. When determining the strength parameters of AIC objects related to structures of medium thickness and thickwalled, it is necessary to use threedimensional prismatic finite elements of the second calculation option.
Acknowledgments
The study was carried out with the financial support of the Russian Foundation for Basic Research and the Administration of the Volgograd Region in the framework of the research project No. 1941340005 r_a.
References
 V.A. Levin, I.S. Manuylovich, V.V. Markov, Numerical Simulation of Spinning Detonation in Circular Section Channels, Computat. Mathem. and Mathem. Phys., 56(6), 1102–1117 (2016) [CrossRef] [Google Scholar]
 E.A. Storozhuk, I.S. Chernyshenko, A.V. Yatsura, StressStrain State Near a Hole in a ShearCompliant Composite Cylindrical Shell with Elliptical CrossSection, Int. Appl. Mechan., 54(5), 559–567 (2018) [CrossRef] [Google Scholar]
 I.B. Badriev, V.N. Paimushin, Refined Models of Contact Interaction of a Thin Plate with Positioned on Both Sides Deformable Foundations, Lobachevskii J. of Mathem., 38(5), 779–793 (2017) [CrossRef] [Google Scholar]
 I.B. Badriev, M.V. Makarov, V.N. Paimushin, Contact Statement of Mechanical Problems of Reinforced on a Contour Sandwich Plates with TransversallySoft Core, Russ. Mathem., 61(1), 69–75 (2017) [CrossRef] [Google Scholar]
 A.S. Solodovnikov, S.V. Sheshenin, Numerical Study of Strength Properties for a Composite Material with Short Reinforcing Fibers, Moscow Univer. Mechan. Bull., 72(4), 94–100 (2017) [CrossRef] [Google Scholar]
 V.N. Shlyannikov, A.P. Zakharov, A.V. Tumanov, Nonlinear fracture resistance parameters for elements of aviation structures under biaxial loading, Russ. Aeronaut., 61(3), 340–346 (2018) [CrossRef] [Google Scholar]
 F. Aldakheel, B. Hudobivnik, P. Wriggers, Virtual selection of phase definitions, Int. J. for Multiscale Computat. Engineer. (2018) [Google Scholar]
 L.U. Sultanov, Analysis of Finite ElastoPlastic Strains. Medium Kinematics and Constitutive Equations, Lobachevskii J. of Mathem., 37(6), 787–793 (2016) [CrossRef] [Google Scholar]
 R.A. Kayumov, Postbuckling Behavior of Compressed Rods in an Elastic Medium, Mechan. of Solids, 52(5), 575–580 (2017) [CrossRef] [Google Scholar]
 K.J. Bathe, Finite Element Procedures (Prentice Hall, Englewood Cliffs, 1996) [Google Scholar]
 O.C. Zienkiewicz, R.L. Taylor, J.Z. Zhu, The Finite Element Method: Its Basis and Fundamentals (ButterworthHeinemann, United Kingdom, 2013) [Google Scholar]
 J.T. Oden, Finite Elements of Nonlinear Continua (McGrawHill, New York, 1971) [Google Scholar]
 A.M. Belostotsky, S.B. Penkovoy, S.V. Scherbina, P.A. Akimov, T.B. Kaytukov, Correct Numerical Methods of Analysis of Structural Strength and Stability of HighRise Panel Buildings, Part 1: Theoretical Foundations of Modelling, Key Engineering Materials, 685, 217–220 (2016) [Google Scholar]
 E. Artioli, L.B.D.D. Veiga, C. Lovadina, E. Sacco, Arbitrary order 2d virtual elements for polygonal meshes: Part I, elastic problem, Computat. Mechan., 60, 355–377 (2017) [CrossRef] [Google Scholar]
 E. Artioli, L.B.D.D. Veiga, C. Lovadina, E. Sacco, Arbitrary order 2d for polygonal meshes: Part II, inelastic problem, Computat. Mechan., 60, 643–657 (2017) [CrossRef] [Google Scholar]
 P. Wriggers, W. Rust, B. Reddy, A virtual element method for contact, Computat. Mechan., 58, 1039–1050 (2016) [CrossRef] [Google Scholar]
 Ren Hui, Fast and robust triangular elements for thin plates/shells, with large deformations and large rotations, Trans. ASME. J. Comput. and Nonlinear Dyn., 10(5), 051018/1051018/13 (2015) [Google Scholar]
 Y.V. Klochkov, A.P. Nikolaev, O.V. Vakhnina, T.A. Kiseleva, Using the Lagrange multiplexes in the triangles of the interregolation of displacements, J. of Applied and Industr. Mathem., 11(4), 535–544 (2017) [CrossRef] [Google Scholar]
 L.P. Zheleznov, V.V. Kabanov, D.V. Boiko, Elliptical Cylindrical Composite Shells under Torsion and Internal Pressure, Russ. Aeronaut., 61(2), 175–182 (2018) [CrossRef] [Google Scholar]
 L.P. Zheleznov, V.V. Kabanov, D.V. Boiko, Nonlinear Deformation and Stability of Discretely Reinforced Elliptical Cylindrical Shells under Transverse Bending and Internal Pressure, Russ. Aeronaut., 57(2), 118–126 (2014) [CrossRef] [Google Scholar]
 V.P. Agapov, R.O. Golovanov, K.R. Aidemirov, Calculation of LoadBearing Capacity of Prestressed Reinforced Concrete Trusses by the Finite Element Method, IOP Conf. Ser. Earth and Environmental Sci., 90, 012018 (2017) [Google Scholar]
 Yu.V. Klochkov, A.P. Nikolaev, T.A. Kiseleva, S. S. Marchenko, Comparative analysis of the results of finite element calculations based on an ellipsoidal shell, J. of Machin. Manufact. and Reliab., 45(4), 328–336 (2016) [CrossRef] [Google Scholar]
 Yamashita Hirok, Valkeapaa Antti I., Jayakumar Paramsothy, Syqiyama Hiroyuki, Continuum mechanics based bilinear shear deformable shell element using absolute nodal coordinate formulation, Trans. ASME. J. Comput. and Nonlinear Dyn., 10(5), 051012/1051012/9 (2015) [Google Scholar]
 Nhung Nguyen, Anthonym Waas, Nonlinear, finite deformation, finite element analysis, ZAMP. Z. Angew. math. and Phys., 67(9), 35/135/24 (2016) [Google Scholar]
 Zhen Lei, F. Gillot, Jezeguel, Development of the mixed grid isogeometric ReissnerMindlin shell: quadrate basis and modified quadrature, Int. J. Mech., 54, 105–119 (2015) [Google Scholar]
 A. Javili, A. Mc Bride, P. Steinmann, B.D. Reddy, A curvilinearcoordinate based finite elemental methodology, Comput. Mech., 54(3), 745–762 (2014) [Google Scholar]
 He Xiaocong, Torsional analysis of free vibration of adhesively bonded singlelap joints, Int. J. Adhes. and Adnes., 48, 59–66 (2014) [CrossRef] [Google Scholar]
 S.N. Yakupov, N.M. Yakupov, ThinLayer Films and Coatings, Mechan. of Solids, 857(1), 012056 (2017) [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 Computat. Civil and Struct. Engineer., 13(2), 9 (2018) [CrossRef] [Google Scholar]
 Yu.V. Klochkov, A.P. Nikolaev, T.A. Kiseleva, S.S. Marchenko, Comparative analysis of the results of finite element calculations based on an ellipsoidal shell, J. of Machin. Manufact. and Reliab., 45(4), 328–336 (2016) [CrossRef] [Google Scholar]
All Figures
Fig. 1. Design scheme of the fragment of the pipeline using twodimensional finite elements 

In the text 
Fig. 2. Twodimensional finite element 

In the text 
Fig. 3. Design scheme of a fragment of a pipeline using threedimensional finite elements 

In the text 
Fig. 4. Threedimensional finite element 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.