OPTIMIZATION OF LAMINATED COMPOSITE PLATES AND SHELLS USING GENETIC ALGORITHMS, NEURAL NETWORKS AND FINITE ELEMENTS

. Structural optimization using computational tools has become a major research ﬁeld in recent years. Methods commonly used in structural analysis and optimization may demand considerable computational cost, depending on the problem complexity. Therefore, many techniques have been evaluated in order to diminish such impact. Among these various techniques, artiﬁcial neural networks may be considered as one of the main alternatives, when combined with classic analysis and optimization methods, to reduce the computational effort without affecting the ﬁnal solution quality. Use of laminated composite structures has been continuously growing in the last decades due to the excellent mechanical properties and low weight characterizing these materials. Taken into account the increasing scientiﬁc effort in the different topics of this area, the aim of the present work is the formulation and implementation of a computational code to optimize manufactured complex laminated structures with a relatively low computational cost by combining the Finite Element Method (FEM) for structural analysis, Genetic Algorithms (GA) for structural optimization and Artiﬁcial Neural Networks (ANN) to approximate the ﬁnite element solutions. The modules for linear and geometrically non-linear static ﬁnite element analysis and for optimize laminated composite plates and shells, using GA, were previously implemented. Here, the ﬁnite element module is extended to analyze dynamic responses to optimize problems based in frequencies and modal criteria, and a module with perceptron ANN is added to approximate ﬁnite element analyses. Several examples are presented to show the effectiveness of ANN to approximate solutions obtained using the FEM and to reduce signiﬁcatively the computational cost.


INTRODUCTION
The structural optimization is not a new field.Galileo in his text "Discorsi e Dimostrazioni Matematiche intorno a due Nuove Scienze" (1638) studied the problem which consists of finding the shape of a beam where every transversal section has the same stress distribution.
In composite materials structures, some experiences has shown that Genetic Algorithms (GA) perform better than traditional gradient based techniques due to the discrete nature of the design variables.
A problem that arises when GAs are used is the high computational cost demanded by this method.For this reason some techniques to reduce this cost have been tested.One of them, consists of replacing the complete Finite Element Analysis (FEA) by some approximation technique.The Artificial Neural Network (ANN) have been shown to be a good alternative to avoid the large number of FEA involved in a GA.
In this work these two techniques are combined to make the process faster and cheaper in terms of computational cost.
This work is based on [2], from which some GA parameters, objective functions and results are taken in order to compare the effectiveness of substituting a complete FEA by ANN.

STRUCTURAL OPTIMIZATION
The structural optimization can be understand as a process to search the configuration which gives the best performance, within some criteria and subjected to certain design constraints.
To model the structural optimization as a mathematical optimization problem, the following concepts are used: • Design variables: they are the characteristics that can be modified by the mathematical optimization algorithm to obtain the best structural performance.
• Design constrains: are the restrictions applied to the structure, such as a limit to avoid material failure, maximum or minimum value of design variables and others that depend on the problem being analyzed.
• Objective function: it is a mathematical expression in which the design variables and constrains are involved.This function represents a number which has to be maximized or minimized during the optimization process.

Structural analysis
The analysis of the composite structures is carried out using the FEM.The element used is a triangular flat plate bending element with 18 degree of freedom called DKT (Discrete Kirchhoff Triangle) combined with the CST (Constant Strain Triangle) to take into account membrane effects.This element was developed by [3] for isotropic materials and it was extended by [1] for laminated composite materials.
For the modal analysis the mass matrix is obtained using the formulation given by [8].
Latin American Journal of Solids and Structures 8(2011) 413 -427 To solve geometrically non linear problems, the generalized displacement control method (GDCM) described by [12] is used.
The Tsai-Wu failure criterion is employed for failure prediction in a ply see [4].

Genetic algorithms for composite materials
Genetic Algorithm (GA) is a computational search tool based on concepts of natural selection and survival of the fittest individual.One aspect of fundamental importance in GAs is the way the solutions are tracked.Instead of using derivatives or gradients, as in deterministic optimization algorithms, GAs work with the objective function based on simple values of individuals.This feature makes the method suitable for problems involving discontinuous functions, and/or non-defined derivatives like in integer programming.Moreover, unlike deterministic optimization methods, which perform the search focusing on a single solution at a time, the GAs work with a population of individuals in each generation.Thus, as several search points are maintained, the convergence or stagnation to local minima, if the starting point is poorly chosen, is prevented.All these aspects result in more chances of finding the optimal solution, even on problems having hard search spaces with multiple local minimum [5].
The design of the optimal sequence of layers in laminated composite materials is a problem of global minimum.Due to the stochastic characteristics of GAs, they are more suitable to optimize than deterministic methods of optimization, which often converge to solutions representing a local minimum.Moreover, in commercial designs, fiber orientation angles and the amount and thickness of layers are discrete variables, a fact which confirms the suitability of GAs for these kinds of problems.
More details related to the use of the method for optimization of composite structures can be found in [2], [9] and [10].
The GA approach adopted here was proposed by [11] and extended by [1].The classical GA is modified to manage composite materials structures.Two genes are used, one for materials and another for the orientation of reinforcement fibers on each ply.The genetic operations used here are crossover, mutation and gene swap.The selection scheme adopted is the Multiple Elitist 1 proposed by [11].In this scheme some of the best individuals from the parents and sons generations are used in the son generation to maintain and to improve the evolution.N e is defined as the number of individuals of each generation.The criterion to stop the GA is defined by two parameters.The first one is the number that limits the number of generations (N LG ), which is used to limit the total number of generations in the GA.The second one, is the number of generation with the same optimum design (N SD ), which indicates that the convergence rate has been zero for a defined number of generations.

ARTIFICIAL NEURAL NETWORKS
Artificial neural networks are mathematical models imitating the human brain.Based on the experience, the ANN can learn complex relationships between inputs and outputs data.The use of ANN to predict structural responses has been growing in recent years.To continue this research line here, neural networks will be applied to learn the structural response of laminated composite structures.The neural network model used here is briefly described below.More details may be found in [6].

Multilayer perceptrons
A general Multilayer Perceptron Network (MPN) architecture consists on a layered network fully connected, i.e., all neurons belonging to a layer are, each one, connected to the previous and the next layer.Such architecture representation is indicated by the number of input vector followed by the number of neurons on each layer and the output vector (e.g.1:12:1 indicates an Artificial Neural Network with one input vector, followed by one hidden layer with 12 neurons and an output layer with one neuron).
The ANN Architecture depends on the complexity of the system being modeled.For a more complex system, a larger number of hidden units are needed.In the training process of neural networks, for an input pattern X pi (where index p means "pattern" and index i means an input neuron), the weight (w pi ) adjustments will take place in the links of the neural network in order to get a desired output O pk close to the function which are being fitted Y po (o means an output neuron).For each input-output pattern, the square of the error and the average system error should be minimized.They may be written, respectively, as where k is the number of neurons in the output layer, and p is the number of training patterns.Employing any algorithm to minimize the error function, the weights can be evaluated and an approximated fit may be obtained.The standard back-propagation algorithm is used to adjust the different weights as well as the derivatives of O pk with respect to the input data which will be fitted.After this first adjustment takes place, the network will pick up another pair of X pi and Y po , and will again adjust weights for this new pair.In a similar way, the process will go on until all the input-output pairs are considered.After some training epochs, the network will have a single set of stabilized weights satisfying all the input-output pairs with an average system error lower than a tolerance (10 −3 ).More details regarding the use of the MPN can be found in [7] and [6].

Artificial neural networks and genetic algorithms
To combine the ANN and GA, the scheme shown in Fig. 1 is adopted.To generate the training set for the ANN in each case, three GA are executed with two generations having large populations each one.This approach is used to make a randomly distributed generation in the first one and some elitist generation in the second one.Three executions of the GA make the generations disperse and they are enough to guarantee that almost every zone of the design space is covered.
So, when the ANNs are trained with the training set generated by the previously described process, every structural response is taken from the ANN.In the case of parameters that can be directly calculated from the genetic codification, such as cost, and weight, for example, they are calculated directly.

NUMERICAL EXAMPLES AND DISCUSSION
To evaluate the quality of the optimization scheme and the time that can be saved using ANN, three examples are shown.In each one, the GA and the training time as well as the error of the ANN (with respect to the Finite Element Method (FEM) solution) are presented, together with a comparison of GA-FEM with respect to processing time and quality of results are also presented.All examples have been executed in a 2.4 Ghz Quad Core 2 computer with 4 Gb of RAM.

