Next Article in Journal
Microstructural Characteristics of High-Pressure Die Casting with High Strength–Ductility Synergy Properties: A Review
Next Article in Special Issue
StressLife: A Short-Time Approach for the Determination of a Trend S-N Curve in and beyond the HCF Regime for the Steels 20MnMoNi5-5 and SAE 1045
Previous Article in Journal
Development of Cemented Paste Backfill with Superfine Tailings: Fluidity, Mechanical Properties, and Microstructure Characteristics
Previous Article in Special Issue
Determination of Dimple Core Design Configurations for Sandwich Panel Enhancement Using Fuzzy-Hybrid Decision Making Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Mechanical Engineering, Sogang University, Seoul 04107, Republic of Korea
2
Polymer R&D Team, GS Caltex R&D Center, Daejeon 34122, Republic of Korea
*
Author to whom correspondence should be addressed.
Materials 2023, 16(5), 1952; https://doi.org/10.3390/ma16051952
Submission received: 31 January 2023 / Revised: 23 February 2023 / Accepted: 25 February 2023 / Published: 27 February 2023
(This article belongs to the Special Issue Fatigue Damage and Fracture Mechanics of Materials)

Abstract

:
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.

1. 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 fiber-reinforced 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.

2. 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.
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.

2.1. 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.
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.

2.2. 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 [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:
D a i j D t = 1 2 ( ω i k a k j a i k ω k j ) + 1 2 λ ( γ ˙ i k a k j + a i k γ ˙ k j 2 [ a i j k l γ ˙ k l + ( 1 κ ) ( L i j k l M i j m n a m n k l ) ] γ ˙ k l ) + 2 κ C I γ ˙ ( δ i j 3 a i j )
Here, a i j is the fiber orientation tensor, 0.5 w i j the vorticity tensor, 0.5 γ ˙ i j 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 mm3. 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.
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 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.
σ e f f = ( α σ 11 β σ 22 ) 2 + ( β σ 22 β σ 33 ) 2 + ( β σ 33 α σ 11 ) 2 + 6 [ ( σ 12 ) 2 + ( σ 23 ) 2 + ( σ 31 ) 2 ] 2
σ = E 1 / n K ( n 1 ) / n ε p , e f f 1 / n
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.
σ ( λ I ) = θ + [ ( α m θ ) ( λ m , I 1 / 2 ) ] ( λ I 1 / 2 ) ,   β ( λ I ) = θ + [ ( β m θ ) ( λ m , I 1 / 2 ) ] ( λ I 1 / 2 )
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.

2.3. 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].
f ( ε ) = ( ε p , max ε p , min ) 2 = Δ ε p = A ( N f ) c
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].
ε a , t = ε a , e + ε a , p = σ f E ( 2 N f ) b + ε f ( 2 N f ) c
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.
f ( ε , σ , η ) = ( ε max ε min 2 ε f , i ) ( σ max , i σ min , i 2 σ f , i ) ln ( σ f , i σ max , i ) ( η i ) ln ( 1 + ε f , i ) = A ( N f ) c
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.

3. Results

3.1. 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 load-displacement 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–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.
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.

3.2. 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:
v ˜ f = ( v f , x v f , y v f , z ) = l ( sin α cos β sin α sin β cos α )
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 90o 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 load-displacement curve obtained from the static test with the results from the structural analysis. The triaxiality calculated through structural analysis is shown in Table 4 and Table 5 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.

3.3. 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 Figure 17 and Figure 18, respectively.

4. 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 (R2) 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.

5. 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:
  • 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.
  • 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.
  • The semi-empirical fatigue life prediction model based on the energy function showed high prediction accuracy of PA6-CF and PP-CF with a correlation coefficient of 0.9813 and 0.9797, respectively.
  • The developed fatigue life prediction methodology allows designers to reduce trial and error when designing products of complex geometries composed of injection molded short CFRPs and effectively predict fatigue life.

Author Contributions

