A Methodology to Predict the Fatigue Life under Multi-Axial Loading of Carbon Fiber-Reinforced Polymer Composites Considering Anisotropic Mechanical Behavior

Carbon fiber-reinforced polymers (CFRP) have been actively employed as lightweight materials; yet, evaluating the material’s reliability under multi-axis stress states is still challenging owing to their anisotropic nature. This paper investigates the fatigue failures of short carbon-fiber reinforced polyamide-6 (PA6-CF) and polypropylene (PP-CF) by analyzing the anisotropic behavior induced by the fiber orientation. The static and fatigue experiment and numerical analysis results of a one-way coupled injection molding structure have been obtained to develop the fatigue life prediction methodology. The maximum deviation between the experimental and calculated tensile results is 3.16%, indicating the accuracy of the numerical analysis model. The obtained data were utilized to develop the semi-empirical model based on the energy function, consisting of stress, strain, and triaxiality terms. Fiber breakage and matrix cracking occurred simultaneously during the fatigue fracture of PA6-CF. The PP-CF fiber was pulled out after matrix cracking due to weak interfacial bonding between the matrix and fiber. The reliability of the proposed model has been confirmed with high correlation coefficients of 98.1% and 97.9% for PA6-CF and PP-CF, respectively. In addition, the prediction percentage errors of the verification set for each material were 38.6% and 14.5%, respectively. Although the results of the verification specimen collected directly from the cross-member were included, the percentage error of PA6-CF was still relatively low at 38.6%. In conclusion, the developed model can predict the fatigue life of CFRPs, considering anisotropy and multi-axial stress states.


Introduction
Carbon fiber-reinforced polymer (CFRP) is widely utilized in various industries, such as automobiles and aerospace, for its high specific strength, corrosion resistance, and superior dimensional precision [1][2][3][4]. In addition, short fiber-reinforced composites manufactured through injection molding have noble productivity [5][6][7][8][9]. The material properties of CFRP can vary depending on the type of plastic. Thermosetting plastics have excellent mechanical properties and chemical resistance but have recycling difficulties and extended production time [10]. On the other hand, thermoplastic is flexible and reusable, but upon reheating, the physical properties can be reduced [11,12]. This study adopted polyamide-6 (PA6) and polypropylene (PP) thermoplastic resins as the CFRP matrix.
Despite CFRP's remarkable advantages, the fiber-reinforced polymer (FRP) has been over-designed in the last decades to compensate for its high possibility of debonding and breakage at the fiber-matrix interface [13]. Furthermore, the complex short fiber orientation formed by the polymer flow in the injection molding process generates severe anisotropy of CFRP, making it challenging to secure design reliability compared to conventionally used metals [14][15][16][17].
Therefore, fatigue life prediction considering the anisotropic behavior induced by the fiber orientation is essential for enhancing the durability of CFRP parts and their economic potential. The alignment of the fibers caused by the polymer flow has a fatal effect on fatigue strength [18,19]. In particular, short carbon fiber compounds have a lower fatigue strength than long fiber-reinforced composites due to their shorter nominal fiber length [20,21]. Moreover, if the fiber orientation and principal stress in a specific region susceptible to breakage are transverse to each other, the fiber can be easily separated from the matrix, resulting in reduced static and cyclic mechanical properties [22]. Nevertheless, the utilization of short fibers is advantageous in the injection molding process, for their ease of application. In addition, short fiber is more likely to be selected for producing recycled fiber. Therefore, there is a need to analyze the fatigue behavior of short fiberreinforced plastics.
Fiber orientation plays a crucial role in determining the anisotropy of injection-molded CFRP [23]. The surface layer near the cavity wall has a fiber orientation aligned in the direction of the polymer flow. In contrast, the fibers in the core region are aligned in the lateral direction of the flow [24][25][26][27]. However, for short FRP, the fibers are randomly distributed in the core region [28]. Still, the mechanical properties of the short FRP can be characterized from a macro-mechanical perspective using specific sample angles [29].
Many studies have utilized microcomputer tomography, experiments, and numerical simulations to analyze the anisotropic behavior of FRP [30][31][32][33][34][35]. For instance, Vincent et al. [30] developed an injection molding analysis model by computing fiber orientation tensors through two-dimensional optical image data. Meanwhile, studies on the durability characteristics of FRPs in complex stress states were also conducted [36][37][38][39]. In particular, Gude et al. proposed failure criteria through S-N curves obtained through multi-axis fatigue experiments of CFRP [36]. Furthermore, Tanaka et al. [39] conducted a study that simultaneously considered the notch effect and anisotropy by using single-edge notched plate specimens in different directions.
The fatigue fracture mechanisms of short fiber reinforced PA6 and PP have been investigated based on experimental and Scanning Electron Microscope (SEM) images [40][41][42][43]. PA6 fiber reinforcement showed local ductile fracture behavior at room temperature with little fiber pullout. In addition, cohesive deformation occured around the separated fibers due to the high interfacial bonding [40,41]. Meanwhile, the PP fiber reinforcement showed crack propagation due to fiber failure or pullout. Furthermore, the ductile fracture behavior is reduced as the fiber volume fraction is lowered [42,43].
As concerns the fatigue life prediction model, traditional analytical models, such as Manson-Coffin [44] and Smith Watson and Topper (SWT) [45], have been adopted to describe the anisotropy and complex stress effects [46][47][48][49][50][51][52][53][54]. Branco et al. [48] suggested the modified Manson-Coffin model capable of predicting fatigue life in complex stress states by adopting the equivalent strain energy density concept. Meanwhile, many semi-empirical models considering anisotropy have been developed by using angles or orientations in the model. In particular, Choi et al. [47] suggested the strain-based model, which accounts for the anisotropic behavior in multi-axial stress states by introducing an angle term of the specimens. In addition, Sonsino and Moosbrugger [55] predicted the fatigue behavior of glass FRP by calculating the mean stress and fiber orientation.
In summary, the durability of FRP can be described by calculating the orientation of individual fibers or by categorizing them into specific angles based on specimen cutting orientation. However, the above studies have limitations in that the fatigue life can be estimated only after the fiber orientation of the region of interest is known. In addition, few prior studies have evaluated anisotropy and various geometric effects as direct stress-strain relationships. Thus, there is a strong need to develop a methodology that can predict the cyclic fracture of CFRP based on numerical analysis and semi-experience models.
For the reasons mentioned above, a one-way coupled injection molding-structure numerical simulation and fatigue life model considering the anisotropy of CFRP in complex stress states have been developed. A Ramberg-Osgood flow stress model was selected to determine the anisotropic behavior. Model parameters were determined utilizing the tensile test results collected from different angles. In addition, the reduced strain closure (RSC) model was used to calculate the fiber orientation induced by the polymer flow, and the coefficient was derived via 3D X-ray CT (XCT) data. The fiber orientation mapping results derived from the finite volume method (FVM) were adopted as the initial material property of the finite element method (FEM). Finally, the FEM calculates the stress-strain and triaxiality. It was confirmed that all proposed procedures could be applied to PA6-CF and PP-CF and are not limited to specific materials.
Notch and complex shape specimens were collected from three angles (0 • = injection direction, 45 • , and 90 • ), and a fatigue test was conducted under various strain conditions. A semi-empirical fatigue life prediction model was developed based on the experimental and simulation results. The proposed model does not utilize fiber orientation angles or tensors in specific regions. Instead, the model developed through a combination of stress and strain terms predicts fatigue life based on differences in mechanical properties caused by the principal stress and fiber direction relationship. Accordingly, durability may be evaluated using only stress and strain generated by various load conditions of essential components in an operating environment. The proposed methodology allows designers to consider the anisotropic effect of CFRP components, assess vulnerability in advance, and pursue the convenience of varying the design without needing a prototype.