Cost and weight minimization of an in-plane loaded composite laminated plate
In this problem the number of plies, the material of each ply and the orientation of reinforcement fibers are the design variables.The right combination of these variables can determine the cheapest and lightest structure.The constraints of the problem are the material failure (λ f ) derived from the Tsai-Wu criterion [4], and the structural elastic stability (λ b ).Both must be greater or equal to 1.0.The cost is proportional to the material consumption in the laminated structure, and each material has its own cost per unit weight denoted by C. The plate model with boundary and load conditions is shown in Fig. 2. The finite element mesh has 3000 elements.The objective function is defined by Eq. 1 , where W t * and C * are the normalized weight and cost respectively; this normalization is shown in Eq. 2 and ϕ is the weighting factor of both objectives; in this example the weighting factor is 0.5.
The elastic constants, strength parameters, specific weight, ply thickness and the parameter of cost per unit weight for the Kevlar-epoxy and Graphite-epoxy are shown in Table 1.The elastic constants are the Young's modulus in the fiber direction (E 1 ) and transverse to the fiber direction (E 2 ), the shear modulus (G 12 ) and the Poisson's ratio (ν 12 ), respectively.Strength parameters for tension and compression for longitudinal and transversal directions are given by F 1t , F 1c , F 2t , and F 2c , respectively.The remainder parameters are the shear strength (F 6 ), the specific weight (ρ) and the thickness (t).
The alphabet used in the GA is shown in Table 2; the code 0 is used to indicate an empty layer.The laminated plate can adopt different number of plies varying from 12 to 24 and the symmetry condition is applied.The size of the design space (SDS) is 55944.
The parameters of the GA are shown in Table 3 where P is the size of population, N e is the elitist scheme parameter, p ma and p mm are the probabilities of angle and material mutation respectively, p gs is the probability of gene swap, p pa and p pd are the probabilities of ply addition and ply deletion, respectively, N LG and N SD are stop criterion parameters.
The training set consists in 2347 designs and the time used to generate this training set was 112.74 hours.The neural network architecture used here is 12:18:18:1.A neural network is used to approximate (λ f ) and another one is used to approximate (λ b ).The cost and the weight are calculated directly from the genetic code of each design.The training times and  the error (with respect to the FEM solution) using neural networks are shown in Table 4.In this table "uc/N" means unit cost per newton and the "time" is referred to the time spent whith trained ANN (time spent in the training process was not considered) The optimal design found using GA-ANN and values of the parameters are shown in Table 5.
Performing a FE analysis, using the design obtained from the GA-ANN optimization, the real differences between parameters λ f and λ b are shown in Table 6.To verify the optimization quality, the design obtained with GA-ANN is compared to the design found using GA-FEM.This last design and parameters values are shown in Table 7. Differences between the designs obtained with GA-FEM and GA-ANN are shown in Table 8.These results show that the design obtained with the GA-ANN is a near optimum design.
To evaluate the processing time saved using ANN, a comparison is shown in Table 9.
In fact this problem can be solved by a mesh having much less elements, but the aim of this example is to show that major gains in processing time occurs when the FE analysis is time consuming.In this case the difference in processing time is 57.37%.

Stiffness maximization of a composite laminated shell with geometrically nonlinear behavior
This optimization problem aims to maximize the stiffness of a composite shallow shell under pressure load.Figure 3 shows the shell with its boundary and load conditions.The mesh has 800 elements.The nonlinear analysis is made using the GDCM method [12], with a load increment parameter λ i = 0.05.The fitness function is defined in Eq. 3, where N C crit is the critical load level(when curve load -displacement reaches its first limit point), U max is the maximum displacement, which is taken at the end of incremental load process or when material failure is observed, N C max is the maximum load level without material failure and V nlc is a penalization for designs having more than 4 plies with same fiber orientation.
In the example glass-epoxy is used and the material properties are shown in Table 10.The genetic alphabet shown in Table 11 is used.In this example only the orientation of reinforcement fibers are the design variables, while the material, thickness and number of layers are fixed.The laminated has 28 plies (with thickness t = 0.45mm for each ply), the symmetry condition is imposed, each gene control 2 plies and the chromosome has 7 genes to control the laminate.Due to the long codification, for the crossover is used a double break point.The size of the design space (SDS) is 2187.The GA parameters used here are P = 18, N e = 3, p ma = 5%, p mm 0%, p gs = 80%, p pa = 0%, p pd = 0%, N LG = 108, N SD = 36 with the same meaning as in the previous example.
In this example, three neural networks are used to approximate each one of the parameters involved in the objective function (N C max , N C crit and U max ); each ANN for the trained parameters, errors and processing time are show in Table 12.The architecture adopted is 7:15:15:1.The training set has 226 designs and the time spent to generate this training set was 117.31 minutes.To evaluate the quality of the approximation using neural networks, the parameter values are compared to the parameter values for the same design using FEM without ANN.These differences are shown in Table 14.Result and processing time used in the optimization process using GA-FEM is show in Table 15.To verify the quality of the optimization process using GA-ANN, parameter values for both designs, obtained with ANN and FEM, are compared in Table 16.A graphical comparison of both designs is shown in Fig. 4, where N C crit and U max are indicated.
Processing time comparison, using GA-ANN and GA-FEM, is shown in Table 17.
In this example the time saved using GA-ANN is 50.79%; to evaluate the optimization quality, the whole design space analysis was used, and in this space the design obtained by the GA-ANN is the 9th.among the near optimum designs (this set includes the optimum design).used to build this plate, and the material properties are given in Table 18.The design obtained here, as well as ω, are shown in Table 21 The difference of ω, obtained using ANN and FEM, is shown in Table 22.Making a GA-FEM optimization, the same design is found.A comparison of the processing time required for each method is shown in Table 23.In this example the processing time saved is 59%.The optimization quality using ANN is high.This occurs because the function that has to be replaced by the neural network is very simple.

