BIO Web Conf.
Volume 17, 2020International Scientific-Practical Conference “Agriculture and Food Security: Technology, Innovation, Markets, Human Resources” (FIES 2019)
|Number of page(s)||5|
|Published online||28 February 2020|
Calculation of thin-walled axisymmetrically loaded structures of the AIC taking into account FEM-based physical nonlinearity
Volgograd State Agrarian University, 400002 Volgograd, Russia
2 Financial University under the Government of the Russian Federation, 125167 Moscow, Russia
* Corresponding author: firstname.lastname@example.org
The article describes an algorithm for calculating an axisymmetrically loaded shell structure with a branching meridian, taking into account elastic-plastic deformations when loading based on the deformation theory of plasticity without assuming that the material is incompressible during plastic deformations. The correct relations which determine the static conjugation conditions of several revolution shells in the joint assembly are used. A comparative analysis of finite element solutions is presented for various options plasticity matrix development at the loading stage.
© 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.
The theory of thin shells is widely used [1–6]. Due to difficulties in obtaining analytical solutions based on the theory of shells for the designed thin-walled elements, numerical methods have become widespread [7–10]. The action of external loads in the presence of shell meridian branching can cause plastic deformations. Taking into account the nonlinear work of the material, one can evaluate the work of the structure which can develop more economical solutions [11, 12]. When determining the relationships between stress and strain increments based on the deformation theory of plasticity, we used the functional dependence between the average stress and average strain obtained from the tension experiments.
where x is the axis coordinate; r the revolution radius, being function of x;i and k - unit vectors of the Cartesian coordinate system.
In the process of shell deformation, the point of its middle surface is considered in three positions: initial (point M0), after loading steps j (point M, displacement vector v), and after the (j + 1)-th loading step (M*, displacement vector Δv). At a distance ζ from the middle surface, the following pointsM0ζ,MζИM*ζ correspond to it. When taking into account the hypothesis of direct normals, the position of these points is determined by the following radius vectors(9)
where are functions of vectors v, ∆v and coefficients of matrix [n], included in (8).
Deformations and their increments in an arbitrary layer of the shell spaced at a distance ζ from the middle surface can be determined on the basis of equations of continuum mechanics (14)
As a result, deformations and curvatures of the arbitrary point of the shell of revolution included in (14) for j loading steps are expressed by components of the displacement vector and their derivatives(16)
(16) and (17) can be written as
where -matrix of differential and algebraic operators.
Under the axisymmetric loading, the use of a one-dimensional finite element is most effective. A fragment of the shell, bounded by two parallel planes perpendicular to the axis of rotation is displayed on the linear element with coordinates of the nodes η = -1 and η = 1, where η – the local coordinate used for the numerical integration (Fig. 1).
where q a meridional υ1 or normal υ component of the displacement vector.
Discretization element at the loading stage
To take into account the physical nonlinearity of the material used, the relations of the deformation theory of plasticity are used (29)
where – average linear strain; – average linear stress; εi - strain intensity; σi - stress intensity; δij – Kronecker symbol.
where μ- Poisson’s ratio.
Coefficient K in (30) was determined in  using the assumption that, as a result of plastic deformations, the volume does not change.
As can be seen, coefficient K_1 is a function of the parameters of the curve of the material deformation diagram.
where [D] - the matrix of the stress-strain state of the material.
When deriving the FE stiffness matrix at the (j + 1) loading stage of loading, we used the functional F written on the basis of the equality of the work of external and internal forces at a possible displacement(35)
where – the matrix-row of strain increments at the loading stage;
- matrix-rows of stresses and their increments in an arbitrary layer of the shell of revolution;
– the matrix-row of increments of the components of the displacement vector at the (j + 1)-th loading stage;
– surface load and its increments at the (j + 1) -th loading stage.
where [Cpl] is determined from з (34).
The increments of the components of the displacement vector at the (j + 1) th loading step are determined according to (27).
where – the stiffness matrix of the final element of the shell of revolution at the loading stage; – the vector of external loads at the loading stage;
– Newton-Raphson correction.
As an example, the problem was solved to determine the stress-strain state of a complex-jointed shell structure shown in Figure 2.
The following data were used: internal pressure q=0.7 MPa; Poisson’s ratioμ = 0.32; elastic modulus E=7·105 MPa; radius of a cylindrical shell Rc=0.9 m; cylinder length Lc=0.8 m; shell thickness t=0.01m; cone lengths LIK=1.1 m. and Luc=0.2 m; angle α=45°, angle β=30°. The left edge of the structure has a hinged bearing (Fig. 2). A sufficient convergence of the computational process took place when each branch of the structure was divided into 12 FEs.
where σ0 = 200 MPa; εT = 0.0023496 m; p = 18087.03 MPa.
When articulating the three shells of revolution at the branching point of the meridian, the correct kinematic and static conjugation conditions were used: the invariance of the displacement vectors, as well as the equality of the angles of rotation of the normals to the median surfaces before and after deformation .
The calculation was performed in two versions. In the first version (Table 1), when deriving the plasticity matrix at the loading step, we used coefficient K, which determines the relationship between the average linear stress σ0 and the average linear strain ε0, without taking into account compressibility at elastic plastic deformations. In the second version (Table 2), coefficient K1 determined in (32) was applied. The number of loading steps in the first and second versions varied from 30 to 100.
The results of finite element solutions are presented in Tables 1 and 2. They show the numerical values of the meridional stresses at the points of the calculated structure: in the reference section (x = 0.0) and in the branch node of the meridian (x = 0.8 m). Values of the parameters are given for internal ( ) and external ( ) fibers, and for the points of the middle surface of the revolution shell . The branch node represents stresses of the elements of adjacent shells: - for the lower cone and - for the upper one.
Revolution shell with a branching meridian
Calculation option 1.
Calculation option 2.
The results were verified according to the convergence of the computing process and the comparison of the values of the controlled parameters of the stress-strain state with the values calculated using well-known formulas obtained from the equilibrium equation.
The analysis showed that in the reference section of the structure, the convergence of the computational process is satisfactory for both calculation options.
At the branching node of the meridian, that is, in the zone of stress concentration, the results of finite-element solutions differ from each other in the first and second versions of the calculation. In the first version, a significant increase in the values of the meridional stresses in the branching node is observed, which is not reliable with the value of the specified load. With an increase in the number of loading steps in the first embodiment, the controlled parameters of the stress-strain state increase in their absolute values so that the convergence of the computational process is not observed. For the lower cone with a number of steps equal to 100, an increase in the meridional stress indices by 50% is observed in comparison with the values when the number of steps n = 70.
In the second version, there are no sharp increase in stresses in the branching node. In addition, in the second version, according to the calculation scheme (Fig. 2), the lower cone works in tension, and the upper one works in compression. For the first version of the calculation, this correspondence is not observed. With an increase in the number of loading steps in the second version, the values of the meridional stresses at the branching node of the meridian change insignificantly, which indicates the satisfactory convergence of the computational process.
In the reference section, the value of the meridional stress at the points of the middle surface can be calculated based on the equilibrium condition (Fig. 2)(40)
As can be seen from the Tables, both options allow us to obtain the exact solution of the reference section. The calculation error does not exceed 1%.
It can be concluded that the algorithm for the formation of a plasticity matrix at the loading stage, when calculating complex shell structures, allows us to obtain finite element solutions acceptable for the design engineering practice. They can be recommended for using in design systems when analyzing the stress-strain state of agricultural facilities. The relationship between the average linear stress σ0 and the average linear strain ε0 during the plasticity matrix formation can be implemented in the form (32).
The study was financially supported by the Russian Fundamental Research Fund and Volgograd Region Administration No. 19-41-340002 r_a.
- R.A. Kayumov, Postbuckling behavior of compressed rods in an elastic medium, Mechanics of Solids, 52(5), 575–580 (2017) [CrossRef] [Google Scholar]
- I.B. Badriev, V.N. Paimushin, Refined models of contact interaction of a thin plate with postioned on both sides deformable foundations, Lobachevskii J. of Mathem., 38(5), 779–793 (2017) [CrossRef] [Google Scholar]
- V.A. Kozlov, Stress-strain of elements of bridge structures with varying thickness of walls along the length, Russ. J. of Build. Construct. and Architect., 1(37), 67–80 (2018) [Google Scholar]
- V.A. Ignatiev, Calculation of plane frames with large displacement of nodes by the finite element method in the form of the classical mixed method, J. of Construct. and reconstruct., 2(58), 12–19 (2015) [Google Scholar]
- K.P. Pyatikrestovsky, B.S. Sokolov, Nonlinear analysis of statically ideterminate wooden structures and optimization of cross section dimensions of dome ribs, Int. J. for Computat. Civil and Structur. Engineer., 14(4), 130–139 (2018) [CrossRef] [Google Scholar]
- V.V. Karpov, O.V. Ignatev, A.A. Semenov, The stress-strain state of ribbed shell structures, Magazine of Civil Engineer., 74(6), 147–160 (2017) [Google Scholar]
- K.Yu. Bate, Numerical methods (Physic. and mathemat. Lit, Moscow, 2010) [Google Scholar]
- V.P. Agapov, R.O. Golovanov, K.R. Aidemirov, Calculation of load bearing capacity of prestressed reinforsced concrete trusses by the finite element method, IOP Conf. Ser. Earth and Environmental Sci., 90, 01 2018 (2017) [CrossRef] [Google Scholar]
- V.P. Agapov, R.O. Golovanov, Comparative analysis of the simplest finite elements of plates inbending, Advances in Intelig. Syst. and Comput., 692, 1009–1016 (2018) [CrossRef] [Google Scholar]
- Y.V. Klochkov, A.P. Nikolaev, A.Sh. Dzhabrailov, Finite element analysis of axisymmetrically loaded shells of rotation with branching meridian under elastic-plastic deformation, J. of Struct. Mechan. of Engineer. Constr. and Buildings, 3, 50–56 (2013) [Google Scholar]
- E.A. Souza Neto, D. Peric, D.R.J. Owen, Computational Methods for Plasticity, Theory and Applications, Chichester: Wiley Mechan. and Engineer., 268, 704–734 (2008) [Google Scholar]
- C. Miehe, F. Welschinger, F. Aldakheel, Variational gradient plasticity at finite strains Computer Methods in Applied Mechanics and Engineering (2014) [Google Scholar]
- L.I. Sedov, Continuum Mechanics, vol. 1 (Nauka, Moscow, 1976) [Google Scholar]
- M.M. Malinin, Applied Theory of Plasticity and Creep (Engineering, Moscow, 1975) [Google Scholar]
- A.Sh. Dzhabrailov, Yu.V. Klochkov, A.P. Nikolaev, The Finite Elements Analysis of Shells of Revolution with a Branching Meridian, J. Russ. Aeronaut., 52(1), 22–29 (2009) [CrossRef] [Google Scholar]
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.