CC BY
22
2021. 17(1). 51-62
Строительная механика инженерных конструкций и сооружений Structural Mechanics of Engineering Constructions and Buildings
HTTP://JOURNALS.RUDN.RU/STRUCTURAL-MECHANICS
Аналитические и численные методы расчета конструкций Analytical and numerical methods of analysis of structures
DOI 10.22363/1815-5235-2021-17-1-51-62 UDC 69:624.074:624.012.4
RESEARCH ARTICLE / НАУЧНАЯ СТАТЬЯ
Stress state analysis of an equal slope shell
under uniformly distributed tangential load by different methods Olga O. Aleshina1*, Vyacheslav N. Ivanov1, David Cajamarca-Zuniga2
1Peoples' Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow, 117198, Russian Federation 2Catholic University of Cuenca, Av. De las Americas & Humboldt, Cuenca, 010101, Republic of Ecuador *[email protected]
Abstract. Nowadays there are various calculation methods for solving a wide range of problems in construction, hydrodynamics, thermal conductivity, aerospace research and many other areas of industry. Analytical methods that make up one class for solving problems, and numerical calculation methods that make up another class, including those implemented in computing complexes, are used for the design and construction of various thin-walled structures such as shells. Due to the fact that thin-walled spatial structures in the form of various shells are widely used in many areas of human activity it is useful to understand and know the capabilities of different calculation methods. Research works on the study of the stress-strain state of the torse shell of equal slope with an ellipse at the base are not widely available at the moment. For the first time the derivation of the differential equations of equilibrium of momentless theory of shells to determine the normal force Nu from the action of uniformly distributed load tangentially directed along rectilinear generatrixes to the middle surface of the torse of equal slope with a directrix ellipse is presented in this article. The parameters of the stress state of the studied torse are also obtained by the finite element method and the variational-difference method. The SCAD software based on the finite element method and the program SHELLVRM written on the basis of the variational-difference method are used. The numerical results of the parameters of the stress state of the studied torse are analyzed, and the advantages and disadvantages of the analytical method and two numerical calculation methods are determined.
Keywords: thin shell theory, analytical method, momentless state, torse shell, surface of equal slope, finite element method, variational-difference method, SCAD Office computing system, Mathcad system
Article history
Received: October 7, 2020 Revised: January 17, 2021 Accepted: February 4, 2021
For citation
Aleshina O.O., Ivanov V.N., Cajamarca-Zuniga D. Stress state analysis of an equal slope shell under uniformly distributed tangential load by different methods. Structural Mechanics of Engineering Constructions and Buildings. 2021 ;17(1):51—62. http://dx.doi.org/10.22363/1815-5235-2021-17-1-51-62
Olga O. Aleshina, teacher-researcher, assistant of the Department of Civil Engineering of the Academy of Engineering; eLIBRARY SPIN-code: 8550-4986. Vyacheslav N. Ivanov, Professor of the Department of Civil Engineering of the Academy of Engineering, Doctor of Technical Sciences; eLIBRARY SPIN-code: 3110-9909, Scopus Author ID: 57193384761, ORCID iD: https://orcid.org/0000-0003-4023-156X.
David Cajamarca-Zuniga, Docent of the Department of Civil Engineering; eLIBRARY SPIN-code: 6178-4383, ORCID iD: https://orcid.org/0000-0001-8796-4635, WoS ResearcherID: AA0-8887-2020.
© Aleshina O.O., Ivanov V.N., Cajamarca-Zuniga D., 2021
/^N 0 I This work is licensed under a Creative Commons Attribution 4.0 International License https://creativec0mm0ns.0rg/licenses/by/4.Q/
Анализ напряженного состояния оболочки одинакового ската при действии равномерно распределенной касательной нагрузки различными методами
О.О. Алёшина1*, В.Н. Иванов1, Д. Кахамарка-Сунига2
1Российский университет дружбы народов, Российская Федерация, 117198, Москва, ул. Миклухо-Маклая, д. 6 2Католический университет г. Куэнка, Республика Эквадор, 010101, Куэнка, Av. De las Americas & Humboldt *[email protected]
История статьи
Поступила в редакцию: 17 октября 2020 г. Доработана: 17 января 2021 г. Принята к публикации: 4 февраля 2021 г.
Для цитирования
Алёшина О.О., Иванов В.Н., Кахамарка-Сунига Д. Анализ напряженного состояния оболочки одинакового ската при действии равномерно распределенной касательной нагрузки различными методами // Строительная механика инженерных конструкций и сооружений. 2021. Т. 17. № 1. С. 51-62. http://dx.doi.org/10.22363/1815-5235-2021-17-1-51-62
Аннотация. На сегодняшний день существуют различные методы расчета для решения широкого спектра задач в строительстве, гидродинамике, теплопроводности, космических исследованиях и других отраслях. Для проектирования и возведения разнообразных тонкостенных конструкций типа оболочек применяются аналитические методы, составляющие один класс для решения задач, и численные методы расчета, составляющие другой класс, в том числе реализованные в вычислительных комплексах. В связи с тем, что тонкостенные пространственные конструкции в форме разнообразных оболочек широко используются во многих сферах деятельности человека, полезно понимать и знать возможности различных методов расчета. Работы по исследованию напряженно-деформированного состояния торсовой оболочки одинакового ската с эллипсом в основании представлены на данный момент в малом объеме. В статье впервые приводится вывод дифференциальных уравнений равновесия безмоментной теории оболочек для определения нормального усилия Nu от действия равномерно-распределенной нагрузки, направленной по касательной вдоль прямолинейных образующих к срединной поверхности торса одинакового ската с направляющим эллипсом. Также получены параметры напряженного состояния исследуемого торса методом конечных элементов и вариационно-разностным методом. Используются вычислительный комплекс SCAD Office на основе метода конечных элементов и программа SHELLVRM, написанная на базе вариационно-разностного метода. Выполнен анализ числовых результатов параметров напряженного состояния исследуемого торса, установлены плюсы и минусы применения аналитического метода и двух численных методов расчета.
Ключевые слова: теория тонких оболочек, аналитическое решение, без-моментное состояние, торсовая оболочка, поверхность одинакового ската, метод конечных элементов, вариационно-разностный метод, вычислительный комплекс SCAD Office, система Mathcad
Introduction
For the design of diverse engineering structures, various calculation methods are used, such as analytical, numerical and numerical-analytical. In the practice, to get the general parameters of the stress-strain state of spatial-structures, engineers use automated numerical calculation methods because analytical calculation methods are quite complex and time consuming.
The most common numerical calculation method is the finite element method (FEM). Originally, FEM was used for solving mathematical problems in a simpler form. The subsequent development of FEM and automated software systems based on this method such as SIMULIA (www.3ds.com), ANSYS (www.ansys.com), SAP2000 (www.csiamerica.com), SCAD (www.scadsoft.com), PROKON (www.prokon.com) and others, made it
Алёшина Ольга Олеговна, преподаватель-исследователь, ассистент департамента строительства Инженерной академии; eLIBRARY SPIN-код: 8550-4986.
Иванов Вячеслав Николаевич, профессор департамента строительства Инженерной академии, доктор технических наук; eLIBRARY SPIN-код: 3110-9909, Scopus Author ID: 57193384761, ORCID iD: https://orcid.org/0000-0003-4023-156X.
Кахамарка-Сунига Давид, доцент инженерного факультета; eLIBRARY SPIN-код: 6178-4383, ORCID iD: https://orcid.org/0000-0001-8796-4635, WoS ResearcherID: AAO-8887-2020.
possible to apply it to solve a wide range of problems in aerospace research, to model and take into account dynamic loads, to solve various problems in thermal conductivity, hydrodynamics, construction and many other areas.
The idea of discretization on which the FEA is based is very old. Before 1922, Courant used the finite element ideas in Dirichlet's principle. The period of 1962-1972 is known as the golden age of FEM [1]. There are five groups of papers (Courant, Argyris, Turner et al., Clough and Zienkiewicz) which may be considered in the development of the FEM and in one of these the name originated [2]. Clough coined the term "finite elements", Turner perfected the direct stiffness method and the works of Huges, Bathe and Zienkiewicz [3] laid the foundation for further progress of the FEM [1]. In [4], a method for calculating bending plates by the finite element method in stresses is proposed, and a comparison with the results of the finite element method in displacements is made. The solution of plane problems of the theory of elasticity based on the approximation of stresses is considered, the calculations of a cantilever beam and a plate with a hole are performed for various finite element meshes, and comparison is made with solutions by the method of finite elements in displacements and with exact solutions in the work [5]. A special issue including 35 papers is devoted to research in the field of development and application of FEM [6].
The finite-difference energy method (FDEM) [7-10] or so called variational-difference method (VDM) [9-14] is also referred to numerical calculation methods [15; 16]. This method takes into account the geometric characteristics of the middle surface of the shell, which allows a more accurate representation of the stress-strain state of thin-walled structures of complex geometry. The VDM (FDEM) is based on the idea put forward by Courant in 1943 [9; 17; 18], which was continued by Houbolt in 1958 [8], who performed static analysis of beams and plates combining finite difference analog of derivatives with a variational formulation [19]. Further developed by Griffin and Varga in 1963 [20] who introduced finite difference into the variational formulation of strain compatibility and boundary conditions for the analysis of plane elasticity problems [19]. Further Bushnell in 1973, and Brush and Almroth in 1975 [21] who extended the approach to the analysis of other type of structures [22].
The successful application of VDM largely depends on how well the system of basic functions allows the qualitative characteristics of the solution. Consequently, it can be expected that the efficient solution of these variational problems will require numerical schemes that differ from traditional techniques based on continuous approximations [23].
In the Department of Construction of the Academy of Engineering of the RUDN University of Russia, the Doctor of Technical Sciences, Professor V.N. Ivanov together with his postgraduate students (currently PhD) Nasr Younis Ahmed Abboushi (Palestine), Muhammad Rizwan (Pakistan), Bock Hyeng Christian Alain (Cameroon), Govind Prasad Lamichhane (Nepal) led the development of SHELLVRM, a new Variational-Difference Method based program. This program allows to determine the stress-strain state of plates and various types of shells with an orthogonal coordinate system, which middle surfaces are described by analytical equations. The program includes such classes of shells as: flat shells on rectangular and curved planes, shells of revolution, shells in the form of Joachimsthal's channel surfaces, shells in the form of Monge surfaces and normal cyclic surfaces. The program includes a system of plane curves, on the basis of which sections of surface classes are formed and coefficients of quadratic forms are calculated. The basics of the VDM and the text of the program for plate calculations are given in [15].
Analytical calculation methods are used for spatial structures in the form of various surfaces [24]. Analytical methods are quite complex and time-consuming. More than 600 analytical surfaces are described in the Encyclopedia of Analytical Surfaces [25]. The geometry of surfaces and automated possibilities of their construction are considered in the monograph [26].
Among an extensive variety of analytical surfaces, the torse shells of equal slope possesses the ability to unfold onto a plane without folds and breaks [27], and this type of shells are widely used in many areas of industry and manufacturing [28-31].
This article is part of a series of research papers devoted to the study of the geometry and stress state of torse shell of equal slope with an ellipse at the base under the action of different loads. In previous works, the authors have performed calculations this shell under the action of a linear load on the upper edge and under the action of self-weight [32; 33] and with a different restraint of the base ellipse [34]. Also, a design of an awning in the shape of a torse of equal slope was proposed and new results were obtained in the field of geometric studies [35; 36]. In this article, we consider the uniformly distributed load directed along rectilinear generatrixes of the torse. The choice of the load type is determined by the possibilities of the momentless shell theory. The main task of this article is to find an analytical solution and determine the parameters of the stress state of the torse by
the momentless theory (MLT), followed by comparison with the results of two numerical methods (the finite element method and the variational-difference method).
The surface of equal slope is a ruled surface generated by a straight line moving in the normal plane of a flat directrix-curve with a constant angle of inclination to the plane of the directrix. If we take an ellipse as a flat directrix-curve, then straight lines of equal inclination to the plane of the ellipse will generate the torse surface of equal slope (Figure 1). The surfaces of equal slope are surfaces of zero Gaussian curvature (K = 0). The papers [37; 38] describe the basic properties of these surfaces. The equal slope surface also belongs to the class of Monge surfaces [24; 27].
Figure 1. Torse shell of equal slope with an ellipse at the base As its shown in [27], the directrix ellipse can be defined by parametric equations (1):
x = x(v) = a cos v, y = y(v) = b sin v,
(1)
where a and b are the dimensions of the semi-axes of the directrix ellipse at the base of the torse, and the parameter v must be in the limits 0< v < 2n.
According to [27], the parametric form of setting the torse surface with a directrix ellipse is:
x = x(u,v) = acosv —
ubcosacosv
Va2sin2v + b2cos2v
y = y(u, v) = ftsinv
uacosasinv
Va2sin2v + b2cos2v
z = z(u) = usina.
(2)
The family of u lines is the rectilinear generatrixes of the torse surface of equal slope, while the coordinate line u = 0 coincides with the ellipse at the base, a is the angle between the principal normal n = —exk directed inwards of the directrix ellipse and the straight generatrix u (Figure 1).
The coefficients of the basic quadratic forms of a given surface and its main curvatures are [27]:
,,,, B aft sin a A = 1; B = a1/2-u-; F = 0; L = M = 0; N = B-;
ki — ku — 0; k2 —kv —
where = a2 sin2 v + b2 cos2 v, fi = ab cos a.
ab sin a
Bv '
(3)
Let us consider the application of the momentless theory of shell calculation, the finite element method and the variational-difference method on the example of a thin torse shell of equal slope with an ellipse at the base under the action of a uniformly distributed load q = 1 kN/m2 tangentially directed along rectilinear ge-neratrixes to the middle surface of the torse (Figure 2). Thus, the external surface load is X = -q, Y = Z = 0. The geometric parameters of the torse are: a = 3 m, b = 2 m, a = 60°, the length of the straight generatrixes is u = 2 m. The boundary conditions at the level of the directrix ellipse (u = 0 m) are simple (movable) supports, and at the top (u = 2 m) the edge is free.
To determine the parameters of the stress state of the considered torse, the momentless theory of shell calculation [24; 27], the SCAD integrated system for finite element structural analysis (FEA), and the SHELLVRM program based on the variational-difference method [15; 16] are used.
The differential equilibrium equations of the momentless theory are obtained from the general equilibrium equations of the moment shell theory [24; 27].
The momentless theory is a simplified version of the general theory of thin elastic shells, which neglects the influence of transverse forces and moments. At the same time, the possibility of existence of the momentless stress state of the shell depends on a number of conditions [24; 27]. The shell should have the form of a smoothly changing continuous surface, also the load on the shell should be continuous and smooth, and the supports of the edges should allow the shell to move freely in the direction normal to the middle surface, normal movements and rotation angles at the edges of the shell should not be restrained.
We obtain differential equations of equilibrium for determining the normal force under the action of a uniformly distributed load acting in the direction of a tangent along rectilinear generatrixes to the middle surface of the considered torse.
General differential equations of equilibrium of the momentless theory [24; 27] have the form:
z
Figure 2. Torse under the action of external distributed surface load
Differential equations of equilibrium of the momentless torse shell
д
du
л 9B 1 д Г , Л
{ВК) {A2S) + АВХ = 0;
Adv
д
dv
(4)
For the considered case of load application (Figure 2), we obtain X = -q and Y = Z = 0. The differential equations of equilibrium (4), taking into account expressions (3), are transformed as following:
d dB dS
-iBNJ--Nv+Tv + XB = °;
d Nv 1 d
— -0 R„
(5)
The resulting system of differential equations (5) is of second order. To solve it, it is sufficient to have one boundary condition at each point of the torse shell contour. Thus, at the top of the shell at u = 2 m the force Nu = 0. Moreover, from the second and third equations of system (5) the forces Nv = 0 and S = 0.
By integrating the first equation of system (5), we obtain the expression for the values of normal force Nu along the rectilinear generatrixes u:
N =-
u B(u ,v )
J q B(u,
v ) d u + ( v )
(6)
Here v) is an arbitrary function of integration. Then, by integrating of (6):
f
B (u, v) du = uv
1/2 _
u2B u s „ v
-2^ = 2 (B(u,v) + v1/2) = 2^{v-B 2 (u, v)).
(7)
To satisfy the boundary condition Nu = 0 on the upper free edge under u = n = 2 m, the arbitrary function of integration v) in (7) must be equal to:
(8)
The equation (6) for the calculation of numerical values of the normal forces Nu along the rectilinear generatrixes taking into account the value v) of the arbitrary integration function (8) takes the following form:
AL =-1-
u B(u,v )
V1/2(u
(9)
To find numerical results of normal force Nu (9) we use the engineering math software Mathcad.
Numerical methods for investigation of the stress state of the shell
The investigation of the stress state of the torse of equal slope was performed by the finite element method and the variational-difference method. The first calculation is performed by using SCAD software. The view of
the 3D computational model when approximating the middle surface by a set of quadrangular planar shell elements is shown in Figure 2. The maximum distance between the nodes of the finite elements of the computational model is 0.228 m. The number of finite elements is 1680 and of nodes is 1760.
For the implementation of simple (movable) supports, which is a necessary condition for the momentless work of the torse, the SCAD program has added short bar elements with hinges (Figure 3). The introduction of hinges in these support rod elements releases linear movements along the normal to the torse middle surface (Figure 3, direction z\), angular movements tangent to the surface (FiFi ure3 Im |cmcnlm„ Surc 3, direction vi) and normal to the surface of the shells (Figure 3, direction z\), of momentless state as well as angular movements in the direction of rectilinear generatrixes u (Figure 3, in SCAD software direction X:).
The second calculation is performed in the program SHELLVRM, based on the variational-difference method. The calculated grid is similar to the grid in FEM. This calculation also takes into account and implements all the necessary conditions for the momentless state of the shell. The calculation is performed for a 1/4 segment of the torse shell, taking into account two planes of symmetry.
Results and discussion
The obtained results of the analytical calculation are compared with the results of numerical methods (by the finite element method and the variational-difference method) for 11 cross-sections (Figure 4).
Figure 4. Cross-sections of the torse to compare the results
Figure 5. Normal stress o(N„) by FEM, kN/m2
The maximum deviations of the analytical results of normal force Nu along the rectilinear generatrixes from the results of two numerical methods are: 7.4% in section 1-1 (Table 1), 5.0% in section 2-2, 1.9% in section 3-3, 3.7% in section 4-4, 4.1% in section 5-5 (Table 2), 3.6% in section 6-6, 2.8% in section 7-7, 2.2% in section 8-8, 2.0% in section 9-9, 1.9% in section 10-10, and 1.9% in section 11-11 (Table 3). At nodes of coordinates u = 2.00 m in FEM and VRM the values are different from zero when compared with MLT.
For an overall picture of the stress state of torse shell under the action of uniformly distributed load q tan-gentially applied along rectilinear generatrixes to the torse middle surface, the contour graph distribution of normal stress a(Nu) obtained in the SCAD software is shown in Figure 5.
Table 1
Results of normal force Nu: cross section 1-1
^-axis coordinate, m Nu, MLT, section 1-1, kN/m Nu, FEM, section 1-1, kN/m Deviation, MLT and FEM, section 1-1, % Nu, VDM, section 1-1, kN/m Deviation, MLT and VDM, section 1-1, %
0.000 -1.2500 -1.2438 0.50 -1.2980 3.70
0.200 -1.1432 -1.1906 3.98 -1.1920 4.09
0.400 -1.0353 -1.0841 4.50 -1.0860 4.67
0.600 -0.9258 -0.9757 5.11 -0.9767 5.21
0.800 -0.8143 -0.8640 5.75 -0.8645 5.81
1.000 -0.7000 -0.7477 6.39 -0.7480 6.42
1.200 -0.5818 -0.6253 6.96 -0.6255 6.99
1.400 -0.4579 -0.4942 7.35 -0.4944 7.38
1.600 -0.3250 -0.3503 7.22 -0.3504 7.25
1.800 -0.1769 -0.1871 5.45 -0.1870 5.40
2.000 0.0000 -0.0563 - -0.0002 -
Table 2
Results of normal force Nu: cross section 5- 5
^-axis coordinate, m Nu, MLT, section 5-5, kN/m Nu, FEM, section 5-5, kN/m Deviation, MLT and FEM, section 5-5, % Nu, VDM, section 5-5, kN/m Deviation, MLT and VDM, section 5-5, %
0.000 -1.5623 -1.5004 4.12 -1.5340 1.84
0.200 -1.4292 -1.4004 2.06 -1.4010 2.01
0.400 -1.2930 -1.2647 2.23 -1.2650 2.21
0.600 -1.1531 -1.1262 2.39 -1.1270 2.32
0.800 -1.0090 -0.9840 2.54 -0.9842 2.52
1.000 -0.8599 -0.8375 2.68 -0.8376 2.66
1.200 -0.7050 -0.6859 2.78 -0.6860 2.77
1.400 -0.5432 -0.5284 2.80 -0.5285 2.78
1.600 -0.3731 -0.3637 2.58 -0.3638 2.56
1.800 -0.1928 -0.1892 1.88 -0.1899 1.53
2.000 0.0000 -0.0572 - -0.0002 -
Table 3
Results of normal force N : cross section 11- 11
^-axis coordinate, m Nu, MLT, section 11-11, kN/m Nu, FEM, section 11-11, kN/m Deviation, MLT and FEM, section 11-11, % Nu, VDM, section 11-11, kN/m Deviation, MLT and VDM, section 11-11, %
0.000 -1.7778 -1.7439 1.95 -1.7870 0.51
0.200 -1.6159 -1.6240 0.50 -1.6240 0.50
0.400 -1.4512 -1.4585 0.50 -1.4590 0.53
0.600 -1.2833 -1.2899 0.51 -1.2900 0.52
0.800 -1.1122 -1.1178 0.50 -1.1180 0.52
1.000 -0.9375 -0.9420 0.47 -0.9420 0.48
1.200 -0.7590 -0.7622 0.42 -0.7622 0.42
1.400 -0.5763 -0.5783 0.35 -0.5784 0.36
1.600 -0.3892 -0.3901 0.23 -0.3902 0.26
1.800 -0.1972 -0.1973 0.07 -0.1974 0.10
2.000 0.0000 -0.0553 - 0.0000 -
Comparison of the obtained results of normal force Nu by three calculation methods shows good convergence. The concentration of the largest deviations of the numerical values of the normal force Nu by the mo-mentless theory from the VDM and FEM is in the region of the shell with the largest change in the radius of curvature along the curvilinear directrices, i.e., in the upper nodes of sections 1-1 and 2-2 (Figure 4).
According to the Theory of Strength of Materials, the numerical values of the normal force Nu under the action of uniformly distributed load tangentially along rectilinear generatrixes to the torse middle surface at the nodes of all sections at coordinate u = 2.00 m must be Nu = 0. However, the values of the normal force Nu in the FEM and VDM are different from zero, and the results of the VDM are more accurate compared to the FEM. It is well known that the accuracy of the results of FEM and VDM calculations depends on the correct choice of the size of the finite elements (mesh). Moreover, it is noted in [15] that a comparison of the results of VDM and FEM calculations at the same mesh shows close accuracy, and in some cases, VDM gives even higher accuracy results.
FEM and VDM allow obtaining numerical values also for normal forces Nv along curved directrices, bending moments Mu, Mv, tangential forces S and shear forces Qu, Qv. The normal forces Nv by VDM range from -0.0246 to 0.0148 kN/m, and by FEM ranges from 0.0490 to 0.0216 kN/m. The shear forces Qu, Qv range from -0.01 to 0.01 kN/m. The tangential forces S ranges from -0.0354 to 0.0354 kN/m by FEM, and from -0.0162 to 0.0067 kN/m by VDM. The values of bending moment Mu range from -0.0261 to 0.4244 Nm/m by VDM, and by FEM from -0.1143 to 0.4733 N m/m. The values of bending moment Mv by VDM range from -0.4140 to 1.4030 Nm/m, and by FEM range from -0.4398 to 1.5562 Nm/m.
The bending moments Mu and Mv are of particular interest, since the values of bending stresses when compared with normal stresses can be used to infer the bending state of the torse shell under the action of the considered load.
The normal stress oM and oN from normal forces Nu,v and moments Mu v are determined as follows:
— (10)
-0.1143 -0.04Ш
0,0326
0;0326 0.106Î
§HSW 0,1755
$1795 шш
"ЙЩЗ о ,зш
0.3Щ- 0,3999
Figure 6. Bending moment Mu by FEM, N m/m The results of the VDM for the maximum ratio of stresses gm,, to gn,. are: in the cross section 1-1 is
156.5% in the node of coordinate u = 2.00 m, 27.2% in the node with the coordinate u = 1.80 m, 13.1% in the node with the coordinate u = 1.60 m, 7.2% in the node with the coordinate u = 1.40 m, in other nodes does not exceed 4.4%. In cross section 2-2 is 61.9% in the node of coordinate u = 2.00 m, 25.0% in the node with coordinate u = 1.80 m, 12.3% in the node with the coordinate u = 1.60 m, 6.8% in the node with the coordinate u = 1.40 m, other nodes do not exceed 4.2%. In section 3-3 is 25.9% in the node with the coordinate u = 2.00 m, 18.5% in the node with the coordinate u = 1.80 m, 9.9% in a node with coordinate u = 1.60 m, and in the other
nodes does not exceed 5.8%. In section 4-4 is 7.0% in node with coordinate u = 2.00 m, 9.3% in the node with coordinate u = 1.80 m, 6.7% in the node with coordinate u = 1.60 m, and in all other nodes does not exceed 4.5%. In section 5-5 is 45.5% in the node with the coordinate u = 2.00 m; in section 6-6 is 13.3% in the node with the coordinate u = 2.00 m; in section 7-7 is 17.4% in the node with the coordinate u = 2.00 m; in section 8-8 is 10.9% in the node with the coordinate u = 2.00 m; in section 10-10 is 8.9% in the node with the coordinate u = 2.00 m; in section 11-11 is 34.8% in the node with the coordinate u = 2.00 m, 5.8% in the node with the coordinate u = 1.80 m. In other nodes of sections 5-5 to 11-11, the stress ratio does not exceed 5.2%.
Figure 7. Bending moment Mv by FEM, Nm/m
The bending stresses oMv arising from a uniformly distributed load directed tangentially along rectilinear generatrixes to the middle surface, in the VDM have an even greater influence on the bending state of the considered torse shell with a directrix ellipse at the base.
The results of studying the influence of bending stresses and in FEM show a similar character. Figures 6 and 7 show the contour graph distribution of bending moments Mu and Mv obtained in the SCAD software.
Conclusion
The research is carried out at the Academy of Engineering of the Peoples' Friendship University of Russia (RUDN University). In the field of geometry and stress-strain state of various shells, in particular torse shells class, works at RUDN University have been carried out since 1960's. An undeniable contribution to modern theory of shells was made by Prof. V.G. Rekach, Prof. S.N. Krivoshapko and Prof. V.N. Ivanov and their postgraduate students (today PhD in Technical Sciences). Currently, S.N. Krivoshapko and V.N. Ivanov continue their research in the field of shell theory [39-41].
This paper for the first time presents the differential equations of equilibrium for a torse shell of equal slope with a directrix ellipse and the expression for the normal force Nu determination under the action of uniformly distributed load tangentially directed along rectilinear generatrixes to the torse middle surface.
Determination of the internal force Nu of the investigated torse shell by the analytical method is a complex and time-consuming task that requires a lot of time and increased concentration of attention on its implementation, since a slight inaccuracy can lead to incorrect results. The comparison of the results of the momentless theory with the results of the finite element method and the variational-difference method shows good convergency, which indicates the correctness of the obtained differential equilibrium equations and the expression for determining the values of the normal force Nu. The use of SHELLVRM and SCAD programs simplifies the solution of this task. However, the calculation in the SHELLVRM program is possible if there is the program text for its implementation, and in the SCAD program it becomes difficult to implement a momentless state (introduction of simple-movable supports). When choosing a method of solving the problem, SCAD program, based on the finite element method, is the simplest and most versatile way for solving the research problem.
The values of the normal force Nu along the rectilinear generatrixes of the shell indicate that the considered shell is working in compression. Thus, it is a big plus when selecting materials for the design and manufacture of torse shells. Considering the property of this class of shells to be flattened on the plane without folds and breaks, this is also an advantage when selecting torse shells among similar shaped.
Due to the results of the FEM and VDM, it was found that the bending stresses oMu and aMv have a significant influence on the torse shell stress state. Therefore, it is necessary to consider the bending moments Mv and Mu when designing different structures in the form of this class of shells. The momentless theory does not allow us to obtain these parameters of the stress state of the torse. Thus, it may be concluded that the momentless theory is not applicable for the considered torse shell of equal slope with ellipse directrix.
References
1. Sabat L., Kundu C.K. History of finite element method: a review. Recent Developments in Sustainable Infrastructure. 2021:395-404. https://doi.org/10.1007/978-981-15-4577-132
2. Gupta K.K., Meek J.L. A brief history of the beginning of the finite element method. International Journal for Numerical Methods in Engineering. 1996;39(22):3761-3774. https://doi.org/10.1002/(SICI)1097-0207(19961130)39:22<3761:: AID-NME22>3.0.C0;2-5
3. Zenkevich O., Morgan K. Konechnye elementy i approksimaciya [Finite elements and approximation]. Moscow: Mir Publ.; 1986. (In Russ.)
4. Tyukalov Yu.Ya. Finite element models in stresses for bending plates. Magazine of Civil Engineering. 2018;6(82): 170-190. https://doi.org/10.18720/MCE.82.16
5. Tyukalov Yu.Ya. Finite element models in stresses for plane elasticity problems. Magazine of Civil Engineering. 2018;1(77):23-37. https://doi.org/10.18720/MCE.773
6. Cen S., Li C., Rajendran S., Hu Z. Advances in finite element method. Mathematical Problems in Engineering. 2014:206369. https://doi.org/10.1155/2014/206369
7. Bushnell D., Almroth B.O., Brogan F. Finite-difference energy method for nonlinear shell analysis. Computers and Structures. 1971;1(3):361-387. https://doi.org/10.1016/0045-7949(71)90020-4
8. Barve V.D., Dey S.S. Isoparametric finite difference energy method for plate bending problems. Computers and Structures. 1983;17(3):459-465. https://doi.org/10.1016/0045-7949(83)90137-2
9. Maksimyuk V.A., Storozhuk E.A., Chernyshenko I.S. Variational finite-difference methods in linear and nonlinear problems of the deformation of metallic and composite shells. International Applied Mechanics. 2012;48(6):613-687. https://doi.org/10.1007/s10778-012-0544-8
10. Trushin S., Goryachkin D. Numerical evaluation of stress-strain state of bending plates based on various models. Procedia Engineering. 2016;153:781-784. https://doi.org/10.1016Zj.proeng.2016.08.242
11. Ivanov V.N., Kushnarenko I. Stiffeners in variational-difference method for calculating shells with complex geometry. Vestnik MGSU. Proceedings of Moscow State University of Civil Engineering. 2014;(5):25-34. (In Russ.)
12. Ivanov V., Rynkovskaya M. Analysis of thin walled wavy shell of monge type surface with parabola and sinusoid curves by variational-difference method. MATEC Web of Conferences. 2017;95:1-5. https://doi.org/10.1051/matecconf/20179512007
13. Govind P.L. Complicated features and their solution in analysis of thin shell and plate structures. Structural Mechanics of Engineering Constructions and Buildings. 2018;14(6):509-515. https://doi.org/10.22363/1815-5235-2018-14-6-509-515
14. Dzhavadyan A.D. Grid selection in the variation-difference method for solving second-order elliptic equations with quasidegenerate quadratic form. USSR Computational Mathematics and Mathematical Physics. 1989;29(6):22-33. https://doi.org/doi:10.1016/s0041-5553(89)80004-7
15. Ivanov V.N. Osnovy metoda konechnyh elementov i variacionno-raznostnogo metoda [Fundamentals of the finite element method and the variational-difference method]. Moscow: RUDN Publ.; 2008. (In Russ.)
16. Ivanov V.N. The variational-difference method and the method of global elements in the calculation of interfaces of shell compartments. Structural Mechanics of Engineering Constructions and Buildings. 2003;12:34-41. (In Russ.)
17. Mikhlin S.G. Variational-difference approximation. Journal of Soviet Mathematics. 1978;10(5):661-787. https://doi.org/https://doi.org/10.1007/BF01083968
18. Courant R. Variational methods for the solution of problems of equilibrium and vibrations. Bulletin of the American Mathematical Society. 1943;49(1):1-23.
19. Zhong H., Yu T. A weak form quadrature element method for plane elasticity problems. Applied Mathematical Modelling. 2009;33(10):3801-3814. https://doi.org/10.1016/j.apm.2008.12.007
20. Griffin D.S., Varga R.S. Numerical solution of plane elasticity problems. Journal of the Society for Industrial and Applied Mathematics. 1963;11(4):1046-1062.
21. Brush D.O., Almroth B.O. Buckling of bars, plates, and shells. New York: McGraw-Hill; 1975.
22. Xing Y., Liu B., Liu G. A differential quadrature finite element method. International Journal of Applied Mechanics. 2010;2(1):207-227. https://doi.org/10.1142/S1758825110000470
23. Repin S.I. A variational-difference method of solving problems with functionals of linear growth. USSR Computational Mathematics and Mathematical Physics. 1989;29(3):35-46. https://doi.org/10.1016/0041-5553(89)90145-6
24. Ivanov V.N., Krivoshapko S.N. Analiticheskie metody rascheta obolochek nekanonicheskoj formy [Analytical methods for calculating shells of non-canonical form]. Moscow; 2010. (In Russ.)
25. Krivoshapko S.N., Ivanov V.N. Encyclopedia of analytical surfaces. Switzerland: Springer International Publishing AG; 2015.
26. Ivanov V.N., Romanova V.A. Konstruktsionnye formy prostranstvennykh konstruktsii. Vizualizatsiya poverkh-nostei v sistemakh MathCad, AutoCad [Constructive forms of space constructions. visualization of the surfaces at the systems "MathCAD"and "AutoCAD"]. Moscow: ASV Publishing House; 2016. (In Russ.)
27. Krivoshapko S.N. Geometriya linejchatyh poverhnostej s rebrom vozvrata i linejnaya teoriya rascheta torsovyh obolochek [Geometry of ruled surfaces with cuspidal edge and linear theory of analysis of torse shells]. Moscow; 2009. (In Russ.)
28. Krivoshapko S.N. The application, geometrical and strength researches of torse shells: the review of works published after 2008. Structural Mechanics and Analysis of Constructions. 2018;2:19-25.
29. Krivoshapko S.N. Perspectives and advantages of tangential developable surfaces in modeling machine-building and building designs. Vestnik Grazhdanskix Inzhenerov [Proceedings of Civil Engineers]. 2019;16(1):20-30. (In Russ.) https://doi.org/10.23968/1999-5571-2019-16-1-20-30
30. Aleshina O.O. New information about the use of shells with tangential developable middle surfaces. Process Management and Scientific Developments. Birmingham: Infinity; 2020. p. 140-146.
31. Chen M., Tang K. A fully geometric approach for developable cloth deformation simulation. Visual Computer. 2010;26(6-8):853-863. https://doi.org/10.1007/s00371-010-0467-5
32. Ivanov V.N., Alyoshina O.O. Comparative Analysis of the stress-strain state's parameters of equal slope shell with the director ellipse using three calculation methods. Structural Mechanics and Analysis of Constructions, 2020;3(290):37-46. (In Russ.) https://doi.org/10.37538/0039-2383.20203.37.46
33. Aleshina O.O., Ivanov V.N., Grinko E.A. Investigation of the equal slope shell stress state by analytical and two numerical methods. Structural Mechanics and Analysis of Constructions. 2020;6:2-13. https://doi.org/10.37538/0039-2383.2020.6.2.13
34. Ivanov V.N., Alyoshina O.O. Comparative analysis of the results of determining the parameters of the stressstrain state of equal slope shell. Structural Mechanics of Engineering Constructions and Buildings. 2019;15(5):374-383. http://dx.doi.org/10.22363/1815-5235-2019-15-5-374-383 (In Russ.)
35. Aleshina O.O. Studies of geometry and calculation of torso shells of an equal slope. Structural Mechanics and Analysis of Constructions. 2019;3:63-70. (In Russ.)
36. Alyoshina O.O. Definition of the law of setting closed curves torso shells of the equal slope. Proceedings of the scientific and practical conference with international participation "Engineering Systems - 2020", dedicated to the 60th anniversary of the Peoples' Friendship University of Russia, Moscow, October 14-16, 2020. 2020;1:22-30. (In Russ.)
37. Kumudini Jayavardena M.K. Geometry and example of strength analysis of thin elastic shell in the form of a torse-helicoid. Questions of the strength of spatial systems: materials of the XXVIII Scientific Conference of the Engineering Faculty. Moscow: RUDN Publ.; 1992. p. 48-51. (In Russ.)
38. Krivoshapko S.N., Krutov A.B. Cuspidal edges, lines of the unit and self-intersectionsof some technological surfaces of slope. Journal of Engineering Researches. 2001;1:98-104. (In Russ.)
39. Ivanov V.N., Lamichane G.P. Compound space constructions. Proceedings of the scientific and practical conference with international participation "Engineering Systems - 2020", dedicated to the 60th anniversary of the Peoples' Friendship University of Russia, Moscow, October 14-16, 2020. 2020;1:31-39. (In Russ.)
40. Krivoshapko S.N. The opportunities of umbrella-type shells. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(4):271-278. http://dx.doi.org/10.22363/1815-5235-2020-16-4-271-278
41. Krivoshapko S.N. Analytical ruled surfaces and their complete classification. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(2):131-138. http://dx.doi.org/10.22363/1815-5235-2020-16-2-131-138 (In Russ.)