Conceptualization, N.K.; methodology, J.C.; software, H.L. (Hyungyil Lee); validation, J.C.; formal analysis, J.C.; investigation, J.C.; resources, H.L. (Hyungyil Lee).; data curation, J.C.; writing—original draft preparation, J.C.; writing—review and editing, J.C.; visualization, J.C.; supervision, N.K.; project administration, N.K.; funding acquisition, H.L. (Hyungtak Lee). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea), grant numbers 10083609, 20004983, and 20013060.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bidaine, B.; Adam, L.; Assaker, R.; Robert, G.; Moulinjeune, O.; Perret, L.; Charrier, P.; Taveau, D.; Chauvet, L.; Haouala, S.; et al. Fiber Reinforced Plastic Durability: Nonlinear Multi-Scale Modeling for Structural Part Life Predictions. Available online: https://saemobilus.sae.org/content/2019-26-0278/ (accessed on 30 January 2023). [CrossRef]
  2. Kim, J.H.; Jung, Y.H.; Lambiase, F.; Moon, Y.H.; Ko, D.C. Novel Approach toward the Forming Process of CFRP Reinforcement with a Hot Stamped Part by Prepreg Compression Molding. Materials 2022, 15, 4743. [Google Scholar] [CrossRef]
  3. Kim, J.H.; Song, B.K.; Min, K.J.; Choi, J.C.; Eun, H.S. Optimizing Heat Treatment Conditions for Measuring CFRP and GFRP Resin Impregnation. Materials 2022, 15, 8182. [Google Scholar] [CrossRef] [PubMed]
  4. Zhuge, P.; Tao, G.; Wang, B.; Jie, Z.; Zhang, Z. Effects of High Temperatures on the Performance of Carbon Fiber Reinforced Polymer (CFRP) Composite Cables Protected with Fire-Retardant Materials. Materials 2022, 15, 4696. [Google Scholar] [CrossRef] [PubMed]
  5. Mandell, J.F. Fatigue Behavior of Short Fiber Composite Materials, 2nd ed.; Elsevier Science Publishers B.V.: Amsterdam, The Netherlands, 1991; Volume 4. [Google Scholar]
  6. Phelps, J. Processing-Microstructure Models for Short- and Long-Fiber Thermoplastic Composites; University of Illinois at Urbana-Champaign: Urbana, IL, USA, 2009; p. 198. [Google Scholar]
  7. Qureshi, H.J.; Saleem, M.U.; Khurram, N.; Ahmad, J.; Amin, M.N.; Khan, K.; Aslam, F.; Al Fuhaid, A.F.; Arifuzzaman, M. Investigation of CFRP Reinforcement Ratio on the Flexural Capacity and Failure Mode of Plain Concrete Prisms. Materials 2022, 15, 7248. [Google Scholar] [CrossRef] [PubMed]
  8. Gillani, F.; Zahid, T.; Bibi, S.; Khan, R.S.U.; Bhutta, M.R.; Ghafoor, U. Parametric Optimization for Quality of Electric Discharge Machined Profile by Using Multi-Shape Electrode. Materials 2022, 15, 2205. [Google Scholar] [CrossRef]
  9. Gillani, F.; Khan, M.Z.; Shah, O.R. Sensitivity Analysis of Reinforced Aluminum Based Metal Matrix Composites. Materials 2022, 15, 4225. [Google Scholar] [CrossRef]
  10. Pan, Y.; Yan, D. Study on the Durability of GFRP Bars and Carbon/Glass Hybrid Fiber Reinforced Polymer (HFRP) Bars Aged in Alkaline Solution. Compos. Struct. 2021, 261, 113285. [Google Scholar] [CrossRef]
  11. Xian, G.; Guo, R.; Li, C.; Wang, Y. Mechanical Performance Evolution and Life Prediction of Prestressed CFRP Plate Exposed to Hygrothermal and Freeze-Thaw Environments. Compos. Struct. 2022, 293, 115719. [Google Scholar] [CrossRef]
  12. Zhou, P.; Tian, J.; Li, C.; Tang, Z. Comparative Study of Durability Behaviors of Thermoplastic Polypropylene and Thermosetting Epoxy Exposed to Elevated Temperature, Water Immersion and Sustained Bending Loading. Polymers 2022, 14, 2953. [Google Scholar] [CrossRef]
  13. Shao, Y.; Okubo, K.; Fujii, T.; Shibata, O.; Fujita, Y. Effect of Matrix Properties on the Fatigue Damage Initiation and Its Growth in Plain Woven Carbon Fabric Vinylester Composites. Compos. Sci. Technol. 2014, 104, 125–135. [Google Scholar] [CrossRef]
  14. Harris, B. Fatigue in Composites. Harris, B., Ed.; Woodhead Publishing Ltd.: Cambridge, UK, 2003; ISBN 9781626239777. [Google Scholar]
  15. Rezaei, F.; Yunus, R.; Ibrahim, N.A. Effect of Fiber Length on Thermomechanical Properties of Short Carbon Fiber Reinforced Polypropylene Composites. Mater. Des. 2009, 30, 260–263. [Google Scholar] [CrossRef]
  16. Shao-Yun, F.; Bernd, L. Effects of Fiber Length and Fiber Orientation Distributions on the Tensile Strength of Short-Fiber-Reinforced Polymers. Compos. Sci. Technol. 1996, 56, 1179–1190. [Google Scholar]
  17. Bernasconi, A.; Davoli, P.; Basile, A.; Filippi, A. Effect of Fibre Orientation on the Fatigue Behaviour of a Short Glass Fibre Reinforced Polyamide-6. Int. J. Fatigue 2007, 29, 199–208. [Google Scholar] [CrossRef]
  18. Fu, S.Y.; Lauke, B.; Mäder, E.; Yue, C.Y.; Hu, X. Tensile Properties of Short-Glass-Fiber- and Short-Carbon-Fiber-Reinforced Polypropylene Composites. Compos. Part A Appl. Sci. Manuf. 2000, 31, 1117–1125. [Google Scholar] [CrossRef]
  19. Noda, K.; Takahara, A.; Kajiyama, T. Fatigue Failure Mechanisms of Short Glass-Fiber Reinforced Nylon 66 Based on Nonlinear Dynamic Viscoelastic Measurement. Polymer 2001, 42, 5803–5811. [Google Scholar] [CrossRef]
  20. Meneghetti, G.; Ricotta, M.; Lucchetta, G.; Carmignato, S. An Hysteresis Energy-Based Synthesis of Fully Reversed Axial Fatigue Behaviour of Different Polypropylene Composites. Compos. Part B Eng. 2014, 65, 17–25. [Google Scholar] [CrossRef]
  21. Panin, S.V.; Bogdanov, A.A.; Eremin, A.V.; Buslovich, D.G.; Alexenko, V.O. Estimating Low-and High-Cyclic Fatigue of Polyimide-CF-PTFE Composite through Variation of Mechanical Hysteresis Loops. Materials 2022, 15, 4656. [Google Scholar] [CrossRef]
  22. Zago, A.; Springer, G.S. Fatigue Lives of Short Fiber Reinforced Thermoplastics Parts. J. Reinf. Plast. Compos. 2001, 20, 606–620. [Google Scholar] [CrossRef]
  23. Brunbauer, J.; Mösenbacher, A.; Guster, C.; Pinter, G. Fundamental Influences on Quasistatic and Cyclic Material Behavior of Short Glass Fiber Reinforced Polyamide Illustrated on Microscopic Scale. J. Appl. Polym. Sci. 2014, 131, 1–14. [Google Scholar] [CrossRef]
  24. Tseng, H.C.; Chang, Y.J.; Wang, T.C.; Hsu, C.H.; Chang, R.Y. Three Dimensional Predictions of Fiber Orientation for Injection Molding of Long Fiber Reinforced Thermoplastics. In Proceedings of the SPE ACCE Conference, Michigang, AT, USA, 11–13 September 2013; Volume 2, pp. 1051–1058. [Google Scholar]
  25. De Monte, M.; Moosbrugger, E.; Quaresimin, M. Influence of Temperature and Thickness on the Off-Axis Behaviour of Short Glass Fibre Reinforced Polyamide 6.6—Cyclic Loading. Compos. Part A Appl. Sci. Manuf. 2010, 41, 1368–1379. [Google Scholar] [CrossRef]
  26. Amorim, L.; Santos, A.; Nunes, J.P.; Viana, J.C. Bioinspired Approaches for Toughening of Fibre Reinforced Polymer Composites. Mater. Des. 2021, 199, 109336. [Google Scholar] [CrossRef]
  27. Klimkeit, B.; Nadot, Y.; Castagnet, S.; Nadot-Martin, C.; Dumas, C.; Bergamo, S.; Sonsino, C.M.; Büter, A. Multi-axial Fatigue Life Assessment for Reinforced Polymers. Int. J. Fatigue 2011, 33, 766–780. [Google Scholar] [CrossRef]
  28. Shaharuddin, S.I.S.; Salit, M.S.; Edi Syams, Z. A Review of the Effect of Moulding Parameters on the Performance of Polymeric Composite Injection Moulding. Turkish J. Med. Sci. 2006, 30, 23–34. [Google Scholar] [CrossRef]
  29. Mortazavian, S.; Fatemi, A. Effects of Fiber Orientation and Anisotropy on Tensile Strength and Elastic Modulus of Short Fiber Reinforced Polymer Composites. Compos. Part B Eng. 2015, 72, 116–129. [Google Scholar] [CrossRef]
  30. Vincent, M.; Giroud, T.; Clarke, A.; Eberhardt, C. Description and Modeling of Fiber Orientation in Injection Molding of Fiber Reinforced Thermoplastics. Polymer 2005, 46, 6719–6725. [Google Scholar] [CrossRef]
  31. Advani, S.G.; Tucker, C.L. The Use of Tensors to Describe and Predict Fiber Orientation in Short Fiber Composites. J. Rheol. 1987, 31, 751–784. [Google Scholar] [CrossRef]
  32. Kanters, M.J.W.; Douven, L.F.A.; Savoyat, P. Fatigue Life Prediction of Injection Moulded Short Glass Fiber Reinforced Plastics. Procedia Struct. Integr. 2019, 19, 698–710. [Google Scholar] [CrossRef]
  33. Magino, N.; Köbler, J.; Andrä, H.; Welschinger, F.; Müller, R.; Schneider, M. A Space-Time Upscaling Technique for Modeling High-Cycle Fatigue-Damage of Short-Fiber Reinforced Composites. Compos. Sci. Technol. 2022, 222, 109340. [Google Scholar] [CrossRef]
  34. Ivan, R.; Sorgato, M.; Zanini, F.; Lucchetta, G. Improving Numerical Modeling Accuracy for Fiber Orientation and Mechanical Properties of Injection Molded Glass Fiber Reinforced Thermoplastics. Materials 2022, 15, 4720. [Google Scholar] [CrossRef]
  35. Mortazavian, S.; Fatemi, A. Fatigue Behavior and Modeling of Short Fiber Reinforced Polymer Composites Including Anisotropy and Temperature Effects. Int. J. Fatigue 2015, 77, 12–27. [Google Scholar] [CrossRef]
  36. Gude, M.; Hufenbach, W.; Koch, I.; Protz, R. Fatigue Failure Criteria and Degradation Rules for Composites under Multi-axial Loading. Mech. Compos. Mater. 2006, 42, 443–450. [Google Scholar] [CrossRef]
  37. Russo, A.; Sellitto, A.; Curatolo, P.; Acanfora, V.; Saputo, S.; Riccio, A. A Robust Numerical Methodology for Fatigue Damage Evolution Simulation in Composites. Materials 2021, 14, 3348. [Google Scholar] [CrossRef]
  38. Zhang, X.; Zhu, H.; Lv, Z.; Zhao, X.; Wang, J.; Wang, Q. Investigation of Biaxial Properties of CFRP with the Novel-Designed Cruciform Specimens. Materials 2022, 15, 7034. [Google Scholar] [CrossRef]
  39. Tanaka, K.; Oharada, K.; Yamada, D.; Shimizu, K. Fatigue Crack Propagation in Short-Carbon-Fiber Reinforced Plastics Evaluated Based on Anisotropic Fracture Mechanics. Int. J. Fatigue 2016, 92, 415–425. [Google Scholar] [CrossRef]
  40. Mortazavian, S.; Fatemi, A. Fatigue of Short Fiber Thermoplastic Composites: A Review of Recent Experimental Results and Analysis. Int. J. Fatigue 2017, 102, 171–183. [Google Scholar] [CrossRef]
  41. Ahmadifar, M.; Benfriha, K.; Shirinbayan, M. Thermal, Tensile and Fatigue Behaviors of the PA6, Short Carbon Fiber-Reinforced PA6, and Continuous Glass Fiber-Reinforced PA6 Materials in Fused Filament Fabrication (FFF). Polymers 2023, 15, 507. [Google Scholar] [CrossRef] [PubMed]
  42. Belmonte, E.; De Monte, M.; Hoffmann, C.J.; Quaresimin, M. Damage Initiation and Evolution in Short Fiber Reinforced Polyamide under Fatigue Loading: Influence of Fiber Volume Fraction. Compos. Part B Eng. 2017, 113, 331–341. [Google Scholar] [CrossRef]
  43. Tanaka, K.; Kitano, T.; Egami, N. Effect of Fiber Orientation on Fatigue Crack Propagation in Short-Fiber Reinforced Plastics. Eng. Fract. Mech. 2014, 123, 44–58. [Google Scholar] [CrossRef]
  44. Coffin, L.F.J. A Study of the Effects of Cyclic Thermal Stresses on a Ductile Metal. TraIns. Am. Soc. Mech. Eng. 1954, 76, 931–950. [Google Scholar] [CrossRef]
  45. Smith, K.N.; Topper, T.H. A Stress-Strain Function for the Fatigue of Metals. J. Mater. 1970, 5, 767–778. [Google Scholar]
  46. Choi, J.; Quagliato, L.; Shin, J.; Kim, N. Investigation on the Static and Cyclic Anisotropic Mechanical Behavior of Polychloroprene Rubber (CR) Reinforced with Tungsten Nano-Particles. Eng. Fract. Mech. 2020, 235, 107183. [Google Scholar] [CrossRef]
  47. Choi, J.; Quagliato, L.; Lee, S.; Shin, J.; Kim, N. Multi-axial Fatigue Life Prediction of Polychloroprene Rubber (CR) Reinforced with Tungsten Nano-Particles Based on Semi-Empirical and Machine Learning Models. Int. J. Fatigue 2021, 145, 106136. [Google Scholar] [CrossRef]
  48. Branco, R.; Costa, J.D.; Antunes, F.V. Fatigue Behaviour and Life Prediction of Lateral Notched Round Bars under Bending-Torsion Loading. Eng. Fract. Mech. 2014, 119, 66–84. [Google Scholar] [CrossRef]
  49. Mi, C.; Huang, Z.; Wang, H.; Zhang, D.; Xiong, T.; Jian, H.; Tang, J.; Yu, J. An Energy-Based Method for Lifetime Assessment on High-Strength-Steel Welded Joints under Different Pre-Strain Levels. Materials 2022, 15, 4558. [Google Scholar] [CrossRef]
  50. Tomažinčič, D.; Klemenc, J. Estimate of Coffin–Manson Curve Shift for the Porous Alloy AlSi9Cu3 Based on Numerical Simulations of a Porous Material Carried Out by Using the Taguchi Array. Materials 2022, 15, 2269. [Google Scholar] [CrossRef] [PubMed]
  51. Lopez-Crespo, P.; Moreno, B.; Lopez-Moreno, A.; Zapatero, J. Study of Crack Orientation and Fatigue Life Prediction in Biaxial Fatigue with Critical Plane Models. Eng. Fract. Mech. 2015, 136, 115–130. [Google Scholar] [CrossRef]
  52. Ghorbani, H.; Chakherlou, T.N.; Taghizadeh, H. On the Estimation of Fatigue Life in Bolt Clamped Al-Alloy 2024-T3 Plates. Eng. Fract. Mech. 2016, 164, 74–92. [Google Scholar] [CrossRef]
  53. Chen, Q.; Guo, H.; Avery, K.; Su, X.; Kang, H. Fatigue Performance and Life Estimation of Automotive Adhesive Joints Using a Fracture Mechanics Approach. Eng. Fract. Mech. 2017, 172, 73–89. [Google Scholar] [CrossRef]
  54. Choi, J.; Choi, J.; Lee, K.; Hur, N.; Kim, N. Fatigue Life Prediction Methodology of Hot Work Tool Steel Dies for High-Pressure Die Casting Based on Thermal Stress Analysis. Metals 2022, 12, 1744. [Google Scholar] [CrossRef]
  55. Sonsino, C.M.; Moosbrugger, E. Fatigue Design of Highly Loaded Short-Glass-Fibre Reinforced Polyamide Parts in Engine Compartments. Int. J. Fatigue 2008, 30, 1279–1288. [Google Scholar] [CrossRef]
  56. Wang, J.; O’Gara, J.F.; Tucker, C.L. An Objective Model for Slow Orientation Kinetics in Concentrated Fiber Suspensions: Theory and Rheological Evidence. J. Rheol. 2008, 52, 1179–1200. [Google Scholar] [CrossRef]
  57. Fröhler, B.; Weissenböck, J.; Schiwarth, M.; Kastner, J.; Heinzl, C. Open_iA: A Tool for Processing and Visual Analysis of Industrial Computed Tomography Datasets. J. Open Source Softw. 2019, 4, 1185. [Google Scholar] [CrossRef]
  58. Manson, S.S. Fatigue: A Complex Subject—Some Simple Approximations. Soc. Exp. Stress Anal. 1964, 15, 1–23. [Google Scholar] [CrossRef] [Green Version]
  59. Basquin, O.H. The Exponential Law of Endurance Tests. In Proceedings of the American Society for Testing and Materials Proceedings, Atlantci City, NJ, USA, 28 June–2 July 1910; pp. 625–630. [Google Scholar]
Figure 1. The procedure of fatigue life prediction for CFRP.
Figure 1. The procedure of fatigue life prediction for CFRP.
Materials 16 01952 g001
Figure 2. Dimensional specifications of test specimens denoted in mm. (a) ASTM-D638 Type IV standard uniaxial specimen. (b) Notched specimens. (c) Complex specimens.
Figure 2. Dimensional specifications of test specimens denoted in mm. (a) ASTM-D638 Type IV standard uniaxial specimen. (b) Notched specimens. (c) Complex specimens.
Materials 16 01952 g002
Figure 3. (a) Direction and location of specimens cut from injection molded plates. (b) Locations and geometries of verification specimens in injection-molded automobile cross-member components.
Figure 3. (a) Direction and location of specimens cut from injection molded plates. (b) Locations and geometries of verification specimens in injection-molded automobile cross-member components.
Materials 16 01952 g003
Figure 4. The one-way coupled injection molding-structure numerical analysis procedure for fatigue life prediction of CFRPs.
Figure 4. The one-way coupled injection molding-structure numerical analysis procedure for fatigue life prediction of CFRPs.
Materials 16 01952 g004
Figure 5. (a) The procedure to determine the RSC parameter by comparison of fiber orientation tensor obtained from XCT data and Moldflow simulation results. (b) The locations of XCT samples in the PA6-CF injection molded plate.
Figure 5. (a) The procedure to determine the RSC parameter by comparison of fiber orientation tensor obtained from XCT data and Moldflow simulation results. (b) The locations of XCT samples in the PA6-CF injection molded plate.
Materials 16 01952 g005
Figure 6. Boundary conditions and simulation models. (a) The injection gate and tetrahedra mesh for the injection molding simulation. (b) Mechanical boundary conditions for structural analysis.
Figure 6. Boundary conditions and simulation models. (a) The injection gate and tetrahedra mesh for the injection molding simulation. (b) Mechanical boundary conditions for structural analysis.
Materials 16 01952 g006
Figure 7. True stress and strain curves of the PA6-CF and PP-CF materials along 0° (injection direction), 45°, and 90°.
Figure 7. True stress and strain curves of the PA6-CF and PP-CF materials along 0° (injection direction), 45°, and 90°.
Materials 16 01952 g007
Figure 8. True stress–time curves of the PA6-CF and PP-CF obtained from the stress relaxation tests.
Figure 8. True stress–time curves of the PA6-CF and PP-CF obtained from the stress relaxation tests.
Materials 16 01952 g008
Figure 9. SEM images on fracture cross-sectional surface of (a) PP-CF and (b) PA6-CF.
Figure 9. SEM images on fracture cross-sectional surface of (a) PP-CF and (b) PA6-CF.
Materials 16 01952 g009
Figure 10. The XCT images from PA6-CF injection molded plate. Fibers and matrices appear as white and grey elements. (a) Full 3D view of a single XCT data. (b) Surface and core layer along thickness and injection direction. (c) Surface layer view from normal to the thickness direction.
Figure 10. The XCT images from PA6-CF injection molded plate. Fibers and matrices appear as white and grey elements. (a) Full 3D view of a single XCT data. (b) Surface and core layer along thickness and injection direction. (c) Surface layer view from normal to the thickness direction.
Materials 16 01952 g010
Figure 11. (a) The analysis points on the injection molded plate and the fiber orientation vectors on point 7. (b) The fiber orientation tensor component results on points 3, 6, and 9.
Figure 11. (a) The analysis points on the injection molded plate and the fiber orientation vectors on point 7. (b) The fiber orientation tensor component results on points 3, 6, and 9.
Materials 16 01952 g011
Figure 12. Fiber orientation distribution of surface and middle layer of injection molded plate, and cross-member.
Figure 12. Fiber orientation distribution of surface and middle layer of injection molded plate, and cross-member.
Materials 16 01952 g012
Figure 13. The mapping results of the (a) uniaxial, (b) notched, (c) complex, and (d) cross-member specimens.
Figure 13. The mapping results of the (a) uniaxial, (b) notched, (c) complex, and (d) cross-member specimens.
Materials 16 01952 g013
Figure 14. Load-displacement curves (line: experimental data, scatter: FEM results). (a) PA6-CF and PP-CF materials for the uniaxial specimens. (b) Cross-member specimens.
Figure 14. Load-displacement curves (line: experimental data, scatter: FEM results). (a) PA6-CF and PP-CF materials for the uniaxial specimens. (b) Cross-member specimens.
Materials 16 01952 g014
Figure 15. Load-displacement curves of PA6-CF and PP-CF materials (line: experimental data, scatter: FEM results). (a) The notched Type I. (b) The notched Type II.
Figure 15. Load-displacement curves of PA6-CF and PP-CF materials (line: experimental data, scatter: FEM results). (a) The notched Type I. (b) The notched Type II.
Materials 16 01952 g015
Figure 16. Load-displacement curves of PA6-CF and PP-CF materials (line: experimental data, scatter: FEM results). (a) The complex Type I. (b) The complex type II.
Figure 16. Load-displacement curves of PA6-CF and PP-CF materials (line: experimental data, scatter: FEM results). (a) The complex Type I. (b) The complex type II.
Materials 16 01952 g016
Figure 17. Energy function vs. fatigue life graph with regression line of the PA6-CF (log-linear scale).
Figure 17. Energy function vs. fatigue life graph with regression line of the PA6-CF (log-linear scale).
Materials 16 01952 g017
Figure 18. Energy function vs. fatigue life graph with regression line of the PP-CF (log-linear scale).
Figure 18. Energy function vs. fatigue life graph with regression line of the PP-CF (log-linear scale).
Materials 16 01952 g018
Table 1. Process parameter setting values in the injection analysis for the plate and cross-member.
Table 1. Process parameter setting values in the injection analysis for the plate and cross-member.
Process ParametersPlateCross-Member
Mold surface temperature [°C]85
Melt temperature [°C]285
Switch-over by %volume filled [%]99
Injection time [s]220
Cooling time [s]20140
Table 2. Model constants for the Ramberg–Osgood flow stress model and a modified Hill’48 yield function.
Table 2. Model constants for the Ramberg–Osgood flow stress model and a modified Hill’48 yield function.
Materials K [MPa] n [−] α m [−] β m [−] λ m , I [−] E m [GPa] E f [GPa]
PA6-CF95.0711.561.771.140.854.92136.57
PP-CF18.2712.602.761.940.851.0182.55
Table 3. Mechanical properties of PA6-CF and PP-CF.
Table 3. Mechanical properties of PA6-CF and PP-CF.
MaterialSpecimen Angle [°]TS [MPa] E [GPa] ε f [−]
PA6-CF0201.914.70.024
45141.67.50.038
90125.16.80.030
PP-CF037.28.30.027
4531.17.10.032
9027.84.30.044
Table 4. Summary of fatigue experiments on the uniaxial, notched, complex, and cross-member specimens of PA6-CF.
Table 4. Summary of fatigue experiments on the uniaxial, notched, complex, and cross-member specimens of PA6-CF.
GeometryTriaxiality
[ η ]
Specimen Angle
[°]
ε max
[−]
ε min
[−]
( σ v M ) max
[MPa]
( σ v M ) min
[MPa]
Nf
[cycles]
Uniaxial0.33300.0120.003150.148.41800
0.0130.007157.7103.846,700
0.0170.007183.1103.8230
450.0110.00475.129.0>106
0.0130.00792.258.4428,500
0.0170.007112.558.416,100
900.0120.00380.925.129,800
0.0130.00786.853.5346,900
0.0170.007103.253.51790
Notched Type I0.40800.1160.087300.2275.883,400
0.1160.058300.2227.538,600
0.475450.1630.108210.4172.9134,000
0.1630.072210.4125.63280
0.481900.0910.046184.3106.441,000
0.0910.057184.3128.3304,000
Notched Type II0.39300.1110.069288.7210.630,900
0.1110.055288.7176.64400
0.469450.1430.067223.5164.24980
0.1430.072223.5138.15900
0.498900.1800.113228.3163.66200
0.1800.090228.3136.5860
Complex Type I0.55500.1320.093266.8228.910,400
0.1320.078226.8214.11720
0.612450.1590.114238.6207.629,600
0.1590.129230.3187.81450
0.641900.2150.158230.3206.542,300
0.2150.129230.3187.84200
Complex Type II0.54400.1070.071268.5233.216,400
0.1070.053268.5209.6980
0.602450.1390.092237.5207.924,000
0.1390.069237.5180.21480
0.655900.1230.082206.7176.814,500
0.1230.061206.7151.21060
Location I0.529-0.1720.094198.9122.488,700
Location II0.601-0.3380.193331.1220.67960
Table 5. Summary of fatigue experiments on the uniaxial, notched, and complex specimens of PP-CF.
Table 5. Summary of fatigue experiments on the uniaxial, notched, and complex specimens of PP-CF.
GeometryTriaxiality
[ η ]
Specimen Angle
[°]
ε max
[−]
ε min
[−]
( σ v M ) max
[MPa]
( σ v M ) min
[MPa]
Nf
[Cycles]
Uniaxial0.33300.0140.00437.329.91780
0.0160.00437.629.91040
0.0170.00637.832.9520
450.0130.00430.723.89780
0.0160.00431.123.81240
0.0170.00631.127.12480
900.0130.00426.817.3123,800
0.0160.00427.617.312,400
0.0170.00728.021.639,800
Notched Type I0.37900.1520.102161.7131.1262,400
0.1520.076161.7114.68700
0.389450.1630.109147.4122.3312,600
0.1630.081147.4104.224,500
0.404900.1260.049115.356.07800
0.1260.078115.382.2176,400
Notched Type II0.37600.1210.073145.3112.6196,000
0.1210.061145.3100.929,400
0.421450.1460.087141.8101.6156,700
0.1460.073141.888.521,600
0.435900.1450.100123.697.1846,200
0.1450.078123.680.036,500
Complex Type I0.46700.1370.082154.3126.389,500
0.1370.064154.3110.416,750
0.481450.1650.110140.0120.0746,500
0.1650.088140.0107.078,400
0.512900.2130.142134.9114.2216,500
0.2130.113134.9102.023,900
Complex Type II0.46700.1340.089153.0132.4314,500
0.1340.071153.0119.434,000
0.479450.1770.118140.8118.2226,500
0.1770.095140.8105.815,400
0.501900.1560.094176.2131.548,600
0.1560.073176.2110.24700
Table 6. The deviation between the experimental and simulation load-displacement curves.
Table 6. The deviation between the experimental and simulation load-displacement curves.
MaterialsGeometry Specimen Angle [°]Deviation [%]
PA6-CFUniaxial01.04
452.05
902.68
Notched Type I02.08
452.79
901.38
Notched Type II02.81
452.37
902.51
Complex Type I02.64
452.87
902.92
Complex Type II02.58
451.73
902.51
Cross-member location 1-2.94
Cross-member location 2-3.16
PP-CFUniaxial01.28
452.78
902.64
Notched Type I01.13
452.41
902.95
Notched Type II02.43
452.55
902.89
Complex Type I02.42
451.68
901.27
Complex Type II02.61
451.42
902.88
Table 7. Model constants for the proposed semi-empirical fatigue life prediction models.
Table 7. Model constants for the proposed semi-empirical fatigue life prediction models.
Materials A [−] c [−]
PA6-CF0.4958−0.183
PP-CF0.5915−0.167
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Choi, J.; Lee, H.; Lee, H.; Kim, N. A Methodology to Predict the Fatigue Life under Multi-Axial Loading of Carbon Fiber-Reinforced Polymer Composites Considering Anisotropic Mechanical Behavior. Materials 2023, 16, 1952. https://doi.org/10.3390/ma16051952

AMA Style

Choi J, Lee H, Lee H, Kim N. A Methodology to Predict the Fatigue Life under Multi-Axial Loading of Carbon Fiber-Reinforced Polymer Composites Considering Anisotropic Mechanical Behavior. Materials. 2023; 16(5):1952. https://doi.org/10.3390/ma16051952

Chicago/Turabian Style

Choi, Joeun, Hyungtak Lee, Hyungyil Lee, and Naksoo Kim. 2023. "A Methodology to Predict the Fatigue Life under Multi-Axial Loading of Carbon Fiber-Reinforced Polymer Composites Considering Anisotropic Mechanical Behavior" Materials 16, no. 5: 1952. https://doi.org/10.3390/ma16051952

APA Style

Choi, J., Lee, H., Lee, H., & Kim, N. (2023). A Methodology to Predict the Fatigue Life under Multi-Axial Loading of Carbon Fiber-Reinforced Polymer Composites Considering Anisotropic Mechanical Behavior. Materials, 16(5), 1952. https://doi.org/10.3390/ma16051952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop