1. Introduction
Numerous industries, including automotive, aeronautics, aerospace, construction, military, and sports equipment, are currently witnessing significant technological advancements in the use of composite materials. Composite materials are preferred due to their ability to reduce weight and, as a result, minimize energy consumption. They are also valued for their resistance to various environmental factors and their flexibility, which allows them to be shaped in different ways compared to metals [
1]. Among the various types of composites, thermoplastic matrix composites have gained popularity and continue to experience rapid growth. Despite their many advantages, however, it is often challenging to obtain parts that meet the expectations of manufacturers. Extreme precision is required during the implementation of composites, and the slightest defect in impregnation or manufacturing can cause significant deformations within the material [
2,
3]. As a result, its mechanical characteristics are weakened, making it imperative to reduce the deformations within the material. Many factors influence the properties of composite materials, such as the type of fibers and matrix used or the processing method. However, it is widely acknowledged that the fiber–matrix interface plays a critical role because stress transfer occurs through this interface, and good compatibility between the two parts is necessary. The quality of the interface is closely linked to the presence of a specific surface treatment known as sizing. The impact of sizing on the properties of glass fiber-reinforced composite materials has been extensively studied, and it is evident that the quality of the interface significantly affects the properties of the material. Currently, short carbon fiber-reinforced thermoplastic matrix composites have the best adhesion properties, while glass fiber composites still require improved adhesion properties.
The use of short fiber reinforced thermoplastic composites has gained significant attention in various industrial fields and applications, including transportation, sports, aerospace, and automotive industries. These materials must have superior mechanical and chemical properties, be malleable to form complex shapes and be lightweight. Many authors have investigated the mechanical behavior of polyamide 66 (PA66GFs) reinforced with glass fibers in previous studies [
4,
5,
6,
7,
8,
9,
10,
11,
12]. The properties of composite material are influenced by various factors, such as fiber type, matrix type, fiber volume fraction and orientation, interface, and fabrication process. Injection molding is a highly efficient and productive method for processing PA66GFs, and its mechanical properties depend heavily on this processing technique.
In addition, it has been shown that the mechanical properties of PA66GFs are primarily governed by three main variables: fiber characteristics (diameter, length, orientation, and content), molding conditions (temperature, pressure, and molding speed), and mechanical testing conditions (deformation speed, test temperature, and hygrothermal condition). Polyamide 66 reinforced by short glass fiber composites (PA66GFs) display anisotropic properties, and it has been found that the material properties such as Young’s modulus, fatigue strength, and tensile strength are significantly higher in the flow direction than perpendicular to the polymer flow direction when considering the material as a three-dimensional solid. This orientation dependency of PA66GFs mechanical properties makes it an anisotropic material, as highlighted by several studies [
13,
14,
15,
16,
17,
18]. Haut du formulaire.
Composite materials are known to experience damage accumulation, with stiffness reduction often beginning early in the fatigue life due to various types of damage (matrix cracking, fiber fracture, fiber pullout, and matrix fiber delamination). The type of failure that occurs is dependent on the material variables and test conditions. Moreover, the damage mechanisms can interact differently and have different predominance depending on these factors. In the case of short fiber reinforced PA66GF subjected to cyclic loading, the damage evolution occurs in three stages [
19,
20,
21,
22,
23]. The first stage corresponds to material softening and damage initiation, characterized by a significant reduction in stiffness during the initial cycles. The second stage involves coalescence and the propagation of micro-cracks that occurred during the first stage. This phase is characterized by interfacial fiber–matrix damage, and a relatively steady reduction in stiffness is observed, corresponding to the behavior accommodation of the material. In the final stage, fiber fractures occur with the appearance of macroscopic cracks, leading to the total failure of the material. Degradation mechanisms result in a rapid reduction in stiffness.
Mechanical components made of PA66GFs are frequently subjected to random or cyclic loading during their service life, which can result in irreversible damage [
24]. To ensure that these components meet the necessary requirements of safety, economy, and reliability, it is critical to evaluate their cyclic behavior and the evolution of fatigue damage. However, these properties are usually determined through expensive, complex, and time-consuming experimental tests. Therefore, numerical simulations should be employed to predict the evolution of damage and to provide a reliable estimate of fatigue life. Nevertheless, to enable numerical simulations, the development of simple and accurate fatigue damage models is essential.
The fatigue damage of composite materials has been extensively investigated and documented in the literature. In their work, Amjadi and Fatemi [
25] introduced a fatigue damage model that focuses on the critical plane approach for predicting the fatigue life of SGFR thermoplastics under tension-tension or tension-compression loading. Their model takes into account the effects of fiber orientation and mean stress, as well as temperature and frequency. They applied this damage model to a general fatigue model in order to incorporate these additional effects.
Several researchers have proposed different fatigue damage models for short glass fiber-reinforced polyamide composites. Toubal et al. [
21] developed an analytical model based on cumulative damage to predict damage evolution. The model used a new damage parameter determined via the final Young’s modulus [
20]. Chebbi et al. [
9,
26] also proposed a 1D fatigue damage model based on young modulus degradation for predicting damage evolution in short glass fiber reinforced PA66GFs composite materials, using a similar approach and modifying the phenomenological model proposed by Chaboche [
27]. Poursartip et al. [
28] developed a phenomenological model where the unidirectional damage growth rate is a power function of mean stress and stress magnitude, independent of damage when cycle loading reaches constant stress magnitude. Recently, Chen et al. [
29] proposed a novel damage model that integrates multiple deformation mechanisms using the modified Mori–Tanaka model and Transformation Field Analysis (MT–TFA) approach for predicting monotonic and cyclic stress–strain responses applied on short fibers reinforced polyamide composites. Liu et al. [
30] and Kawai [
31] used an expression for damage that varies with the compression and tension regime and is related to the damage itself and the amount of strain applied. Van Paepegem et al. [
32] proposed a fatigue damage model that uses the amplitude of stress instead of the strain amplitude, which was validated by experimental results prepared on a flat woven glass/epoxy composite. Taheri-Behrooz et al. [
33] presented multiple damage models, including four previously studied models for unidirectional stiffness degradation, upon which they based their own proposed models.
The use of Continuum Damage Mechanics (CDM) in conjunction with a fatigue predictor has been a common approach in the development of fatigue-damage formalism. Kachanov [
34] was the first to introduce the CDM theory to explain the process of material degradation, and it was later extended to include fatigue damage due to its effectiveness in modeling the behavior of composite materials.
In the formulation of fatigue damage models, the small strain assumption is often used instead of the hyperelastic approach. The assumptions can be integrated with single or multiple damage variables corresponding to isotropic or anisotropic damage models, respectively. Ladevèze et al. [
35] formulated a new phenomenological fatigue damage model using an anisotropic damage formulation and the CDM theory for laminate composite sheets in the plane stress assumptions. This model was later reformulated by Sedrakian et al. [
36] to account for the strain energy and facilitate finite element implementation.
Nouri et al. [
37] extended this method for predicting nonlinear cumulative fatigue damage in short glass fiber reinforced polyamide 66 under mechanical loading. The damage rate was assumed to be the sum of an exponential law and a power function to capture the three stages of damage typically observed in this material. Avanzini et al. [
38] used a simplified version of the model to predict the fatigue behavior of a PEEK-based short fiber reinforced composite, considering isotropic fatigue damage with one damage variable. In related work, the authors [
39] used a similar approach to model fatigue damage in short fiber reinforced PEEK, modifying the damage state variable to accurately describe the damage behavior. The fatigue damage model was implemented into a finite element code. Payan et al. [
40] modeled the behavior of laminated carbon/epoxy composites under static and cyclic loading with plane stress assumptions and an orthotropic elastic potential. This model was later used by Hochard et al. [
41] to develop a nonlinear cumulative damage model for woven-ply laminates in the framework of CDM subjected to cyclic and static loadings.
Wang et al. [
42] employed the CDM theory in combination with a hyperelastic approach to describe the fatigue life of carbon-filled natural rubbers as a function of nominal strain amplitude under cyclic loading. They derived the damage evolution using Ogden strain energy density (SED) [
43]. Ayoub et al. [
44,
45] improved upon the Continuum Damage Mechanics model proposed by Wang et al. [
42] by utilizing a generalized form of Ogden SED and developing the theoretical framework for multiaxial loading. They tested the model using fatigue tests on a styrene–butadiene rubber (SBR) material under both constant and variable amplitude loading conditions to validate its applicability. The proposed model requires three damage parameters in addition to the constitutive law parameters.
Calvo et al. [
46] proposed an uncoupled directional continuum damage model for biological soft tissues with fiber reinforcement using a hyperelastic approach. Their model considers two damage variables, one for the matrix and one for the fibers. Several damage models have been proposed recently [
47,
48,
49,
50] for metallic materials that can be extended for thermoplastic-based composites.
The objective of this paper is to investigate the fatigue damage behavior of glass fiber-reinforced polyamide 66 under both small and large deformations. To this end, a 3D anisotropic damage model is proposed for short glass fiber-reinforced polyamide based on the CDM theory. The model is developed using a hyperelastic finite deformation description of the elastic response and is an extension of the results reported by Chebbi et al. [
8] for transversely anisotropic hyperelasticity where only the quasi-static uniaxial tensile and bending material behavior is considered. In addition, it is also an extension of the results illustrated in Chebbi et al. [
9,
26] for 3D anisotropic fatigue damage, where only a 1D elastic fatigue damage model is considered based on the degradation of the elastic modulus to predict the damage growth in composite materials.
The proposed 3D anisotropic fatigue damage model for short glass fiber-reinforced polyamide is developed in the framework of CDM theory and includes internal variables that introduce scalar variables for the fibers and matrix damages. The model takes into account the orientation of fibers due to the injection process of composite materials and is used to study the orthotropic material behavior under fatigue loading. To validate the model, material parameters are first determined through tensile tests on PA66GF in longitudinal and transverse directions. The fatigue damage of the PA66GF composite is then studied through both numerical simulations and experimental tests. The numerical simulations are performed via the developed finite element code and are used for sensitivity analysis of the effect of model parameters on damage evolution. The fatigue experimental tests are used to validate the proposed fatigue damage model.
3. Experimental Study
The experimental investigation involved the use of polyamide 66, which was reinforced with short glass fibers at weight percentages of 10%, 20%, and 30%. The molding injection process was utilized to manufacture four grades of composites, namely, PA66GF00, PA66GF10, PA66GF20, and PA66GF30, with glass contents of 0%, 10%, 20%, and 30% based on polyamide 66. Rectangular plates with specific dimensions and ISO 527-2 tensile specimens were produced, as detailed in Chebbi et al. [
9]. The mechanical properties were characterized via tensile tests in accordance with the ISO 527-2 standard specification. The specimens were cut along both the longitudinal and transverse directions, which are parallel and perpendicular, respectively, to the direction of mold flow. Room temperature tensile tests were conducted using a 10 kN INSTRON testing machine equipped with a RUDLPH laser extensometer, with a crosshead speed of 1 mm/min.
The specimens were subjected to cyclic loading using a fatigue-testing machine that was constructed in-house, as described in Chebbi et al. [
9]. The adjustable crank-linkage mechanism was utilized to apply the alternating displacement (
) to one end of the specimen (point B) while the other end was clamped (point A). Consequently, the specimen was loaded as a clamped–clamped beam with one moving clamp (as shown in
Figure 1). The amplitude of the displacement (
) and the displacement ratio (
) were adjustable parameters. The force acting on the composite sample was measured using a force gauge connected to the sample fixture.
The fatigue tests were conducted at room temperature with a constant alternating displacement, using a cyclic frequency of 2 Hz for load cycles. The aim was to minimize the increase in temperature due to the self-heating phenomenon, as proposed by Handa et al. [
59]. The variation of the displacement applied to the specimen during the test, as well as the resultant force, were recorded instantaneously. Three minimum tests were performed for each variation of displacement between 8 mm and 12.5 mm.
4. Identification of the Anisotropic Hyperelastic Behavior of PA66GFs
The parameters of the energy potential (defined in Equations (16) and (20)), including
,
,
,
, and
, were identified via an inverse method with three experimental tests. Two tensile tests were conducted in the longitudinal and transverse directions for PA66GFs, and a tensile test was conducted for the pure matrix. A finite element (FE) model was developed based on the specimens used in the tensile tests. Load control was employed, and an error function was expressed as follows:
This error function is the sum of different contributions from loaded nodes estimated at different load levels.
,
are, respectively, the nodal displacements from the Finite Element simulation and from the experiments. The non-linear least squares minimization problems Equation (50) is solved via the Levenberg–Marquardt algorithm. The FE model uses a four-node quadrilateral shell element proposed by Dammak et al. [
60].
The material properties for the matrix, relative to the Yeoh model, are obtained in Ref. [
8] as
,
, and
.
With
,
, and
obtained, the experimental results of uniaxial tensile tests of PA66GFs along the longitudinal direction are used to obtain the material parameter
. The predicted composite responses are plotted in
Figure 2. A good agreement with the experimental data is observed under moderate deformation.
Finally, the material parameter
for the studied material was obtained by fitting the transverse uniaxial tensile test data as plotted in
Figure 3. The experimental data are in good agreement with predicted results within the applicable domain of the presented model.
The other two material parameters,
and
, in the anisotropic hyperelastic model cannot be determined without additional experimental tests, such as shear tests. As these tests have not yet been conducted, the material parameters are currently considered to be valueless.
The material properties for the longitudinal and transverse directions,
and
, were determined via a nonlinear least squares method and are summarized in
Table 1 for three different volume fractions.
The evolution of
and
can be calculated by second-order polynomial functions with a good correlation factor.
The values of the material properties
and
identified for the three different materials are plotted as a function of the fiber volume fraction,
, in
Figure 4. It is worthwhile to notice the difference between the longitudinal and transverse directions.
The final whole model parameters are summarized in
Table 2.
Young modulus in longitudinal and transverse directions can be expressed as follows:
These moduli are given in
Table 3 for PA66, PA66GF10, PA66GF 20, and PA66GF 30. The results obtained are in good agreement with those obtained through experimental characterization [
8]. It is interesting to note the difference between Young’s modulus in the longitudinal and transverse directions. This difference is induced by the orientation of the fibers in these directions.