Optimization of Hybrid Laminates with Extension-Shear Coupling

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 extension-shear-coupled laminates with immunity to hygrothermal shear distortion (HTSD) are further derived. The extension-shear-coupled 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 extension-shear-coupled effect are simulated and verified, meanwhile the robustness of hybrid laminates is analyzed by Monte Carlo method.


Introduction
Laminated composites are playing an important and irreplaceable role in designing structures with coupling effect.For example, the bending-twisting-coupled wing structure can be designed by using composite extension-shearcoupled laminates [1] and the bending-twisting-coupled wind turbine blades structure can be designed by using composite extension-twisting-coupled laminates [2].However, with the large-scale use of these bending-twisting 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 bendingtwisting-coupled 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 in-layer composites and between-layer composites [5].The in-layer hybrid composites consist of two or more fibers which uniformly dispersed in the same matrix of the lamina; the between-layer hybrid composites are composed of two or more different laminae which consist of different single-fiber composites.In this paper, the bending-twisting-coupled structure is designed by between-layer hybrid laminates.
At present, the widely used materials in the study of the bending-twisting coupled structure are the single-fiber composite laminates [6][7][8][9][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 single-fiber composite laminates with immunity to hygrothermal shear distortion (HTSD) with only extension-shear-coupled effect-the A F B 0 D S 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 A F B 0 D S 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 single-objective optimization problems such as integer problems, real problems, and mixed integer-real problems [12].The single-objective 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 A F B 0 D S laminates with immunity to HTSD is a typical nonlinear mixed integer-real single-objective 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 extension-shear-coupled laminates with immunity to HTSD are derived.Thirdly, the DE_CMSBHS algorithm combined with the penalty function is used to optimize the extensionshear-coupled effect of hybrid laminates with immunity to HTSD.Finally, the mechanical properties of optimized laminates are verified.

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, z k is the position of the k-ply in the entire laminates, n is the number of plies of the laminates, and H 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.

Stiffness Coefficient.
The off-axis stiffness coefficients of the k-ply in the entire hybrid laminates can be defined as Wherein, θ k is the paving angle of the k-ply of the laminates.The superscript "q" is defined to mean which kind of lamina corresponds to: q = ① implies that the lamina corresponds to the type ① lamina and q = ② implies that the lamina corresponds to the type ② lamina.Moreover, U q i i = 1, 2, … , 5 are the material constants of hybrid laminates, which are only related to the material parameters of the laminae, as shown in In which, Q q ij are the stiffness coefficients of two types of lamina.The geometric factors ξ q j j = 1, 2, … , 15, q = ①, ② of hybrid laminates, which are convenient to express the stiffness matrices, are defined as where k = q 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.

Extension Stiffness Matrix
A. 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 2 International Journal of Aerospace Engineering are only related to the material constants and geometric factors, as shown in 2.1.2.Coupling Stiffness Matrix B. 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 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 k-ply in the entire hybrid laminates are wherein α q 1 and α q 2 q = ①, ② are the thermal expansion coefficients of two different types of lamina.The variation of temperature can be expressed by ΔT, and the thermal stress and thermal moment of laminates are expressed as 4 International Journal of Aerospace Engineering 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, U Tq 1 and U Tq 2 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 , 19  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 A F B 0 D S laminates with immunity to HTSD will be derived in this section.The stiffness coefficients of hybrid A F B 0 D S laminates should meet the following relationships [2].
In order to ensure that the hybrid A F B 0 D S laminates will not cause hygrothermal shear distortion, its thermal shear strain should meet following relationship.International Journal of Aerospace Engineering where According to (30) can we find out that Inserting (31) into (10) gives Then inserting (32) into (29) further gives Inserting ( 10), (35), and (37) into (36) gives Therefore, the thermal strains of two main directions of the hybrid A F B 0 D S 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 ξ q 2 , ξ q 4 , and ξ q 13 .

Optimized Design of Laminates
Considering that the extension-shear-coupled effect of hybrid A F B 0 D S laminates with immunity to HTSD is the main performance index, the maximum extension-shearcoupled 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 A F B 0 D S laminates with immunity to HTSD.Take the flexibility coefficient a 16 as the parameter to test the extension-shear-coupled effect of laminates, which can be obtained by inverting the stiffness matrix of laminates, as shown in (39).

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, q k = 1, 2 k = 1, 2, … , n represents the values of "q" for the k-ply, q k = 1 shows that the lamina corresponds to the type ① lamina, and q k = 2 shows that the lamina corresponds to the type ② lamina.If ∑ n k=1 q k = n, it means that all lamina are made of type ① material, and similarly, if ∑ n k=1 q k = 2n, it means that all lamina are made of type ② material.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-20-ply hybrid A F B 0 D S laminates with immunity to HTSD which are optimized by using DE_CMSBHS algorithm to maximize the extension-shearcoupled 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 A F B 0 D S laminates with immunity to HTSD for 1-11-ply laminates; (2) for 12-20-ply hybrid A F B 0 D S 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 extension-shear-coupled 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 extension-shear-coupled effect of laminates more intuitively, this paper further optimizes the extension-shear-coupled    8 International Journal of Aerospace Engineering effect of single-material laminates and laminates with different proportions of two materials.Table 3 shows the optimization results of two kinds of single-material laminates.Taking 16-ply laminates and 17-ply 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 extension-shear-coupled 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 single-glass fiber composite laminates of Table 3.Furthermore, the data in Tables 2 and  3 suggests that the maximum extension-shear-coupled effects of the 14-ply, 17-ply, and 20-ply hybrid laminates are greater than those of single-carbon fiber composite laminates, which indicates that the introduction of glass fiber composite materials into carbon fiber composite laminates can improve the extension-shear-coupled 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 single-carbon material composite laminates, some hybrid laminates can not only significantly reduce costs but also increase the maximum extensionshear-coupled effect.

Verification of Mechanical
Properties of Laminates 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 A F B 0 D S laminates with immunity to HTSD are shown in Table 4.In which, ε T x and ε T y are the thermal strain of two main directions, γ T xy is the thermal shear strain of the laminates, and κ T x , κ T y , and κ T xy 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 high-temperature 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.

Verification of Extension-Shear-Coupled Effect.
In order to verify the extension-shear-coupled effect of two kinds of hybrid A F B 0 D S 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 F = 1000 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 A F B 0 D S 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 A F B 0 D S laminates not only have axial distortion but also shear distortion and Table 5 also shows the theoretical calculation results of two kinds of hybrid A F B 0 D S 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 extension-shear-coupled 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.

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 extension-shear-coupled effect, in order to ensure the practicability of the composite laminates, the slight angle deviation should not have an obvious influence on the extension-shear-coupled effect of the hybrid A F B 0 D S laminates with immunity to HTSD.Now, the extension-shearcoupled 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 k-ply for the hybrid A F B 0 D S laminates with immunity to HTSD is θ k ± Δθ k , in which θ k is the theoretical paving angle of the k-ply and Δθ k is the angle deviation and with the value of Δθ k = 2 ° [17].Figure 6 shows the error distribution of the extensionshear-coupled effect for the 16-ply laminates and the 17-ply 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 extension-shearcoupled effect is in accordance with the normal distribution law and can be controlled within 2%.

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 extension-shear-coupled 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 extensionshear-coupled effect of the hybrid A F B 0 D S laminates with immunity to HTSD.The hygrothermal effect, extensionshear-coupled 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 A F B 0 D S 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 extension-shear-coupled effect of the single-glass fiber composite laminates is greater than that of single-carbon 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 extension-shear-coupled 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 extension-shear-coupled 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.

3 .
The Hybrid A F B 0 D S Laminates with Immunity to HTSD 3.1.Necessary and Sufficient Conditions.From (

A 11 A 12 A 16 A 12 A 22 A 26 A 16 A
connection between thermal stress N T and thermal shear strain γ T xy of A F B 0 D S hybrid laminates is

a 11 a 12 a 16 b 11 b 12 b 16 a
12 a 22 a 26 b 12 b 22 b 26 a 16 a 26 a 66 b 16 b 26 b 66 b 11 b 12 b 16 d 11 d 12 d 16 b 12 b 22 b 26 d 12 d 22 d 26 b 16 b 26 b

θ
k represents the paving angle of the k-ply, whose initial value is selected randomly from the range of [−90 °, 90 °].

Figure 2 :
Figure 2: Maximal extension-shear-coupled effect of hybrid laminates with different proportions of two materials.

Figure 3 :
Figure 3: Finite element model of the laminates for hygrothermal effect verification.

Figure 4 :
Figure 4: Finite element model of the laminates for extension-shear-coupled effect verification.

Figure 5 :
Figure 5: Deformation of two kinds of hybrid A F B 0 D S laminates due to extension load.
2 26 A 16 A 26 − A 12 A 66 A 12 A 26 − A 16 A 22 A 16 A 26 − A 12 A 66 A 11 A 66 − A 2 16 A 16 A 16 − A 11 A 26 A 12 A 26 − A 16 A 22 A 16 A 16 − A 11 A 26 A 11 A 22 − A 2 12 25Substituting (25) into (24) can be used to obtain the expression of thermal shear strain of the hybrid A F B 0 D S laminates.
26 − A 16 A 22 N T x + A 12 A 16 − A 11 A 26 N T y + A 11 A 22 − A 2 12 N T xy 26 Inserting (26) into (23) gives A 12 A 26 − A 16 A 22 N T x + A 12 A 16 − A 11 A 26 N T y + A 11 A 22 − A A 12 A 26 − A 16 A 22 + A 12 A 16 − A 11 A 26 〠 q 13 + A 12 A 26 − A 16 A 22 − A 12 A 16 − A 11 A 26 〠 A 12 A 26 − A 16 A 22 + A 12 A 16 − A 11 A 26 ξ q 13

11 A 12 A 16 B 11 B 12 B 16 A 12 A 22 A 26 B 12 B 22 B 26 A 16 A 26 A 66 B 16 B 26 B 66 B 11 B 12 B 16 D 11 D 12 D 16 B 12 B 22 B 26 D 12 D 22 D 26 B 16 B 26 B 66 D 16 D 26 s 66
66 d 16 d 26 d 66 = A

Table 1 :
Material properties of carbon fiber lamina and glass fiber lamina.

Table 2 :
Hybrid A F B 0 D S laminates with immunity to HTSD.

Table 3 :
Maximum extension-shear-coupled effect of two single-composite 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 A F B 0 D S laminates with immunity to HTSD.In this paper, the hygrothermal effect of laminates in Table 2 is verified by the example of 16-ply and 17-ply

Table 5 :
Simulation results of hybrid A F B 0 D S laminates on extension-shear distortion.Type of laminates Simulation result ε x Theoretical result ε x Error of ε x Simulation result γ xy Theoretical result γ xy Error of γ xy