Modelling and Optimizing Structural Behavior of Advanced Materials for Aerospace
View this Special IssueResearch Article  Open Access
Da Cui, Daokui Li, "Optimization of Hybrid Laminates with ExtensionShear Coupling", International Journal of Aerospace Engineering, vol. 2018, Article ID 9869105, 12 pages, 2018. https://doi.org/10.1155/2018/9869105
Optimization of Hybrid Laminates with ExtensionShear Coupling
Abstract
The introduction of hybrid composites into the structure with coupling effect can greatly reduce the cost of materials. The expressions of stiffness coefficient, thermal stress, and thermal moment for hybrid laminates are derived based on the geometrical factors of laminates, and the necessary and sufficient conditions for the hybrid extensionshearcoupled laminates with immunity to hygrothermal shear distortion (HTSD) are further derived. The extensionshearcoupled effect of hybrid laminates is optimized with improved differential evolution algorithm. Results are presented for the hybrid laminates that consist of carbon fiber and glass fiber composite materials. The hygrothermal effect and extensionshearcoupled effect are simulated and verified, meanwhile the robustness of hybrid laminates is analyzed by Monte Carlo method.
1. Introduction
Laminated composites are playing an important and irreplaceable role in designing structures with coupling effect. For example, the bendingtwistingcoupled wing structure can be designed by using composite extensionshearcoupled laminates [1] and the bendingtwistingcoupled wind turbine blades structure can be designed by using composite extensiontwistingcoupled laminates [2]. However, with the largescale use of these bendingtwisting coupling structures, the common glass fiber composites have difficulties in meeting the requirements of structural reliability. On the other hand, the carbon fiber composites with good comprehensive properties [3, 4] cost about ten times more than glass fiber composites, which restricts its wide range of applications. Therefore, it is necessary to introduce hybrid fiber composites into the design of the bendingtwistingcoupled structure to achieve the purpose that greatly reduces material cost under meeting the structural reliability requirements.
According to different hybrid modes, hybrid fiber composites are mainly divided into two types of composites, the inlayer composites and betweenlayer composites [5]. The inlayer hybrid composites consist of two or more fibers which uniformly dispersed in the same matrix of the lamina; the betweenlayer hybrid composites are composed of two or more different laminae which consist of different singlefiber composites. In this paper, the bendingtwistingcoupled structure is designed by betweenlayer hybrid laminates.
At present, the widely used materials in the study of the bendingtwisting coupled structure are the singlefiber composite laminates [6–10], whose important design parameters contain paving angles merely. However, the important design parameters for hybrid laminates also include the paving materials [11]. Once the paving materials become variable, the design of laminates becomes more complicated and the number of optimization constraints even multiplies at the same time.
J. Li and D. Li [2] have designed a kind of singlefiber composite laminates with immunity to hygrothermal shear distortion (HTSD) with only extensionshearcoupled effect—the laminates (refer to the nomenclature of references [1, 2])—and the sequential quadratic programming (SQP) is used to optimize its coupled effect. However, we found that feasible solutions cannot be found when using the SQP algorithm to optimize the hybrid laminates with immunity to HTSD, for this algorithm cannot meet the strong constraints of hybrid laminates. Therefore, this paper takes the improved differential evolution algorithm DE_CMSBHS to optimize this problem.
The DE_CMSBHS algorithm is an efficient global optimization algorithm, which has the characteristics of simple structure, easy realization, fast convergence, and strong robustness, and can effectively solve the singleobjective optimization problems such as integer problems, real problems, and mixed integerreal problems [12]. The singleobjective nonlinear optimization problems are able to be solved well by combining the penalty function’s ability to handle constraints with the optimal performance of DE_CMSBHS algorithm [13, 14]. The optimization problem of hybrid laminates with immunity to HTSD is a typical nonlinear mixed integerreal singleobjective optimization problem.
In this paper, the model of hybrid laminates is firstly established. Then the expressions of stiffness coefficient, thermal stress, and thermal moment for laminates are derived with important parameter geometric factor [15, 16]. Secondly, the necessary and sufficient conditions for the hybrid extensionshearcoupled laminates with immunity to HTSD are derived. Thirdly, the DE_CMSBHS algorithm combined with the penalty function is used to optimize the extensionshearcoupled effect of hybrid laminates with immunity to HTSD. Finally, the mechanical properties of optimized laminates are verified.
2. Stiffness Coefficient, Thermal Stress, and Thermal Moment of Hybrid Laminates
Introducing the geometrical factors into the design of laminates can effectively improve the efficiency of this progress. In this section, the model of hybrid laminates will be established. Based on the geometric factors, the expressions of stiffness coefficient, thermal stress, and thermal moment of hybrid laminates are derived.
The research object is set as a kind of hybrid laminates which are composed of two different types of lamina, and each lamina has the same thickness, as shown in Figure 1. Wherein, is the position of the kply in the entire laminates, is the number of plies of the laminates, and is the entire thickness of the laminates. and represent two kinds of laminae with different material properties, respectively, and the number and layer order of each kind of lamina are variables.
2.1. Stiffness Coefficient
The offaxis stiffness coefficients of the kply in the entire hybrid laminates can be defined as
Wherein, is the paving angle of the kply of the laminates. The superscript “” is defined to mean which kind of lamina corresponds to: implies that the lamina corresponds to the type lamina and implies that the lamina corresponds to the type lamina. Moreover, are the material constants of hybrid laminates, which are only related to the material parameters of the laminae, as shown in
In which, are the stiffness coefficients of two types of lamina. The geometric factors of hybrid laminates, which are convenient to express the stiffness matrices, are defined as where means the sum of all the type or laminae on the right side of the equation. Obviously, the geometric factors are only related to the paving angle and location of each lamina.
2.1.1. Extension Stiffness Matrix
According to the classical theory of laminates, the extension stiffness matrix of laminates can be expressed as
Substituting (1), (2), and (3) into (4) can be used to obtain the expressions of the extension stiffness coefficients, which are only related to the material constants and geometric factors, as shown in
2.1.2. Coupling Stiffness Matrix
Similarly, the coupling stiffness matrix of the hybrid laminates can be expressed as
Substituting (1), (2), and (3) into (6) can be used to obtain the expression of the coupling stiffness coefficients, which are only related to the material constants and geometric factors, as shown in
2.1.3. Bending Stiffness Matrix
The bending stiffness matrix of the hybrid laminates can be expressed as
Substituting (1), (2), and (3) into (8) can be used to obtain the expressions of the bending stiffness coefficients, which are only related to the material constants and geometric factors, as shown in
To sum up, the stiffness coefficients of the hybrid laminates, which are only related to the material constants and geometric factors, can be expressed as
2.2. Thermal Stress and Thermal Moment
The thermal expansion coefficients of the kply in the entire hybrid laminates are wherein and are the thermal expansion coefficients of two different types of lamina. The variation of temperature can be expressed by , and the thermal stress and thermal moment of laminates are expressed as
Substituting (1) and (13) into (14) can be used to obtain the expressions of thermal stress, which are only related to the material constants and geometric factors, as shown in
In which, and are defined as the thermal invariants of the hybrid laminates, which can be calculated from the thermal expansion coefficients and the invariants of laminates, as shown in
Substituting (1) and (13) into (15) can be used to obtain the expressions of the thermal moment, which are only related to the material constants and geometric factors, as shown in
To sum up, the thermal stress and thermal moment of the hybrid laminates, which are only related to the thermal invariants and geometric factors, can be expressed as
3. The Hybrid Laminates with Immunity to HTSD
3.1. Necessary and Sufficient Conditions
From (10), (11), and (12), for matrices A, B, and D, respectively, can we see that the stiffness coefficients of laminates are only related to geometrical factors and material constants. In order to make the designed laminates suitable for all materials, the necessary and sufficient conditions of geometric factors for the hybrid laminates with immunity to HTSD will be derived in this section. The stiffness coefficients of hybrid laminates should meet the following relationships [2].
Substituting (11) and (12) into Eq. (21) can be used to obtain the necessary and sufficient conditions of geometric factors for the hybrid laminates.
In order to ensure that the hybrid laminates will not cause hygrothermal shear distortion, its thermal shear strain should meet following relationship.
Moreover, the connection between thermal stress and thermal shear strain of hybrid laminates is where
Substituting (25) into (24) can be used to obtain the expression of thermal shear strain of the hybrid laminates.
Inserting (26) into (23) gives
Inserting (19) into (27) and simplifying gives
If (28) is constantly established, the previous coefficient of thermal material constant and must be zero.
According to (30) can we find out that
Inserting (31) into (10) gives
Then inserting (32) into (29) further gives
Furthermore, because of and , (29) is also established. Considering that if , and of laminates are both zero and the laminates will not have extensionshearcoupled effect. Therefore, geometric factor cannot be zero; thus, the necessary and sufficient conditions for “with immunity to HTSD” can be expressed as
In summary, the necessary and sufficient conditions of the hybrid laminates with immunity to HTSD can be expressed as
3.2. Curing Deformation
Considering that the curing deformation phenomenon will occur during the process of molding, which may bring about deformation difference between the expected design shape and the free shape after taking off the model at room temperature, now, the thermal strain of the hybrid laminates with immunity to HTSD will be derived by the necessary and sufficient conditions during the process of curing deformation.
Due to the similarity of influence between humidity changes and temperature changes on composite materials, it only needs to replace the thermal expansion coefficients by humidity expansion coefficients. In order to simplify the analysis and derivation process, the analysis for hygrothermal distortion contains the thermal effect merely.
Substituting (25) into (24) can be used to obtain the expressions of thermal strain of the hybrid laminates with immunity to HTSD.
Substituting (35) into (19) can give the expressions of thermal stress of the hybrid laminates with immunity to HTSD.
Inserting (10), (35), and (37) into (36) gives
Therefore, the thermal strains of two main directions of the hybrid laminates with immunity to HTSD are equal to each other and the value of the thermal strain is not only related to the material constants and temperature variation but also changing with geometric factors , , and .
4. Optimized Design of Laminates
Considering that the extensionshearcoupled effect of hybrid laminates with immunity to HTSD is the main performance index, the maximum extensionshearcoupled effect is therefore the major objective. Regarding the ply materials and ply angles of each lamina as variables, take the ply materials and ply angles as the optimized design variables. The optimized constraint conditions can be obtained by (35), which ensures that the final optimized laminates are hybrid laminates with immunity to HTSD. Take the flexibility coefficient as the parameter to test the extensionshearcoupled effect of laminates, which can be obtained by inverting the stiffness matrix of laminates, as shown in (39).
However, taking into account the actual engineering demands, which should be combined with the good comprehensive properties of carbon fiber composites and the low cost of glass fiber composites, the proportion of two kinds of material must have a minimum requirement. On the one hand, if the proportion of carbon fiber composites is too small, the structural stability of the hybrid laminates cannot be guaranteed. On the other hand, if the proportion of glass fiber composites is too small, the cost reduction effect is not obvious. Therefore, in this paper, the minimum proportions of two types of material are both set as 30% and the mathematical model of optimization problem can be formulated as
Wherein, represents the values of “q” for the kply, shows that the lamina corresponds to the type lamina, and shows that the lamina corresponds to the type lamina. If , it means that all lamina are made of type material, and similarly, if , it means that all lamina are made of type material. represents the paving angle of the kply, whose initial value is selected randomly from the range of [−90°, 90°].
The improved differential evolution algorithm DE_CMSBHS is used to optimize this problem. Because of the strong constraints of optimization problems, which have 18 equality constraints and 2 inequality constraints, the penalty function is applied to deal with the constraints. The type lamina and type lamina are made up of carbon fiber composites and glass fiber composites, respectively, and the corresponding material parameters are shown in Table 1.

