Acessibilidade / Reportar erro

COMPARATIVE ANALYSIS OF THE MECHANICAL PERFORMANCE OF TIMBER TRUSSES STRUCTURAL TYPOLOGIES APPLYING COMPUTATIONAL INTELLIGENCE

ANÁLISE COMPARATIVA DO DESEMPENHO MECÂNICO DE TIPOLOGIAS ESTRUTURAIS DE TRELIÇAS DE MADEIRA APLICANDO INTELIGÊNCIA COMPUTACIONAL

ABSTRACT

Timber application is viable in constructive systems because of its mechanical properties, suitable for structural applications in engineering. Timber is even more interesting to use since it is a renewable source. Among its several applications, timber is widely used in roofing structures, with several typologies of trusses. Therefore, it is necessary to understand the behavior of plane truss in different loading conditions, including dead load, service load, and wind suction load. The mechanical performance of two trusses (Pratt and Scissor) was analyzed and compared, according to the Brazilian standard of timber structures ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997. (methods and calculus premises), with the finite element method, an algorithm of swarm intelligence optimization (structure weight minimization), and a parametric study. Based on minimum weight and maximum displacement as a function of span variation, Pratt typology presented lower weight (3-19%) when compared with Scissor, under the same span and loading conditions. Regarding maximum displacements, Pratt typology presented lower displacement values than the Scissor typology. The difference between these values ranged from two to seven times, indicating that Scissor typology can better distribute normal loads (maximum displacement closer to the normative limit displacement). Variance analysis (5% of significance) confirmed these results.

Keywords:
Pratt; Scissor; Optimization

RESUMO

Aplicação da madeira é viável em sistemas construtivos devido a suas propriedades mecânicas, que são compatíveis com as aplicações estruturais na engenharia A madeira é ainda mais interessante de se usar, uma vez que é uma fonte renovável. Dentre as diversas aplicações, a madeira é amplamente utilizada em estruturas de cobertura, com várias tipologias de treliças. Portanto, deve-se entender o comportamento das treliças planas em diferentes condições de carregamento, incluindo carga de peso próprio, carga de serviço e carga de sucção de vento. O desempenho mecânico de duas treliças (Pratt e Scissor) foi analisado e comparado, de acordo com a norma Brasileira de estruturas de madeira ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997. (métodos e premissas de cálculo), com o método dos elementos finitos, um algoritmo de otimização de inteligência de enxame (minimização do peso da estrutura) e um estudo paramétrico. Com base no peso mínimo e no deslocamento máximo em função da variação do vão, a tipologia Pratt apresentou um peso menor (3-19%) quando comparada com a Scissor, sob a mesma condição de carregamento e vãos. Em relação aos deslocamentos máximos, a tipologia Pratt apresentou valores de deslocamento menores do que a tipologia Scissor. A diferença entre estes valores variou de duas a sete vezes, indicando que a tipologia Scissor distribui melhor os esforços normais (deslocamentos máximos mais próximos do deslocamento limite normativo). Esses resultados foram confirmados pela análise de variância (5% de significância).

Palavras-Chave:
Pratt; Scissor; Otimização

1. INTRODUCTION

Timber has been used as a structural material since ancient times, for it is not only renewable (Bita and Tannert, 2018Bita HM, Tannert T. Numerical optimisation of novel connection for cross-laminated timber buildings. Engineering Structures. 2018;175:273–83. doi: 10.1016/j.engstruct.2018.08.020
https://doi.org/10.1016/j.engstruct.2018...
), but it also has good mechanical properties and is easy to handle with simple tools. This material has been widely used in Structural Engineering (Mam et al., 2020Mam K, Douthe C, Le Roy R, Consigny F. Shape optimization of braced frames for tall timber buildings: Influence of semi-rigid connections on design and optimization process. Engineering Structures. 2020;216:110692. doi: 10.1016/j.engstruct.2020.110692
https://doi.org/10.1016/j.engstruct.2020...
) and mainly in projects of rural buildings (roofing structures).

Regarding sustainability representation, the American agency Energy Information Administration – EIA (2011)Energy Information Administration - EIA. Emission of Greenhouse Gases in the United States 2009. Washington, DC: U.S. Energy Information Administration - EIA; 2011. states that, until 2017, construction sectors responded for 40% of greenhouse gas emission (Mayencourt and Mueller, 2020Mayencourt P, Mueller C. Hybrid analytical and computational optimization methodology for structural shaping: Material-efficient mass timber beams. Engineering Structures. 2020;215:110532. doi: 10.1016/j.engstruct.2020.110532
https://doi.org/10.1016/j.engstruct.2020...
).

Kumar et al. (2016)Kumar JP, Kumar JDC, Kumar MP. Quantitative Study of Howe Truss (A- Type) and Parallel Chord Scissor Truss (B- Type) by Applying External Pre-Stressing. International Journal of Engineering Research & Technology. 2016:5. doi: 10.17577/IJERTV5IS010409
https://doi.org/10.17577/IJERTV5IS010409...
claim that the truss structural system is widely used in roofs of residential and industrial buildings because it can overcome large spans. Timber truss is prevalent in roofing systems, with several geometric typologies. Therefore, we must understand the behavior of typologies to use them in engineering projects.

