Reinforced Concrete Slab Optimization with Simulated Annealing

Featured Application: An optimization method, based on Simulated Annealing, is developed for the design of reinforced concrete slab. Abstract: Flat slabs have several advantages such as a reduced and simpler formwork, versatility, and easier space partitioning, thus making them an economical and e ﬃ cient structural system. When producing structural components in series, every detail can lead to signiﬁcant cost di ﬀ erences. In these cases, structural optimization is of paramount relevance. This paper reports on the structural optimization of reinforced concrete slabs, presenting the case of a rectangular slab with two clamped adjacent edges and two simply supported edges. Using the yield lines method and the principle of virtual work, a cost function can be formulated and optimized using simulated annealing (SA). Thus, the optimal distribution of reinforcing bars and slab thickness can be found considering the ﬂexural ultimate limit state and market materials costs. The optimum result was deﬁned by the orthotropic coe ﬃ cient k = 8, anisotropic coe ﬃ cient g = 1.4, and slab thickness H = 11.8 cm. A sensitivity analysis of the solution was developed considering di ﬀ erent material costs.


Introduction
The earliest studies on structural optimization can be dated back to 1500-1600 when, first Leonardo da Vinci, and later Galileo, carried out experimental tests on scale models specifically designed to enhance certain characteristics [1]. Maxwell's early work in 1869 [2] is commonly considered as the forefather of classic structural optimization, while the first computer-based optimization study was carried out by Schmit in 1960 [3]. Following the improvement of computational power, there was a great development in numerical methods applied to structural optimization.
In the 1980s, important heuristic algorithms were developed, thus establishing the benchmark for structural optimization in the following years. A representative example is simulated annealing [4], which uses numerical techniques to simulate the processes of evolution and selection. The specific language of the procedure is borrowed from the Darwinian evolutionary system: a possible solution to a problem is called an "individual", while a set of solutions is called a "population".
Over the years, these four categories have been merged and updated, giving rise to many hybrid optimization algorithms that share an evolutionary origin. Harmony search was conceptualized by using the musical process of searching for a perfect state of harmony [9]. Evolutionary structural optimization (ESO) algorithms are based on the idea that the optimal structure can be reached by gradually removing ineffectively used materials from the design domain (see [10]). In [11], genetic algorithms (GAs) were integrated with ESO to obtain genetic evolutionary structural optimization (GESO), which uses GA for global optimum searches, based on the ESO approach. The constrained evolution method searches for the optimum design as a single-objective optimization problem, therefore enforcing various constraints. Interesting applications can be found in [12] for pre-stressed concrete beams and in [13] for force-limiting floor anchorage systems. Another example is the particle swarm optimization algorithm, which was developed in the last years of the 1990s. A note-worthy application of this algorithm to structural optimization can be found in [14], while in [15], it was used for the identification of Van der Pol-Duffing non-linear oscillators.
Structural optimization problems can be described with different approaches. Often, the objective of the analysis is to obtain a construction with the minimum cost, or minimum mass. The well-known paper by Sarma and Adeli [16] highlighted that the construction of concrete structures involved at least three different materials: concrete, steel, and the formwork. Thus, depending on the scope, the design of these structures can be based on cost rather than weight minimization. This aspect is particularly relevant in structural retrofitting. Examples dealing with reinforced concrete structures under seismic loads are reported in [17,18], while masonry structure cases are reviewed in [19][20][21].
The research aimed at reducing vibrations and improving the dynamic structural behavior of reinforced concrete structures based on optimization algorithms deserves special attention. In [22], the slab track optimization problem was discussed with a review of the available solutions. A genetic algorithm was employed for the multi-objective optimization of structural passive control systems in [23], while the optimization of a concrete cable stayed bridge under seismic action was discussed in [24]. The optimization of a tuned slab dumper for a metro system was presented in [25]. The optimization of vibration control is of particular relevance in the case of tall buildings (see [26] for a comprehensive analysis), but also in the case of concrete slabs [27,28]).
Structural systems with reinforced concrete slabs are a common structural solution. Their structural behavior is complex and is currently being studied. Fatigue [29], punching shear [30][31][32] resistance, aging [33], and cracking properties [34,35] have been discussed in the current literature.
Flat slabs have several advantages such as a reduced and simple formwork as well as versatility and easy space partitioning, thus making them an economical and efficient structural system. Slab optimization has been thoroughly discussed in the literature. The case of flat slabs with an arbitrary configuration is presented in [36], where materials and construction costs were taken into account. The yield lines method is usually adopted to assess the slab collapse load with low computational cost. The discontinuity layout optimization (DLO) procedure [37] can be used to systematically automate the method to be applied to a wide range of practical slab problems. A hybrid genetic algorithm was proposed for slab optimization in [38]. A very interesting example of design optimization of simply supported slabs using finite elements was described in [39].
An emerging trend in this kind of problem is the analysis of the whole life-cycle of the structure, and consequently the optimization of the total costs of the building from construction to demolition. Various papers have been published on this specific topic in recent years. In particular, it is worth highlighting the results of the work dealing with flat slabs as presented in [40] as well as the sensitivity analysis of the life-cycle assessment in [41], which proves how transport and material waste are significant parameters. CO 2 emissions were the focus of [42,43], which dealt with slab life-cycle optimization. More recently, in [44], deep learning neural networks were used to optimize slab design. Finally, in [45,46], criteria to optimize structural and energy retrofitting in order to optimize greenhouse gas emissions and economic cost were discussed.
Simulated annealing (SA) applications to slab optimization are not common in the literature. In this work, the cost optimization of rectangular anisotropic RC slabs was developed in an innovative way using limit analysis, based on yield lines theory and the SA algorithm.
The main target of this research was to demonstrate the effectiveness of this new optimization method for a flat slab with a given set of boundary conditions. The rest of this paper is organized as follows. Section 2 presents the structural model based on limit analysis and yield lines theory and the equations that will be the target of the optimization. The main characteristics of the SA algorithm are presented in Section 3. The results of the specific case study are shown in Section 4, which also reports a comparison with finite element analysis. Finally, conclusive remarks are made in Section 5.