Materials and Methods
The fatigue life prediction of CFRP regarding anisotropy and multi-axial stress was performed through the procedure shown in Figure 1. First, static and cyclic experiments were conducted to characterize mechanical properties. In addition, the fiber orientation at 27 injection-molded plate regions was measured through XCT. The obtained image data were post-processed to calculate the short fiber orientation according to the polymer flow. For the reasons mentioned above, a one-way coupled injection molding-structure numerical simulation and fatigue life model considering the anisotropy of CFRP in complex stress states have been developed. A Ramberg-Osgood flow stress model was selected to determine the anisotropic behavior. Model parameters were determined utilizing the tensile test results collected from different angles. In addition, the reduced strain closure (RSC) model was used to calculate the fiber orientation induced by the polymer flow, and the coefficient was derived via 3D X-ray CT (XCT) data. The fiber orientation mapping results derived from the finite volume method (FVM) were adopted as the initial material property of the finite element method (FEM). Finally, the FEM calculates the stress-strain and triaxiality. It was confirmed that all proposed procedures could be applied to PA6-CF and PP-CF and are not limited to specific materials.
Notch and complex shape specimens were collected from three angles (0° = injection direction, 45°, and 90°), and a fatigue test was conducted under various strain conditions. A semi-empirical fatigue life prediction model was developed based on the experimental and simulation results. The proposed model does not utilize fiber orientation angles or tensors in specific regions. Instead, the model developed through a combination of stress and strain terms predicts fatigue life based on differences in mechanical properties caused by the principal stress and fiber direction relationship. Accordingly, durability may be evaluated using only stress and strain generated by various load conditions of essential components in an operating environment. The proposed methodology allows designers to consider the anisotropic effect of CFRP components, assess vulnerability in advance, and pursue the convenience of varying the design without needing a prototype.