Regardless of the material selected for the structural system, designing a structural part is a hard and time-consuming task that allows for the design engineer to determine ideal conditions (mass reduction combined with normative requirements for mechanical performance) of a structural system. This task is complex, since these objectives so-called “ideal goals” are conflicting. As an example, reducing the mass or volume of a structure implies reducing its stiffness (Lemonge et al., 2021Lemonge ACC, Carvalho JPG, Hallak PH, Vargas DênisEC. Multi-objective truss structural optimization considering natural frequencies of vibration and global stability. Expert Systems with Applications. 2021;165:113777. doi: 10.1016/j.eswa.2020.113777
https://doi.org/10.1016/j.eswa.2020.1137...
) and result in greater displacements.

In structural projects, computational science was essential for designers to act quickly regarding development, scale reproduction, and improving test capacity and refinement.

With this evolution, structural optimization became a viable tool, and nowadays it is applied in several structural systems. Regarding timber structures, it is possible to observe some applications in timber beams (Mayencourt et al., 2017Mayencourt PL, Giraldo JS, Wong E, Mueller CT. Computational Structural Optimization and Digital Fabrication of Timber Beams. Proceedings of IASS Annual Symposia 2017. 2017:1–10.; Pech et al., 2019Pech S, Kandler G, Lukacevic M, Füssl J. Metamodel assisted optimization of glued laminated timber beams by using metaheuristic algorithms. Engineering Applications of Artificial Intelligence. 2019;79:129–41. doi: 10.1016/j.engappai.2018.12.010
https://doi.org/10.1016/j.engappai.2018....
; Schietzold et al., 2021Schietzold FN, Graf W, Kaliske M. Multi-objective Optimization of Tree Trunk Axes in Glulam Beam Design Considering Fuzzy Probability Based Random Fields. ASCE-ASME J. Risk and Uncert in Engrg. Sys. Part B Mech. Engrg. 2021. doi: 10.1115/1.4050370
https://doi.org/10.1115/1.4050370...
), braced frames in timber (Mam et al., 2020Mam K, Douthe C, Le Roy R, Consigny F. Shape optimization of braced frames for tall timber buildings: Influence of semi-rigid connections on design and optimization process. Engineering Structures. 2020;216:110692. doi: 10.1016/j.engstruct.2020.110692
https://doi.org/10.1016/j.engstruct.2020...
), and timber trusses (Villar et al., 2016Villar JR, Vidal P, Fernández MS, Guaita M. Genetic algorithm optimisation of heavy timber trusses with dowel joints according to Eurocode 5. Biosystems Engineering. 2016;144:115–32. doi: 10.1016/j.biosystemseng.2016.02.011
https://doi.org/10.1016/j.biosystemseng....
; Villar-García et al., 2019Villar-García JR, Vidal-López P, Rodríguez-Robles D, Guaita M. Cost optimisation of glued laminated timber roof structures using genetic algorithms. Biosystems Engineering. 2019;187:258–77. doi: 10.1016/j.biosystemseng.2019.09.008
https://doi.org/10.1016/j.biosystemseng....
; Ching, 2020Ching HYE. Truss topology optimization of steel-timber structures for embodied carbon objectives [thesis]. Massachusetts Institute of Technology -MIT; 2020.). Few studies have addressed the optimization of trussed structures (or other structural systems) that follow the fulfillment normative performance requirements to compare typologies which are generally used in structural projects.

In Optimization, there are two types of algorithms: probabilistic and deterministic. The probabilistic model incorporates aspects of random variations in its formulations (Amaran et al., 2016Amaran S, Sahinidis NV, Sharda B, Bury SJ. Simulation optimization: a review of algorithms and applications. Ann Oper Res. 2016;240:351–80. doi: 10.1007/s10479-015-2019-x
https://doi.org/10.1007/s10479-015-2019-...
). Given a specific input dataset, the procedure executed by the algorithm of a deterministic method will always reproduce the same output (Lin et al., 2012Lin M-H, Tsai J-F, Yu C-S. A Review of Deterministic Optimization Methods in Engineering and Management. Mathematical Problems in Engineering. 2012;2012:e756023. doi: 10.1155/2012/756023
https://doi.org/10.1155/2012/756023...
). Deterministic methods can conduct or approach the global minimum values of the Objective Function (OF), but since derivates of the OF are hard to obtain, they have a high computational cost (Lin et al., 2012Lin M-H, Tsai J-F, Yu C-S. A Review of Deterministic Optimization Methods in Engineering and Management. Mathematical Problems in Engineering. 2012;2012:e756023. doi: 10.1155/2012/756023
https://doi.org/10.1155/2012/756023...
; Segovia-Hernández et al., 2015Segovia-Hernández JG, Hernández S, Bonilla Petriciolet A. Reactive distillation: A review of optimal design using deterministic and stochastic techniques. Chemical Engineering and Processing: Process Intensification. 2015;97:134–43. doi: 10.1016/j.cep.2015.09.004
https://doi.org/10.1016/j.cep.2015.09.00...
). Regardless, deterministic methods achieve the global minimum of the function with less iterations than probabilistic methods (Huang et al., 2017Huang C, Hami AE, Radi B. Metamodel-based inverse method for parameter identification: elastic–plastic damage model. Engineering Optimization. 2017;49:633–53. doi: 10.1080/0305215X.2016.1206537
https://doi.org/10.1080/0305215X.2016.12...
).

Probabilistic optimization technique, has a subdivision of methods called bio-inspired, as their mathematical formulation is inspired by the nature. This subgroup includes the Firefly Algorithm (FA), which is prevalent in structural optimization processes (Gandomi et al., 2011Gandomi AH, Yang X-S, Alavi AH. Mixed variable structural optimization using Firefly Algorithm. Computers & Structures. 2011;89:2325–36. doi: 10.1016/j.compstruc.2011.08.002.
https://doi.org/10.1016/j.compstruc.2011...
; Lieu et al., 2018Lieu QX, Do DTT, Lee J. An adaptive hybrid evolutionary firefly algorithm for shape and size optimization of truss structures with frequency constraints. Computers & Structures. 2018;195:99–112. doi: 10.1016/j.compstruc.2017.06.016
https://doi.org/10.1016/j.compstruc.2017...
; Pereira et al., 2020Pereira LLM, Santos DC, Moraes MHM, Gonçalves Filho GM, Ancioto Junior EM, Pereira Junior WM, et al. Estudo de Sensibilidade do Algoritmo de Colônia de Vagalumes para um Problema de Engenharia Envolvendo Dimensionamento de Treliças. Tend. Mat. Apl. Comput. 2020;21:583. doi: 10.5540/tema.2020.021.03.583
https://doi.org/10.5540/tema.2020.021.03...
; Talatahari et al., 2014Talatahari S, Gandomi AH, Yu n GJ. Optimum design of tower structures using Firefly Algorithm. The Structural Design of Tall and Special Buildings. 2014;23:350–61. doi: 10.1002/tal.1043
https://doi.org/10.1002/tal.1043...
).

From the viable optimal design perspective, optimization models are divided into three basic strands: a) Section Optimization; b) Shape Optimization; and c) Topology Optimization (Zhou et al., 2004Zhou M, Pagaldipti N, Thomas HL, Shyy YK. An integrated approach to topology, sizing, and shape optimization. Structural and Multidisciplinary Optimization. 2004;26:308–17. doi: 10.1007/s00158-003-0351-2
https://doi.org/10.1007/s00158-003-0351-...
).

Several typologies are used to design timber roofing structures, with emphasis on Pratt and Scissor (Gomes, 2020Gomes AFF. Influência do valor médio do módulo de elasticidade à compressão paralela às fibras no projeto de treliças em madeira [thesis]. Universidade Federal de São Carlos - UFSCar; 2020.). However, it is hard to accurately compare mechanical performance with the usual design strategy and to choose the most suitable method for a design. Therefore, this research intended to assess and to compare the mechanical performances of Pratt and Scissor typologies, with the aid of the Finite Element Method (FEM), the standard premises and the calculation methods from ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997., the Firefly Algorithm (FA) (minimization of structural weight), and a parametric study [variation of the span (L) and of the horizontal distance (B) and vertical distance (H) between nodes].

2. MATERIAL AND METHODS

To design timber trusses, precepts were used according to ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.. The FA, the penalization method and the calculation of the OF were used alike and are presented below, including details of the parametric study and description of the statistic tool used to compare mechanical behavior between the two typologies considered.

Notably, the classic formulation of the FEM was considered in the deduction of stiffness matrix of the truss bar element (two nodes and four degrees of freedom in all). Based on the solution of the system of equilibrium equation (equilibrium between external and internal nodal loads) of the FEM, nodal displacements, specific strains, stresses, and normal loads on each truss bar were determined. Stresses and normal loads were essential for the structure design.

2.1.Algorithm used for truss bars optimization

The Firefly Algorithm (FA) was used, a bioinspired probabilistic optimization proposed by Yang (2008)Yang X-S. Nature-Inspired Metaheuristic Algorithms. 1o. Frome: Luniver Press; 2008.. FA is a method of population characteristic, in which one particle walks through the sample space searching for the optimal viable solution. In this method, concepts of random variables are used to generate the initial population, which is a random event limited by the problem (Yang, 2008Yang X-S. Nature-Inspired Metaheuristic Algorithms. 1o. Frome: Luniver Press; 2008.).

This algorithm design was inspired by the bioluminescence and the influence of iteration between fireflies at mating season. Therefore, the FA optimization method is based on how fireflies can emit light and be perceived by other individuals of the same population.

To design the algorithm, Yang (2008)Yang X-S. Nature-Inspired Metaheuristic Algorithms. 1o. Frome: Luniver Press; 2008. defined some precepts to aid in the development, including: all fireflies have a single genus and having a single gender, they are attracted to each other; the attraction capacity of each firefly is proportional to their own brightness, but larger distances between individuals decrease such capacity.

When the initial population is created, the firefly (design variable) begins a random walk (Equation 1), so that x “moves” according to an update function of the design variable (ω), in which x is the vector of design variables, ω is the vector of update function of x, and t is the number of iterations.

Eq.1 x t + 1 = x t + ω t

Based on this new direction, new positions and possible solutions for the optimal point of design are originated (Wang et al., 2017Wang H, Wang W, Zhou X, Sun H, Zhao J, Yu X, et al. Firefly algorithm with neighborhood attraction. Information Sciences. 2017;382–383:374–87. doi: 10.1016/j.ins.2016.12.024
https://doi.org/10.1016/j.ins.2016.12.02...
). Thus, fireflies move at each step of the iterative process (Equation 2).

Eq.2 ω t = β ( x j t x i t ) + α ( η 0 , 5 ε )

Where β is a term of attractiveness between fireflies i and j; xi is the firefly i; xj is the firefly j; η is the vector of random numbers between 0 and 1; α is the randomness factor (Equation 3); and ε is a unit vector.