FINAL REMARKS
Some important remarks can be outlined from the examples, such as following: • An important amount of processing time can be saved using GA-ANN instead of GA-Latin American Journal of Solids and Structures 8(2011) 413 -427 FEM.This is particularly true when Finite Elements Analyses are applied for structures requiring refined meshes and for cases where the structure has a nonlinear behavior.
• The optimum design is very little affected if well trained ANNs are used substituting complete Finite Element Analyses, or, in other words, accuracy of ANNs can be improved with a longer and more elitist training set and increasing the number of GA application to generate samples of the training set.
• In cases where the cost function is simple, final design using GA-ANN will be very similar to that obtained using GA-FEM.
In future works other kind of ANNs, such as ANN with Radial Basis, (substituting Multilayer Perceptrons) will be tested.These tools (GAs and ANNs) may be employed also in Reliability Based Optimization Problems.

Figure 1
Figure 1 Flowchart of GA-ANN.

Figure 2
Figure 2 Composite laminated plate with its boundary and load conditions.

Figure 3
Figure 3 Composite laminated shell whit boundary and load conditions.

Figure 5
Figure 5 Composite laminated plate with its boundary conditions and the FE mesh.

Table 2
Genetic codification alphabet and possible values -Example 1.

Table 4
Neural networks training time and errors -Example 1.

Table 6
Differences between ANN and FEM -Example 1.

Table 9
Processing time comparison (in hours) -Example 1.

Table 12 Neural
Networks training time and errors -Example 2.

Table 14
Differences between ANN and FEM -Example 2.

Natural frequency maximization of a laminated plate
The objective of this optimization is to maximize the first natural frequency of vibration (ω) of a simply supported composite laminated plate.No constraints are imposed to the problem.To obtain ω the eingenvalues problem is solved.The GA modifies only the reinforcement fibers orientation in each ply.The material and number of plies are fixed.The number of plies is fixed in 8, with thickness equal to 2 mm each one.Then the total height is 16 mm.Figure5shows the plate, boundary conditions and the mesh (with 256 elements).Graphite-epoxy is Cardozo et al / Optimization of plates and shells using genetic algorithms, neural networks and finite elements 425

Table 18
Properties of Graphite-epoxy.The genetic alphabet used in this example is shown in Table19.The chromosome length is 8. Double break point is used in the crossover operation.The size of the design space (SDS) is 6561.

Table 19
Genetic codification alphabet -Example 3.The GA parameters used here are P = 30, N e = 3, p ma = 5%, p mm 0%, p gs = 80%, p pa = 0%, p pd = 0%, N LG = 200, N SD = 100, with the same meaning as in the first example.The neural network architecture used in this example is 8:17:17:1.The training time and error are shown in Table20.The training set has 420 designs and the time to generate the training set was 26.75 minutes.

Table 20 Neural
Network training time and error for Example 3.Table 21 Result using GA-ANN -Example 3.
Table 22 Difference between ANN and FEM -Example 3.

Table 23
Processing time comparison (in minutes) -Example 3.