Materials and Methods
The fatigue life prediction of CFRP regarding anisotropy and multi-axial stress was performed through the procedure shown in Figure 1. First, static and cyclic experiments were conducted to characterize mechanical properties. In addition, the fiber orientation at 27 injection-molded plate regions was measured through XCT. The obtained image data were post-processed to calculate the short fiber orientation according to the polymer flow. Secondly, the flow stress model coefficients for calculating the elastoplastic behavior were determined by reverse engineering the true stress-strain curve obtained in the tensile experiments. Furthermore, the RSC fiber orientation distribution parameter was derived from XCT data. Next, one-way coupled numerical analysis was modeled by combining the injection molding process and structural analysis utilizing the FVM and FEM. The injection molding parameters have been adopted based on the actual manufacturing process. Secondly, the flow stress model coefficients for calculating the elastoplastic behavior were determined by reverse engineering the true stress-strain curve obtained in the tensile experiments. Furthermore, the RSC fiber orientation distribution parameter was derived from XCT data. Next, one-way coupled numerical analysis was modeled by combining the injection molding process and structural analysis utilizing the FVM and FEM. The injection molding parameters have been adopted based on the actual manufacturing process.
Finally, a semi-experiential fatigue life prediction model in the form of an energy function was developed based on the stress-strain-triaxiality relationship by evaluating the anisotropy and the multi-axial stress state effects. The value of each term used in the model was calculated from the one-way coupled numerical analysis results. The stress and strain used in the model were calculated in the selected stabilized hysteresis loop. As a result, the fatigue life of CFRP can be successfully predicted through the developed procedure.

Materials and Experimental Methods
Short fiber injection molded PA6-CF and PP-CF materials were selected to evaluate the mechanical properties of CFRP under a multi-axial stress state. Tensile and displacement control fatigue tests are conducted to analyze static and cyclic behavior. The uniaxial and notched specimens were machined based on the ASTM D 638 02a TYPE IV standard [50]. The fracture behavior due to stress and strain concentration occurring in the notch can be analyzed using two types of notched geometry. The dimensions of complex specimens were specified based on the numerical analysis results such that the specimen covers a wide range of triaxialities. Detailed information on the dimensions of the specimens used is reported in Figure 2. As illustrated, the notch and complex specimens were designated as "Type I" and "Type II" depending on their geometry. The experiment and numerical analysis results were described using the same denotator.
Finally, a semi-experiential fatigue life prediction model in the form of an energy function was developed based on the stress-strain-triaxiality relationship by evaluating the anisotropy and the multi-axial stress state effects. The value of each term used in the model was calculated from the one-way coupled numerical analysis results. The stress and strain used in the model were calculated in the selected stabilized hysteresis loop. As a result, the fatigue life of CFRP can be successfully predicted through the developed procedure.

Materials and Experimental Methods
Short fiber injection molded PA6-CF and PP-CF materials were selected to evaluate the mechanical properties of CFRP under a multi-axial stress state. Tensile and displacement control fatigue tests are conducted to analyze static and cyclic behavior. The uniaxial and notched specimens were machined based on the ASTM D 638 02a TYPE IV standard [50]. The fracture behavior due to stress and strain concentration occurring in the notch can be analyzed using two types of notched geometry. The dimensions of complex specimens were specified based on the numerical analysis results such that the specimen covers a wide range of triaxialities. Detailed information on the dimensions of the specimens used is reported in Figure 2. As illustrated, the notch and complex specimens were designated as "Type I" and "Type II" depending on their geometry. The experiment and numerical analysis results were described using the same denotator. The verification specimen was obtained by partially cutting the injection-molded automobile cross-member. All specimens except the verification were machined on the injection molded plate at 0° (injection direction), 45°, and 90°. Furthermore, specimens were gathered 80 mm from the injection gate to exclude the effect of random fiber orientation due to unstable polymer flow near the injection gate. Figure 3 shows the direction and location of the machined specimen and the details of the verification specimen. The verification specimen was obtained by partially cutting the injection-molded automobile cross-member. All specimens except the verification were machined on the injection molded plate at 0 • (injection direction), 45 • , and 90 • . Furthermore, specimens were gathered 80 mm from the injection gate to exclude the effect of random fiber orientation due to unstable polymer flow near the injection gate. Figure 3 shows the direction and location of the machined specimen and the details of the verification specimen.
Tensile and fatigue tests were performed at room temperature (20 • C) and with a relative humidity of 50% utilizing the Instron 3367 and 8871 test systems, respectively. In addition, all tests were conducted after drying at 120 • C for 2 h in an oven dryer and cooling at room temperature for 30 min to remove moisture from the sample. In the case of ASTM D 638 and notched specimens, displacement was measured by an extensometer with a 25 mm gauge length. However, the complex specimens could not attach to the extensometer, so machine displacements have been employed. Tensile tests were carried out at a 10 mm/min tensile speed, and the parameters of the Ramberg-Osgood flow stress model described in Section 2.2 were derived using the obtained stress-strain curve. In addition, stress relaxation experiments were conducted to analyze the mechanical properties depending on the relaxing behavior. Fatigue tests were controlled by mean displacement and amplitude at 1 Hz frequency. Materials 2023, 16, x FOR PEER REVIEW 5 of 23 Tensile and fatigue tests were performed at room temperature (20 °C) and with a relative humidity of 50% utilizing the Instron 3367 and 8871 test systems, respectively. In addition, all tests were conducted after drying at 120 °C for 2 h in an oven dryer and cooling at room temperature for 30 min to remove moisture from the sample. In the case of ASTM D 638 and notched specimens, displacement was measured by an extensometer with a 25 mm gauge length. However, the complex specimens could not attach to the extensometer, so machine displacements have been employed. Tensile tests were carried out at a 10 mm/min tensile speed, and the parameters of the Ramberg-Osgood flow stress model described in Section 2.2 were derived using the obtained stress-strain curve. In addition, stress relaxation experiments were conducted to analyze the mechanical properties depending on the relaxing behavior. Fatigue tests were controlled by mean displacement and amplitude at 1 Hz frequency.

