Open access

Introductory Chapter: An Introduction to the Optimization of Composite Structures

Written By

Karam Maalawi

Submitted: 06 February 2018 Published: 05 November 2018

DOI: 10.5772/intechopen.81165

From the Edited Volume

Optimum Composite Structures

Edited by Karam Y. Maalawi

Chapter metrics overview

1,361 Chapter Downloads

View Full Metrics

1. Introduction

Structural applications of composite materials are increasing in several engineering areas where high stiffness and strength-to-weight ratios, long fatigue life, superior thermal properties, and corrosive resistance are most beneficial [1, 2, 3, 4]. Common types include laminated composites [5], functionally graded material (FGM) structures, and nanocomposites as well as smart composite structures [6]. In fact composite structures are usually tailored, depending upon the specific objectives, by choosing the individual constituent materials and their volume fractions, fiber orientation angles, and laminas thickness and number, as well as the fabrication procedure. To attain the best results, adequate optimization models have to be implemented to find practical optimal solutions satisfying a given set of design constraints.

This introductory chapter provides a brief review on the optimum design of composite structures and the relevant optimization techniques that are capable of finding the needed optimal solutions. Several problems can be addressed, including the structural design for maximum stability, maximum natural frequencies, and minimum mass or maximum stiffness subject to limits on strength, deflections, and side constraints. The relevant design variables include geometrical dimensions and material properties as well. A numerical example is given at the end of this chapter to demonstrate a real and practical application of the optimum composite structures.

Advertisement

2. The optimal design problem

Several research papers and text books exist in the field of optimal design of composite structures with a variety of valuable applications in civil, mechanical, ocean, and aerospace engineering. An important stage has now been reached at which an investigation of such developments and their practical possibilities should be made and presented. Two distinct review papers have been published covering the development of the optimum design of composites over more than 40 years. The first paper by Sonmez [7] presented a comprehensive survey for more than 1000 journal papers, conference papers, textbooks, and web links from the year 1969 to 2009. Sonmez classified the papers according to the type of the composite structure, loading conditions, optimization model, failure criteria, and the utilized search algorithm. The second paper by Ganguli [8] covered a historical review from 1973 to 2013. It provides the growth of the field by including more than 90 references dealing with a variety of optimization methods utilized for tailoring composites to achieve certain design objectives. Applications of several optimization techniques were presented, including feasible direction methods, sequential quadratic programming, and stochastic optimization such as particle swarm and ant colony algorithms. Ganguli classified the published work into five categories named pioneering research for the work published in the 1970s, early research in the 1980s, moving toward design in the 1990s, the new century in the 2010s, and the current research for papers published after 2010.

In general, design optimization seeks the best values of design variables, Xnx1, to achieve, within certain constraints, Gmx1(X) placed on the system behavior, allowable stresses, geometry, or other factors; its goal of optimality is defined by the a vector of objective functions, Fkx1(X), for specified environmental conditions. Mathematically, design optimization may be cast in the following standard form [9]:

Find the set of design variables Xnx1 that will

minimizeFX¯=i=1kwfiFiX¯E1
subject toGjX¯0,j=1,2,IE2
GjX¯=0,j=I+1,I+2,mE3

where wfi is the weighting factors measuring the relative importance of Fi(x) with respect to the overall design goal:

0wfi1i=1kwfi=1E4

Figure 1 shows the overall structure of an optimization approach to design. Major objectives in mechanical and structural engineering involve minimum fabrication cost, maximum product reliability, maximum stiffness/weight ratio, minimum aerodynamic drag, maximum natural frequencies, maximum critical shaft speeds, etc. Design variables describe configuration, dimensions and sizes of elements, and material properties as well. In the design of structural components, such as those of an automobile structure, the main design variables represent the thickness of the covering skin panels and the spacing, size, and shape of the transverse and longitudinal stiffeners. The sizes of the constituent elements of the system are measured by such properties as the cross-sectional dimensions, section areas, area moments of inertia, torsional constants, plate’s thicknesses, etc. If the skin and/or stiffeners are made of layered composites, the orientation of the fibers and their proportion can become additional variables. If one optimizes for configuration, the design variables will include spatial coordinates. Also, in dynamic problems, the location of nonstructural masses and their magnitudes can be additional design variables.