Problem Definition
In general, structural optimization problems can always be expressed using an objective function F(x 1 , x 2 , . . . , x n ) and constraints g n (x 1 , x 2 , . . . , x n ), which depend on the characteristic parameters of the problem.
where the objective function F(x 1 , x 2 , . . . , x n ) is expressed as the sum of the products between the unit prices p i and the corresponding quantities m i (x 1 , x 2 , . . . , x n ) of the construction components.
The constraint system represents the conditions necessary to satisfy the limit states expressed through the problem variables.

Yield Line Method
The yield line method makes it possible to model the structural behavior of a structure in the ultimate limit state (ULS). To better explain the physics of the problem, it is useful to describe the limit behavior of reinforced concrete (RC) slabs: under increasing load, the structure has linear elastic behavior, followed by non-linear plastic behavior characterized by tensile reinforcement yielding and plasticization of the cross section. In this condition, a redistribution of internal stresses is expected and reinforcement yielding is spread over a wide set of points until a collapse mechanism is reached. Such behavior can be modeled with some simplifications:

•
The materials present a plastic constitutive law. The yield line method can be divided into two distinct steps: (1) Taking into account the geometrical, constraint and load characteristics of the slab, one or more families of possible failure mechanisms are determined. Each family is defined by a certain number of geometric parameters that can express the ultimate load as a scalar function of the same parameters p = p(x 1 , x 2 , . . . , x n ). (2) Enforcing the kinematic limit analysis theorem, the failure mechanism corresponding to the smallest limit load is identified for each mechanism family. The smallest value among the loads of each family will be the collapse load of the slab. The core of the method is the expression of the collapse load as a function of the mechanism parameters p = p(x 1 , x 2 , . . . , x n ). Typically, this expression is non-linear and its minimization is a complex operation.

Internal Work
Consider a rectangular plate, whose dimensions are L x and L y , with both top and bottom reinforcements arranged in the direction of the axes of the reference system x, y. A x + represents the area, per unit length, of the reinforcement located in the bottom part of the section parallel to the x-axis, while A x − denotes the area, per unit length, of the top reinforcements in the x-direction.
Similarly, in the y direction, A y + and A y + represent the area, per unit length, of the bottom and top reinforcements, respectively. It is necessary to introduce the following coefficients of k, characterizing the orthotropy, and g, the anisotropy of the slab: The ultimate limit bending moment of the cross section can be expressed as a function of its geometrical and mechanical characteristics. In a conservative approach, the role of the compressive reinforcements is neglected, thus, considering the bottom reinforcements along the x-axis area, A + x , is characterized by a tensile strength, f yd , it is possible to define the ultimate bending moment M + x : where z represents the distance between the tensile force and the compressive force.
Considering that RC slabs have relatively low levels of reinforcement compared to beams and the values of z are almost constant, the ultimate bending moments can be considered proportional to the reinforcement areas. This approximation makes it possible to express the other bending moments as a function of the M + x as defined in Equation (4): where k is the orthotropic coefficient and g denotes the anisotropic coefficient. Given the k-th yield line characterized by length ∆l parallel to the x-axis, the internal work D i can be expressed by Equation (8).
where M n is the ultimate bending moment per unit length corresponding to the yield line while θ is the relative rotation between the rigid parts separated by the given yield line. Considering a generic yield line, characterized by an inclination angle ϕ with respect to the x-axis, the ultimate bending moments corresponding to the top M − p and to the bottom M + p reinforcements can be expressed following Johansen's criterion [47] (see Figure 1a).

Yield Line Internal Work
By applying Equations (9) and (10), it is possible to have a more general expression of the internal work corresponding to a generic yield line, whose inclination with respect to the x-axis is : where the characteristic length of the slab , = Δ / , and = • are given. Thus, the internal work of the whole set of n yield lines is:

Fan Internal Work
In orthotropic RC slabs, the internal work dissipated in the j-th fan can be expressed with Equation (15) (see Figure 1b).
where is the virtual displacement of the fan vertex; is the angle at the fan center, while and are the angles between the two outer fan rays and the x axis. Given = • , the internal work dissipated in the whole fan set is:

External Work
The external load work can be obtained by integrating the product between the specific force per surface unit and the corresponding displacements over the whole slab surface A: In this case, the principle of virtual work is defined by Equation (11).
where W i is the internal work corresponding to the yield lines; W ϕ is the fan internal work; and W e is the external work.

Yield Line Internal Work
By applying Equations (9) and (10), it is possible to have a more general expression of the internal work corresponding to a generic yield line, whose inclination with respect to the x-axis is ϕ: where the characteristic length of the slab L x , l = ∆l/L x , and w k = D i L x ·M + x are given. Thus, the internal work of the whole set of n yield lines is:

Fan Internal Work
In orthotropic RC slabs, the internal work dissipated in the j-th fan can be expressed with Equation (15) (see Figure 1b).
where δ vj is the virtual displacement of the fan vertex; Ψ j is the angle at the fan center, while α j and β j are the angles between the two outer fan rays and the x axis. Given w ϕ j = , the internal work dissipated in the whole fan set is:

External Work
The external load work can be obtained by integrating the product between the specific force per surface unit and the corresponding displacements over the whole slab surface A: where p(x, y) = P(x, y)L 2 x / 6 × M + x is the non-dimensional surface load and v(x, y) = V(x, y)/L x is the non-dimensional displacements for each slab point.
In the case of constant surface load p, the external load work becomes: where

Reinforcement Arrangement Optimization
Now, by applying Equations (10), (14), (16)-(18), it is possible to obtain the ultimate slab load as a function of the selected collapse mechanism: Actually, the problem can be formulated in order to find the minimum of the ratio between the non-dimensional reinforcements volume v a = V s /(A + x L x ) and the non-dimensional slab collapse load p (see Equation (19)). In this way, it is possible to find the optimal reinforcement arrangement. Indeed, the optimization of Equation (20) Equation (20) is typically a very complex non-linear expression. It makes it possible to express the ultimate load of the slab as a function of the geometric parameters of the considered collapse mechanism.
Applying the kinematic theorem of the limit analysis gives the definition of the slab collapse load: it is the load corresponding to the combination of parameters that minimizes the value of Equation (20). This optimization was developed using the simulated annealing algorithm described in the next paragraph.

Slab Cost Optimization
In order to optimize the total cost of the slab, it is necessary to take into account the different prices of concrete p c and steel p s . For the sake of simplicity, the cost of formworks, casting operations etc. will be included in these prices. Equation (21) presents the cost equation given the volume of concrete V c and steel V s .
In the following, a slab cross section with unitary width will be considered. The materials' constitutive laws are presented in Figure 2, showing the concrete stress-block model on the left and the steel ideal elasto-plastic model on the right. In the same figure, y n denotes the neutral axis depth, f cd is the concrete design compressive strength, and f yd is the steel design tensile strength. Now, by applying the equilibrium conditions of the slab cross section, it is possible to express the total cost as a function of the prices of materials and of the bottom reinforcements area . Indeed, the slab thickness, H, can be obtained by considering a certain concrete cover c, bar diameter , and applying equilibrium equations given the material mechanical characteristics and the reinforcements: Thus, the total cost can be obtained by applying the definition of H in Equation (22) into Equation (21). The latter equation is a function of given the value of k and g, obtained by the optimization of Equation (20).
In order to have sufficient ductility for the collapse mechanism development, the ratio between the neutral axis, yn, and the effective depth, d, should be limited [48]: Equation (23) has been minimized by means of the simulated annealing (SA) algorithm. Thus, the optimization procedure can be synthetized as follows: 1. Input of geometrical characteristics and boundary conditions. 2. Expression of the collapse mechanism and definition of collapse load (Equation (20)) as a function of the mechanism parameters: ( , , . . , ). 3. Optimization of ( , , . . , ) varying anystropic and orthotropy coefficients (g, k) in order to find the reinforcement distribution that minimizes Equation (20)    Now, by applying the equilibrium conditions of the slab cross section, it is possible to express the total cost as a function of the prices of materials and of the bottom reinforcements area A + x . Indeed, the slab thickness, H, can be obtained by considering a certain concrete cover c, bar diameter φ, and applying equilibrium equations given the material mechanical characteristics and the reinforcements: Thus, the total cost can be obtained by applying the definition of H in Equation (22) into Equation (21). The latter equation is a function of A + x given the value of k and g, obtained by the optimization of Equation (20).
In order to have sufficient ductility for the collapse mechanism development, the ratio between the neutral axis, y n , and the effective depth, d, should be limited [48]: Equation (23) has been minimized by means of the simulated annealing (SA) algorithm. Thus, the optimization procedure can be synthetized as follows: 1. Input of geometrical characteristics and boundary conditions. 2.
Expression of the collapse mechanism and definition of collapse load (Equation (20)) as a function of the mechanism parameters: p(x 1 , x 2 , . . . , x n ).

3.
Optimization of p(x 1 , x 2 , . . . , x n ) varying anystropic and orthotropy coefficients (g, k) in order to find the reinforcement distribution that minimizes Equation (20) by means of the SA algorithm.

4.
Optimization of slab thickness (Equation (23)) given the reinforcement distribution by means of the SA algorithm.

Simulated Annealing
The simulated annealing algorithm is often defined as a general technique for solving combinatorial optimization problems. It is a heuristic algorithm that finds the solution through a stochastic and iterative procedure. In its original form, the SA was proposed by Kirkpatrik, Gelatt, and Vecchi [4]. The word "annealing" denotes a physical process where a solid is progressively heated until fusion and then cooled down slowly. If a sufficiently high temperature has been reached and the cooling process has been slow, the molecules are arranged in the minimum potential energy configuration.
In this work, a MATLAB™ version of the SA algorithm was adopted, see [49].
The input values of the algorithm parameters are: • c 0 = 1 initial value of the control parameter; • c min = 10 −8 final value of the control parameter; and • k B = 1 Boltzmann constant.
The maximum number of iterations allowed before decreasing the control parameter is fixed at 150. The decrement procedure of the control parameter is expressed by Equation (25):

Test Case
Consider the rectangular slab shown in Figure 3 with two simply supported edges and two clamped edges subjected to a uniform surface load. The collapse mechanism that produces the lowest collapse load is depicted in Figure 3 and is characterized by three fans and five geometrical parameters x 1 , x 2 , x 3 , x 4 , and x 5 . For the sake of clarity, these parameters are non-dimensional because all lengths refer to the largest slab size, L x .

Simulated Annealing
The simulated annealing algorithm is often defined as a general technique for solving combinatorial optimization problems. It is a heuristic algorithm that finds the solution through a stochastic and iterative procedure. In its original form, the SA was proposed by Kirkpatrik, Gelatt, and Vecchi [4]. The word "annealing" denotes a physical process where a solid is progressively heated until fusion and then cooled down slowly. If a sufficiently high temperature has been reached and the cooling process has been slow, the molecules are arranged in the minimum potential energy configuration.
In this work, a MATLAB™ version of the SA algorithm was adopted, see [49].
The input values of the algorithm parameters are: • c0 = 1 initial value of the control parameter; • cmin = 10 −8 final value of the control parameter; and • kB = 1 Boltzmann constant.
The maximum number of iterations allowed before decreasing the control parameter is fixed at 150. The decrement procedure of the control parameter is expressed by Equation (25):

Test Case
Consider the rectangular slab shown in Figure 3 with two simply supported edges and two clamped edges subjected to a uniform surface load. The collapse mechanism that produces the lowest collapse load is depicted in Figure 3 and is characterized by three fans and five geometrical parameters , , , , and . For the sake of clarity, these parameters are non-dimensional because all lengths refer to the largest slab size, Lx. The reinforcement distribution is presented in Figure 4. The bottom reinforcements were uniformly distributed over the slab surface while the top ones were compliant with the collapse mechanism shown in Figure 4.
The material characteristics and prices are shown in Table 1. The reinforcement distribution is presented in Figure 4. The bottom reinforcements were uniformly distributed over the slab surface while the top ones were compliant with the collapse mechanism shown in Figure 4.
The material characteristics and prices are shown in Table 1.

Reinforcement Arrangement Optimization
The non-dimensional parameter, , defined in Equation (20) was minimized for various combinations of k (from 0 to 10) and g (from 1 to 2) using the SA algorithm.
The adopted constraints enforce the geometrical limits that characterize the mechanism: The values are presented in Figure 5 where a red dot indicates the minimum value: = 0.294 corresponding to k = 8 and g = 1.4. The collapse mechanism parameter values for this case are shown in Table 2. It is important to highlight that this result was obtained for a non-dimensional problem. In order to find the optimal slab thickness, it is necessary to introduce material characteristics, load conditions, and material prices.

Reinforcement Arrangement Optimization
The non-dimensional parameter, ρ, defined in Equation (20) was minimized for various combinations of k (from 0 to 10) and g (from 1 to 2) using the SA algorithm. The adopted constraints enforce the geometrical limits that characterize the mechanism: The ρ values are presented in Figure 5 where a red dot indicates the minimum value: ρ = 0.294 corresponding to k = 8 and g = 1.4. The collapse mechanism parameter values for this case are shown in Table 2. It is important to highlight that this result was obtained for a non-dimensional problem. In order to find the optimal slab thickness, it is necessary to introduce material characteristics, load conditions, and material prices.

Reinforcement Arrangement Optimization
The non-dimensional parameter, , defined in Equation (20) was minimized for various combinations of k (from 0 to 10) and g (from 1 to 2) using the SA algorithm. The adopted constraints enforce the geometrical limits that characterize the mechanism: The values are presented in Figure 5 where a red dot indicates the minimum value: = 0.294 corresponding to k = 8 and g = 1.4. The collapse mechanism parameter values for this case are shown in Table 2. It is important to highlight that this result was obtained for a non-dimensional problem. In order to find the optimal slab thickness, it is necessary to introduce material characteristics, load conditions, and material prices.

Slab Thickness Optimization
The optimal reinforcement arrangement obtained in Section 4.2 (k = 8 and g = 1.4) was the starting point for the slab thickness optimization. Indeed, these values of orthotropic and anisotropic coefficients represent a constraint in this section.
The slab thickness that minimizes the slab material cost and satisfies the ultimate limit state can be found considering Lx = 7 m and Ly = 4.9 m. The service load applied to the slab surface was 3 kN/m 2 . The material characteristics and prices are shown in Table 1, where the concrete and reinforcement steel prices came from the current Italian construction market.
The optimization of Equation (21) gives A + x = 38 mm 2 /m and H = 11.8 cm. The authors studied the stability of the solution for different material prices and found no significant differences.

Finite Element Method Comparison
In order to test the accuracy of the obtained results, the authors developed a nonlinear finite element in ANSYS™. Concrete was modeled with 900 SOLID65 brick elements while reinforcements were represented by 830 LINK180 elements.
The reinforcement distribution, slab thickness, material characteristics, and models were equal to those adopted for the test case.
A nonlinear static analysis was obtained and the results are presented in Figure 6, which shows the FE load-displacement curve. The results are non-dimensionalized with respect to the SA algorithm load. It can be seen how the ratio between the SA ultimate load and the FE element load was similar to one, and how the collapse of the slab corresponded to a mid-span displacement of about 0.06 m.

Slab Thickness Optimization
The optimal reinforcement arrangement obtained in Section 5.2 (k = 8 and g = 1.4) was the starting point for the slab thickness optimization. Indeed, these values of orthotropic and anisotropic coefficients represent a constraint in this section.
The slab thickness that minimizes the slab material cost and satisfies the ultimate limit state can be found considering Lx = 7 m and Ly = 4.9 m. The service load applied to the slab surface was 3 kN/m 2 . The material characteristics and prices are shown in Table 1, where the concrete and reinforcement steel prices came from the current Italian construction market.
The optimization of Equation (21) gives = 38 mm 2 /m and H = 11.8 cm. The authors studied the stability of the solution for different material prices and found no significant differences.

