Magazine of Civil Engineering. 2023. 119(3). Article No. 11903
2 ¡1 ^ Magazine of Civil Engineering issn
PMViV 2712-8172
5 r h ni.
ok ©33 1»' journal homepage: http://engstroy.spbstu.ru/
Research article UDC 624.35:531.391.3 DOI: 10.34910/MCE.119.3
Hexagonal rod pyramid: deformations and natural oscillation
frequency
M.N. Kirsanov B
National Research University "Moscow Power Engineering Institute", Moscow, Russia Mc216@ya.ru
Keywords: spatial truss, hexagonal dome, natural vibrations, lower frequency bound, upper frequency bound, Dunkerley method, Rayleigh method, Maple, induction
Abstract. A new scheme of a statically determinate dome truss is proposed. The purpose of the study is to obtain exact formulas for structural deflections under a uniform load and to find upper and lower analytical estimates of the first frequency of natural oscillations depending on the number of panels, sizes, and masses concentrated in the truss nodes. Calculation of forces in the truss rods is performed by cutting nodes. The system of equations in projections on the coordinate axes, compiled in the Maple software, includes the forces in the rods and the reactions of vertical supports located along two contours of the structure at the base. The amount of deflection and stiffness of the entire truss is calculated using the Mohr integral. To determine the lower estimate of the first frequency an approximate Dunkerley method is used. The formula for the upper limit of the first frequency is derived by the Rayleigh energy method. In the Rayleigh method, the shape of the deflection from the action of a uniformly distributed load is taken as the deflection of the truss. Displacements of loads are assumed to be only vertical. The overall dependence of the solution on the number of panels is obtained by induction on a series of solutions for trusses with a successively increasing number of panels. The operators of the Maple system of symbolic mathematics are used. Based on the calculation results, it was concluded that the distribution of forces over the structure rods does not depend on the number of panels. Asymptotes were found on the graphs of the obtained analytical dependences of the deflection on the number of panels for different truss heights. The estimates of the first natural frequency are compared with the numerical solution obtained from the analysis of the natural frequency spectrum. The coefficients of the frequency equation are found using the eigenvalue search operators in the Maple system. It is shown that the lower analytical estimate based on the calculation of partial frequencies differs from the numerical solution by no more than 54 %, and the upper estimate by the Rayleigh method has an error of about 2 %. The formula for the lower Dunkerley frequency estimate is simpler than the Rayleigh estimate.
Funding: This work was financially supported by the Russian Science Foundation 22-21-00473.
Citation: Kirsanov, M.N. Hexagonal rod pyramid: deformations and natural oscillation frequency. Magazine of Civil Engineering. 2023. 119(3). Article no. 11903. DOI: 10.34910/MCE.119.3
1. Introduction
The study of the stress state, deformations, and stability of spatial multifaceted dome structures are of both practical and theoretical interest [1, 2]. In the calculations of building structures of this type, as a rule, the finite element method is used [3, 4]. Regular (with periodic structure elements) statically determinate truss schemes are quite rare. The search for such schemes was even called "hunting" by R.G. Hutchinson and N.A. Fleck [3, 4]. One of the advantages of regular schemes [5] is that for them it is possible to analytically derive the dependences of the characteristics of the stress-strain state on their order (the number of periodicity elements, for example, the number of panels). To search for analytical solutions in the form of compact formulas, operators of symbolic mathematics, such as Maple, Mathematica [6], and
© Kirsanov, M.N., 2023. Published by Peter the Great St. Petersburg Polytechnic University.
others, are applicable. To obtain analytical solutions in the form of closed formulas, when modeling structures, it is necessary to make some simplifications. The construction must be statically determined. If at the same time it is regular, then for such a design it is possible to obtain calculation formulas for an arbitrary number of repeating elements. Regular trusses are, for example, planar or spatial trusses with identical panels or groups of panels. In this case, the analytical solution has a great advantage over the numerical one, not only due to the saving of computation time but also due to the fundamental possibility to calculate the truss with a very large number of panels without a loss of accuracy. A larger number of elements causes the inevitable effect of accumulation of rounding errors in the numerical calculation. Analytical solutions are especially effective for preliminary draft calculations, for assessing the accuracy of numerical solutions, and in truss optimization problems.
Analytical solutions for building structures do not always lead to the final compact design formula. In such works, an algorithm for calculations in the system of symbolic mathematics is given [7, 8]. The purpose of this work is to obtain formulas for calculating deflections and estimates of the first natural frequency of a three-dimensional truss. Most often, to solve such problems for planar [9-13], and three-dimensional [14] trusses, the induction method is used. Solutions of problems on the deformation of some planar arch trusses [15] with an arbitrary number of panels are obtained inductively and some problems on natural frequencies of regular structures are solved [16-19]. Calculation formulas for deflections and oscillation frequencies of spatial regular trusses were obtained in [20].
In the analytical form, it is impossible to obtain a solution directly from the analysis of the entire frequency spectrum in the general case.
Therefore, to estimate its lower limit of the first (lowest) frequency of natural oscillations, we will use the Dunkerley method [21-23], and for the upper one, the Rayleigh method [24, 25]. These approximate methods are based on the calculation of partial frequencies, for which it is not necessary to compose highorder characteristic polynomials (in terms of the number of degrees of freedom).
Pyramidal trusses are studied in connection with the design of structural panels (composite) in which the trusses act as a kind of reinforcement. Hexagonal trusses are also used in mesh coatings [26, 27]. Eccentrically braced frames and beam-type spatial trusses were studied numerically and experimentally in [28]. An overview of analytical solutions for planar statically determinate regular trusses is given in [29].
In this paper, we propose a new scheme of a statically determinate dome-type hexagonal spatial truss. The truss has architectural expressiveness and can be used in public buildings (circus, airport building, railway station, etc.).
The construction can be used as a basis for complicated statically indeterminate systems of this type.
2.. Methods
2.1. The truss scheme
The truss in the form of a regular pyramid 2hi + ^ high with a hexagonal base side na contains
ns = 36n -15 rods, including 6n vertical support posts h1 high located along the outer contour of the
structure and 6(n-2) posts 2h1 high supporting the upper contour (Fig. 1, 2). There are no posts in the
corner nodes D of the upper contour. A similar hexagonal cover, but with a dome attached to the corner points of the outer (upper) contour, is considered in [30].
Figure 1. Truss scheme, n = 3, vertical load.
C
Figure 2. The truss dimensions, n = 3.
The lower horizontal rod contour consists of 6n rods of length a, the upper one consists of
6(n -1) similar rods. The braces connecting the contours have a length of c = \Ja2 + h2. The lengths
of six identical braces emanating from vertex C depend on the number of panels: (n -1) c. Corner node
A rests on a spherical support hinge, modeled by three mutually perpendicular rods, one of which is a vertical post. Node B is a cylindrical hinge corresponding to two support bars. The following ratios of sizes are chosen: h = h, h2 =(n-1)h. All connections of the truss rods are hinged. An analytical calculation
of the deformations of a spatial coating with a similar structure was performed in [20]. In [31] an optimization problem is solved for an irregular spatial truss of 25 rods.
Calculation of forces in rods in symbolic form is performed based on a program written in the language of computer mathematics Maple [32]. To do this, we introduce the coordinates of the nodes (Fig. 3), using the annular periodic structure of the truss:
16 15 15 14 M 13 13
Figure 3. Numbers of nodes and rods of contours n = 3.
The coordinates of the nodes of the lower contour look like this:
xi+jn = L cos ^ - a(i -1) cos P, yi+jn =L sin ^ + a(i -1) sin P,
zi+jn = ai=n j=0,..A
where L = na, ^ = jn/3, P = n/3
The coordinates of the nodes of the smaller (upper) contour:
+j(n-1)+6n = (L -a) cos ^ - a (i - 1) cos P y+j(n-1)+6n = (L -a) sin ^ + a (i - 1) sin P, zi+j( n-1)+6 n = ^ i = 1,..., n - 1, j = 0,...A Vertex C coordinates^ x12n-5 = y12n-5 = ^ z12n-5 = h1 + h2.
The order of connecting the bars of the lattice is entered into the program using ordered lists Of, i = 1,..,ns. of numbers of the ends of the corresponding bars, similar to how graphs are given in discrete mathematics. For example, the bars of the lower chord are encoded with the following vertex lists
Oi =[i, i +1], i = 1,..,6n-1,O3n =[6n,l].
Upper chord bar code:
0i+6n =[i + 6n, i + 6n +1], i = 1, ..., 6n-7, 012n_6 =[12n-6, 6n +1].
The numbers of ends and other rods are set similarly.
2.2. Calculation of forces in bars
Let us represent the system of equilibrium equations of nodes in the projection on the coordinate axes in matrix form GS = T, where S is the vector of unknown forces, including the reactions of the
supports, G is the matrix of coefficients (projections of unit forces in the rods), T is the vector of loads on the nodes. For each node in the matrix, three rows are assigned, corresponding to projections onto three coordinate axes. Similarly, in the elements of the load vector of the form Y3^2, where i is the number
of the node, the loads on node i in the projection on the x-axis are written. Elements ¥3^1 contain projections of external forces in projection onto the y -axis. Vertical loads on nodes are recorded in elements ^3i.
Matrix G elements are calculated according to the data on the structure of the connection of the bars and the coordinates of the nodes
gxi = ( x0„, - x0„ )/'< • gy' =(y0„ - y0,,2 )/'< •
/2 2 2
where = Jlx i + lyt + lz i is the length of the rod i. The number of rods also includes three horizontal
support rods at angles A and B. The matrix of coefficients of equilibrium equations in projections is filled in rows. Every three lines correspond to the projection equations on the x, y, and z axes, respectively:
G3<0 ! - 2, i = gxi, G30i,1 -1,i = gyi, G30i,1,i = gzi, G30i 2-2,i = -gxi, G30i,2-1,i = -gyi, G30i,2,i = -gzi.
If a uniform vertical load is applied to the truss nodes (Fig. 1), then the non-zero elements of the load vector have the form: ^3i = P, i = 1,..., 12n - 5. Numerical calculation of forces for a structure with n = 3, a = 5.0m, h = 1.0m gives a picture of the distribution of forces shown in Fig. 4. Compressed rods are highlighted in blue, tension rods are highlighted in red. The thickness of the line is greater, the greater the modulus of force in the corresponding rod. The force value is related to the value of the nodal load P with an accuracy of two significant digits.
Figure 4. Distribution of forces in the rods, n = 3, a = 5.0 m, h = 1.0 m.
The upper contour of the truss under such a load is compressed, the lower one is stretched. The braces connecting the contours have zero forces. An interesting feature of the stressed state of the truss
was noticed. The patterns of distribution of forces in the rods from a uniform load for trusses of different orders are similar. The order of the truss is equal to the number of rods in the side edge of the lower contour. The compressive forces in the rods of the upper contour for any n are equal to
Si = -aP/h, i = 6n + 1,...,12n-6. (1)
The greatest tensile forces are observed in the lower belt
Si = 7aP/(6h), i = 1,...,6n. (2)
The six lower corner ribs are most compressed
Si = -7Pc/(6h), i = 12n-5,...,17n-5. (3)
The forces in the six upper rods of the dome, connected at the top C, are equal to - PcI ( 6h ).
The forces in the six corner support posts do not depend on the dimensions of the structure: -13P/ 6. The reactions of the supports of the intermediate posts along the outer (lower) and inner (upper) contours are equal to P.
2.3. Deflection
The formula for the dependence of the deflection of the top C on the dimensions of the structure, the load, and the number of panels will be obtained by induction. Under the deflection An we mean the
vertical displacement of the node C of the truss of order n. To calculate the deflection value, we use the Mohr integral
n±SJsJlJ
a n = I-E^' (4)
j=1 EF
where lj is the length of the rod, Sj is the force in theyth rod from the action of the load, Sj is the force
in the rod from the action of a single vertical force applied to the vertex C, E is the modulus of elasticity
of the rods, F is the cross-sectional area. The summation is carried out over all the bars of the structure. The elastic moduli and cross-sectional areas are the same for all rods. Sequential calculation of the deflection of a series of trusses with an increasing number of panels gives the following results
A2 = P (14a3 + 8c3 + 13h3 )/(6h2EF),
A 3 = P ( 21a 3 + 9c 3 + 13h3 )/( 6h2 EF ),
A 4 = P ( 28a 3 + 10c3 + 13h3 ( 6h 2 EF ),
A5 = P (35a3 + 11c 3 + 13h3 )/( 6h2 EF ),
A 6 = P (42a3 + 12c3 + 13h3 )/( 6h2 EF ).
In the general case, we have the form of a formula for the deflection:
A n = P (C1a3 + C2c3 + C3h3)/(h 2EF). (5)
The coefficients in this expression are functions of the number of panels n. The common members of the sequences they form can be found using the special operators rsolve and rgf_findrecur from the Maple system. Equally effective in finding common members of sequences are the operators of the Mathematica computer mathematics system. The common terms of the sequences of coefficients at n = 2,3,...,7 are linear concerning n
C1 = 7n/6, C2 =( n + 6)/6, C3 = 13/3. (6)
Compared to similar well-known solutions for planar trusses, which have a shape that is non-linear in terms of the number of panels, the solution turned out to be much more compact. In part, this can be explained by the observed feature of the stress state of the structure, which does not depend on the number of panels. The same simple solution is obtained in the problem of the deflection of node D at the corner (not supported) point of the upper contour. The coefficients in (5) in this case have the form
C1 = (13n - 6)/27, C2 = 7/6, C3 = 13/6. (7)
Let us plot the solution graphs (5), (6). Let us denote the total load on the truss psum = P (12 n - 5)
and the length of the outer side of the cover L = na. Let us introduce the designation for the
dimensionless deflection: A'= A nEF / (PsumL). The graphs of the curves of solutions (6) and (7) at
L = 50 m show that in this setting the dimensionless deflection at points C and D decreases monotonically with an increase in the number of panels (Fig. 5).
\
\
\ \
\
\ V \ \
- \ \
\ \ \ \
\ \
\ \
s h=0.5n
K.
=lru"
Figure 5. Dependence of dimensionless deflection of top C and node D on the number of panels.
For small n, the deflection of the non-supported node D is one and a half times greater than the
deflection of the vertex C. Horizontal asymptotes of the solutions (ultimate deflection) are noted. Using the analytical form of the solution, using the operators of the Maple system, we obtain the lower limits of the relative deflections: lim AC = h/(72L), lim A'D = 0. It follows that curves A'D and AC must
n n
intersect. Calculations show that the intersection of the curves for the height h = 1.0m occurs at n = 76 and for h = 0.5m at n = 151.
In practice, there are problems about the deflection of a structure under the action of a load on only part of its surface. This corresponds, for example, to a snow load applied to one half of the roof (Fig. 6). The derivation of the formula in this case is no different from the previous tasks. The form of the solution obtained by induction on eight girders with a successively increasing number of panels on one edge coincides with solution (5). The coefficients look like:
q = 5n/ 6, C2 =( 4 + n )/6, C3 = 3/2.
Figure 6. Vertical load on half of the truss surface, n = 4.
2.4. Natural oscillation frequency
The calculation of the first (lowest) frequency of natural oscillations is included in most dynamic calculations of the structure and is of independent interest. In addition to the direct calculation of the spectrum of natural frequencies, which is performed numerically [33-35], there are also known approximate methods for obtaining its upper and lower estimates of the first natural oscillation frequency [21-25]. These methods are based on the calculation of partial frequencies, the values of which can be found analytically. For regular constructions, analytical estimates can be generalized to an arbitrary number of panels using the [16] induction method.
The inertial properties of the truss are modeled by concentrated masses in the nodes. In the simplest setting, the masses of loads m are the same. Only vertical vibrations of nodes are considered. The number of degrees of freedom of the truss weight system of order n is equal to the number of nodes K = 12n - 5.
The dynamics of the system is described by a system of differential equations for the movement of goods in matrix form:
M k Z
D K Z = 0,
(8)
where Dk is the structural stiffness matrix, Z is the vector of vertical displacements of masses 1,..., K, Mk is the inertia matrix of size K x K, Z is the acceleration vector. The inertia matrix is proportional to the identity matrix Mk = mlK if the masses are the same. The elements of the compliance matrix Bk, which is the inverse of the stiffness matrix Dk , can be found using the Mohr integral
h j =isH ) S/ \/( EF ),
a=1 /
(9)
where
sO' )
is the force in the rod a from the action of a unit vertical force at node i. The problem can
be reduced to the problem of matrix eigenvalues Bk. To do this, we multiply (8) from the left by B
'K
and, taking into account the replacement Z = -q Z, which is valid for harmonic oscillations of the form
zt =ui sin (ot + ^o )
we obtain BKZ = AZ, where A = lj(mo2 ) is the eigenvalue of the matrix B frequency of oscillations.
From here we obtain the calculation formula for the frequency q = a//(mA).
K>
(10)
q is the natural
The forces S„ ) in the truss rods included in the elements of the matrix BK are determined by solving the system of equations of the truss nodes, which also includes the reactions of the supports. Consider approximate methods that give upper and lower estimates of the first frequency.
2.5. Energy method. Top rating
The Rayleigh formula, which follows from the equality of the maximum values of the kinetic and potential energies, has the form:
T = n (11)
The kinetic energy of a system consisting of K identical masses m located at the nodes of the
K 2 /
structure has the form: T = X mvi /2.
i=1 /
The vertical velocity vi of the mass i according to (7) has the form: vt = zi = raui sin (rat + 90 ). The maximum kinetic energy corresponds to equality max (sin (rat + 9 )) = 1. From here we get:
Tmax =ra2m Xu2 /2, (12)
i=1 /
where the amplitude of the vertical displacement is calculated using the Mohr integral:
ut = X Sf)S£\ / (EF) =P X sa)S«/a / (EF) = Put .
a=1 a=1
The previous designations are used: S^) is force in the rod a = 1,...,ns from the action of the load P, uniformly distributed over the nodes, S^) is force in the same rod from a single (dimensionless)
load applied to the mass with number i, S^) = S^)/P. The choice of a uniformly distributed load is
determined by the proximity of the shape corresponding deflection to the form of vibrations of the system of weights with the first frequency. Thus, (12) takes the form:
Tmax = P 2ra2 Xmu2/2, (13)
i=1
ns /
where ui = ui/P = X S^S$/a (EF) is amplitude of displacements of the mass with number i
a=1 /
under the action of a distributed load, referred to the value of P. The potential energy of deformation of the truss rods has the form:
n max
: X SiPWa 2 = XX (SiP))2 /a/(2EF).
a=1 / a=1 'I
N
Due to the linearity of the problem concerning loads, we have S^ = P X S^). From here we get:
i=1
nmax = P2 XSf )XSt')/a/(2EF) =
a=1 i=1 /
K ns ( ) ( ) / N /
= P2 XX SiP)Sahl{2EF) = P2 X u~ /2.
i =1 a / i=1 /
From (11), (13), (14) we obtain a formula for the upper estimate of the first oscillation frequency of the truss (the Rayleigh formula):
k ik 2 = X ui X mUi. (15)
i=1 / i=1
We find the displacements ui as functions of n. We generalize the solution obtained for a different
K K
number of panels concerning n. Consider the sums X u^ and X U separately.
i=1 i=1
The calculation of displacement for trusses with the different number of panels shows that the form
K
of the solution X ui does not depend on n. Note that irregular trusses do not have this property. The
i=1
numerator in (15) can be represented as:
f u = ( Caa3 + Cce3 + Chh3 )/( h2 EF ), i=1
or in a more compact form
k
f Ui = f mCaa3/(h2EF), (16)
i=1 a=[a, c,h] /
where the coefficients Ca, Cc, Ch are obtained by induction, generalizing a series of solutions for different n :
k
n = 2, f ut = (l34a3 + 50c3 + 205h3 6h2 EF ), i=1
K
n = 3, f Ui =(219a3 + 51c3 + 313h3 6h2EF), i =1
k
n = 4, f U =(304a3 + 52c3 + 421h36h2EF), i =1
k
n = 5, f ut = (389a3 + 53c3 + 529h3 )/( 6h2 EF ), i =1
k
n = 6, f U =(474a3 + 54c3 + 637h3 6h 2EF ),... i =1
As a result, we have the coefficients
Ca =(85n-36), Cc =(n + 48^6, Ch =(108n-11)/6. (17)
The denominator (15) has a more complex form:
fTmU2k = f mCapa3p3 /(h4E2F2 ), (18)
k=1 a,P=[a, c, h] /
where
Caa = (1063n2 -936n + 216)/36, Ccc =( n2 + 12n + 330)136,
(19)
Chh =(1080n + 253)/ 36,
Cac = ( 7n2 + 588n - 252)/36,
Cah =(1105n - 468V 36, Cch =(13n + 624)/36.
Thus, the upper estimate of the first frequency of the truss, depending on the number of panels, can be obtained by the formula:
= h
EF f
a=[a,c,h]
Caa
3
303
m f Capa3P
a,p=[a,c,h]
(20)
with coefficients (17), (19) depending only on the construction order n.
2.6. Dunkerley score We obtain a lower estimate of the first frequency of oscillations using the Dunkerley formula:
-2 K -9
®D =f®-2, (21)
i=1
where is the oscillation frequency of one mass m located at node i. To calculate the partial frequencies , we compose equation (8) in the scalar form:
mzi + Dizi = 0
where Di is the scalar stiffness coefficient (i is the mass number). The frequency of vibrations of the
load is rai =^Di/m. The stiffness coefficient, the reciprocal of the compliance coefficient, is determined by the Mohr integral (4):
8i = VDi = f (sa))2 la (eEF).
Arguing in the same way as when calculating the frequencies of a system with many degrees of freedom, we obtain
2 K K 1 K
®D =f®i = m f — = m fSf = i=1 i=1Di i=1
= m f f (Sia))2 la {EF ) = mf n /(h2EF ).
i=1 a=1 /
2 K u2
Let us successively calculate the sums f n = h f f ( Sa M la for n = 2, 3, 4,
i=1 a=1 '
(22)
_ 94a3 + 154c3 + 113k3 12 _ 12
_ 219a3 + 231c3 + 601k3 13 _ 18 1124a3 + 5
14 _ 48
_ 4115a3 + 1139c3 + 10555k3 15 _ 150
1124a3 + 584c3 + 2495k3
Let us calculate the common terms of the sequences of coefficients in these expressions.
: I >
a_[a,c, h]
3
We get In _ I raa , where
ra _(49n2 - 60n +18)I( 6n), rc _(n3 + 36n2 + 186n - 216)/(6n2 ), (23)
(108n3 -101 n2 - 60n + 30y(6n2 ).
rh = (108« -107« - 60n + 30 ]
When deriving expressions for the coefficients (23), the operators for compiling and solving the recursive equations of the Maple system were used. Some complication was the dependence on n not only of the numerators of the sequence members but also of their denominators. Maple system operators are adapted to define common members of such sequences. The result in the form of expressions (23) was obtained only because it was possible to guess the form of the denominators.
Thus, we obtain the lower estimate for the first frequency according to Dunkerley:
raD = h
EF
(24)
m X raa3
a=[ a, c, h]
The form of the Dunkerley estimate (24) almost coincides with formula (20) obtained by the Rayleigh method, but formula (24) is much simpler. The desired coefficients are contained here only in the denominator.
3. Results and Discussion
To estimate the error of the estimates found, consider an example of a truss with n panels for a = 6.0 m and h = 1.0 m . The stiffness of the steel rods of the truss will be taken EF = 1.8-105 kN. Fig. 7 plots the dependences of the upper estimate of the first frequency Qr of natural oscillations of the
truss (20), obtained by the Rayleigh energy method, ®d the Dunkerley estimate (24), and the numerical solution found as the minimum frequency of the entire frequency spectrum.
The numerical value of the lowest natural frequency of a system with K = 12« - 5 degrees of freedom is so close to the Rayleigh estimate that for n > 3 the curves merge. To refine the error estimates, we introduce the relative values Sd = |®d , SR = |®R .
CO,S1 20
15
10
\ \
\A
(!)
s -.1
--- ---- ---- —--
n
23456789 10
Figure. 7. The first oscillation frequency obtained in three ways.
Depending on the number of panels, the error of the Dunkerley solution is relatively large (Fig. 8),
but it changes little, increasing sharply only at n = 3 The Rayleigh estimate error (Fig. 9) is very small and, in principle, it can be neglected, considering the analytical solution (20) to be exact.
0.52
0.51
/ \ \
i i \ \
1 i \
I j \ \ V
i / 1 / ■s t h-2 r n
ll ¡1 1= In
i
0.10
0.06
0.04
10
\ i i
i i i
\ i
\ i
\ i \ i \ i
\ i \ i \ i h =2 m
______ -----
1 i=l ni
n
10
Figure 8. Dunkerley's estimation. Figure 9. Rayleigh estimation.
With an increase in the number of panels up to a certain number n, the Dunkerley estimate error slightly decreases, and then, starting from n = 12 at h = 2m , it grows a little.
The error in the Rayleigh estimate increases slightly while remaining very small. The error of both methods increases with the height of the truss.
3.1. Discussion
A new scheme of a statically determined truss of a three-dimensional hexagonal cover is proposed. The truss can only partly be considered regular. It has only a regular base with two hexagonal horizontal contours, while the upper part has the shape of a regular hexagonal pyramid of rods connected at the vertex C. Despite this, the inductive method using a computer mathematics system made it possible to obtain analytical exact solutions for the deflections of its characteristic peaks and a two-sided estimate of the first frequency.
The design under consideration can be used in coverings of public buildings and structures of a large span, for example, arenas, buildings of stations and airports, circuses.
The external static indeterminacy, which means that the reactions of the supports can only be found from the solution of the joint system of equilibrium equations for all nodes simultaneously with the forces in the rods, did not complicate the solution. On the contrary, some obtained formulas turned out to be linear
in the number of panels, which distinguishes this solution from similar formulas even for planar trusses. This is largely due to the property of the truss itself. An important property of similarity of the stressed states of trusses of various orders is noticed. The calculation showed that the pattern of force distribution over the structure rods does not depend on the number of panels. When n changes, the forces in the chords do not change, the formulas for the forces (1) - (3) and the reactions of the supports do not contain the number of panels n.
Comparison with the numerical solution of the complete problem of oscillations of a mass system with many degrees of freedom confirmed the well-known fact [16] that the Rayleigh formula for the upper estimate gives much greater accuracy than the Dunkerley method for the lower estimate of the first frequency.
4. Conclusion
Main results of the work:
1. A scheme of a spatial statically determinate rod structure of the dome type has been developed. Analytical solutions for deflection are obtained for both symmetrical and asymmetric loads. The truss design does not contain a central support post and can be used in large-span structures without central supports.
2. The calculation showed that the pattern of force distribution over the structure rods does not depend on the number of panels.
3. Formulas are derived both for the deflections of the truss and the boundaries of its first frequency of natural oscillations for an arbitrary number of panels. The obtained estimate of the first frequency from above has very high accuracy. Formulas can be used to evaluate numerical solutions for a very large number of panels, that is, precisely in those cases where the accumulation of calculation errors in numerical form is most likely.
4. The closed analytical form of the obtained formulas allows the use of mathematical analysis tools to identify their features and to search for combinations of design parameters that are optimal in terms of strength, rigidity, or stability.
References
1. Kobielak, S., Zamiar, Z. Oval concrete domes. Archives of Civil and Mechanical Engineering. 2017. 17 (3). Pp. 486-501. DOI: 10.1016/J.ACME.2016.11.009
2. Rezaiee-Pajand, M., Rajabzadeh-Safaei, N. Exact post-buckling analysis of planar and space trusses. Engineering Structures. 2020. 223. Pp. 111146. DOI: 10.1016/J.ENGSTRUCT.2020.111146
3. Hutchinson, R.G., Fleck, N.A. Microarchitectured cellular solids - The hunt for statically determinate periodic trusses. ZAMM Zeitschrift fur Angewandte Mathematik und Mechanik. 2005. 85 (9). Pp. 607-617. DOI: 10.1002/zamm.200410208
4. Hutchinson, R.G., Fleck, N.A. The structural performance of the periodic truss. Journal of the Mechanics and Physics of Solids. 2006. 54 (4). Pp. 756-782. DOI: 10.1016/j.jmps.2005.10.008
5. Zok, F.W., Latture, R.M., Begley, M.R. Periodic truss structures. Journal of the Mechanics and Physics of Solids. 2016. 96. Pp. 184-203. DOI: 10.1016/j.jmps.2016.07.007.
6. Zotos, K. Performance comparison of Maple and Mathematica. Applied Mathematics and Computation. 2007. 188 (2). Pp. 14261429. DOI: 10.1016/j.amc.2006.11.008
7. Goloskokov, D.P., Matrosov, A.V. A Superposition Method in the Analysis of an Isotropic Rectangle. Applied Mathematical Sciences. 2016. 10 (54). DOI: 10.12988/ams.2016.67211. (date of application: 17.06.2020).
8. Goloskokov, D.P., Matrosov, A.V. Comparison of two analytical approaches to the analysis of grillages. 2015 International Conference on "Stability and Control Processes" in Memory of V.I. Zubov, SCP 2015 - Proceedings. 2015. Pp. 382-385. DOI: 10.1109/SCP.2015.7342169
9. Rakhmatulina, A.R., Smirnova, A.A. The dependence of the deflection of the arched truss loaded on the upper belt, on the number of panels. Science Almanace. 2017. 28 (2-3). Pp. 268-271. DOI: 10.17117/na.2017.02.03.268. (date of application: 9.05.2021).
10. Kitaev, S.S. Derivation of the formula for the deflection of a cantilevered truss with a rectangular diagonal grid in the computer mathematics system Maple. Postulat. 2018. 5-1. Pp. 43. URL: http://e-postulat.ru/index.php/Postulat/article/view/1477 (date of application: 3.03.2021).
11. Ilyushin, A. The formula for calculating the deflection of a compound externally statically indeterminate frame. Structural mechanics and structures. 2019. 3 (22). Pp. 29-38. URL: https://www.elibrary.ru/download/elibrary_41201106_54181191.pdf
12. Arutyunyan, V.B. Calculation of the deflection of a statically indeterminate beam truss. Postulat. 2018. 6(6). URL: http://vuz.exponenta.ru/1/ar18.pdf (date of application: 26.02.2021).
13. Dai, Q. Analytical Dependence of Planar Truss Deformations on the Number of Panels. AlfaBuild. 2021. 17. Pp. 1701. DOI: 10.34910/ALF.17.1
14. Kirsanov, M.N., Zaborskaya, N. Deformations of the periodic truss with diagonal lattice. Magazine of Civil Engineering. 2017. 71(3). DOI: 10.18720/MCE.71.7
15. Voropay, R., Domanov, E. Analytical solution of the problem of shifting a movable support of a truss of arch type in the Maple system. Postulat. 2019. 1. URL: http://vuz.exponenta.ru/1Zvd.pdf (date of application: 27.02.2021).
16. Vorobev, O.V. Bilateral Analytical Estimation of the First Frequency of a Plane Truss. Construction of Unique Buildings and Structures. 2020. 92 (7). Pp. 9204-9204. DOI: 10.18720/CUBS.92.4. URL: https://unistroy.spbstu.ru/article/2020.92.4 (date of application: 27.02.2021).
17. Petrenko, V.F. The natural frequency of a two-span truss. AlfaBuild. 2021. (20). Pp. 2001. DOI: 10.34910/ALF.20.1
18. Vorobyev, O. About methods of obtaining analytical solution for eigenfrequencies problem of trusses. Structural mechanics and structures. 2020. 1 (24). Pp. 25-38. URL: http://vuz.exponenta.ru/PDF/NAUKA/elibrary_42591122_21834695.pdf
19. Kirsanov, M., Safronov, V. Analytical estimation of the first natural frequency and analysis of a planar regular truss oscillation spectrum. Magazine of Civil Engineering. 2022. 111 (3). DOI: 10.34910/MCE.111.14
20. Kirsanov, M. Deformations And Spatial Structure Vibrations Frequency of The Rectangular Contour Type Cover: Analytical Solutions. AlfaBuild. 2021. 98. Pp. 9805. DOI: 10.4123/CUBS.98.5
21. Low, K.H. A modified Dunkerley formula for eigenfrequencies of beams carrying concentrated masses. International Journal of Mechanical Sciences. 2000. 42(7). Pp. 1287-1305. DOI: 10.1016/S0020-7403(99)00049-1
22. Trainor, P.G.S., Shah, A.H., Popplewell, N. Estimating the fundamental natural frequency of towers by Dunkerley's method. Journal of Sound and Vibration. 1986. 109(2). Pp. 285-292. DOI: 10.1016/S0022-460X(86)80009-8
23. Levy, C. An iterative technique based on the Dunkerley method for determining the natural frequencies of vibrating systems. Journal of Sound and Vibration. 1991. 150(1). Pp. 111-118. DOI: 10.1016/0022-460X(91)90405-9
24. Mochida, Y., Ilanko, S. On the Rayleigh-Ritz Method, Gorman's Superposition Method and the exact Dynamic Stiffness Method for vibration and stability analysis of continuous systems. Thin-Walled Structures. 2021. 161. Pp. 107470. DOI: 10.1016/j.tws.2021.107470
25. Low, K.H. Natural frequencies of a beam-mass system in transverse vibration: Rayleigh estimation versus eigenanalysis solutions. International Journal of Mechanical Sciences. 2003. 45 (6-7). Pp. 981-993. DOI: 10.1016/j.ijmecsci.2003.09.009
26. Clough, E.C., Ensberg, J., Eckel, Z.C., Ro, C.J., Schaedler, T.A. Mechanical performance of hollow tetrahedral truss cores. International Journal of Solids and Structures. 2016. 91. Pp. 115-126. DOI: 10.1016/J.IJSOLSTR.2016.04.006
27. Kang, K.J. Wire-woven cellular metals: The present and future. Progress in Materials Science. 2015. 69. Pp. 213-307. DOI: 10.1016/J.PMATSCI.2014.11.003
28. Haji, M., Azarhomayun, F., Azad, A.R.G. Numerical investigation of truss-shaped braces in eccentrically braced steel frames. Magazine of Civil Engineering. 2021. 102(2). DOI: 10.34910/MCE.102.8
29. Tinkov, D.V. The optimum geometry of the flat diagonal truss taking into account the linear creep. Magazine of Civil Engineering. 2016. 61(1). Pp. 25-32. DOI: 10.5862/MCE.61.3
30. Kirsanov, M. Model of a spatial dome cover. Deformations and oscillation frequency. Construction of Unique Buildings and Structures. 2022. 99 Article No 9904. DOI: 10.4123/CUBS.99.4
31. Bekdaç, G., Yucel, M., Nigdeli, S.M. Evaluation of metaheuristic-based methods for optimization of truss structures via various algorithms and lèvy flight modification. Buildings. 2021. 11(2). Pp. 1-25. DOI: 10.3390/BUILDINGS11020049
32. Buka-Vaivade, K., Kirsanov, M.N., Serdjuks, D.O. Calculation of deformations of a cantilever-frame planar truss model with an arbitrary number of panels. Vestnik MGSU. 2020. (4). Pp. 510-517. DOI: 10.22227/1997-0935.2020.4.510-517
33. Siriguleng, B., Zhang, W., Liu, T., Liu, Y.Z. Vibration modal experiments and modal interactions of a large space deployable antenna with carbon fiber material and ring-truss structure. Engineering Structures. 2020. 207. Pp. 109932. DOI: 10.1016/J.ENGSTRUCT.2019.109932
34. Liang, L., Li, X., Sun, Y., Gong, Z., Bi, R. Measurement research on vibro-acoustic characteristics of large-span plate-truss composite bridge in urban rail transit. Applied Acoustics. 2022. 187. Pp. 108518. DOI: 10.1016/J.APACOUST.2021.108518
35. Vo, N.H., Pham, T.M., Bi, K., Hao, H. Model for analytical investigation on meta-lattice truss for low-frequency spatial wave manipulation. Wave Motion. 2021. 103. Pp. 102735. DOI: 10.1016/J.WAVEMOTI.2021.102735
Information about author:
Mikhail Kirsanov, Doctor of Physics and Mathematics ORCID: https://orcid.org/0000-0002-8588-3871 E-mail: c216@ya.ru
Received 18.02.2022. Approved after reviewing 07.12.2022. Accepted 07.12.2022.