Homogenization Method to Calculate the Stiffness Matrix of Laminated Composites

New analytical models have been developed for predicting equivalent Young’s and shear moduli of laminate composites. Sets of procedures and calculations are presented in order to obtain equivalent properties in all levels, lamina and laminate. An ultimate path to predict the mechanical properties of laminated composites on the perspective of the material, orientation, and thickness has been developed. By calculating the mechanical properties using Chamis model then an Objectif function with five norms, these norms allow the mechanical properties to be examined and the ultimate answer to be predicted. Another model discusses an alternative concept of equivalent lamina elements (ELEs) by first using Chamis model for hybrid composites. Next, the ELEs are laminated in the direction and integrated into the compliance matrix for each ply. In addition, four new Generalization models for equivalence in thickness and in angle are presented in this paper. The analytical results are validated against other developed models in published articles as well as experimental results. Numerical case studies were conducted to assess the precision of results from the suggested models. The results demonstrated the capability and efficiency of the presented models for predicting the mechanical properties of multi-layer/multi-material laminate composites under different orientation conditions.


Introduction
Composite materials consist of two or more components such as fiber and matrix and may or may not be fillers in case of polymer composite materials. They have a wide range of applications, such as industrial, automotive engineering medical aerospace electrical, transportation, and sports [1].
The integration of two materials gives better properties compared with individual components [2].
Composites structures are often classified as four basic structures or levels, i.e., lamina, laminate, sandwich, and woven composites. A lamina (also called a ply or layer) is a single flat layer of unidirectional fibers or woven fibers arranged in a matrix, a thin layer of a composite material, while a laminate is a stack of plies of composites. Each layer can be laid at various orientations and can be made up of different material systems.
In recent years, conventional materials are continuously being replaced by a variety of composite materials. This fact has caused an increasing interest on the modeling of composites. Several approaches have been developed, but there is still a strong need of predicting models that can be used for the stiffness and strength assessment of this type of materials in practical situations. [1] Constitutive models for fiber-reinforced composite laminates can be classified according to the scale in which they are defined [3].

Problematic
In composite materials, the properties are different depending on the fiber orientation in the matrix, so knowing the properties or obtaining the correct and more accurate properties is a very important key when it comes to analyzing the structures [4]. Using composite materials in engineering industries requires complicated analysis and modeling which in most cases computer software runs. It is difficult and time-consuming to calculate the properties of a very thick laminate with each layer having different stacking sequences. It would be simpler and less time-consuming to use average properties that would represent the desired model, instead of going through calculations and obtaining stiffness matrices for each layer, i.e., ply-by-ply approach.
While experimental tests made to determine the effective elastic properties of composites, i.e., Young's and shear moduli (Ex, Ey, Ez, Gxy, Gxz, and Gyz) and Poisson's ratios (νxy, νxz, and νyz), are expensive in terms of money and time, analytical and numerical modeling could present a good alternative solution for industries and the obtained results are supposed to be more accurate. These methods have the disadvantages of a great deal of time consumed on the tasks of defining the geometry of the composite, meshing, and calculation procedures, besides the need of super computers to make the computation. Analytical modeling presents a more flexible and easier tool in terms of geometric modeling and the time saved in calculation, without the need of super computers. Therefore, in order to overcome these difficulties, equivalent properties of lamina or laminate are often used. Equivalent properties or effective stiffness properties are an average measure of the stiffness of the desired material. The actual averaging process must be done carefully and delicately in order to obtain equivalent properties.
Up to now, many models have been proposed to analyze the microstructures and mechanical performance of laminates composites. Many scholars have paid great efforts to the prediction of mechanical properties of laminate composite, but there is no physical model or method that can predict the mechanical properties and behavior of the composite materials with different sets of parameters in terms of orientation, material properties, and thicknesses. The novelty of this paper is to provide the most accurate results of the mechanical properties using an effective algorithm that can be used as an alternative to software or experimental studies. Different analytical models predict different elastic properties of composites. Additional problems involve estimating model adequacy and choosing the most accurate model. Different models are more accurate, depending on the composite content, the shape of the inclusions, and the relationship between the properties of the inclusions and the matrix [2].
The works on the mechanical analysis and optimization design of laminate composites are extremely limited. To the author's knowledge, there is no literature on applying the classical laminate theory to analyze the mechanical properties of laminate composites. In addition, there is no literature on using the ABD matrix (The ABD matrix is a 6 × 6 matrix that serves as a connection between the applied loads and the associated strains in the laminate. It essentially defines the elastic properties of the entire laminate) to find and predict the mechanical/engineering properties of laminated composites, such as compliance or stiffness matrices. Therefore, it is valuable to propose a theoretically analytical model to investigate the mechanical properties of laminate composites due to its advantages in low modeling cost and high computational efficiency.
It is important to calculate strain and stress to know the effectiveness of compound combinations. Combined applications are widely available in the field of automotive, airports, and sports.