Finite Element Method Comparison
In order to test the accuracy of the obtained results, the authors developed a nonlinear finite element in ANSYS™. Concrete was modeled with 900 SOLID65 brick elements while reinforcements were represented by 830 LINK180 elements.
The reinforcement distribution, slab thickness, material characteristics, and models were equal to those adopted for the test case.
A nonlinear static analysis was obtained and the results are presented in Figure 6, which shows the FE load-displacement curve. The results are non-dimensionalized with respect to the SA algorithm load. It can be seen how the ratio between the SA ultimate load and the FE element load was similar to one, and how the collapse of the slab corresponded to a mid-span displacement of about 0.06 m. Good agreement between the FE and SA results was found considering the crack pattern distribution reported in Figure 7. Indeed, by comparing Figure 7 with Figure 3, it can be seen how the assumed yield lines distribution was very similar to the crack pattern obtained with finite elements in the collapse condition.  Good agreement between the FE and SA results was found considering the crack pattern distribution reported in Figure 7. Indeed, by comparing Figure 7 with Figure 3, it can be seen how the assumed yield lines distribution was very similar to the crack pattern obtained with finite elements in the collapse condition. Appl. Sci. 2019, 9,

Conclusions
In this work, the optimization of RC slabs was developed in an innovative way by using limit analysis and the SA algorithm. The procedure can be divided in two main steps: first, the arrangement of the reinforcements, followed by the optimization of the slab thickness. While in the first step the problem can be fully non-dimensional, in the second, it is necessary to introduce geometrical, mechanical, loading, and economic characteristics.
The case of a rectangular RC slab with two simply supported edges and two clamped edges was discussed and the slab characteristics that minimize the structural cost were found. The collapse mechanism was identified with five geometrical parameters and can be characterized by five yielding lines and three fans. The slab thickness was optimized by using a minimum cost approach.
The limit load analysis approach is reliable if the considered structural element has the ductility necessary to develop the considered mechanism. For this reason, it is important to ensure that the slab characteristics will prevent any brittle failures that would impede the development of the assumed collapse mechanism.
In addition, ultimate limit analysis with yield lines makes it possible to design plates with adequate strength, but it is clear that it does not take into account any serviceability limit states such as deformation or vibration. Thus, a complete design must add serviceability limit state checks to the optimization process. This, however, was beyond the scope of this paper.
It is important to highlight that the kinematic theorem provides an upper limit to the collapse load. For this reason, all the possible mechanisms that can be developed in the structure must be examined, making the presence of the fans indispensable. Furthermore, the comparison between the experimental tests and the theoretical results of the limit analysis with the yield lines method

Conclusions
In this work, the optimization of RC slabs was developed in an innovative way by using limit analysis and the SA algorithm. The procedure can be divided in two main steps: first, the arrangement of the reinforcements, followed by the optimization of the slab thickness. While in the first step the problem can be fully non-dimensional, in the second, it is necessary to introduce geometrical, mechanical, loading, and economic characteristics.
The case of a rectangular RC slab with two simply supported edges and two clamped edges was discussed and the slab characteristics that minimize the structural cost were found. The collapse mechanism was identified with five geometrical parameters and can be characterized by five yielding lines and three fans. The slab thickness was optimized by using a minimum cost approach.
The limit load analysis approach is reliable if the considered structural element has the ductility necessary to develop the considered mechanism. For this reason, it is important to ensure that the slab characteristics will prevent any brittle failures that would impede the development of the assumed collapse mechanism.
In addition, ultimate limit analysis with yield lines makes it possible to design plates with adequate strength, but it is clear that it does not take into account any serviceability limit states such as deformation or vibration. Thus, a complete design must add serviceability limit state checks to the optimization process. This, however, was beyond the scope of this paper.
It is important to highlight that the kinematic theorem provides an upper limit to the collapse load. For this reason, all the possible mechanisms that can be developed in the structure must be examined, making the presence of the fans indispensable. Furthermore, the comparison between the experimental tests and the theoretical results of the limit analysis with the yield lines method confirmed that the approach is safe and well-founded [50,51], and that the comparison between finite elements and limit analysis showed very good agreement.
This method can be easily applied to very different structural systems with different loading and boundary conditions. The main difficulties arise in the definition of the collapse mechanism that produces the lowest collapse load.
The simulated annealing algorithm proved to be efficient and reliable in this context. A comparison with other heuristic optimization algorithms is projected in future developments of this research.
In conclusion, this method represents an effective procedure for the early design of concrete slabs.