Numerical Analysis Methods
The fibers and matrix of the injection molded CFRPs are continuously damaged from exposure to repeated multi-axial stress. This phenomenon is intensified by complex stress states in the actual operating environment. Therefore, it is essential to accurately calculate the fiber orientation induced by the injection process and the stress-strain field in the fracture risk area. For this reason, a one-way coupled simulation of the injection process and the FEM was developed. The numerical analysis procedure is summarized in Figure 4.

Numerical Analysis Methods
The fibers and matrix of the injection molded CFRPs are continuously damaged from exposure to repeated multi-axial stress. This phenomenon is intensified by complex stress states in the actual operating environment. Therefore, it is essential to accurately calculate the fiber orientation induced by the injection process and the stress-strain field in the fracture risk area. For this reason, a one-way coupled simulation of the injection process and the FEM was developed. The numerical analysis procedure is summarized in Figure 4.  First, an input file of Abaqus 2017, without material property data, is generated containing geometry, solid mesh, and boundary conditions. Simultaneously, Python scripts are used to calculate Lamberg-Osgood flow stress parameters from stress-strain data of tensile experiments. The RSC coefficient is also determined using XCT data. After that, the injection process analysis through FVM is performed by Autodesk Moldflow Insight (AMI). The material properties due to fiber directions are then mapped to the input file via Advanced Material Exchange Helius (AME). Finally, the stress-strain field is calculated through structural analysis by FEM, and the fatigue life of the CFRP part is predicted. The entire procedure utilizing the three commercial software is verified by comparing the calculated fatigue life and analysis results with the actual experimental results.
This research applied the RSC model to account for the increased tensile strength of the CFRP when the directions of carbon fiber and principal stress are the same. The RSC model considers the difference in orientation kinetics according to the location of the injection gate and accurately calculates the short-fiber orientation in complex geometries First, an input file of Abaqus 2017, without material property data, is generated containing geometry, solid mesh, and boundary conditions. Simultaneously, Python scripts are used to calculate Lamberg-Osgood flow stress parameters from stress-strain data of tensile experiments. The RSC coefficient is also determined using XCT data. After that, the injection process analysis through FVM is performed by Autodesk Moldflow Insight (AMI). The material properties due to fiber directions are then mapped to the input file via Advanced Material Exchange Helius (AME). Finally, the stress-strain field is calculated through structural analysis by FEM, and the fatigue life of the CFRP part is predicted. The entire procedure utilizing the three commercial software is verified by comparing the calculated fatigue life and analysis results with the actual experimental results.
This research applied the RSC model to account for the increased tensile strength of the CFRP when the directions of carbon fiber and principal stress are the same. The RSC model considers the difference in orientation kinetics according to the location of the injection gate and accurately calculates the short-fiber orientation in complex geometries [56]. This model applies empirical factors and calculates fiber orientation by reducing the growth rate of the directional tensor without changing the rotation vector as follows: Here, a ij is the fiber orientation tensor, 0.5w ij the vorticity tensor, 0.5 . γ ij the deformation rate tensor, λ the particle shape coefficient, κ the scalar factor, L and M the fourth-order tensor, and C I the fiber interaction coefficient. The RSC factor between 0 and 1 is an empirical parameter, and the fiber alignment becomes similar to the polymer flow direction as the factor approaches 1.
The RSC parameters were optimized through reverse engineering by comparing fiber orientation calculated from XCT data with fiber orientation results from numerical simulations. Three steps are required to utilize the fiber orientation of the XCT data: (i) X-ray computed tomography, (ii) image data preprocessing through visualization, and (iii) fiber orientation calculation. First, XCT was conducted in three locations of the injection molded plate with dimensions of 2.5 × 2.5 × 1.96 mm 3 . Three different locations were selected to consider the diversity of fiber orientation. Secondly, for the preprocessing of the image data, the collected raw images have been reconstructed in 3D with lighting and transparency control and sectionalization with Dragonfly 2021. With the open_iA open-source tool for XCT datasets [57], the 3D data are converted to a combination of pixels. Finally, the fiber orientation tensor from the XCT data is calculated and compared with the Moldflow analysis results. The RSC parameters are updated until the difference in fiber orientation becomes less than 10% by the suggested procedure. The developed RSC parameter-determining procedure and the locations of the XCT photographing samples selected from the PA6-CF injection molded plates are described in Figure 5. With the open_iA opensource tool for XCT datasets [57], the 3D data are converted to a combination of pixels. Finally, the fiber orientation tensor from the XCT data is calculated and compared with the Moldflow analysis results. The RSC parameters are updated until the difference in fiber orientation becomes less than 10% by the suggested procedure. The developed RSC parameter-determining procedure and the locations of the XCT photographing samples selected from the PA6-CF injection molded plates are described in Figure 5. Injection process analyses were conducted using Moldflow commercial software with specified process variables and injection locations. Two models were developed using injection plates and automobile cross-member 3D-CAD step files. The simulations consisted of filling and packing processes, and the tetrahedra mesh was generated with 809,708 and 2,658,379 elements, respectively. A total of nine injection gate locations were defined concerning the location of the actual injection gates: one for the plate and eight for the cross-member. In addition, each process parameter used for the plate and the crossmember was designated in the process setting step. The variables related to temperature, Injection process analyses were conducted using Moldflow commercial software with specified process variables and injection locations. Two models were developed using injection plates and automobile cross-member 3D-CAD step files. The simulations consisted Materials 2023, 16,1952 7 of 21 of filling and packing processes, and the tetrahedra mesh was generated with 809,708 and 2,658,379 elements, respectively. A total of nine injection gate locations were defined concerning the location of the actual injection gates: one for the plate and eight for the cross-member. In addition, each process parameter used for the plate and the cross-member was designated in the process setting step. The variables related to temperature, filling, and packing used in the analysis are summarized in Table 1. Material properties were interpolated into the Abaqus 2017 input file considering the fiber orientation calculated in the injection process analysis using the Helius commercial program. The Ramberg-Osgood flow stress model and the modified Hill'48 yield function were used to evaluate the specimens' fiber orientation. (2) Here, E is an elastic modulus, K the strength coefficient, and n the hardening exponent. In addition, α and β are dependent parameters to λ I which define the first eigenvalue of the fiber orientation tensor. α m , β m and λ m,I are the reference values calculated in the Helius of the material.
The material parameters have been determined by minimizing the differences between experimental and numerical load-displacement curves, as shown in Table 2. Each specimen was modeled with 10,000 to 120,000 elements using the types C3D8R and C3D10, which are hexahedral and tetrahedral meshes. Fixed and displacement boundary conditions were assigned to apply experiment conditions. In addition, the rigid body constraint was given to the specimen part attached to the jig of the experimental machine. The boundary conditions adopted in injection and structure analysis are shown in Figure 6.
Each specimen was modeled with 10,000 to 120,000 elements using the types C3D8R and C3D10, which are hexahedral and tetrahedral meshes. Fixed and displacement boundary conditions were assigned to apply experiment conditions. In addition, the rigid body constraint was given to the specimen part attached to the jig of the experimental machine. The boundary conditions adopted in injection and structure analysis are shown in Figure 6.