Eq.3 α = α min + ( α max α min ) θ t

Randomness factor α (Equation 3) follows an exponential decay behavior according to the number of iterations t, and θ is the constant decay of 0.98.

β is the attractiveness between the swarm of fireflies (Equation 4), in which βo is attractiveness for a distance r=0; rij is a Euclidean distance between fireflies i and j (Equation 5); and γ is the parameter of light absorption (Equation 6).

Eq.4 β = β 0 e γ r i j 2 β 0 / ( 1 + γ r i j 2 )

Eq.5 r i j = x i x j = k = 1 d ( x i , k x j , k ) 2

Eq.6 γ = 1 / ( x max x min ) 2

Where k is the k-th component of the design vector variable x; d is the amount of design variables; xmax is the maximum of design variables; and xmin is the minimum of design variables.

FA or any other probabilistic method of optimization requires attention in defining algorithm parameters (attractiveness: β and γ; randomness: α).

Parameter γ is the attractiveness variation {γ ∈ [0, ∞)}, essential to determine convergence speed and algorithm behavior. It mostly ranges between 0.1 and 10 (Yang, 2008Yang X-S. Nature-Inspired Metaheuristic Algorithms. 1o. Frome: Luniver Press; 2008.).

2.1.1. Optimization algorithm parameters

Input parameters of FA were based on the study by Pereira et al. (2020)Pereira LLM, Santos DC, Moraes MHM, Gonçalves Filho GM, Ancioto Junior EM, Pereira Junior WM, et al. Estudo de Sensibilidade do Algoritmo de Colônia de Vagalumes para um Problema de Engenharia Envolvendo Dimensionamento de Treliças. Tend. Mat. Apl. Comput. 2020;21:583. doi: 10.5540/tema.2020.021.03.583
https://doi.org/10.5540/tema.2020.021.03...
. In our study, a total population (Npop) of 15 individuals, randomly generated, were considered. All parameters were based on the sensitivity study by Pereira et al. (2020)Pereira LLM, Santos DC, Moraes MHM, Gonçalves Filho GM, Ancioto Junior EM, Pereira Junior WM, et al. Estudo de Sensibilidade do Algoritmo de Colônia de Vagalumes para um Problema de Engenharia Envolvendo Dimensionamento de Treliças. Tend. Mat. Apl. Comput. 2020;21:583. doi: 10.5540/tema.2020.021.03.583
https://doi.org/10.5540/tema.2020.021.03...
, including: 500 generations (Ngen), Attractiveness between fireflies (βo) of 0.90, Minimum randomness factor mim) of 0.20, Maximum randomness factor (αmax) of 1.00, and Penalty factor (Rp) of 105.

2.1.2. OF

Optimization intends to minimize the total weight of the structural system, considering constraints of nodal displacements, mechanical resistance of the bars, and geometric criteria that can cause lateral instability. Objective Function (OF) is the total weight of a truss (Equation 7), in which Ai. is the cross-section area of a bar i; ρi. is the density of the material; x0i. is the length of truss bar i; and n is the number of truss bars.

Eq.7 O F ( A i , ρ i , L 0 i ) = i = 1 n A i ρ i L 0 i

Other factors besides weight must be considered to optimize truss manufacture. For example, in industrial manufacturing, it is often more interesting to produce similar elements (same sections) than several sections. This and other aspects can affect the cost of the structure (Šilih et al., 2010Šilih S, Kravanja S, Premrov M. Shape and discrete sizing optimization of timber trusses by considering of joint flexibility. Advances in Engineering Software. 2010;41:286–94. doi: 10.1016/j.advengsoft.2009.07.002
https://doi.org/10.1016/j.advengsoft.200...
; Jelusic and Kravanja, 2017Jelusic P, Kravanja S. Optimal design of timber-concrete composite floors based on the multi-parametric MINLP optimization. Composite Structures. 2017;179:285–93. doi: 10.1016/j.compstruct.2017.07.062
https://doi.org/10.1016/j.compstruct.201...
). However, our study considers only the weight for optimization, focusing on the final weight of the truss.

For optimization problems, constraint restrictions must be treated. Therefore, the external penalty technique was used (Kuri-Morales and Gutiérrez-García, 2002Kuri-Morales AF, Gutiérrez-García J. Penalty Function Methods for Constrained Optimization with Genetic Algorithms: A Statistical Analysis. In: Coello Coello CA, de Albornoz A, Sucar LE, Battistutti OC, editors. MICAI 2002: Advances in Artificial Intelligence, Berlin, Heidelberg: Springer; 2002, p. 108–17. doi: 10.1007/3-540-46016-0_12
https://doi.org/10.1007/3-540-46016-0_12...
; Yeniay, 2005Yeniay Ö. Penalty Function Methods for Constrained Optimization with Genetic Algorithms. Mathematical and Computational Applications. 2005;10:45–56. doi: 10.3390/mca10010045
https://doi.org/10.3390/mca10010045...
), in which the OF is modified to become a pseudo OF, where gj. represents inequality constraints and hk restrictions on equality. For penalty, the OF used the complete form of the penalty method (Equation 8), resulting in Penalized OF (W) (Equation 9).

Eq.8 P ( x ) = j = 1 m max [ 0 , g j ( x ) ] 2 + k = 1 n [ h k ( x ) ] 2

