Layup Configuration Effect on Notch Residual Strength in Composite Laminates

The current trend shows an increasing demand for composites due to their high stiffness to weight ratio and the recent progress in manufacturing and cost reduction of composites. To combine high strength and stiffness in a cost-effective way, composites are often joined with steel or aluminum. However, joining of thermoset composite materials is challenging because circular holes are often used to join them with their metal counterparts. These design based circular holes induce high stress concentration around the hole. The purpose of this paper is to focus on layup configuration and its impact on notch stress distribution. To ensure high quality and uniformity, the holes were machined by a 5 kW continuous wave (cw) CO2 laser. The stress distribution was evaluated and compared by using finite element analysis and Lekhnitskii’s equations. For further understanding, the notch strength of the laminates was compared and strain distributions were analyzed using the digital image correlation technique.


Introduction
Fiber reinforced polymer (FRP) laminates are well known for their excellent stiffness and strength to weight ratio. As a result, FRP are increasingly used to replace metals in numerous fields of engineering where lightweight materials with good mechanical properties are required, such as automotive, aeronautics, astronautics, wind energy, and sports industries.
The laminates are often equipped with circular holes in order to join them with other components. Figure 1 shows design-based holes within composite structures used in the BMW 7 series for joining with metals. However, these design-based holes create high stress concentrations close to the notched area due to the interruption of the force flux in the fiber direction. These stress concentrations are often failure initiation regions in structural parts. Due to their good mechanical performance, FRP composites are increasingly used in the design of structural parts. Therefore, intensive study on the fatigue behavior of notched FRP laminates is [1][2][3][4] required. Many researchers have studied the effect of holes within FRP composites [5,6] but only a few with regard to the layup arrangement and its effects on the strength of the material. Woven fabric composites have gained their importance due to their bidirectional reinforcement in a single layer of fabric which provides more balanced properties in the fabric plane than unidirectional laminae. As well as their excellent drapability over complex geometries, low manufacturing costs and impact resistance quantify them as preferred structural materials. Thermoset composites offer excellent dimensional stability, receptiveness to paint, a lower thermal expansion coefficient compared to steel, high impact strength, while relatively low weights quantify them as functional components for automotive, aeronautics, astronautics, and in wind energy.
Many authors have investigated the performance on the notch stress concentration in FRP laminates. In ref. [7], it was identified that there is 30% strength reduction in notched composite, compared with unnotched composite. In ref. [8], the Waddoups, Eisenmann and Kaminski (WEK) model assumes that there is a high stress concentration around the hole in the FRP composite. In ref. [1,2,[4][5][6][7][8], many authors identified experimentally and numerically there is a high stress concentration around the edge of a hole. In ref. [9,10], notched stress on quasi-isotropic FRP laminates were performed and it was found that there is a development of an intense damaged zone near the hole edge. The presence of this stress concentration in a notched composite causes a wide range of damage and a failure mechanism when compared with an unnotched composite specimen. Thus, there is a necessity to identify the extensive full-field strain data to identify the failure behavior, rather than the limited experimental strength reduction values. However, very few studies have been made to identify the failure behavior using full-field strain data. In addition, there have been limited studies that explained the layup sequence and its influence in notch stress concentration. Thus, the aim of this paper is to analyze the strain development with respect to various layup architectures in notched glass/epoxy woven composites. The notch strength of the four woven fabric FRP laminates f [0°]4, f [45,0°]S, f [0,45°]S and f [0°]3 was investigated. By using Lekhnitskii's equation and finite element analysis, the stress concentration was analyzed. For further understanding, digital image correlation (DIC) technique was employed to obtain full-field surface strain measurements in differently notched laminates under tensile loading.

Material
Commercially available 220 g/m² non-crimp E-glass biaxial reinforcement and standard epoxy resin were used as matrix material for the FRP. The construction of biaxial non-crimp fabric (-45°/+45°) is shown in Figure 2. f [0°]4, f [45,0°]S, f [0,45°]S, and f [0°]3 laminates were manufactured using the vacuum assisted resin transfer molding (VARTM) process. The VARTM process was preferred because it allows the use of low-cost tooling for the production of high quality composite parts. For curing, the samples were kept under vacuum at ambient temperature of 24 °C for 24 h. Woven fabric composites have gained their importance due to their bidirectional reinforcement in a single layer of fabric which provides more balanced properties in the fabric plane than unidirectional laminae. As well as their excellent drapability over complex geometries, low manufacturing costs and impact resistance quantify them as preferred structural materials. Thermoset composites offer excellent dimensional stability, receptiveness to paint, a lower thermal expansion coefficient compared to steel, high impact strength, while relatively low weights quantify them as functional components for automotive, aeronautics, astronautics, and in wind energy.
Many authors have investigated the performance on the notch stress concentration in FRP laminates. In ref. [7], it was identified that there is 30% strength reduction in notched composite, compared with unnotched composite. In ref. [8], the Waddoups, Eisenmann and Kaminski (WEK) model assumes that there is a high stress concentration around the hole in the FRP composite. In ref. [1,2,[4][5][6][7][8], many authors identified experimentally and numerically there is a high stress concentration around the edge of a hole. In ref. [9,10], notched stress on quasi-isotropic FRP laminates were performed and it was found that there is a development of an intense damaged zone near the hole edge. The presence of this stress concentration in a notched composite causes a wide range of damage and a failure mechanism when compared with an unnotched composite specimen. Thus, there is a necessity to identify the extensive full-field strain data to identify the failure behavior, rather than the limited experimental strength reduction values. However, very few studies have been made to identify the failure behavior using full-field strain data. In addition, there have been limited studies that explained the layup sequence and its influence in notch stress concentration. Thus, the aim of this paper is to analyze the strain development with respect to various layup architectures in notched glass/epoxy woven composites.

Material
Commercially available 220 g/m 2 non-crimp E-glass biaxial reinforcement and standard epoxy resin were used as matrix material for the FRP. The construction of biaxial non-crimp fabric (-45 • /+45 • ) is shown in Figure 2  Plane strain ( ε 11 , ε 22 ) was calculated by performing tensile fatigue tests on non-notched specimens according to the standard ASTM D3039 [11]. Elastic modulus (E11 and E22) results presented in this section were obtained by testing five specimens for each layup and standard deviation (SD) between the samples is compiled in Table 1. Since the woven fabric was used to make FRP laminates, laminates were considered to be transversely isotropic. Thus, the elastic moduli in the x-y plane directions are equal (E11 = E22) and the laminate will only have five independent elastic constants in the elasticity tensor. Shear modulus was calculated experimentally using the V-Notched Rail Shear Method according to the standard ASTM D 7078/D 7078M-05 [12]. The dimensions used for tensile and shear tests, respectively, are shown in Figure 3. The tests were carried out at room temperature with a linear torsion all-electric dynamic test instrument (Instron ElectroPuls E10000, Darmstadt, Germany). The crosshead speed was maintained Plane strain (ε 11 , ε 22 ) was calculated by performing tensile fatigue tests on non-notched specimens according to the standard ASTM D3039 [11]. Elastic modulus (E 11 and E 22 ) results presented in this section were obtained by testing five specimens for each layup and standard deviation (SD) between the samples is compiled in Table 1. Since the woven fabric was used to make FRP laminates, laminates were considered to be transversely isotropic. Thus, the elastic moduli in the x-y plane directions are equal (E 11 = E 22 ) and the laminate will only have five independent elastic constants in the elasticity tensor. Shear modulus was calculated experimentally using the V-Notched Rail Shear Method according to the standard ASTM D 7078/D 7078M-05 [12]. The dimensions used for tensile and shear tests, respectively, are shown in Figure 3. Plane strain ( ε 11 , ε 22 ) was calculated by performing tensile fatigue tests on non-notched specimens according to the standard ASTM D3039 [11]. Elastic modulus (E11 and E22) results presented in this section were obtained by testing five specimens for each layup and standard deviation (SD) between the samples is compiled in Table 1. Since the woven fabric was used to make FRP laminates, laminates were considered to be transversely isotropic. Thus, the elastic moduli in the x-y plane directions are equal (E11 = E22) and the laminate will only have five independent elastic constants in the elasticity tensor. Shear modulus was calculated experimentally using the V-Notched Rail Shear Method according to the standard ASTM D 7078/D 7078M-05 [12]. The dimensions used for tensile and shear tests, respectively, are shown in Figure 3. The tests were carried out at room temperature with a linear torsion all-electric dynamic test instrument (Instron ElectroPuls E10000, Darmstadt, Germany). The crosshead speed was maintained The tests were carried out at room temperature with a linear torsion all-electric dynamic test instrument (Instron ElectroPuls E10000, Darmstadt, Germany). The crosshead speed was maintained at 2 mm/min, shear modulus (G 12 ), uniform slope between shear stress and shear strain (G 12 = ∆τ/∆ε). Engineering shear strain ε and shear stress τ are calculated using Equation (1).
where A is the area of the specimen and ε +45 , ε −45 are measured by digital image correlation (DIC) technique and the technique is explained in Section 3.4. GOM Correlate Evaluation Software was utilized for post processing. Three specimens were tested in the V-Notched Rail Shear Method and all the specimens were spray painted and a video was recorded. The average maximum force of three specimens was calculated. To reduce the post processing experiments, the specimen which had a close value to the average force was used for post processing. The corresponding shear force value (P) was identified from the experimental data. By measuring the uniform slope between shear stress and shear strain, the shear modulus (G 12 ) was calculated and is listed in the Table 1.

Preliminary Study on the Quality of Circular Holes
A preliminary study was conducted with quasi isotropic ([45 • , -45 • , 0 • , 90 • ] S ) carbon-epoxy composite laminate to choose the best method to produce circular holes with good quality in FRP laminates [2]. The holes were produced by drilling processes or laser methods. Introduction of the holes by cw-CO 2 lasers caused a reduction in maximum force of 27.94% compared to the unnotched specimen and drilling results a reduction of 26.06%, respectively. So, there were only slightly improved characteristic values obtained for drilling compared to cw-CO 2 laser process. The laser system can offer high cutting speeds, better flexibility, less tool wear, and less noise and dust compared to drilling. The laser system allows a high degree of automation compared to the milling machine because changing the drill bit for various hole diameters increases the processing time. In addition, tool wear increases delamination during drilling.
Hence a 3D robot assisted TRUMPF TruLaser Cell 7040 5 kW cw-CO 2 laser system was used for hole machining in this study ( Figure 4).

Lekhnitskii's Equation
The strength of the laminates (σ ∞ /σ O ) was calculated for various layup configurations with varying hole diameters (4, 6, and 8 mm). The infinite width notch tensile strength σ ∞ is calculated by multiplying the finite width notch tensile strength (σ ) by the finite width correction factor ( ∞ / ). σ O is the tensile strength of the corresponding laminates [13,14] without notch. Infinite width notch tensile strength and tensile strength of the unnotched laminates in the failure region are explained in Figure 5 and the corresponding values are calculated using the Lekhnitskii's equation. For an infinite width orthotropic plate containing a circular hole of radius R subjected to a uniform stress σ ∞ in the y direction, the normal stress σ , along the x-axis in front of the hole can be approximated by Equation (2) [ 7,[15][16]. Please check all the variables mentioned in the main text. The format (italic or not) should be unified through the main text.
where ∞ is the orthotropic stress concentration factor for a circular hole in an infinite width plate as expressed by Equation (2).
, , and are axial, transversal, and shear modules respectively. is the major Poisson's ratio. To modify equation 1 for a finite width plate, a correction factor is implemented and

Lekhnitskii's Equation
The strength of the laminates (σ N ∞ /σ O ) was calculated for various layup configurations with varying hole diameters (4, 6, and 8 mm). The infinite width notch tensile strength σ N ∞ is calculated by multiplying the finite width notch tensile strength (σ N ) by the finite width correction factor (K T ∞ /K T ).
σ O is the tensile strength of the corresponding laminates [13,14] without notch. Infinite width notch tensile strength and tensile strength of the unnotched laminates in the failure region are explained in Figure 5 and the corresponding values are calculated using the Lekhnitskii's equation.

Lekhnitskii's Equation
The strength of the laminates (σ ∞ /σ O ) was calculated for various layup configurations with varying hole diameters (4, 6, and 8 mm). The infinite width notch tensile strength σ ∞ is calculated by multiplying the finite width notch tensile strength (σ ) by the finite width correction factor ( ∞ / ). σ O is the tensile strength of the corresponding laminates [13,14] without notch. Infinite width notch tensile strength and tensile strength of the unnotched laminates in the failure region are explained in Figure 5 and the corresponding values are calculated using the Lekhnitskii's equation. For an infinite width orthotropic plate containing a circular hole of radius R subjected to a uniform stress σ ∞ in the y direction, the normal stress σ , along the x-axis in front of the hole can be approximated by Equation (2) [ 7,[15][16]. Please check all the variables mentioned in the main text. The format (italic or not) should be unified through the main text.
where ∞ is the orthotropic stress concentration factor for a circular hole in an infinite width plate as expressed by Equation (2).
, , and are axial, transversal, and shear modules respectively. is the major Poisson's ratio. To modify equation 1 for a finite width plate, a correction factor is implemented and For an infinite width orthotropic plate containing a circular hole of radius R subjected to a uniform stress σ ∞ in the y direction, the normal stress σ y , along the x-axis in front of the hole can be approximated by Equation (2) [7,15,16]. Please check all the variables mentioned in the main text. The format (italic or not) should be unified through the main text.
where K T ∞ is the orthotropic stress concentration factor for a circular hole in an infinite width plate as expressed by Equation (2).
E xx , E yy , and G xy are axial, transversal, and shear modules respectively. ϑ xy is the major Poisson's ratio. To modify Equation (1) for a finite width plate, a correction factor is implemented and for the composite laminates, equivalent modules for the laminates are replaced. The finite width correction (FWC) factor for orthotropic laminates with a circular hole is given in Equation (3). K T and K T ∞ are the stress concentration factors for finite and infinite width orthotropic plates, respectively.
Equation (5) shows the sample dimension related expression to calculate M, where D (=2R) is the diameter of the holes, W the width of the sample, and K T the orthotropic stress concentration factor for a circular hole in a finite width plate.

Finite Element Analysis (FEA)
FEA can improve the understanding of structural behavior and local failure mechanisms due to notches inside the plate. The normal stress distribution (NSD) depending on the geometry of the hole was executed using the commercially available ANSYS R18.1 software. NSD is the ratio of normal stress σ y (x = 0) along the x-axis in front of the hole to the stress applied on the laminate [17,18]. NSD was calculated using FEA and the results were compared with the values obtained from Lekhnitskii's equation. This helps to understand damage distribution and evolution of stress at the notch edge.

Digital Image Correlation (DIC)
DIC technique is a full field non-contact optical technique used to measure the surface deformation. DIC registers two (or more) images of the same structure subjected to various loadings (e.g., mechanical loading or thermal loading) to obtain full-field surface strain measurements of the entire specimen [17,[19][20][21]. Optical full-field measurement techniques are increasingly used in material testing, fracture mechanics, and vibration analysis. Samples are coated with black/white spray paints to create a random speckle pattern to obtain suitable surface contrast for DIC. Videos of the deformation were recorded with a digital single-lens reflex camera and the videos were processed with GOM Correlate Evaluation Software (2017) to analyze the strain distribution.

Strength of the Laminates
The strength of the laminates (σ N ∞ /σ O ) was calculated as explained in Section 3.2. The strength of the notched glass/epoxy laminates with different layup configurations are listed in Table 2.
All tensile tests were carried out according to ASTM D 5766/D 5766M standard [22]. For all calculations, the width of the laminates (W) was maintained at 36 mm and three different sizes of the hole (2R), 2R/W ratio were varied uniformly in all laminates.
Representative force-displacement curves of the unnotched laminates are shown in Figure 6.  Representative force-displacement curves of the unnotched laminates are shown in Figure 6. The results show that mixing the layup orientation ( f [0,45°]S and f [45,0°]S) result in a reduction in elastic modulus but at the same time in an improvement in the strength of the material.

Normal Stress Distribution
Finite element analysis was conducted under static structural analysis. A four noded element with six degrees of freedom at each node (SHELL181) was used because it allows the layered section with different orientations to be defined. A rectangular plate with a circular hole in the center as shown in Figure 1 was taken for calculating stresses induced by the applied load. The values of inplane Young's moduli (i.e., E1 and E2), Poisson's ratio (ʋ 12 ), and shear moduli (G12) were calculated using the experimental results and are shown in Table 1. The current investigation is not sensitive to the values in thickness direction (values are considered as identical in the 90° direction, ʋ 12 = ʋ 13 and G12 = G13). A plate with identical dimension asin the experiments was used in FEA. A typical finite element coarse mesh division was used to model the laminates.
To measure the stress along the x-axis in front of the hole, a path geometry was constructed to obtain the stress along that path. The calculated values were divided by the applied stress to calculate the NSD. The effect of a hole in the maximum principal stress distribution for a 4mm notched f [0°]4 laminate is shown in Figure 8. There is a substantial stress concentration at the hole edges in the X-Y direction which drops low at a certain distance from the hole.

Normal Stress Distribution
Finite element analysis was conducted under static structural analysis. A four noded element with six degrees of freedom at each node (SHELL181) was used because it allows the layered section with different orientations to be defined. A rectangular plate with a circular hole in the center as shown in Figure 1 was taken for calculating stresses induced by the applied load. The values of in-plane Young's moduli (i.e., E 1 and E 2 ), Poisson's ratio (υ 12 ), and shear moduli (G 12 ) were calculated using the experimental results and are shown in Table 1. The current investigation is not sensitive to the values in thickness direction (values are considered as identical in the 90 • direction, υ 12 = υ 13 and G 12 = G 13 ). A plate with identical dimension asin the experiments was used in FEA. A typical finite element coarse mesh division was used to model the laminates.
To measure the stress along the x-axis in front of the hole, a path geometry was constructed to obtain the stress along that path. The calculated values were divided by the applied stress to calculate the NSD. The effect of a hole in the maximum principal stress distribution for a 4mm notched f [0 • ] 4 laminate is shown in Figure 8. There is a substantial stress concentration at the hole edges in the X-Y direction which drops low at a certain distance from the hole.

Normal Stress Distribution
Finite element analysis was conducted under static structural analysis. A four noded element with six degrees of freedom at each node (SHELL181) was used because it allows the layered section with different orientations to be defined. A rectangular plate with a circular hole in the center as shown in Figure 1 was taken for calculating stresses induced by the applied load. The values of inplane Young's moduli (i.e., E1 and E2), Poisson's ratio (ʋ 12 ), and shear moduli (G12) were calculated using the experimental results and are shown in Table 1. The current investigation is not sensitive to the values in thickness direction (values are considered as identical in the 90° direction, ʋ 12 = ʋ 13 and G12 = G13). A plate with identical dimension asin the experiments was used in FEA. A typical finite element coarse mesh division was used to model the laminates.
To measure the stress along the x-axis in front of the hole, a path geometry was constructed to obtain the stress along that path. The calculated values were divided by the applied stress to calculate the NSD. The effect of a hole in the maximum principal stress distribution for a 4mm notched f [0°]4 laminate is shown in Figure 8. There is a substantial stress concentration at the hole edges in the X-Y direction which drops low at a certain distance from the hole. The results of Lekhnitskii could be validated with the results of FEA by calculating the stress concentration around the circular hole with respect to various layup arrangements. The results will give in-depth understanding of the stress concentration. As shown in Figure 9, Lekhnitskii and FEA The results of Lekhnitskii could be validated with the results of FEA by calculating the stress concentration around the circular hole with respect to various layup arrangements. The results will give in-depth understanding of the stress concentration. As shown in Figure 9, Lekhnitskii and FEA results have slight variations at the hole edges from 2.1 mm to 2.25 mm but in sequel both graphs are similar. This behavior is obvious in all these graphs. According to both techniques, the highest stress occurs at the edge of the hole. Similar to Section 3, in unsymmetric laminate [0] 3 the highest stress concentration occurs compared to all unsymmetric laminates. As explained in Section 3 there is a very high strength difference between the laminates [45,0] S and [0,45] S , therefore these two laminates were compared in more detail. As expected, [0,45] S laminate (strength reduction −14%) has low stress concentration compared to [45,0] S laminate (strength reduction −25%). But the difference in NSD was not proportional to the strength reduction. Further investigations on this strength reduction were required for a proper understanding, so digital image correlation (DIC) technique was used to analyze the damage mechanisms inside the laminates.

Strain Evolution
As explained in Section 3.4, all samples were prepared with a random speckle pattern using a black/white spray and the longitudinal strain (ε 22 ) field close to the notch was calculated from the DIC contrast. The strain fields were compared with the ones of other laminates. All strains were calculated around 1% prior to the ultimate failure of the laminates. In DIC, the subset size is a user defined parameter, which is critical for the measured accuracy of the samples [14]. To have a uniform result in all samples, a uniform subset size of 19 × 16 pixels was used for processing. As Figure 10 shows, the [45,0]S laminate (strength reduction −25%) has the highest strain concentration of all laminates. The strain within this laminate is in the range of 4%-5.5%, whereas in all other laminates strains in the range of 2%-2.5% occur. This could be the reason for the higher strength reduction in the laminate [45,0]S. In unsymmetric laminate, even though it has a lower strain range compared to

Strain Evolution
As explained in Section 3.4, all samples were prepared with a random speckle pattern using a black/white spray and the longitudinal strain (ε 22 ) field close to the notch was calculated from the DIC contrast. The strain fields were compared with the ones of other laminates. All strains were calculated around 1% prior to the ultimate failure of the laminates. In DIC, the subset size is a user defined parameter, which is critical for the measured accuracy of the samples [14]. To have a uniform result in all samples, a uniform subset size of 19 × 16 pixels was used for processing. As Figure 10 shows, the [45,0] S laminate (strength reduction −25%) has the highest strain concentration of all laminates. The strain within this laminate is in the range of 4%-5.5%, whereas in all other laminates strains in the range of 2%-2.5% occur. This could be the reason for the higher strength reduction in the laminate [45,0] S . In unsymmetric laminate, even though it has a lower strain range compared to [45,0] S , the strain is significantly higher than in all other laminates. So, this could be the reason for the lowest strength of this laminate compared to the symmetric laminates.

Conclusion
The effect of notch stress concentration in glass/epoxy woven composite with four different layup arrangements was investigated and compared. Analytical, numerical, and experimental test investigations were performed and analyzed. The results confirmed that the layup architecture has a significant influence on the notch strength. 5%.  Unsymmetric laminate shows high stress distribution at the edge of the hole compared to symmetric laminates.  All laminates show very high stress concentration near the edge of hole. Stress distribution is linear after 2 mm from the edge of hole. Until 2 mm, normal stress distribution increases exponentially towards the edge of the hole.

Acknowledgments:
The paper was created in the framework of the "LocPro-Local production by smart value chains cluster" in the Graduate Research School established by Brandenburg University of Technology (BTU) Cottbus-Senftenberg. The authors would like to thank BTU for their continuous support throughout the project.
Author Contributions: Venkateswaran Santhanakrishnan Balakrishnan, and Holger Seidlitz contributed to the design of the experiments, paper idea, experiments, and simulation. Sabine Weiß contributed to drafting and revising the paper and providing valuable help in structuring. I would like to assure that this review article has not been published previously and is not under consideration for publication elsewhere. I also confirm that the manuscript has been read and approved for submission by all the named authors.

Conflicts of Interest:
The authors declare no conflict of interest

Conclusions
The effect of notch stress concentration in glass/epoxy woven composite with four different layup arrangements was investigated and compared. Analytical, numerical, and experimental test investigations were performed and analyzed. The results confirmed that the layup architecture has a significant influence on the notch strength. The following behaviors were observed in this study, All laminates show very high stress concentration near the edge of hole. Stress distribution is linear after 2 mm from the edge of hole. Until 2 mm, normal stress distribution increases exponentially towards the edge of the hole.

Acknowledgments:
The paper was created in the framework of the "LocPro-Local production by smart value chains cluster" in the Graduate Research School established by Brandenburg University of Technology (BTU) Cottbus-Senftenberg. The authors would like to thank BTU for their continuous support throughout the project.
Author Contributions: Venkateswaran Santhanakrishnan Balakrishnan, and Holger Seidlitz contributed to the design of the experiments, paper idea, experiments, and simulation. Sabine Weiß contributed to drafting and revising the paper and providing valuable help in structuring. I would like to assure that this review article has not been published previously and is not under consideration for publication elsewhere. I also confirm that the manuscript has been read and approved for submission by all the named authors.

Conflicts of Interest:
The authors declare no conflict of interest