Fatigue Life Prediction Model
Predicting the fatigue life of short CFRPs that undergo complex stress conditions is challenging. In the case of short FRPs, fibers are randomly distributed at the center in the

Fatigue Life Prediction Model
Predicting the fatigue life of short CFRPs that undergo complex stress conditions is challenging. In the case of short FRPs, fibers are randomly distributed at the center in the thickness direction of the product. Conversely, fibers are arranged with directionality on the product's surface, resulting in considerable anisotropy. In addition, fibers exhibit excellent strength in the longitudinal direction but low mechanical properties in the transverse direction.
Many previous studies have been conducted using the concept of local fiber orientation, geometry, or deformation concentration in a fatigue life model to predict the fatigue life of short CFRPs. In this study, a semi-empirical model based on the strain-stress-based energy function has been proposed. The interaction between the fiber orientation and the principal stress determines the durability of the fragile part. Therefore, the fatigue life prediction model has been developed based on strain-stress relationships induced by the macroscopic fiber orientation, even if fiber orientation itself is not included in the model.
The CFRP's semi-empirical fatigue life prediction model is developed based on the Manson-Coffin model, which relates the strain amplitude and the failure cycle [58].
This model investigated the relationship between the plastic strain amplitude ∆ε p and the fatigue fracture cycle N f . The plastic strain amplitude is calculated by the max and min plastic strain, ε p,max and ∆ε p,min . Furthermore, two material constants, A and c, have been utilized to predict the failure cycle.
In addition, the Manson-Coffin-Basquin model has introduced elastic deformation and failure stress to the Manson-Coffin Model [59].
This model investigated the relationship between the total strain amplitude ε a,t and the fatigue fracture cycle N f by dividing it into elasticity and plastic strain amplitude, ε a,e and ε a,p . The elastic strain amplitude can be figured with the fracture stress σ f divided elastic modulus E term with b, and the plastic strain amplitude can be calculated with the fracture strain ε f with c. b and c are the material constants.
The proposed model in this research has been developed by combining the above models and adding stress amplitudes and triaxiality terms.
The significant difference of the proposed model is that it consists of the energy function as the product of stress and strain amplitude. The fracture values at the static experiments have been utilized for the failure stress σ f ,i and strain ε f ,i in the model. The first term allows the adoption of the model for CFRPs with various matrices. CFRP shows a difference in strain depending on the strength and rigidity of the base material. This study confirmed whether the life model could consider the difference in physical properties according to the matrix, using PA6-CF and PP-CF.
The second term is the stress amplitude term, which is calculated as half the difference between the maximum and minimum equivalent stresses, σ max,i and σ min,i . In addition, by providing the maximum and failure stress to the exponential term in logarithmic form, the stress term is multiplied twice when the natural logarithm is taken on both sides. By multiplying the stress term twice, the varying stress differences induced by the directionality of the specimens can be effectively calculated. Thus, allowing consideration of the effect of anisotropy on fatigue life. Furthermore, the third term accounts for the relationship between the fracture strain and the triaxiality η i , enabling consideration of the complex stress state. All stress, strain, and triaxiality were extracted from the numerical analysis model. Fatigue model constants A and c were determined by regression.