Eq.9 W ( A i , ρ i , L 0 i , x i ) = O F ( A i , ρ i , L 0 i ) + R p P ( x )

Note that, P(x) (Equation 8) is the static outer penalty function j; k is j-th inequality constraints and k-th equality constraints, respectively; m, n are the total number of inequality and equality constraints, respectively; x is the solution vector (random population); g, h are the set of inequality and equality constraints; and W(x) is the penalized OF.

2.2. Description of trusses and design variables

This study analyzed two plane trusses: Pratt and Scissor. Trusses were divided into seven spans with a total length of 5, 7.5, 10, 12.5, 15, 17.5, and 20 meters, based on Gomes (2020)Gomes AFF. Influência do valor médio do módulo de elasticidade à compressão paralela às fibras no projeto de treliças em madeira [thesis]. Universidade Federal de São Carlos - UFSCar; 2020.. For each of the trusses, the optimization algorithm was performed 30 times to obtain a dispersion of the results. Nodal distances (B=L/4 and H=L/8) were used for both typologies, and design variables (bar position) considered x1, x2, x3, x4 and x5 (Figure 1).

Figure 1
Generic representation of the nodal distances of trusses.
Figura 1
Representação genérica das distâncias nodais das treliças.

xi is the generic design variable (Equation 10) in which xi is the variable vector of design i; bi is the thickness of the cross section of variable i; and hi is the height of the cross section of variable i. For each type of bar, a design variable was used: x1 for the bottom chord, x2 for the top chord, x3 for secondary vertical bars, x4 for the main vertical bars, and x5 for diagonal bars.

Eq.10 x i = ( b i ; h i )

For optimization problems, a set of distinct variables was considered, in which each variable had an ordered set of values (Hsu, 1985Hsu CS. A discrete method of optimal control based upon the cell state space concept. J Optim Theory Appl. 1985;46:547–69. doi: 10.1007/BF00939159
https://doi.org/10.1007/BF00939159...
). Therefore, we considered nominal values for coniferous sawn timber according to ABNT NBR ISO 3179 (2011)Associação Brasileira de Normas Técnicas - ABNT. Madeira serrada de coníferas – Dimensões nominais - NBR ISO 3179. Rio de Janeiro – RJ: ABNT; 2011.. Based on the standard nominal values, variations of the design variables were: thickness of the cross-section of a generic bar i (bi) assuming values of [16, 19, 22, 25, 32, 38, 50, 63, 75, 100, 125, 150, 175, 200, 250, and 300 mm] and cross-section height of a generic bar i (hi) assuming values of [75, 100, 115, 125, 150, 160, 175, 200, 225, 250, 275, and 300 mm].

For optimization process, coniferous timber was considered and belonging to the resistance class C30, whose properties are presented in the ABNT NBR 7190 standard (1997). Sawn timber was classified as permanent loading class, with moisture class I, 12% equilibrium moisture, and 1st category timber.

2.3. Loads and limit states

The vertical loads to which the trusses are submitted consist of roofing load for a sandwich tile (G) of 150 N/m²; structure dead load (DL) of 150 N/m²; service load (S) of 25 N/m²; and wind suction load (W) of 600 N/m².

To verify Limit States, three combinations were adopted based on the normative requirements of ABN NBR 7190 (1997), two in Ultimate Limit State (ULS): combination 1 (1.3 DL + 1.3 G + 1.4 S + 0.7 W); combination 2 (1.0 DL + 1.0 G + 0.56 S + 1.4 W); and one for the Serviceability Limit State (SLS): combination 3 (1.0 DL + 1.0 G + 0.3 S + 0.2 W). Soon Afterwards, such combinations were checked in all OF constraints.

Based on the area of influence and considering a distance of 5 meters between gantries, distributed loads were transformed into nodal loads. The IDE (Integrated Development Environment) of the MATLAB® software was used for structural analysis and optimization.

2.4. Design constraints

Constraints x (Equations 11, 12, 13, 14, 15) where i is a generic bar (i=1,...,,nbars), nbars is the number of analyzed truss bars , n is a generic node (n=1,...,,nnodes), and nnodes is the number of analyzed truss joints.

Eq.11 g j ( x ) = σ i / σ lim 1 0

Eq.12 g j ( x ) = u n / u lim 1 0

Eq.13 g j ( x ) = λ i / λ lim 1 0

Eq.14 g j ( x ) = A min / A i 1 0

Eq.15 g j ( x ) = b min / b i 1 0

ULS (Equation 11) was checked for normal stresses (σi) (tensile and compression), in which σlim is the normal stress limit. SLS (Equation 12) was also checked, in which un is a nodal displacement and ulim is the displacement limit, with ulim = L/200 (ABNT NBR 7190, 1997Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.), in which L is the span of the analyzed structural system. The geometric constraint (Equation 13) of bars based on the slenderness index (λi) was also analyzed, in which λlim is the limit slenderness index, with 140 as maximum value according to ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.. Geometric constraints generated by the ratio between the cross-sectional area (Ai) and minimum area (Amin) of a bar (Equation 14), with a minimum of 50 cm² for the cross-sectional area (ABNT NBR 7190, 1997Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.). Finally, we analyzed geometric constraints (Equation 15) that associate the thickness (bi) and minimum thickness (bmin) of a given bar. ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997. limits the minimum thickness to 5 cm for main bars.

2.5. Analysis of Variance (ANOVA)

The analysis of variance (ANOVA), with 5% significance level, was used to verify potential differences in the mechanical performance of the two typologies. For this, the ANOVA was considered on the results of constraints with normal loads and displacements. For each span (5.00, 7.50, 10.00, ..., 20.00 m) considering a type of constraint (stress or displacement), we verified whether both typologies (Pratt and Scissor) were statistically equivalent to each other regarding the mean restriction. According to ANOVA, p-value (probability p) greater than or equal to the significance level implies an equivalence of restriction means, while lower significance implies a non-equivalence.

Following ANOVA, the Anderson-Darling test was also used with 5% significance level to verify the normality in the distribution of ANOVA residues. According to the test, p-value greater than or equal to the significance level implies normality of residue distribution, therefore validating the ANOVA results.

3. RESULTS

This section shows the results obtained from the optimization process of trusses. Results (Table 1) of the 30 optimization processes for different trusses can be observed considered, in which Wmax and Wmin represent the maximum and minimum values of the penalized OF, respectively; amplitude (A), median (μ), mean (x), standard deviation (σ), and feasibility rate (FR) represent the ratio between tests with satisfied constraints and total number of tests (30). To sum the results, trusses of type x-y will be identified, where “x” is the typology (P for Pratt and S truss for Scissor truss) and “y” is the truss span in meters (5, 7.5, 10, 12.5, 15, 17.5, and 20).

Table 1
Summary of results obtained from truss optimization process.
Tabela 1
Resumo dos resultados obtidos do processo de otimização das treliças.

With Box plot, we observed the weight variation of Pratt and Scissor trusses as a function of the span (L) (Figure 2).

Figure 2
Box plot of Minimum weight (W) of trusses as a function of the variation of span L.
Figura 2
Box plot do Peso mínimo (W) das treliças em função da variação do vão L.

After optimization analyses, we created curves (via quadratic regression models) of the OF (W -minimum weight) as a function of the span (L) for conditions delimited in this study (Figure 3a). Such curves (Figure 3a) allow estimating the minimum weight of the trusses as a function of span, even for L values not included here.

Figure 3
Curves of: a) Minimum weight (W) of trusses as a function of span (L); and b) Maximum displacement as a function of span (L).
Figura 3
Curvas do: a) Peso mínimo (W) das treliças em função do vão (L); e b) Máximo deslocamento (u) em função do vão (L).

Once we had the total weight of the trusses, an analysis of the maximum displacement of each truss was required to ensure the structure would follow normative requirements. Similarly, truss weight curves were created to predict the maximum displacement (Figure 3b).

We obtained the convergence curve from the best response among 30 repetitions of the weight of the trusses of 5 and 20 meters (Figure 4).

Figure 4
Convergence curves for minimum weight (W) for P-5, P-20, S-5, and S-20.
Figura 4
Curvas de convergência para o peso mínimo (W) para P-5, P-20, S-5 e S-20.

3.1. Assessment of stress constraints

Constraints showed how the tension, displacement, slenderness, area, and dimensions of cross-sections are close to the established limit. Restrictions (g values) with values closer to 0 indicate more loaded bars (stress). Stress constraints were considered the combination 1 (ULS combination, with the service load as the main accidental load) and combination 2 (ULS combination, with the wind load as the main accidental load).

Based on the ANOVA of normal stress constraints, we verified major differences between typologies. Normal stress restrictions had p-values between 0.004 and 0.046, lower than the significance level (α = 0.05), indicating great differences between means. P-values of the Anderson-Darling normality test on ANOVA residues were greater than 0.05 (5%), thus validating the ANOVA model. Moreover, Pratt truss had mean values between −0.9087 and −0.8333, while Scissor truss had values between −0.7989 and −0.6646.

3.2. Assessment of displacement constraints

Based on the analysis of variance (ANOVA) of displacement constraints, we verified major differences between typologies. P-values of the Anderson-Darling normality test were higher than the significance level, which validates the ANOVA model. All trusses had p-values between 0.002 and 0.004. Displacement restrictions presented p-values lower than the significance level (α=0.05), indicating significant differences between the mean values of such restrictions. The Pratt truss had mean values between −0.9875 and −0.9380, while the Scissor truss values between −0.7457 and −0.6545.

4. DISCUSSIONS

Our results help to understand the mechanical behavior of Pratt and Scissor trusses. From the Box plot (Figure 2) of weight variation as a function of span (L), we observed that the dispersion of the results increases according to the increase of the L parameter. This variation occurs because when the span increases, loads also increase, and it is harder to find viable minimum values through optimization.

A graph of the minimum weights of trusses (Figure 3a) shows that trusses differ according to the spans. As a result, the Pratt truss presented OF values (OF - minimum weight) between 3% and 19% lower in relation to the Scissor truss.

Regarding displacement (Figure 3b), both trusses are within the ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997. standard. Thus, based on resulting curves, maximum displacements of the Pratt truss are two to seven times lower than those of the Scissor truss.

Regarding convergence of the best result (Figure 4), the Pratt truss was the first to converge to the minimum weight for the 5 m span truss. For the 20 m span truss, Scissor converged firstly. P-5 (Pratt with 5 m span) converged in 160 iterations, with 456 iterations for S-5, 433 iterations for P-20, and 278 iterations for S-20, based on a tolerance rate of 10-2.