State of Art
There have been a number of publications regarding predicting composite materials properties and equivalent or effective properties of composite materials. Richardson [5] predicted stiffness using rule of mixtures (ROM), ROM with efficiency factors, classical laminate analysis, and empirical formulas. To discuss effective properties of composites, Lydzba [6] used the homogenization method to define the heterogeneous medium that possesses the property of statistical homogeneity as an equivalent homogeneous material that would have the same average properties. Chan [7] represented the lumping procedures, in which a group of plies lumps into a single layer to reduce the size of the problem and to increase the efficient computation, which is often used in thick composite structure modelling. On the other hand, Kumar [8] presented a new analytical model to replace the original laminate ABD by much simpler ABD-equivalent virtual material models. These virtual material models, under the usual assumptions made in the lamination theory, have the same constitutive relationship as the corresponding two-dimensional (2D) plate model of the original laminate but require only a small fraction of computational integration costs in three dimensions.
It is noted that all the numerical approaches consist in building a representative volume element (RVE) able to statistically represent the overall material [9]. Prescribed loads or displacements are applied at the RVE, and the resulting stresses and strains are correlated with the whole structure at the macroscale by means of the averaging principle [10]. However, 2D RVEs may not give accurate results for composite materials [5], especially for out-of-plane properties, demanding more laborious and time-consuming three-dimensional (3D) models [6]. As an alternative, analytical models are used for composite homogenization, usually being derived from the strength of materials or semiempirical relations. Numerical and analytical models usually share important assumptions in the homogenization process, namely neglecting the fiber/matrix interface and, homogeneous fibers and matrixes and perfect fiber alignment [9]. In addition, in the evaluation process, the properties of the individual constituents are required, and these may be hard to obtain, especially for fibers [9].
Alternatively, as proposed by Chamis [13], HC (Hybrid Composite) properties can be estimated by splitting it on the primary composite (PC) and the secondary composite (SC), where the first is reinforced by fiber 1 and the second is reinforced by fiber2, with the matrix proportionally distributed. The PC and SC properties are calculated analogously to single-fiber composites [8]. Then, the ROM is applied to compute E2 from PC and SC. This same methodology is adopted by Chamis to compute G12 and G23. ν23 can be computed under the transversally isotropic consideration [9].
The model presented by Chou [11] applies the same division and equations proposed by Chamis [8] to compute PC and SC properties. However, the final properties are evaluated using an inverse ROM.
The relations presented above are well accepted in the literature for laminate composites [14]. However, many new models with greater complexity were developed, avoiding assumptions related to long fibers [15] and straight fibers [16] or focusing on other reinforcement geometries, such as particles or nanotubes [17]. The semi-analytical model proposed by Aboudi [18] was also applied for HC [19], yielding good mechanical predictions. The tensile loadings of specimens produce different failure types, and researchers have studied the inhomogeneity of the matrix [20]. Research interest has been focused on towards fiber-reinforced polymeric composites, as these composite materials exhibit excellent mechanical properties. Parvanesh et al. [21] studied the mechanical behavior of PVC (poly(vinyl chloride)) nanocomposites. They have studied the Young's modulus and tensile strength of polymer composites, and the result suggested that good stress transfer can be obtained at an amorphous interface, depending on the polymer.
The local/global stiffness matrix approach is similar to the transfer and flexibility matrix formulation presented by Bufler [22] for isotropic layered media with constant elastic properties. This approach was later reformulated in terms of the local stiffness matrix by Rowe and Booker [23] and applied to nonhomogeneous isotropic layered soils. Chatterjee et al. [24] extended the method to anisotropic layered media for generalized plane deformation problems and used it to solve the problem of an interlaminar crack in an arbitrarily laminated composite plate under three-point loading.