Static and Fatigue Mechanical Properties
The tensile strength (TS), elastic modulus (E) and elongation at fracture (ε f ) have been summarized in Table 3. The true stress-strain curves obtained from uniaxial specimen tensile tests of PA6-CF and PP-CF materials are shown in Figure 7. In addition, the loaddisplacement curves of all specimens with FEM results have been reported in Section 3.2. Static and fatigue tests were conducted on both CFRPs by cutting specimens at 0 • (injection direction), 45 • , and 90 • . The static experiment results can confirm the effect of anisotropy and triaxiality on the material properties of the CFRPs. PA6-CF has higher tensile strength than PP-CF but breaks at a lower strain. All CFRPs showed high tensile strength at 0 • , followed by 45 • and 90 • . As the triaxiality increased, the fracture stress and strain decreased. All static test results were repeated three times, and then the median value was used. The true stress-time curves obtained from the stress relaxation tests of PA6-CF and PP-CF are shown in Figure 8. Stress losses occurred more rapidly to PA6-CF than to PP-CF. This rapid reduction can be inferred from the fact that the interfacial bonding between the matrix and the fiber is higher in PA6-CF, resulting in quicker stress dispersion. The strength loss was relatively low when both CFRP were 0 • than 90 • in the results according to the three angles.  The true stress-time curves obtained from the stress relaxation tests of PA6-CF and PP-CF are shown in Figure 8. Stress losses occurred more rapidly to PA6-CF than to PP-CF. This rapid reduction can be inferred from the fact that the interfacial bonding between the matrix and the fiber is higher in PA6-CF, resulting in quicker stress dispersion. The  The true stress-strain values of the uniaxial specimen were calculated directly using the experiment results, and for other specimens, the stress-strain was obtained from FEM. The fatigue test conditions and results are presented in Table 4 and Table 5, respectively. All values were extracted from a stabilized hysteresis loop with a difference of less than 5% from the failure stress-strain hysteresis loop. The true stress-strain values of the uniaxial specimen were calculated directly using the experiment results, and for other specimens, the stress-strain was obtained from FEM. The fatigue test conditions and results are presented in Tables 4 and 5, respectively. All values were extracted from a stabilized hysteresis loop with a difference of less than 5% from the failure stress-strain hysteresis loop.   SEM images were taken at PP-CF and PA6-CF fracture cross sections, as shown in Figure 9. Although both materials had carbon fibers with a volume fraction of 20%, they showed a clear difference in fracture characteristics. In PP-CF, many fibers were pulled out due to the weak interfacial bonding between the matrix and the fiber. Conversely, PA6-CF had a relatively high interfacial bonding showing cohesive failure without many fiber pullouts. The fracture behavior of PP-CF occurred in the order of matrix cracking, fiber breakage, and fiber pullout, while fiber breakage and matrix cracking occurred simultaneously for PA6-CF.

Numerical Analysis Results
The XCT data were obtained to optimize the RSC parameter, and the fiber orientations have been visually confirmed. XCT data of the PA6-CF injection molded plate is shown in Figure 10. Black, gray, and white represent the CFRP's void, matrix, and fiber, respectively. Figure 10a shows a full 3D image of single XCT data, and Figure 10b reveals the fiber distribution along the thickness and injection directions. The fibers in the surface layer are aligned with the injection direction. In contrast, the fibers are randomly distributed and appear as dots in the core layer. Figure 10c is an enlarged image of the surface layer fiber orientation viewed normal to the thickness direction.

Numerical Analysis Results
The XCT data were obtained to optimize the RSC parameter, and the fiber orientations have been visually confirmed. XCT data of the PA6-CF injection molded plate is shown in Figure 10. Black, gray, and white represent the CFRP's void, matrix, and fiber, respectively. Figure 10a shows a full 3D image of single XCT data, and Figure 10b reveals the fiber distribution along the thickness and injection directions. The fibers in the surface layer are Materials 2023, 16,1952 13 of 21 aligned with the injection direction. In contrast, the fibers are randomly distributed and appear as dots in the core layer. Figure 10c is an enlarged image of the surface layer fiber orientation viewed normal to the thickness direction. Figure 9. SEM images on fracture cross-sectional surface of (a) PP-CF and (b) PA6-CF.