Table 2 shows the 12–20ply hybrid laminates with immunity to HTSD which are optimized by using DE_CMSBHS algorithm to maximize the extensionshearcoupled effect. The subscripts “c” and “gl” in the table indicate that the lamina is made of carbon fiber composites and glass fiber composites, respectively. As can be seen from the table, (1) there is no hybrid laminates with immunity to HTSD for 1–11ply laminates; (2) for 12–20ply hybrid laminates with immunity to HTSD, the proportions of two kinds of lamina are both not less than 30%; and (3) for optimized laminates with the maximum extensionshearcoupled effect, paving materials are symmetrical about the geometric middle plane and have a relatively uniform distribution, which is the one in the middle and the other on both sides.

In order to reflect the influence of the hybrid form on the extensionshearcoupled effect of laminates more intuitively, this paper further optimizes the extensionshearcoupled effect of singlematerial laminates and laminates with different proportions of two materials. Table 3 shows the optimization results of two kinds of singlematerial laminates. Taking 16ply laminates and 17ply laminates as an example, the optimized results of laminates with different proportions of two materials are exhibited in Figure 2.

It can be obtained from Table 3 that the maximum extensionshearcoupled effect of glass fiber composite laminates is greater than that of carbon fiber composite laminates in the case of the same number of plies. In the process of optimization by DE_CMSBHS algorithm, if the minimum proportion limits of the two materials are removed, which means the paving material is free, the final optimal results are totally the same as those of singleglass fiber composite laminates of Table 3. Furthermore, the data in Tables 2 and 3 suggests that the maximum extensionshearcoupled effects of the 14ply, 17ply, and 20ply hybrid laminates are greater than those of singlecarbon fiber composite laminates, which indicates that the introduction of glass fiber composite materials into carbon fiber composite laminates can improve the extensionshearcoupled effect of partial numbers of plies laminates.
The abscissa “n_{gl}” in Figure 2 represents the number of glass fiber lamina. As can be seen from Figure 2, (1) for hybrid laminates with different proportions of two materials, there is no feasible solution when n_{gl} takes some values, for which the geometrical factors of those kinds of laminates cannot satisfy the necessary and sufficient condition of (35) and (2) compared with singlecarbon material composite laminates, some hybrid laminates can not only significantly reduce costs but also increase the maximum extensionshearcoupled effect.
5. Verification of Mechanical Properties of Laminates
5.1. Verification of Hygrothermal Effect
The finite element method is used to verify the distortion caused by the temperature change of the hybrid laminates with immunity to HTSD. In this paper, the hygrothermal effect of laminates in Table 2 is verified by the example of 16ply and 17ply laminates and the conclusion of other laminates is the same as that of these laminates.
Based on the finite element software MSC.Patran, the 10 m × 1 m finite element model is established and 360 shell units are divided. In order to simulate the displacement boundary condition of the composite laminates, the geometric center of the finite element model is fixed, as shown in Figure 3. The typical temperature difference of the hightemperature curing process is −180°C to this finite element model. Then the finite element software MSC.Nastran is used to compute with the linear statics calculation function.
The calculated results of two kinds of hybrid laminates with immunity to HTSD are shown in Table 4. In which, and are the thermal strain of two main directions, is the thermal shear strain of the laminates, and , , and are expressed as the surface curvature and distortion of the laminates caused by temperature changes. It can be seen from the table that the shear strain of these laminates are all zero during the hightemperature curing process, which means two kinds of laminates will not cause hygrothermal distortion. The thermal strains of two directions are equal to each other, but this regulation is not suitable for laminates with different numbers of plies. Furthermore, their bending curvature and twist rate are all zero, which indicates that two kinds of laminates will not cause hygrothermal warping distortion.