Paper Organization
The first section introduces basic ideas about composites including fibers and matrices and basic definitions used in the mechanics of composites. This section starts with a review of basic topics of stress, strain, and elastic moduli. Then, the mechanical behaviors of a single lamina were discussed, including concepts about the stress-strain relationship for a lamina and the stiffness and strength of a lamina, and finally, the equivalence in orientation was discussed by studying the case of rotating fibers and rotating the whole ply, i.e., fiber and matrix in a laminate composite.
Secondly, two ways to find the equivalence in thickness were presented using the greatest common divisor (GCD) and another mathematical method based on the stressstrain relation and Hook's Law.
Thirdly, two new analytical models were developed for predicting equivalent Young's and Shear moduli of a laminate composite. The mechanical properties was calculated using Chamis model and an Objectif function with five norms, which allowed the mechanical properties to be examined and ultimate answer to be predicted. Moreover, another model, an alternative concept of equivalent lamina elements (ELEs), was first established using Chamis model for hybrid composites. Next, the ELEs were laminated in the direction and integrated into the compliance matrix for each ply. Finally, the results obtained with these models were compared to experimental results from other papers and journals, and some conclusions are shown.
The objective of the paper is to present new analytical models for analyzing the mechanical properties of laminate composites in all levels, i.e., lamina, laminate, and structure, from the perspectives of the material, orientation, and thickness. The contribution of all unit cells to the elastic properties of specimen was considered in the analytical model. In addition, the effects of the fiber orientation and thickness on the elastic properties were discussed in detail. Finally, some valuable conclusions are drawn.
An efficient recursive algorithm, the stiffness matrix method, has been developed for generally anisotropic multi-layered media. This algorithm has the computational efficiency and simplicity of the standard transfer matrix method and is unconditionally computationally stable for high frequency and layer thickness. In this algorithm, the stiffness (compliance) matrix is calculated for each layer and recursively applied to generate a stiffness (compliance) matrix for a layered system. This algorithm allows researchers and end users to calculate the stiffness and compliance matrix of any composite material, with any form, orientation, or thickness, by saving time and cost while finding the optimal solution.

Methodology
In this paper, we calculated the mechanical properties of laminate composites on three different levels and on three different categories. In one category, the variety of materials used and their homogenization that consisted of two different methods were discussed. In the other categories or paths, a variety of thicknesses and orientation, in which we are able to provide researchers the ultimate path and algorithm to calculate the stiffness matrix of laminated composites and to predict the mechanical properties, was discussed.
The ultimate path is called the K-PY algorithm, and it is shown in the conclusion part and based on the results found from comparing with experimental data from the literature ( Figure 1).

UD (unidirectional)-Homogenization:
Two new analytical models were proposed to predict the mechanical performance of laminate composites based on Chamis model [25]. The multi-scale modeling strategy was adopted to predict the elastic properties of the composites ( Figure 2). The KY analytical model consists totally of three basic steps. Firstly, the elastic performance of constituent components in lamina (or layer), i.e., the fiber and the pure resin matrix, is computed. Secondly, the stiffness matrix of each layer can be obtained by the stiffness-volume-averaging theory. Thirdly, an objective function is implemented to optimize the calculation and to find the optimal solution. Another model was investigated, which consists of finding equivalence of fibers in angle to be in the same direction 0° and to adopt the new fibers into Chamis and Sinclair's hybrid method [8].

UD (Unidirectional)-Homogenization:
Two new analytical models were proposed to predict the mechanical performance of laminate composites based on Chamis model [25]. The multi-scale modeling strategy was adopted to predict the elastic properties of the composites ( Figure 2). The KY analytical model consists totally of three basic steps. Firstly, the elastic performance of constituent components in lamina (or layer), i.e., the fiber and the pure resin matrix, is computed. Secondly, the stiffness matrix of each layer can be obtained by the stiffness-volume-averaging theory. Thirdly, an objective function is implemented to optimize the calculation and to find the optimal solution. Another model was investigated, which consists of finding equivalence of fibers in angle to be in the same direction 0 • and to adopt the new fibers into Chamis and Sinclair's hybrid method [8].