Numerical Analysis Results
The XCT data were obtained to optimize the RSC parameter, and the fiber orientations have been visually confirmed. XCT data of the PA6-CF injection molded plate is shown in Figure 10. Black, gray, and white represent the CFRP's void, matrix, and fiber, respectively. Figure 10a shows a full 3D image of single XCT data, and Figure 10b reveals the fiber distribution along the thickness and injection directions. The fibers in the surface layer are aligned with the injection direction. In contrast, the fibers are randomly distributed and appear as dots in the core layer. Figure 10c is an enlarged image of the surface layer fiber orientation viewed normal to the thickness direction. The image is converted into coordinate data in pixel units through the Open_Ai open source tool. Three-dimensional data are stored as a three-dimensional coordinate system with a unique ID per pixel. The data of each pixel ID consist of Node ID and fiber vector. Additionally, fiber orientation vectors of Txx, Tyy, Tzz, Txy, Txz, and Tyz are calculated for each Node. The fiber orientation vector is defined with the fiber length l with orientation angles as follows: The image is converted into coordinate data in pixel units through the Open_Ai open source tool. Three-dimensional data are stored as a three-dimensional coordinate system with a unique ID per pixel. The data of each pixel ID consist of Node ID and fiber vector. Additionally, fiber orientation vectors of T xx , T yy , T zz , T xy , T xz , and T yz are calculated for each Node. The fiber orientation vector is defined with the fiber length l with orientation angles as follows: Figure 11a shows the 9 data extraction locations and compares the fiber orientation vector between the image and injection analysis results. Fiber orientation results show high accuracy with a 3.96%, 5.27%, and 4.38% difference in the principal direction. Fiber orientation tensor components were automatically calculated using fiber vectors at each point of the injection plate. The optimal RSC parameter was selected as 0.472, and the resulting fiber orientation tensor components are shown in Figure 11b. A total of 61 iterative calculations were made to generate the results.
The injection process simulations were performed using the obtained optimal RSC parameter and the proposed one-way coupled analysis method. The fiber orientation distribution results of the injection molding process analysis through FVM are shown in Figure 12. For the injection plate, the fiber distribution of the surface layers of both CFRPs is shown at the top. In addition, the results in the middle layer of the injection plate are shown at the bottom. The results of the injection cross-member are displayed by dividing them into top and bottom views on the right.
The fiber orientation of the middle layer is more consistent with the injection direction when compared to the surface layer. Since the center region is filled first and thus exposed to molten polymer flow the longest, the middle layer fiber orientation is closer to 1. In the surface layer, the overall fiber orientation value is 0.73. Meanwhile, PA6-CF showed a larger fiber orientation tensor than PP-CF. This phenomenon explains why the anisotropy of PA6-CF was more significant in static and fatigue tests. In addition, it was found that the rib segments of the cross-member were more susceptible to fracture due to their considerable difference in fiber orientation. Fiber orientation mapping was performed using AME to consider the anisotropy in structural analysis. Figure 13 shows the mapping results of the 90 o cases of each specimen and the cross-member.
The FEM model calculated stress, strain, and triaxiality based on the applied FVM fi-ber orientation results. The one-way coupled analysis was verified by comparing the loaddisplacement curve obtained from the static test with the results from the structural analysis. The triaxiality calculated through structural analysis is shown in Tables 4 and 5 with fatigue experiment results. Furthermore, the comparison results for all specimens are reported in Figures 14-16, respectively. The minimum and maximum average deviation were calculated to be 1.04% and 3.16% by comparing the area integration of experiment-analysis curves, supporting the reliability of the analysis model. In addition, it was confirmed that each specimen's most prominent stress position and fracture position were identical. In conclusion, it was demonstrated that failure of CFRP parts due to stress concentration could be predicted through injection-structure one-way couped numerical analysis method.  Figure 11a shows the 9 data extraction locations and compares the fiber orientation vector between the image and injection analysis results. Fiber orientation results show high accuracy with a 3.96%, 5.27%, and 4.38% difference in the principal direction. Fiber orientation tensor components were automatically calculated using fiber vectors at each point of the injection plate. The optimal RSC parameter was selected as 0.472, and the resulting fiber orientation tensor components are shown in Figure 11b. A total of 61 iterative calculations were made to generate the results.    Figure 11a shows the 9 data extraction locations and compares the fiber orientation vector between the image and injection analysis results. Fiber orientation results show high accuracy with a 3.96%, 5.27%, and 4.38% difference in the principal direction. Fiber orientation tensor components were automatically calculated using fiber vectors at each point of the injection plate. The optimal RSC parameter was selected as 0.472, and the resulting fiber orientation tensor components are shown in Figure 11b. A total of 61 iterative calculations were made to generate the results.  1. In the surface layer, the overall fiber orientation value is 0.73. Meanwhile, PA6-CF showed a larger fiber orientation tensor than PP-CF. This phenomenon explains why the anisotropy of PA6-CF was more significant in static and fatigue tests. In addition, it was found that the rib segments of the cross-member were more susceptible to fracture due to their considerable difference in fiber orientation. Fiber orientation mapping was performed using AME to consider the anisotropy in structural analysis. Figure 13 shows the mapping results of the 90 o cases of each specimen and the cross-member. The FEM model calculated stress, strain, and triaxiality based on the applied FVM fiber orientation results. The one-way coupled analysis was verified by comparing the loaddisplacement curve obtained from the static test with the results from the structural analysis. The triaxiality calculated through structural analysis is shown in this paper with fatigue experiment results. Furthermore, the comparison results for all specimens are reported in Figure 14, Figure 15, and Figure 16, respectively. The minimum and maximum average deviation were calculated to be 1.04% and 3.16% by comparing the area integration of experiment-analysis curves, supporting the reliability of the analysis model. In addition, it was confirmed that each specimen's most prominent stress position and fracture position were identical. In conclusion, it was demonstrated that failure of CFRP parts due to stress concentration could be predicted through injection-structure one-way couped numerical analysis method.        The experiment and numerical analysis results showed very high agreement. The differences in area integral beneath the experiment-simulation curves by material, geometry, and angle are summarized in Table 6.