Figure 1.

Design optimization process.

Advertisement

3. Optimization techniques

The class of optimization problems described by Eqs. (1)(3) may be thought of as a search in an n-dimensional space for a point corresponding to the minimum value of the overall objective function and such that it lies within the region bounded by the subspaces representing the constraint functions. Iterative techniques are usually used for solving such optimization problems in which a series of directed design changes (moves) are made between successive points in the design space. Several optimization techniques are classified according to the way of selecting the search direction [9]. The most commonly used approaches are the random search, conjugate directions, and conjugate gradients methods. Other algorithms for solving global optimization problems may be classified into heuristic methods that find the global optimum only with high probability and methods that guarantee to find a global optimum with some accuracy. The simulated annealing technique and the genetic algorithms (GAs) belong to the former type, where analogies to physics and biology to approach the global optimum are utilized. The simulated annealing technique is an iterative search method based on the simulation of thermal annealing of critically heated solids. Hasancebi et al. [10] applied it to find the optimum design of fiber composite structures as an efficient method to solve multi-objective optimization models. On the other hand, the GAs [11, 12] are based on the principles of natural genetics and natural selection. GAs do not utilize any gradient information during the searching process. Narayana Naik et al. [12] used GA and various failure mechanisms based on different failure criteria to reach an optimal composite structure. Another robust algorithm in solving complex problems of optimal structural design is named particle swarm optimization algorithm (PSOA). This algorithm is based on the behavior of a colony of living things, such as a swarm of insects like ants, bees, and wasps, a folk of birds, or a school of fish. Omkar et al. [13] applied PSOA to achieve a specified strength with minimizing weight and total cost of a composite structure under different failure criteria. To the author’s knowledge, GA has been the most efficient stochastic method for obtaining the global optimum design of composite structures.

Advertisement

4. Application: buckling optimization of anisotropic cylindrical shells

Structural buckling failure due to high external hydrostatic pressure is a major consideration in designing cylindrical shell-type structures. This section presents a direct approach for enhancing buckling stability limits of thin-walled long cylinders that are fabricated from multi-angle fibrous laminated composite lay-ups. The mathematical formulation employs the classical lamination theory for calculating the critical buckling pressure, where an analytical solution that accounts for the effective axial and flexural stiffness separately as well as the inclusion of the coupling stiffness terms is presented. The associated design optimization problem of maximizing the critical buckling pressure has been formulated in a standard nonlinear mathematical programming problem with the design variables encompassing the fiber orientation angles and the ply thicknesses as well. The physical and mechanical properties of the composite material are taken as preassigned parameters. The proposed model deals with dimensionless quantities in order to be valid for thin shells having different thickness-to-radius ratios. Results have been obtained for cases of filament wound cylinders fabricated from different types of composite materials.

The basic analysis and analytical formulation presented in this chapter are based on the work given by Maalawi [14], which provides good sensitivity to lamination parameters and allows the search for the needed optimal stacking sequences in a reasonable computational time. Referring to the structural model depicted in Figure 2, the significant strain components are the hoop strain (εss0) and the circumferential curvature (Kss) of the mid-surface. The reduced form of the stress-strain relationships in matrix form is

NssMss=A22B22B22D22εssoκssE5

where Nss and Mss are the resultant distributed force and moment and (Aij, Bij, Dij) are the extensional, coupling, and bending stiffness coefficients, respectively [1].

Figure 2.

Laminated composite cylindrical shell under external pressure (u displacement in the axial direction x, v in the tangential direction s, w in the radial direction z).

4.1. Analytical buckling model

The governing differential equations of anisotropic long cylinders subjected to external pressure are cast in the following:

Mss+RNssβNss=βpR2E6.1
MssRNss+βNss+pwo+vo=pR2E6.2

where uo, vo, and wo are the displacements of a generic point (x, s) on the shell middle surface (z = 0) in x, s, and z directions, respectively. The prime denotes differentiation with respect to the angular position φ and β=vowo/R. For the case of thin cylinders with thickness-to-radius ratio (h/R≤ 0.1, the critical buckling pressure can be determined using the mathematical expression [14]:

pcr=3D22R31ψ2/α1+α+2ψE7.1
ψ=1RB22A22E7.2
α=1R2D22A22E7.3

4.2. Definition of the baseline design

It is convenient first to normalize all variables and parameters with respect to a baseline design, which has been selected to be a unidirectional orthotropic laminated cylinder with the fibers parallel to the shell axis x. Optimized designs shall have the same material properties, mean radius R, and total shell thickness h of the baseline design. Expressions for calculating the critical buckling pressure (Pcro) of the baseline design are defined in Table 1, which depend upon the type of composite material utilized and the shell thickness-to-radius ratio (h/R) as well.

Material typeOrthotropic mechanical properties* (GPa)Pcro × (h/R)3 (GPa)
E11E22G12ν12
E-Glass/vinyl ester41.066.732.50.2991.708
Graphite/epoxy130.07.06.00.281.757
S-Glass/epoxy57.014.05.70.2773.567

Table 1.

Material properties and critical buckling pressure of the baseline design (Pcro).

E11 = longitudinal modulus, E22 = hoop modulus, ν12 = Poisson’s ratio for axial load, ν21 = ν12E22/E11.


4.3. Optimization model

The search for the optimized lamination can be performed by coupling the analytical buckling shell model to a standard nonlinear mathematical programming procedure. The resulting optimization problem may be cast in the following standard form to

minimizep̂crE8.1
subject tohLĥkhU,E8.2
θLθkθUk=1,2,….nE8.3
k=1nĥk=1E8.4

where p̂cr=pcr/pcro is the dimensionless critical buckling pressure, and (hL, hU) are the lower and upper bounds imposed on the individual dimensionless ply thicknesses ĥk=hk/h.

According to the filament-winding manufacturing process, each ply is characterized by its angle θk with respect to the cylinder axis x. The stacking sequence is denoted by [θ1/θ2//θn], where the angles are given in degrees, starting from the outer surface of the shell. In addition, in a real-world manufacturing process, the filament-winding angles θk must be chosen from a limited range of allowable lower (θL) and upper (θU) values according to technology references. It is important to mention here that the volume fractions of the constituent materials of the composite structure is assumed to not significantly change during optimization, so that the total structural mass remains constant at its reference value of the baseline design.

4.4. Optimal solutions

The functional behavior of the candidate objective function, as represented by maximization of the dimensionless buckling pressure p̂cr, is thoroughly investigated in order to see how it is changed with the optimization variables in the selected design space. The final optimum designs recommended by the model will directly depend on the mathematical form and behavior of the objective function.

4.4.1. Two-layer anisotropic long cylinder

The first case study to be considered herein is a long thin-walled cylindrical shell fabricated from E-glass/vinyl ester composites with the lay-up composed of only two plies (n = 2) having equal thicknesses (ĥ1=ĥ2=0.5) and different fiber orientation angles. Considering the case of ±63° angle ply, the present model gives p̂cr= 4.23, i.e., Pcr = 4.23 × 1.708 × (h/R)3 GPa, depending on the shell thickness-to-radius ratio. The actual dimensional values of the critical buckling pressure for the different thickness ratios are given in Table 2 for the cases of baseline design [0°], helically wound [±63°], and [±90°] hoop layers. The unconstrained maximum value of p̂cr= 6.1 occurs at the design points [θ1/θ2= [±90, ±90].

Baseline [0°]Helically wound [±63°]Hoop plies [±90°]
p̂cr = 1.004.236.10
(h/R)
(1/15)506.072140.693087.05
(1/20)213.50903.111302.35
(1/25)109.31462.39666.80
(1/50)13.6657.8083.35

Table 2.

Critical buckling pressure for E-glass/vinyl ester cylinder with different lay-ups.

[Pcr = p̂cr × 1.708 × 106 (h/R)3 KPa].

For a two-ply long cylinder fabricated from graphite/epoxy composites, Figure 3 shows the developed level curves of the dimensionless buckling pressure, p̂cr (also named isomerits or isobars) in the (θ1θ2) design space. As seen in the figure, the maximum value of p̂cr reaches a value of 18.57 for a hoop wound construction. Table 3 presents the solutions for the [±45°] angle-ply and the [90°] cross-ply constructions for different thickness-to-radius ratios. These solutions are also valid for lay-ups [03°]s, [903°]s, [452°/−452°]s, and [45°/−45°/45°/−45°]s. The case of a helically wound lay-up construction [/−θ] with unequal play thicknesses ĥ1 and ĥ2, such that their sum is held fixed at a value of unity, has also been investigated. Computer solutions have shown that no significant change in the resulting values of the critical buckling pressure can be remarked in spite of the wide change in the ply thicknesses. This is a natural expected result since the stiffness coefficients A22, B22, and D22 remain unchanged for such lay-up construction.

Figure 3.

p̂cr-isomerits for a graphite/epoxy, two-layer cylinder in [θ1 /θ2] design space (ĥ1=ĥ2=0.5).

Baseline []Helically wound [±45°]Hoop plies [±90°]
p̂cr = 1.005.918.57
(h/R)
(1/15)520.593071.509667.40
(1/50)14.0682.93261.02
(1/120)1.025.9918.88

Table 3.

Critical buckling pressure, Pcr for graphite/epoxy cylinder with different lay-ups.

[Pcr = p̂cr × 1.757 × 106 (h/R)3 KPa].

4.4.2. Three-layer anisotropic long cylinder

Results for a cylinder constructed from three, equal-thickness layers with stacking sequence denoted by [θ1/θ2/θ1] are given in Table 4. The same behavior can be observed as before but with slight change in the attained values. It was found that for the range −30° > θ1 > 30° the critical buckling pressure is not much affected by variation in the ply angle θ2. A substantial increase in the critical buckling pressure by changing the ply angles can be observed. Similar solutions were obtained for the stacking sequences [0°2/90°]s and [90°2/0°]s.

Baseline [0°3][0°/90°/0°][90°/0°/90°]
p̂cr = 1.001.65117.92
(h/R)
(1/15)520.59859.579331.19
(1/50)14.0623.21251.94
(1/120)1.021.6818.23

Table 4.

Critical buckling pressure, Pcr for graphite/epoxy cylinder [θ1/θ2/θ1].

[Pcr = p̂cr × 1.757 × 106 (h/R)3 KPa].

4.4.3. Four-layer sandwiched anisotropic cylinder

The same graphite/epoxy cylinder is reconsidered here with changing the stacking sequence to become ±20° equal-thickness layers sandwiched in between outer and inner 90° hoop layers with unequal thicknesses, i.e., (ĥ2=ĥ3) and (ĥ1ĥ4), such that the thickness equality constraint k=14ĥk=1 is always satisfied. Figure 4 shows the developed p̂cr -isomerits in the (ĥ1,ĥ2) design space. The contours inside the feasible domain, which is bounded by the three lines ĥ1=0 and ĥ2=0 and ĥ1+2ĥ2=1 (i.e., ĥ4=0), are obliged to turn sharply to be asymptotes to the line ĥ4=0, in order not to violate the thickness equality constraint. This is why they appear in the figure as zigzagged lines. At the design point (ĥ1,ĥ2) = (0.25, 0.25), the dimensionless buckling pressure p̂cr = 16.43 (see Figure 4 and Table 5). As a general observation, as the thickness of the hoop layers increases, a substantial increase in the critical buckling pressure will be achieved, e.g., at (ĥ1,ĥ2) = (0.33, 0.17), p̂cr= 17.92 representing a percentage increase of (17.9216.43)/16.43 = 9.1%.

Figure 4.

Design space for a sandwich lay-up graphite/epoxy cylinder [90/±20/90].

(h/R)(1/15)(1/20)(1/25)(1/50)
Pcr8553.03609.31847.5231.0

Table 5.

Critical buckling pressure, Pcr (KPa), for graphite/epoxy cylinder [90/±20/90].

[p̂cr = 16.43, Pcro = 1.757 × 106 (h/R)3 KPa, Pcr = p̂cr × Pcro].

Finally, the obtained results have indicated that the optimized laminations induce significant increases, always exceeding several tens of percent, of the buckling pressures with respect to the reference or baseline design. It is assumed that the volume fractions of the composite material constituents do not significantly change during optimization, so that the total structural mass remains constant. It has been shown that the overall stability level of the laminated composite shell structures under considerations can be substantially improved by finding the optimal stacking sequence without violating any imposed side constraints. The stability limits of the optimized shells have been substantially enhanced as compared with those of the reference or baseline designs.

References

  1. 1. Daniel I, Ishai O. Engineering Mechanics of Composite Materials. 2nd ed. Oxford University Press, New York; 2006. 432p. ISBN: 978-0195150971
  2. 2. Kollar LP, Springer GS. Mechanics of Composite Structures. 1st ed. United Kingdom: Cambridge University Press; 2009. 500p. ISBN: 978- 0521126908
  3. 3. Niu MCY. Composite Airframe Structures. 3rd ed. Hong Kong: Hong Kong Conmilit Press Ltd; 2011. 500p. ISBN: 978-9627128069
  4. 4. Kassapoglou C. Design and Analysis of Composite Structures: With Applications to Aerospace Structures. 2nd ed. New Jersey: Wiley; 2013. 410p. ISBN: 978-1118401606
  5. 5. Ye J. Laminated Composite Plates and Shells, 3D Modelling. 1st ed. London: Springer-Verlag; 2003. DOI: 10.1007/978-1447100959. 273p. ISBN: 978-1447100959
  6. 6. Kumar SS. Smart Composite Structures. 1st ed. United States: Lambert Academic Publishing; 2013. 180p. ISBN: 978-3659322532
  7. 7. Sonmez FO. Optimum design of composite structures: A literature survey (1969–2009). Journal of Reinforced Plastics and Composites. 2017;36(1):3-39. DOI: 10.1177/0731684416668262
  8. 8. Ganguli R. Optimal design of composite structures: A historical review. Journal of the Indian Institute of Science. 2013;93(4):557-570. ISSN: 0970-4140 Coden-JIISAD
  9. 9. Maalawi K, Badr M. Design optimization of mechanical elements and structures: A review with application. Journal of Applied Sciences Research. 2009;5(2):221-231
  10. 10. Hasancebi O, Carba S, Saka M. Improving the performance of simulated annealing in structural optimization. Structural and Multidisciplinary Optimization. 2010;41:189-203
  11. 11. Walker M, Smith RE. A technique for the multi-objective optimization of laminated composite structures using genetic algorithms and finite element analysis. Composite Structures. 2003;62:123-128
  12. 12. Narayana Naik G, Gopalakrishnan S, Ganguli R. Design optimization of composites using genetic algorithms and failure mechanism based failure criterion. Composite Structures. 2008;83:354-367
  13. 13. Omkar SN, Senthilnath J, Khandelwal R, Narayana Naik G, Gopalakrishnan S. Artificial bee colony (ABC) for multi-objective design optimization of composite structures. Applied Soft Computing. 2011;11(1):489-499
  14. 14. Maalawi KY. Optimal buckling design of anisotropic rings/long cylinders under external pressure. Journal of Mechanics of Materials and Structures. 2008;3(4):775-793

Written By

Karam Maalawi

Submitted: 06 February 2018 Published: 05 November 2018