Architecture and Engineering
Volume 1 Issue 3
AN INDUCTIVE METHOD OF CALCULATION OF THE DEFLECTION OF THE TRUSS REGULAR TYPE
Mikhail Kirsanov
Moscow Power Engineering Institute (National Research University) 14 Krasnokazarmennya str., Moscow, 111250, Russia
Abstract
In this paper, the author proposes a regular option of statically determinate flat truss with a lattice of Shukhov's type working under beam scheme. By induction, using Maple computer algebra system, the exact formula for the mid-span deflection under the action of the uniform load is obtained. The deflection of the elastic structure is defined by the formula of Maxwell - Mohr. Forces in the rods are found using cut nodes method with the composition of the management of the matrix of equations in vector form. The analysis of the coefficients in the truss evaluation with a different number of panels revealed the recurrence equation satisfied by these coefficients, and then the determined operators of a specialized package genfunc give a solution of the resulting equation that determines the required general formula. The case of a uniform load on the top horizontal zone is considered. The author obtained a simple asymptotic evaluation of this solution, proving the cubic nature of the growth of deflection for a fixed span length and a given nodal load. The solution found for an even number of panels and rods of the same cross-section, has a polynomial form; the dependence is non-monotonic and detects extremum, which allows optimizing the design.
Keywords: truss, deflection, induction, exact solution, Maple.
Introduction
The number of diagrams of statically determinate trusses of the regular type, which used in mechanical engineering as elements of bearing systems, lifting mechanisms of lightweight truss manipulators is very limited. Based on them more complex (regular or irregular) statically indefinable constructions are created. Solutions for basic systems can be used as a part of the solution for the real system by the method of forces or just as a test, and debug for various kinds of numerical methods or software packages based on them. The problem of finding diagrams of statically determinate trusses is called "hunting", and this "hunting" is keep going (Hutchinson and Fleck, 2005). In this paper, author proposes a new scheme of the regular type, where the cell performs the regular Shukhov' truss. This type of truss doesn't allow the direct computation of forces in the rods (regardless of the current load). In this paper we consider regular truss (Figure 1). Denote n is the number of panels (cell periodicity in the form of Shukhov's truss).
Figure 1. Truss, k=1
The method of cutting out of knots in the form in which it is usually applied (calculation from node to node with a consistent definition of the forces of the "chain") are not suitable. There are no such nodes which are attached to only two rods with unknown forces and with which to be-
Mikhail Kirsanov — Pages 14-17 AN INDUCTIVE METHOD OF CALCULATION OF THE DEFLECTION OF THE TRUSS REGULAR TYPE
gin the calculation. The method of sections (Ritter's method) is also not suitable since the respective cross sections for n>1 does not exist. Only the calculation of the complete system of equilibrium equations of all components and definition forces across all terminals, including three rods, simulating the support is possible. This approach will be used below. Moreover, the aim is not just to obtain a formula for a specific truss, but to deduce the General formula for all trusses of this scheme, i.e. to obtain the solution for an arbitrary number of panels. The solution will be obtained by the methods of symbolic mathematics using the method of induction. Previously this method was used for deriving the formulas of deflection in a plane and spatial structures (Kirsanov, 2014, 2015a,b; Tinkov, 2015).
The geometry of the truss and the system of equations of equilibrium
We will create the system of equilibrium equations of all the nodes in the projections on the coordinate axes to compute the forces in the rods using cut nodes. The origin of coordinates we choose the left fixed support. Calculation program is compiled in the Maple language (Monagan et al., 2005). Assume that the truss contains an even number of panels: n = 2k . The number of the rod is ns = 8n + 4 , including three support rods. The number of nods along with three nods mounted on the base m = 4n + 5 . We will set coordinates of nodes on which it will be possible to determine the directing cosines of the forces included in the equation of equilibrium:
xi = 4a- ^ y = 0
xi+n+1 = xi, yi+n+1 = b, i = 1,-; n + 1,
Xi+2n+2 = a + 2a(i - 1)' yi+2n+2 = b + C 1 =
Coordinates of the supports:
xm-2 = ym-2 = 0,
Xm-1 = Xl, ym-l =
Xm = Xn+1, ym = — 1
The order of connection of nodes and cores in the truss we will introduce by special vectors Vi, i = 1,...,ns . Here the first component is the number of the joint at the beginning of the corresponding rod, the second is the number of the hinge at the end. Of course, the choice of special direction vectors of the rods is not associated with the sign of force in them. We have the following vectors:
V = [i,2(i + n +1)], Vi+n = [i,i + n + 2], Vi+2n = [i +1,i + n +1],
V+3n = [i + 1,2(i + n) + 1], V+4n = [i + n +1,2(i + n) + 1],
Vi+5n = [i + n + 2,2(i + n +1)], i = 1,...,n,
Vi+6n = [i + 2(n +1),i + 2n + 3],i = 1,...,2n -1.
The following vectors present the vertical side rods of the truss:
V8„ = [1,n + 2], r8n+1 = [1 + n,2n + 2].
Three rods of the supports correspond to the following vectors:
V„s-2 = [1,m-2], V„s- = [1,m-1], V„s = [n +1,m].
The length of the rods and the projection of the vector representations of these rods need to calculate the direction cosines:
h + h,i , h,i = xv2J - xvu, 12,i = yv2i - yvu, i = !,•••,ns ■
The first index in the number Vjtakes on the values 1 or 2 and corresponds to the number of components of the vector Vj; the second is a number of the rod. The matrix O guides of the cosines has the following elements
®k,i =-lj,t /li, k = 2Vi,i -2 + j, k < ns, j = 1,2, i = 1,...,ns,
®k i = /1, k = 2Vn -2 + j, k < ns, j = 1,2, i = 1,...,ns.
The forces can be determined as:
OS = B, (1)
where S = {S1,...,Sns} is the vector of forces in rods, B = {Px,\,Py,i'...'Px,«s,Py,ns} is the vector of the right parts (external loads applied to nodes). In order to calculate the deflection of the truss by the formula of Maxwell - Mohr:
ns -3
A = P£ S(P)S,(1)lt /(EF), it is necessary to determine the
i=1
force from the action of a single force applied at the middle of the lower zone (unit k+1, Figure 2). S(p) is the force in the rods from a given load and calculates the length of the rods . The rigidity of the rods EF accepts the same.
Figure 2. The numbering of the nodes in the truss, n=2k=4
Solution
We consider the case of loading of the upper belt truss when b = c. In this case, the non-zero elements of the vector of the right part of the system of equilibrium equations of the nodes will have the following: B2i =-P, i = 2n + 3,..., 4n + 2. Solution (1) with such right-hand side gives the vector of forces S{P). The vector of forces from a single power at node k+1 selected at B2( k+1) = 1 (the other elements are zero). The calculation is produced in the analytical form consistently for a truss with two panels, then four, three, etc. In every solution it is possible to identify specific components and their ratios. The deflection of the truss with k panels in half of the span takes the form
EF A t = P( Aka a + Bkb3 + Gkg3 + Dkd3 + Hkrf)/ (25b2), (2)
where d = V 16a2 + b1, g = Va2 + b1, h = J9a2 + b2 is the length of the braces.
Architecture and Engineering
Volume 1 Issue 3
The regularity in the sequence of coefficients can be found by the methods of Maple computer algebra system. With the command rgf_findrecur of package genfunc for next sequence 225, 4400, 21925, 69600, 169325, ... 17753600 of Ak of length 16 we obtain the recurrence equation as follows:
A = 2 A - + 2 A-2 - 6 A-3 + 6 A-5 - 2 A-6 - 2 A-7 + A-8
The operator rsolve gives the solution of this equation. The following general member of the sequence has the form:
Ak = 25(65k4 / 6 + (3(-1)k + 5)k2 /12 + (-1)k -1).
In a similar way the next coefficients are obtained:
Bk = (25/2)k (1 - (-1)k ),
Gk = (10/ 3)k4 + ((-1)k + 23 / 3)k2 + 3(1 - (-1)k ),
Dk = (10/ 3)k4 + ((-1)k +1 / 6)k2 - 3(1 - (-1)k ) / 4,
Hk = (5 / 6)k4 + (3(-1)k + 53)k2 /12 - 3(1 - (-1)k).
Analysis
Graph of (2) deflection A = ! Psum from the number of panels with constant span length L = 4an = 40 m and a given total magnitude of the vertical load Psum = 2nP (Figure 3, the dimensions in meters) shows that this dependence is non-monotonic for small heights of the truss. The kinks of the curves appear for small k and define "flashing" terms (-1f in the coefficients. Of course, the initial part of the curve with b = 1 m has no practical value, because the corresponding panel length is too big for this altitude truss a = 40/8 = 5 m, but the fact of crossing curves can be used in a practical sense in case of choosing the optimal size of a truss. The dependence of the deflection height (Figure 4) shows a clear minimum, which can also be used to optimize the structures. For a larger number of panels, extrem-um is stronger. It is interesting to perform the marginal analysis of the obtained formula. In the Maple system, using the limit operator (Delta/k*3,k=infinity) we find that the increase
of deflection is cubic lim A / k = 2b/15 and it is determined
k —
only by the height of the truss.
1000 -
500 -
b=6
7>=4
6=3
Figure 3. The dependence of the deflection on the number of panels
Figure 4. The dependence of deflection from the height
Conclusions
We proposed and calculated analytically statically determinate trusses in Maple regular scheme. Although the machines and mechanisms, in which the use of such schemes is assumed, undergo mainly dynamic loads, the static solution makes sense as a base, or estimated for more complex tasks, loads and types of supports. The solution has a polynomial form; the dependence is non-monotonic and detects extremum, which allows optimizing the design. The author provided the plots of deflection against a number of panels and height of the truss, showing the link between the deflection number of panels and height of the truss. Calculation and analysis are performed for a particular case when
b = c . Based on the proposed algorithm, this solution can be easily generalized to any instance. However, the final formula, in this case, will be more complex. Marginal analysis revealed the cubic nature of the change of the deflection depending on the number of panels fixed on the length of the passage. Without any changes, the whole calculation can be performed in other symbolic mathematics packages. Similar design with precise formulas can be embedded in building structures to give them architectural expression and optimum rigidity. The analysis and optimization of beam structures was examined in a number of related researches (Farshi and Alinia-ziazi, 2010; Heyman, 2010; Kaveh and Talataharib, 2009).
Mikhail Kirsanov — Pages 14-17 AN INDUCTIVE METHOD OF CALCULATION OF THE DEFLECTION OF THE TRUSS REGULAR TYPE
References
Farshi, B., Alinia-ziazi, A. (2010). Sizing optimization of truss structures by method of centers and force formulation. International Journal of Solids and Structures, 47, pp. 2508 - 2524. doi:10.1016/j.ijsolstr.2010.05.009
Heyman, J. (2010). Design of a simple steel truss. Proceedings of the Institution of Civil Engineers: Structures and Buildings, 163(1), pp. 53 - 56. doi: 10.1680/stbu.2009.163.1.53
Hutchinson, R.G., Fleck, N.A. (2005). Microarchitectured cellular solids - the hunt for statically determinate periodic trusses. ZAMMZ. Angew. Math. Mech., 85(9), pp.: 607 - 617. doi: 10.1002/zamm.200410208
Kaveh, A., Talataharib, S. (2009). Size optimization of space trusses using Big Bang-Big Crunch algorithm. Computers and Structures, 87, pp. 1129 - 1140. doi:10.1016/j.compstruc.2009.04.011
Kirsanov, M.N. (2014). Bending, torsion and the asymptotic analysis of the spatial core console [Izgib, kruchenie i asimptotich-eskij analiz prostranstvennoj sterzhnevoj konsoli]. Magazine of Civil Engineering, 5(49), 37 - 43. doi: 10.5862/MCE.49.4 (in Russian).
Kirsanov, M.N. (2015a). Deflection analysis of rectangular spatial coverage truss [Analiz progiba fermy pryamougol'nogo pros-transtvennogo pokrytiya]. Magazine of Civil Engineering, 1(53), pp. 32-38. doi: 10.5862/MCE.53.4 (in Russian).
Kirsanov, M.N. (2015b) Analysis of the deflection of a strut-type lattice girder truss [Analiz progiba reshetchatoj balochnoj fermy raspornogo tipa]. Magazine of Civil Engineering, 5(57), pp. 58 - 65. doi: 10.5862/MCE.57.5 (in Russian).
Monagan, M.B., Geddes, K.O., Heal, K.M., Labahn, G., Vorkoetter, S.M., McCarron, J., DeMarco, P. (2005). Maple 10 Advanced Programming Guide, Waterloo Maple Inc.
Tinkov, D.V. (2015). Comparative analysis of analytical solutions to the problem of truss structure deflection [Sravnitel'nyj analiz analiticheskih reshenij zadachi o progibe fermennyh konstrukcij], Magazine of Civil Engineering, 5(57), pp. 66 - 73. doi: 10.5862/ MCE.57.6 (in Russian).