Fatigue Life Prediction
A semi-empirical fatigue life prediction model is developed based on the fatigue test results at varying angles of five different geometries. The model constants of the life model were derived from 29 datasets for both PA6-CF and PP-CF, and a verification set was constructed using the four test results of the complex specimen. In the case of PA6-CF, the fatigue life predictability of the actual injection product was reviewed through two types of cross-member parts. The model coefficients A and c were derived using the least squares method and are reported in Table 7. The regression line and stress-strain-triaxiality-based energy function-to-failure cycle graphs calculated by the obtained model constants are shown in Figures 17 and 18, respectively.

Discussion
To develop a fatigue life prediction model for short CFRP, experimental and numerical analyses of PA6-CF and PP-CF were conducted. Anisotropy and exposure to complex stress state significantly affected CFRP's static and cyclic mechanical properties. In particular, it was found that when a load was applied in a direction consistent with the fiber arrangement, the tensile strength significantly increased. In addition, higher triaxiality resulted in a more complicated fracture pattern and shorter fatigue life, making it difficult to predict the life of injection molded CFRPs. This study employed a triaxiality range of 0.333 to 0.655 to investigate the multi-axial stress state. Furthermore, it was confirmed that the triaxiality was higher when the load was applied in the opposite direction to the fiber arrangement with the same geometry specimen.
The effect of anisotropy caused by fiber orientation can be explained through the experimental results of the directionally machined specimens. In the case of static tests, the

Discussion
To develop a fatigue life prediction model for short CFRP, experimental and numerical analyses of PA6-CF and PP-CF were conducted. Anisotropy and exposure to complex stress state significantly affected CFRP's static and cyclic mechanical properties. In particular, it was found that when a load was applied in a direction consistent with the fiber arrangement, the tensile strength significantly increased. In addition, higher triaxiality resulted in a more complicated fracture pattern and shorter fatigue life, making it difficult to predict the life of injection molded CFRPs. This study employed a triaxiality range of 0.333 to 0.655 to investigate the multi-axial stress state. Furthermore, it was confirmed that the triaxiality was higher when the load was applied in the opposite direction to the fiber arrangement with the same geometry specimen.
The effect of anisotropy caused by fiber orientation can be explained through the experimental results of the directionally machined specimens. In the case of static tests, the highest tensile strength was shown when the load was applied similar to the injection direction. PA6-CF showed a more severe anisotropy than PP-CF. This phenomenon occurs because PA6-CF's fiber orientation is more consistent with the injection direction, as shown in the numerical analysis results.
One-way coupled injection-structure analysis was conducted using FVM and FEM. The FVM analysis model was developed using the RSC parameter optimized from XCT data to consider anisotropy accurately. A widely used Ramberg-Osgood model in CFRP was adopted to reflect the anisotropic effects of fiber orientation. In addition, the FVM results were mapped to the structural analytical model to provide the material property data for stress-strain calculation.
The FEM model was finally utilized to investigate stress, strain, and triaxiality and to evaluate stress concentrations. The load-displacement curves of the structural analysis for all specimens were compared with the experiments, and the maximum mean deviation was 3.16%, demonstrating the analytical model's reliability. Furthermore, the actual fracture locations and the calculated stress concentration points of the cross-member specimens were consistent, thus, indicating the possibility of applying the proposed one-way analytical model for complex-shaped injection products.
A developed fatigue life prediction model introduces a strain-stress-triaxiality-based energy function considering both the anisotropy and complex stress states. Although the model does not contain a direct fiber orientation term, fatigue life can be predicted using the fiber orientation-induced stress difference. This is because anisotropy due to fiber orientation causes a difference in generated stress.
Therefore, the proposed model multiplies the stress term twice by taking the natural logarithm on both sides of the equation. This stress term allows the influence of the fiber to be implied in the model. In addition, damage due to strain and stress concentration in the complex stress state was predicted by adding a term consisting of a triaxiality and fracture strain. The developed semi-empirical fatigue life prediction results showed correlation coefficients (R 2 ) of 0.9813 and 0.9797, respectively, showing PA6-CF and PP-CF, respectively, proving the model's reliability.
The proposed fatigue life prediction procedure is highly accurate and can be widely used in designing CFRP injection molding products of complex shapes. The proposed methodology's advantage lies in the fact that the micro-local fiber orientation does not have to be repeatedly analyzed when predicting the mechanical properties of injection molded products. In addition, this research presents the potential for designers to predict fatigue life without using prototypes through one-way coupled numerical analysis and lab-scale experiments. However, future research is still required to evaluate the proposed life prediction model further. By applying this procedure to design actual injection products and to a broader range of CFRP materials, a more universal and validated form of the proposed model could be derived.

Concluding Remarks
This study predicted the fatigue life considering the anisotropy and complex stress states of two CFRPs, PA6-CF and PP-CF. The concluding remarks are as follows: 1.
The tensile strength of a material increases as a load is applied in the same direction as the fiber orientation. In addition, the PA6-CF shows a more significant anisotropy than the PP-CF.

2.
A one-way coupled injection-structural analysis model was developed, and fiber orientation along the polymer flow was considered through XCT photography and the Ramberg-Osgood model.