The results of the stress constraints evaluation showed that the Scissor truss presented better distribution of normal loads (mean value of the constraint closer to 0), further using the load capacity of the trussed structures. Analyzing the results of the displacement constraints evaluation showed that Scissor is more displaceable than Pratt, and that both trusses respected the limit (L/200) established by the Brazilian standard ABNT NBR 7190 (1997)Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.. This fact allows choosing the best truss for roofing structures.

Normal stress and displacement constraints results showed that, the normal stresses guided the design for Pratt truss, whereas regarding the Scissor truss, displacements were the main guide.

5. CONCLUSION

Our results showed that the Scissor truss had greater weight, since the curve (relation between minimum weight and span) of the Pratt truss provided OF values (minimum weight) between 3% and 19% lower than the Scissor truss. Regarding displacement-span curves, the Pratt truss presented displacements between two and seven times lower than Scissor truss, and nodal displacements of Scissor trusses were closer to the L/200. Therefore, the Pratt truss was less displaceable.

This establishes parameters for choosing between the two typologies. Since both trusses satisfy the design criteria, and the designer can choose between better mechanical performance, displaceability and material consumption.

For broader results, other common typologies for performance analysis must be studied, focusing on properties of hardwood timber (dicotyledons) and other optimization methods.

6. ACKNOWLEDGEMENTS

This study was partially supported by the Coordination for the Improvement of Higher Education Personnel (CAPES) – Brazil – Finance Code 001.