UD (unidirectional)-Homogenization:
Two new analytical models were proposed to predict the mechanical performance of laminate composites based on Chamis model [25]. The multi-scale modeling strategy was adopted to predict the elastic properties of the composites ( Figure 2). The KY analytical model consists totally of three basic steps. Firstly, the elastic performance of constituent components in lamina (or layer), i.e., the fiber and the pure resin matrix, is computed. Secondly, the stiffness matrix of each layer can be obtained by the stiffness-volume-averaging theory. Thirdly, an objective function is implemented to optimize the calculation and to find the optimal solution. Another model was investigated, which consists of finding equivalence of fibers in angle to be in the same direction 0° and to adopt the new fibers into Chamis and Sinclair's hybrid method [8].

Algorithm Objectif Function
For laminate composites, every kind of unit cell actually includes two types of constituent components, i.e., the fiber and the pure resin matrix. First, the compliance matrix of the lamina was calculated. Similar to many micro-mechanical models, the fibers and the Eng 2021, 2 421 resin matrix assume transversely isotropic, and both of them assume to be linearly elastic in the model. To obtain the elastic properties of the laminate composites, the fiber-matrix ROM proposed by Chamis was chosen to compute the engineering elastic constants. The Chamis micro-mechanical model is the most widely used and trusted model, which gives an equation for all five independent elastic properties as follows [25]: where V f is the fiber volume fraction; E f 11 is the Young's elastic modulus of the fiber in principle axis 1; E f 22 is the Young's elastic modulus of the fiber in principle axis 2; G f 12 is the longitudinal shear modulus of the fiber; G f 23 is the transverse shear modulus of the fiber; ν f 12 is the primary Poisson's ratio of the fiber; E m , ν m , and G m represent the Young's elastic modulus, Poisson's ratio, and shear modulus of the matrix, respectively.
The compliance matrix of the fiber [S] can be calculated easily. For the pure resin matrix pocket, it is generally regarded as isotropic materials. Once the elastic modulus and Poisson's ratio of the resin matrix are given, it is easy to determine its constitutive relationship.
It is easy to obtain their stiffness matrix in the local material coordinate systems (1, 2, 3) to (x, y, z) by inverting their compliance matrix [Sij] [7].
Then, we can obtain the stiffness in the global coordinate system by transforming the stresses and strains via the generalized transformation matrix as the following form: where l i , m i , and n i (i = 1, 2, 3) are the direction cosines and are defined as the cosines of the angle between the axes of the local and global coordinate systems before and after rotation: Each of the matrices is evaluated in ply k, and then the contributions of all the plies n are added.

• Objective function
The objective function indicates how much each variable contributes to the mechanical properties to be optimized. It takes the following general form: Norm 0: Norm ∞: Norm Cartesian: Norm 3: The second-order equation is shown as: ∆ < 0 is represented by the Cauchy-Schwarz inequality [26], so Norm 3 is not taken into consideration.

Hybrid Method
This section presents a review of the current status of hybrid composite materials technology, in terms of materials available and properties, and an outline of some of the trends, obvious and speculative, with emphasis on analytical models (Figure 3) including some details of smart hybrid composites. The equations are a simple extension of those applied for unidirectional single-fiber composites [8]. HC properties can be estimated by splitting it on PCs and SCs [13].

•
Chamis model for Hybrid Composites: As mentioned earlier, for the Young's modulus (E1) and the in-plane Poisson ratio (ν12), the ROM is able to provide accurate results, even for hybrid composites [7][8][9][10][11]. The equations are a simple extension of those applied for unidirectional single-fiber composites.
For the other engineering constants, Banerjee [7] extended the semi-empirical models proposed by Halpin [12].
The compliance matrix [S] for each fiber in its local coordinate system (1,2,3) was calculated. Using transformation matrix [T], the equivalent matrix for each fiber in the global coordinate system (x,y,z) was calculated [27]. The mechanical properties of equivalent fibers were obtained. Using Chamis model for hybrid composites, the mechanical properties of the whole lamina were calculated.
Eng 2021, 2 where subscripts m and fi designate the matrix and the fibers, respectively; E 2 is the transverse Young's modulus, G 12 is the in-plane shear modulus, V is the volume fraction, and ξ is a parameter associated with the fiber geometry (for circular fibers, it is 1.165 and 1.01). It is noted that the fiber volume fraction should be divided into a number of plies so that their sum should be equal with the matrix volume fraction (the percentage of fiber volume in the entire volume of a fiber-reinforced composite material) up to 100%.