5.2. Verification of ExtensionShearCoupled Effect
In order to verify the extensionshearcoupled effect of two kinds of hybrid laminates in Table 2, the finite element method is also used based on the finite element software MSC.Patran and the finite element model of a rectangular plate with 12 m × 1 m is established. A total of 800 shell units and 891 nodes (six degrees of freedom) are divided, and a multipoint constraint element (RBE2) is used to connect the nodes in the 1.2 m × 1 m region at both ends of the model. The axial tension N is applied to the multipoint confinement unit to ensure that the intermediate 9.6 m × 1 m area bears a uniform load. The geometric center of the finite element model is fixed, as shown in Figure 4. The finite element software MSC.Nastran is used to compute with the linear statics calculation function.
The calculated displacement nephogram of two kinds of hybrid laminates with immunity to HTSD under axial extension force is shown in Figure 5, and the concrete results of distortion are shown in Table 5. It can be seen from the table that under the axial extension force of 1000 N, the two kinds of hybrid laminates not only have axial distortion but also shear distortion and the extensionshearcoupled effect of 16ply laminates is obviously more than that of 17ply laminates.
(a) 16ply laminates
(b) 17ply laminates

Table 5 also shows the theoretical calculation results of two kinds of hybrid laminates subjected to axial extension distortion, and through comparison, it is found that the results of finite element analysis agree well with the theoretical results. The error is controlled within 2%, and the extensionshearcoupled effect of laminates is validated. The reason for the error is that the finite element simulation is loaded on the unit node and cannot be completely equivalent to the linear loading method.
5.3. Robustness Analysis
Considering that there may be human error and equipment error in the actual process of paving the laminates, which may have a negative effect on the extensionshearcoupled effect, in order to ensure the practicability of the composite laminates, the slight angle deviation should not have an obvious influence on the extensionshearcoupled effect of the hybrid laminates with immunity to HTSD. Now, the extensionshearcoupled effect of laminates is analyzed under the ply angle existing deviation, which is based on the Monte Carlo composite laminates robustness analysis method.
Assume that the paving angle of the kply for the hybrid laminates with immunity to HTSD is , in which is the theoretical paving angle of the kply and is the angle deviation and with the value of [17]. Figure 6 shows the error distribution of the extensionshearcoupled effect for the 16ply laminates and the 17ply laminates in Table 2 under 10,000 random samples. The figure suggests that when there exists a random error in the paving angle of the laminates, the error of extensionshearcoupled effect is in accordance with the normal distribution law and can be controlled within 2%.
(a) 16ply laminates
(b) 17ply laminates
6. Summary
In this paper, the model of hybrid laminates is established and the expressions of stiffness coefficient, thermal stress, and thermal moment of laminates are derived based on the geometric factor and the necessary and sufficient conditions for hybrid extensionshearcoupled laminates with immunity to HTSD are further derived. The method of combing the improved differential evolution algorithm DE_CMSBHS with the penalty function is used to optimize the extensionshearcoupled effect of the hybrid laminates with immunity to HTSD. The hygrothermal effect, extensionshearcoupled effect, and robustness of laminates are verified. Some conclusions of optimization and simulation can be obtained as follows. (1)Although the optimized constraints is enhanced, which is mainly due to the introduction of variable parameter paving materials, the feasible solutions of the hybrid laminates with immunity to HTSD can still be obtained by using the DE_CMSBHS algorithm combined with the penalty function. The optimal solution can be obtained according to the different requirements of different proportions of two materials, which achieved the synchronous optimization of paving angle and paving materials.(2)The maximum extensionshearcoupled effect of the singleglass fiber composite laminates is greater than that of singlecarbon fiber composite laminates in the case of the same number of plies, and the introduction of glass fiber composite material into singlecarbon fiber composite laminates can improve the extensionshearcoupled effect of the partial number of plies laminates.(3)The optimized laminates are able to meet the conditions of no hygrothermal shearing distortion and no hygrothermal warping distortion. When there exists a random error in the paving angle of the laminates, the error of extensionshearcoupled effect is in accordance with the normal distribution law and can be controlled within 2%.
The expressions for parameters of the hybrid laminates can also be applied to laminates with other coupling effects. The DE_CMSBHS algorithm combined with the penalty function can also optimize the different objectives of other types of laminates.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of the National Natural Science Foundation of China (Grant no. 11472003).
References
 C. B. York, “On extension–shearing coupled laminates,” Composite Structures, vol. 120, pp. 472–482, 2015. View at: Publisher Site  Google Scholar
 J. Li and D. Li, “Extensionshear coupled laminates with immunity to hygrothermal shearing distortion,” Composite Structures, vol. 123, pp. 401–407, 2015. View at: Publisher Site  Google Scholar
 D. Berry and T. Ashwill, “Design of 9meter carbonfiberglass prototype blades: CX100 and TX100,” Tech. Rep., SAND20070201, Sandia National Laboratories, Warren, RI, USA, 2007. View at: Google Scholar
 C. J. Wang and Y. F. Dai, “Application of carbon fiber composite in aerospace,” Development & Innovation of Machinery & Electrical Products, vol. 2, 2010. View at: Google Scholar
 R. M. Jones, Mechanics of Composite Materials, Scripta Book Company, Washington, DC, USA, 1975.
 M. Moore and D. Frei, “X29 forward swept wing aerodynamic overview,” in Applied Aerodynamics Conference, Danvers, MA, USA, 1983. View at: Publisher Site  Google Scholar
 N. J. Krone Jr., “Divergence elimination with advanced composites,” in Aircraft Systems and Technology Meeting, Los Angeles, CA, USA, 1974. View at: Publisher Site  Google Scholar
 V. Tischer and V. Venkayya, “Plyorientation as a variable in multidisciplinary optimization,” in 4th Symposium on Multidisciplinary Analysis and Optimization, Cleveland, OH, USA, 1992. View at: Publisher Site  Google Scholar
 M. Patil and M. Patil, “Aeroelastic tailoring of composite box beams,” in 35th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 1997. View at: Publisher Site  Google Scholar
 S. Guo, “Aeroelastic optimization of an aerobatic aircraft wing structure,” Aerospace Science and Technology, vol. 11, no. 5, pp. 396–404, 2007. View at: Publisher Site  Google Scholar
 S. Adali and V. E. Verijenko, “Optimum stacking sequence design of symmetric hybrid laminates undergoing free vibrations,” Composite Structures, vol. 54, no. 23, pp. 131–138, 2001. View at: Publisher Site  Google Scholar
 Y. Zhu, H. Wang, and J. Zhang, “Spacecraft multipleimpulse trajectory optimization using differential evolution algorithm with combined mutation strategies and boundaryhandling schemes,” Mathematical Problems in Engineering, vol. 2015, Article ID 949480, 13 pages, 2015. View at: Publisher Site  Google Scholar
 J. A. Snyman, N. Stander, and W. J. Roux, “A dynamic penalty function method for the solution of structural optimization problems,” Applied Mathematical Modelling, vol. 18, no. 8, pp. 453–460, 1994. View at: Publisher Site  Google Scholar
 G. Di Pillo and L. Grippo, “An exact penalty function method with global convergence properties for nonlinear programming problems,” Mathematical Programming, vol. 36, no. 1, pp. 1–18, 1986. View at: Publisher Site  Google Scholar
 C. York, “Unified approach to the characterization of coupled composite laminates: hygrothermally curvaturestable configurations,” International Journal of Structural Integrity, vol. 2, no. 4, pp. 406–436, 2011. View at: Publisher Site  Google Scholar
 C. B. York, “Tapered hygrothermally curvaturestable laminates with nonstandard ply orientations,” Composites Part A: Applied Science and Manufacturing, vol. 44, pp. 140–148, 2013. View at: Publisher Site  Google Scholar
 R. J. Cross, R. A. Haynes, and E. A. Armanios, “Families of hygrothermally stable asymmetric laminated composites,” Journal of Composite Materials, vol. 42, no. 7, pp. 697–716, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Da Cui and Daokui Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.