7. REFERENCES

  • Amaran S, Sahinidis NV, Sharda B, Bury SJ. Simulation optimization: a review of algorithms and applications. Ann Oper Res. 2016;240:351–80. doi: 10.1007/s10479-015-2019-x
    » https://doi.org/10.1007/s10479-015-2019-x
  • Associação Brasileira de Normas Técnicas - ABNT. Madeira serrada de coníferas – Dimensões nominais - NBR ISO 3179. Rio de Janeiro – RJ: ABNT; 2011.
  • Associação Brasileira de Normas Técnicas - ABNT. Projeto de estruturas de madeira - NBR 7190. Rio de Janeiro – RJ: ABNT; 1997.
  • Bita HM, Tannert T. Numerical optimisation of novel connection for cross-laminated timber buildings. Engineering Structures. 2018;175:273–83. doi: 10.1016/j.engstruct.2018.08.020
    » https://doi.org/10.1016/j.engstruct.2018.08.020
  • Ching HYE. Truss topology optimization of steel-timber structures for embodied carbon objectives [thesis]. Massachusetts Institute of Technology -MIT; 2020.
  • Energy Information Administration - EIA. Emission of Greenhouse Gases in the United States 2009. Washington, DC: U.S. Energy Information Administration - EIA; 2011.
  • Gandomi AH, Yang X-S, Alavi AH. Mixed variable structural optimization using Firefly Algorithm. Computers & Structures. 2011;89:2325–36. doi: 10.1016/j.compstruc.2011.08.002.
    » https://doi.org/10.1016/j.compstruc.2011.08.002.
  • Gomes AFF. Influência do valor médio do módulo de elasticidade à compressão paralela às fibras no projeto de treliças em madeira [thesis]. Universidade Federal de São Carlos - UFSCar; 2020.
  • Hsu CS. A discrete method of optimal control based upon the cell state space concept. J Optim Theory Appl. 1985;46:547–69. doi: 10.1007/BF00939159
    » https://doi.org/10.1007/BF00939159
  • Huang C, Hami AE, Radi B. Metamodel-based inverse method for parameter identification: elastic–plastic damage model. Engineering Optimization. 2017;49:633–53. doi: 10.1080/0305215X.2016.1206537
    » https://doi.org/10.1080/0305215X.2016.1206537
  • Jelusic P, Kravanja S. Optimal design of timber-concrete composite floors based on the multi-parametric MINLP optimization. Composite Structures. 2017;179:285–93. doi: 10.1016/j.compstruct.2017.07.062
    » https://doi.org/10.1016/j.compstruct.2017.07.062
  • Kumar JP, Kumar JDC, Kumar MP. Quantitative Study of Howe Truss (A- Type) and Parallel Chord Scissor Truss (B- Type) by Applying External Pre-Stressing. International Journal of Engineering Research & Technology. 2016:5. doi: 10.17577/IJERTV5IS010409
    » https://doi.org/10.17577/IJERTV5IS010409
  • Kuri-Morales AF, Gutiérrez-García J. Penalty Function Methods for Constrained Optimization with Genetic Algorithms: A Statistical Analysis. In: Coello Coello CA, de Albornoz A, Sucar LE, Battistutti OC, editors. MICAI 2002: Advances in Artificial Intelligence, Berlin, Heidelberg: Springer; 2002, p. 108–17. doi: 10.1007/3-540-46016-0_12
    » https://doi.org/10.1007/3-540-46016-0_12
  • Lemonge ACC, Carvalho JPG, Hallak PH, Vargas DênisEC. Multi-objective truss structural optimization considering natural frequencies of vibration and global stability. Expert Systems with Applications. 2021;165:113777. doi: 10.1016/j.eswa.2020.113777
    » https://doi.org/10.1016/j.eswa.2020.113777
  • Lieu QX, Do DTT, Lee J. An adaptive hybrid evolutionary firefly algorithm for shape and size optimization of truss structures with frequency constraints. Computers & Structures. 2018;195:99–112. doi: 10.1016/j.compstruc.2017.06.016
    » https://doi.org/10.1016/j.compstruc.2017.06.016
  • Lin M-H, Tsai J-F, Yu C-S. A Review of Deterministic Optimization Methods in Engineering and Management. Mathematical Problems in Engineering. 2012;2012:e756023. doi: 10.1155/2012/756023
    » https://doi.org/10.1155/2012/756023
  • Mam K, Douthe C, Le Roy R, Consigny F. Shape optimization of braced frames for tall timber buildings: Influence of semi-rigid connections on design and optimization process. Engineering Structures. 2020;216:110692. doi: 10.1016/j.engstruct.2020.110692
    » https://doi.org/10.1016/j.engstruct.2020.110692
  • Mayencourt P, Mueller C. Hybrid analytical and computational optimization methodology for structural shaping: Material-efficient mass timber beams. Engineering Structures. 2020;215:110532. doi: 10.1016/j.engstruct.2020.110532
    » https://doi.org/10.1016/j.engstruct.2020.110532
  • Mayencourt PL, Giraldo JS, Wong E, Mueller CT. Computational Structural Optimization and Digital Fabrication of Timber Beams. Proceedings of IASS Annual Symposia 2017. 2017:1–10.
  • Pech S, Kandler G, Lukacevic M, Füssl J. Metamodel assisted optimization of glued laminated timber beams by using metaheuristic algorithms. Engineering Applications of Artificial Intelligence. 2019;79:129–41. doi: 10.1016/j.engappai.2018.12.010
    » https://doi.org/10.1016/j.engappai.2018.12.010
  • Pereira LLM, Santos DC, Moraes MHM, Gonçalves Filho GM, Ancioto Junior EM, Pereira Junior WM, et al. Estudo de Sensibilidade do Algoritmo de Colônia de Vagalumes para um Problema de Engenharia Envolvendo Dimensionamento de Treliças. Tend. Mat. Apl. Comput. 2020;21:583. doi: 10.5540/tema.2020.021.03.583
    » https://doi.org/10.5540/tema.2020.021.03.583
  • Schietzold FN, Graf W, Kaliske M. Multi-objective Optimization of Tree Trunk Axes in Glulam Beam Design Considering Fuzzy Probability Based Random Fields. ASCE-ASME J. Risk and Uncert in Engrg. Sys. Part B Mech. Engrg. 2021. doi: 10.1115/1.4050370
    » https://doi.org/10.1115/1.4050370
  • Segovia-Hernández JG, Hernández S, Bonilla Petriciolet A. Reactive distillation: A review of optimal design using deterministic and stochastic techniques. Chemical Engineering and Processing: Process Intensification. 2015;97:134–43. doi: 10.1016/j.cep.2015.09.004
    » https://doi.org/10.1016/j.cep.2015.09.004
  • Šilih S, Kravanja S, Premrov M. Shape and discrete sizing optimization of timber trusses by considering of joint flexibility. Advances in Engineering Software. 2010;41:286–94. doi: 10.1016/j.advengsoft.2009.07.002
    » https://doi.org/10.1016/j.advengsoft.2009.07.002
  • Talatahari S, Gandomi AH, Yu n GJ. Optimum design of tower structures using Firefly Algorithm. The Structural Design of Tall and Special Buildings. 2014;23:350–61. doi: 10.1002/tal.1043
    » https://doi.org/10.1002/tal.1043
  • Villar JR, Vidal P, Fernández MS, Guaita M. Genetic algorithm optimisation of heavy timber trusses with dowel joints according to Eurocode 5. Biosystems Engineering. 2016;144:115–32. doi: 10.1016/j.biosystemseng.2016.02.011
    » https://doi.org/10.1016/j.biosystemseng.2016.02.011
  • Villar-García JR, Vidal-López P, Rodríguez-Robles D, Guaita M. Cost optimisation of glued laminated timber roof structures using genetic algorithms. Biosystems Engineering. 2019;187:258–77. doi: 10.1016/j.biosystemseng.2019.09.008
    » https://doi.org/10.1016/j.biosystemseng.2019.09.008
  • Wang H, Wang W, Zhou X, Sun H, Zhao J, Yu X, et al. Firefly algorithm with neighborhood attraction. Information Sciences. 2017;382–383:374–87. doi: 10.1016/j.ins.2016.12.024
    » https://doi.org/10.1016/j.ins.2016.12.024
  • Yang X-S. Nature-Inspired Metaheuristic Algorithms. 1o. Frome: Luniver Press; 2008.
  • Yeniay Ö. Penalty Function Methods for Constrained Optimization with Genetic Algorithms. Mathematical and Computational Applications. 2005;10:45–56. doi: 10.3390/mca10010045
    » https://doi.org/10.3390/mca10010045
  • Zhou M, Pagaldipti N, Thomas HL, Shyy YK. An integrated approach to topology, sizing, and shape optimization. Structural and Multidisciplinary Optimization. 2004;26:308–17. doi: 10.1007/s00158-003-0351-2
    » https://doi.org/10.1007/s00158-003-0351-2

Publication Dates

  • Publication in this collection
    08 July 2022
  • Date of issue
    2022

History

  • Received
    06 July 2021
  • Accepted
    19 Oct 2021
Sociedade de Investigações Florestais Universidade Federal de Viçosa, CEP: 36570-900 - Viçosa - Minas Gerais - Brazil, Tel: (55 31) 3612-3959 - Viçosa - MG - Brazil
E-mail: rarvore@sif.org.br