Principal Globalization
Generally, a laminate does not consist only of unidirectional laminae because of their low stiffness and strength properties in the transverse direction. Therefore, in most laminates, some laminae are placed at an angle ( Figure 4) It is thus necessary to develop the stress-strain relationship for an angle lamina. The coordinate system is used for showing an angle lamina. The axes in the 1-2 coordinate system are called the local axes or the material axes. The direction 1 is parallel to the fibers, and the direction 2 is perpendicular to the fibers. In some references, direction 1 is also called the longitudinal direction L, and the direction 2 is called the transverse direction T. The axes in the x-y coordinate system are called the global axes or the off-axes. The angle between the two axes is denoted by angle θ. chanical properties of the whole lamina were calculated. ₁ = ∑ ᶠⁱ ᵢ ᶠⁱ + ᵐVᵐ , ₂ = ₃ = ᵐ[ where subscripts m and fi designate the matrix and the fibers, respectively; E2 is the transverse Young's modulus, G12 is the in-plane shear modulus, V is the volume fraction, and ξ is a parameter associated with the fiber geometry (for circular fibers, it is 1.165 and 1.01) It is noted that the fiber volume fraction should be divided into a number of plies so that their sum should be equal with the matrix volume fraction (the percentage of fiber volume in the entire volume of a fiber-reinforced composite material) up to 100%.

Principal Globalization
Generally, a laminate does not consist only of unidirectional laminae because of their low stiffness and strength properties in the transverse direction. Therefore, in most laminates, some laminae are placed at an angle ( Figure 4) It is thus necessary to develop the stress-strain relationship for an angle lamina. The coordinate system is used for showing an angle lamina. The axes in the 1-2 coordinate system are called the local axes or the material axes. The direction 1 is parallel to the fibers, and the direction 2 is perpendicular to the fibers. In some references direction 1 is also called the longitudinal direction L, and the direction 2 is called the transverse direction T. The axes in the x-y coordinate system are called the global axes or the off-axes. The angle between the two axes is denoted by angle θ where [T] is called the transformation matrix.
The mechanical properties of a typical unidirectional lamina are severely limited in the transverse direction. If one stacks several unidirectional layers ( Figure 5) this may be an optimum laminate for unidirectional loads. However, for complex loading and stiffness requirements, this would not be desirable. This problem can be overcome by making a laminate with layers stacked at different angles for given loading and stiffness requirements. This approach increases the cost and weight of the laminate, and thus, it is necessary to optimize the ply angles. Moreover, layers of different composite material systems may be used to develop a more optimum laminate. 5) this may be an optimum laminate for unidirectional loads. However, f plex loading and stiffness requirements, this would not be desirable. Th lem can be overcome by making a laminate with layers stacked at differen for given loading and stiffness requirements. This approach increases the c weight of the laminate, and thus, it is necessary to optimize the ply angles over, layers of different composite material systems may be used to de more optimum laminate. The idea is simply to find the equivalent lamina (or fiber) in order to to have a stack of laminae in the same direction. This equivalent can be f two different ways ( Figure 6). Algorithm:  The idea is simply to find the equivalent lamina (or fiber) in order to be able to have a stack of laminae in the same direction. This equivalent can be found in two different ways ( Figure 6).
Algorithm: lem can be overcome by making a laminate with layers stacked at different ang for given loading and stiffness requirements. This approach increases the cost a weight of the laminate, and thus, it is necessary to optimize the ply angles. Mo over, layers of different composite material systems may be used to develop more optimum laminate. The idea is simply to find the equivalent lamina (or fiber) in order to be ab to have a stack of laminae in the same direction. This equivalent can be found two different ways ( Figure 6). Algorithm:

Thickness
Size effects influence not only the strength of individual fibers, but also the longitudinal tensile strength and failure process of composite specimens and structures [28,29]. Most authors agree that the statistics of fiber strength are essential for establishing the relationship between the composite longitudinal tensile strength and size; in addition, other sources of size effects [30][31][32][33][34][35] are often considered as well. Size effects pose a challenge for the design of large composite structures based on experimental data measured in small coupons, and predictive models are therefore of utmost importance. This section suggests a generalized method for equivalence in thickness of laminate composites (Figure 7). relationship between the composite longitudinal tensile strength and size; in addition, other sources of size effects [30][31][32][33][34][35] are often considered as well. Size effects pose a challenge for the design of large composite structures based on experimental data measured in small coupons, and predictive models are therefore of utmost importance. This section suggests a generalized method for equivalence in thickness of laminate composites (Figure 7).

GCD of Plies
This idea is simply based on finding the GCD for the thicknesses in composite materials in order to unify the thicknesses provided. The GCD with two or more integers, which are not all zero, is the largest positive integer that divides each of the integers . By calculating the GCD, the number of plies is equal to the initial thickness, and thus, we replaced the initial thickness with a number of plies with a uniform thickness, even if the composite materials have different types of fibers or matrices.

Stress-Strain Relation (Hook's Law Method)
A mechanical structure takes external forces, which act upon a body as surface forces (for example, bending a stick) and body forces (for example, the weight of a standing vertical telephone pole on itself). These forces result in internal forces inside the body. Knowledge of the internal forces at all points in the body is essential because these forces need to be less than the strength of the material used in the structure. Stress, which is defined as the intensity of the load per unit area, determines this knowledge because the strengths of a material are intrinsically known in terms of stress. The determination of stress distributions of beams in necessary for determining the level of performance for the component. In particular, stress-based failure theories require determination of the maximum combined stresses in which the complete stress state must be either measured or calculated. Consider a bar of cross-section A and length L (Figure 8). A uniform tensile load P is applied to the two ends of the rod; by finding the state of stress, and strain energy per unit volume of the body. Assume that the rod is made of a homogeneous isotropic material of Young's modulus, E. Normal Stress: Having derived the proportionality relation for strain, εₓ, in the x-direction, the variation of stress, σₓ

GCD of Plies
This idea is simply based on finding the GCD for the thicknesses in composite materials in order to unify the thicknesses provided. The GCD with two or more integers, which are not all zero, is the largest positive integer that divides each of the integers. By calculating the GCD, the number of plies is equal to the initial thickness, and thus, we replaced the initial thickness with a number of plies with a uniform thickness, even if the composite materials have different types of fibers or matrices.

Stress-Strain Relation (Hook's Law Method)
A mechanical structure takes external forces, which act upon a body as surface forces (for example, bending a stick) and body forces (for example, the weight of a standing vertical telephone pole on itself). These forces result in internal forces inside the body. Knowledge of the internal forces at all points in the body is essential because these forces need to be less than the strength of the material used in the structure. Stress, which is defined as the intensity of the load per unit area, determines this knowledge because the strengths of a material are intrinsically known in terms of stress. The determination of stress distributions of beams in necessary for determining the level of performance for the component. In particular, stress-based failure theories require determination of the maximum combined stresses in which the complete stress state must be either measured or calculated. Consider a bar of cross-section A and length L (Figure 8). A uniform tensile load P is applied to the two ends of the rod; by finding the state of stress, and strain energy per unit volume of the body. Assume that the rod is made of a homogeneous isotropic material of Young's modulus, E. Normal Stress: Having derived the proportionality relation for strain, ε x , in the x-direction, the variation of stress, σ x ng 2021, 2, FOR PEER REVIEW Using generalized Hook's law: Since the whole idea is to find equivalent where: Since ε 11 = ε 11 (eq), ε 22 = ε 22 (eq),Y 12 =Y 12 (eq) (43) And From Equations (38)-(43), we have: Assuming that: Eng 2021, 2

Materials
In order to verify the applicability of this analytical model based on the equivalence in thickness, a model with two different ply is used in this work. A comparison between the results of present modeling, experimental data and previously developed multi-scale modeling. Four different examples are generated. As case studies, different mechanical properties are considered and the results are extracted.

UD-Homogenization
In order to verify the applicability of the analytical models, three types of examples were chosen from the specimens studied by Toledo and Menezes and their co-workers [32] to calculate the engineering elastic constants of laminate composites.
The laminates studied in this section are the following: Elastic properties of fibers and matrices are shown in [32,33].

Principal Globalization
A three-ply [0/90]s graphite/epoxy laminate with elastic properties as shown in Tables 1 and 2:

Equivalence in Thickness
The elastic properties of the materials are derived from Tables 3 and 4 [36]. Table 3. Elastic properties of fibers and matrices.

Results and Discussion
With the purpose of validating the developed modeling technique, the results were compared with experimental data published in open literature. A comparison between the results of present modeling, experimental data, and previously developed multi-scale modeling was conducted. As case studies, different mechanical properties were considered, and the results are shown in Tables 5-11, each of which shows a clear comparison between the experimental results and the presented algorithm.  For example, Tables 5-11 show the comparison of longitudinal and transversal Young's moduli and shear modulus between four different materials and fibers with different orientations and the two methods suggested in the UD-homogenization section of the K-PY algorithm, as well as a comparison implied between the two suggested methods, which are the 5 norms and the hybrid method.  Table 9. Ex Young's modulus comparison.  Table 11. Gxy Shear modulus comparison. On the other hand, Table 8 shows a clear comparison for the section of principal globalization and compares the results not only with experimental results with the literature but also those between the two suggested algorithms, i.e., the fiber method and the PLY method.

Type of Composite Thickness (mm) CLT Experimental GCD Method Hook's Law Method
Finally, Tables 9-11 show comparisons related to the thicknesses study methods with the classical laminate theory and the two suggested methods by the author, i.e., the GCD and the hook's law method.

UD-Homogenization
For all case studies, the fiber volume fraction was considered constant in all layers. The obtained results are presented in Tables 5-7. As it was predicted, norm 0 shows lower results among the five presented models, while norm ∞ was the highest. The Cartesian norm was the most accurate between all presented methods due to its direct relation to Chamis model, as discussed earlier. Norm 3 is not dependable, since its solution is always imaginary at ply number that exceeds two due to the Cauchy-Schwarz inequality [19], but it is notable that the real part is similar to the Cartesian norm because it is a second-order equation. Norm 4 and the hybrid method were not as efficient as the Cartesian norm.
It was observed that the Young's modulus in direction 1 is the most accurate parameter up to 0.1% (in comparison with that in the second experiment) for the longitudinal Young's modulus Ex, and that the maximum average error is shown for the Cartesian Norm model with 4.55%. For the transverse Young's modulus the percentage of error is much higher than for the longitudinal Young's modulus, with an average of 10.2% (the minimum percentage of error was for the second experiment with 0.38%). In addition, the average percentage of error for shear modulus in the Cartesian norm was 13.2% (with 2.24% as the lowest value from experiment 4).
Moreover, the percentage of errors with the hybrid method was higher, so their results are shown and not dependent on the ultimate path of the K-PY algorithm.
In addition to the high accuracy level of the current modeling technique, the very short runtime required of the modeling renders the developed multi-scale modeling as a cost-effective computational tool for estimating the Young's moduli of laminate composites.

Principal Globalization
Since the Cartesian norm was the most accurate among all presented methods, a clear comparison between the experimental values presented by Kaw [26] and the predicted results from the K-PY Equivalence model was implemented.
As shown in Table 8, the comparison was focused on Young's modulus Ex, Ey and Shear modulus Gxy. The fiber method present more accurate and precise values with the percentages of error varied from 0.128% for longitudinal Young's modulus Ex, 1.23% for transversal Young's modulus Ey, and 5% for Shear modulus Gxy, while the PLY method shows higher percentages of errors, such as 2.6% for Ex, 2.62% for Ey, and 6.15% for Gxy.

Equivalence in Thickness
The obtained results are presented in Tables 9-11.  Tables 9-11 give a clear comparison between the experimental values and the predicted results from the present K-PY Equivalence in thickness model in both GCD and Hook's law methods. The comparison was focused on Young's moduli Ex, Ey, and shear modulus Gxy. From Tables 9 and 10, good agreement was obtained between the predicted results and the available experimental values of the GCD method. The relative errors for the Young's moduli were 5.18% and 6.02%, respectively, while for Hook's law method they were around 5.52% and 4.03%, respectively. Moreover, the relative error for shear modulus Gxy was 2.1% for the GCD method, while for Hook's law method it was around 18.7%. As illustrated in Tables 9-11, as for the comparison of Young's modulus, it can be seen that the result were greater than the experimental values by about 5-6% for the GCD method while for Hook's law method it was greater than 4-5%. Similarly, the result for shear modulus Gxy was 2.1% for the GCD method while for Hook's law method it was greater than18%. From the aforementioned analysis, the contribution of K-PY Equivalence in thickness model in the GCD method can be taken into account for precisely predicting the elastic properties. The results showed that the present analytical model based on finding equivalent and with the application of thickness uniformity can be used to effectively evaluate the elastic properties for laminate composites.

Conclusions
The constitutive models presented in this paper are able to reproduce the behavior of laminated composites formed by fiber-reinforced laminae with different orientations or Eng 2021, 2 432 different types of fibers. The suggested hypotheses are the simplest that can be applied to isotropic fibers. As it was shown in the methodology section, applying the algorithm is easy and simple. Other models were also presented to discuss the equivalence of fibers or ply from the perspective of angle or thickness.
Good agreement was observed among numerical results, with smaller differences in the transversal strains, as shown in the results section that the Young's modulus in direction 1 was the most accurate parameter up to 0.1% (in comparison with that from the second experiment) for the longitudinal Young's modulus Ex and that the maximum average error is shown for the Cartesian norm model with 4.55% but with an average percentage of error of transversal strains of 10.2% (the minimum percentage of error was 0.38% for the second experiment); this can be due basically to the lack of precision in the mechanical properties of the components, fundamentally of the fibers, used as input data.
In general, application examples and comparisons with experimental results showed that the models are able to properly reproduce the mechanical behaviors of laminae and composite laminates. The models also describe the failure of the composite, taking into account what happens in each component. It also allows identifying the failure mode of the composite, produced by the failure of one or more of its components. It is able to reproduce complex failure modes that change from the matrix to the fibers as shown for equivalence in thicknesses and orientation where the percentages of errors were only up to 2% and less from results section in all predicted Young's and shear moduli, since only in equivalence in thicknesses the highest value of error was 5% for Gxy, depending on the type of stress state. The model works better for some types of laminates, and further analysis and model calibration are needed.
The models were covered in the paper, based directly on Chamis model with an objective function for optimization and hybrid-based model that discuss the equivalence of fiber in angle. The equivalence in orientation and in thickness was also discussed in this paper, and a new Generalization model was presented. The framework presented is potentially valid for other fiber conditions such as fiber volume fraction, but these cases would require further development.
Ultimately, it is concluded that the most accurate path to predict and reproduce the behavior of laminated composites in the orientation case is to apply the rotation/equivalence on fiber and then to apply K-PY model using the Objectif function with the Cartesian norm. While for equivalence in thickness, it was concluded that the GCD method is the most accurate method to use. As shown in the results, all values from Tables 5-11 had relatively low percentages of error, differing from 0.128% as shown in Table 8. The percentage is related to the lack of the exact and precise database structure of the used material, but it does show a great capability to use the algorithm to find the ultimate and optimal solution. This algorithm (Figure 9) is easily adaptable to laminates with periodicity, such as multiangle lay-up composites. This algorithm allows researchers and end users to calculate the stiffness and compliance matrix of any composite material, with any form, orientation, or thickness, saving so much time and cost while finding the optimal solution.
Optimal Path Algorithm: Eng 2021, 2, FOR PEER REVIEW 18 The models were covered in the paper, based directly on Chamis model with an objective function for optimization and hybrid-based model that discuss the equivalence of fiber in angle. The equivalence in orientation and in thickness was also discussed in this paper, and a new Generalization model was presented. The framework presented is potentially valid for other fiber conditions such as fiber volume fraction, but these cases would require further development.
Ultimately, it is concluded that the most accurate path to predict and reproduce the behavior of laminated composites in the orientation case is to apply the rotation/equivalence on fiber and then to apply K-PY model using the Objectif function with the Cartesian norm. While for equivalence in thickness, it was concluded that the GCD method is the most accurate method to use. As shown in the results, all values from tables 5-11 had relatively low percentages of error, differing from 0.128% as shown in Table 8. The percentage is related to the lack of the exact and precise database structure of the used material, but it does show a great capability to use the algorithm to find the ultimate and optimal solution. This algorithm ( Figure 9) is easily adaptable to laminates with periodicity, such as multi-angle lay-up composites. This algorithm allows researchers and end users to calculate the stiffness and compliance matrix of any composite material, with any form, orientation, or thickness, saving so much time and cost while finding the optimal solution.
Optimal Path Algorithm: