Next Article in Journal
Temperature/Reduction Dual Response Nanogel Is Formed by In Situ Stereocomplexation of Poly (Lactic Acid)
Next Article in Special Issue
Novel Amphiphilic Polyfluorene-Graft-(Polymethacrylic Acid) Brushes: Synthesis, Conformation, and Self-Assembly
Previous Article in Journal
Surface Treatment of Composites with Bismaleimide Resin-Based Wet Peel Ply for Enhanced Adhesive Bonding Performance
Previous Article in Special Issue
The Mechanical and Thermal Properties of Poly(ethylene-co-vinyl acetate) (PECoVA) Composites with Pristine Dolomite and Organophilic Microcrystalline Dolomite (OMCD)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of the Ultimate Strength of Notched and Unnotched IM7/977-3 Laminated Composites Using a Micromechanics Approach

by
Manzar Masud
1,2,†,
Moosa S. M. Al Kharusi
3,*,
Muhammad Umair Ali
4,†,
Aamir Mubashar
1,*,
Shaik Javeed Hussain
5,*,
Adnan Tariq
6,
Gulfam Ul Rehman
6,
Mahmood Hassan Akhtar
7 and
Shama Javeed
8
1
Computational Mechanics Group, Department of Mechanical Engineering, School of Mechanical and Manufacturing Engineering (SMME), National University of Sciences and Technology (NUST), H-12, Islamabad 44000, Pakistan
2
Department of Mechanical Engineering, Capital University of Science and Technology (CUST), Islamabad 44000, Pakistan
3
Department of Mechanical and Industrial Engineering, College of Engineering, Sultan Qaboos University, 123 Al Khoud, Muscat 112, Oman
4
Department of Unmanned Vehicle Engineering, Sejong University, Seoul 05006, Korea
5
Department of Electrical and Electronics, Global College of Engineering and Technology, Muscat 112, Oman
6
Department of Mechanical Engineering, University of Wah, Wah Cantt 47040, Pakistan
7
Department of Chemistry, National University of Technology (NUTech), Islamabad 44000, Pakistan
8
Department of Mathematics, Pusan National University, Busan 46241, Korea
*
Authors to whom correspondence should be addressed.
M.M. and M.U.A. have equally contributed as the first author.
Polymers 2021, 13(20), 3491; https://doi.org/10.3390/polym13203491
Submission received: 24 August 2021 / Revised: 2 October 2021 / Accepted: 6 October 2021 / Published: 11 October 2021
(This article belongs to the Special Issue Frontiers in Polymer Composites)

Abstract

:
This paper proposes a multi-scale analysis technique based on the micromechanics of failure (MMF) to predict and investigate the damage progression and ultimate strength at failure of laminated composites. A lamina’s representative volume element (RVE) is developed to predict and calculate constituent stresses. Damages that occurred in the constituents are calculated using separate failure criteria for both fiber and matrix. Subsequently, the volume-based damage homogenization technique is utilized to prevent the localization of damage throughout the total matrix zone. The proposed multiscale analysis procedure is then used to investigate the notched and unnotched behavior of three multi-directional composite layups, [30, 60, 90, −60, 30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S, subjected to static tension and compression loading. The specimen is fabricated from unidirectionally reinforced composite (IM7/977-3). The prediction of ultimate strength at failure and equivalent stiffness are then benchmarked against the experimental test data. The comparative analysis with various failure models is also carried out to validate the proposed model. MMF demonstrated the capability to correctly predict the ultimate strength at failure for a range of multidirectional composites laminates under tensile and compressive load. The numerically predicted findings revealed a good agreement with the experimental test data. Out of the three investigated composite layups, the simulated results for the quasi-isotropic [0, 45, 90, −45]2S layup agreed extremely well with the experimental results with all the percentage errors within 10% of the measured failure loads.

1. Introduction

Composite laminate-based structures are extensively used in aerospace applications [1]. They are characterized as a significant configuration of composites that are critical for different categories of exceedingly loaded structures. Due to the elevated strength–weight ratio, high specific modulus, and the ability to be customized for a particular application, composite materials offer a variety of advantages compared to other conventional materials [2]. In view of the recent developments in the advanced aerospace industry and progressively high demand for increased performance of laminated composites structures, more sophisticated design and failure prediction models are required [3].
Reliable failure theories and progressive damage models are needed to precisely predict the complex failure phenomenon in the structures made of composites. Therefore, predictive tools that require a reduced number of essentials tests are becoming more important because of the extremely expensive tests on composite structures [4]. Computational multi-scale modelling and simulation tools for the prediction of damage mechanisms, progressive damage, and residual strengths can be used to achieve this ambitious goal [5,6]. However, the concerns like an-isotropy, non-linear stress–strain response, complex failure processes contribute to the difficulty of modeling the structural and damage behavior of composites, especially the failure’s initiation and propagation, subsequently leading to the prediction of ultimate strength [7,8].
Progressive Damage Analysis (PDA) and ultimate strength at failure’s prediction in composite laminates is a difficult task, even after many years of publication of ground-breaking works on failure theories by Rosen [9] and Tsai and Wu [10]. Even with composite structures subjected to constant in-plane loading, ultimate strength prediction has shown serious issues, as mentioned in different World-Wide Failure exercises [11,12]. The main problems and intricacy lie in the interaction and mixture of various damage modes in laminates, such as fiber-failure, matrix-cracking, and delamination, which ultimately result in residual strength and structural integrity loss.
Progressive failure simulation and modeling approaches focus on the growth of damage and the interaction of different phenomena in composite laminates. To evaluate the progressive damage approaches, a benchmarking exercise [13] organized by “Air Force Research Laboratory” (AFRL) [13], titled “Damage Tolerance Design Principles (DTDP),” was arranged to investigate the capability of different PDA models against AFRL provided test data. The main aim of the exercise was to predict the ultimate strength of several IM7/977-3 carbon fiber reinforced laminates under tensile and compressive loading. During this study, nine different analysis teams forecasted the ultimate strength of the unnotched and open-hole specimen. These predictions were made for three multidirectional laminate groups [14], which are [30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S and [0, 60, 60]3S.
The analysis codes presented in the DTDP exercise can be divided into two categories depending upon the type of methodology, i.e., macromechanics and micromechanics [13]. From the presented analysis codes, the b-spline analysis method with mesh independent cracking (BSAM with MIC) [13] and discrete crack network (DCN) [15] was based on the macromechanics level. In macromechanics, damage in the matrix and fibers is computed using composite level macro stresses. In contrast, the micromechanics approach utilizes stress and strain information from each of the composite constituents.
Using the micromechanics approach to perform the PDA of composite laminates is of great significance because of non-linearity before final failure, which most first ply failure criteria cannot capture. Furthermore, the micromechanics approach facilitates accounting precisely for changes in the properties of constituents in composite materials and the impact of the microstructure, such as the fiber-orientation and fiber volume ratio, which makes it a powerful tool to estimate progressive damage in composite based structures.
In the DTDP exercise [13], the analysis codes, namely generalized method of cells (MAC/GMC) [16], helius progressive failure analysis (Helius PFA) [17], enhanced schapery theory (EST) [13], multi-scale design system for linking continuum scales (MDS-C) [18], general optimization analyzer (GENOA) [19], eigendeformation-based reduced order homogenization (EHM) [20], and n-phase cylindrical model (NCYL) [21], were based on the micromechanics approach.
The failure theories based on micromechanics have become increasingly attractive for researchers in modeling the failure patterns of composite-based structures due to the benefit in describing the failure processes at the constituent level [22]. By combining the representative unit cells of the meso level and micro level, the multiscale analysis technique takes both the accuracy of micromechanics-based failure technique and efficiency of macro-based analysis. Gosse et al. [23] suggested a theory based on the stress invariant to predict the commencement of damage in composites by utilizing the strain-based enhancement factors.
A micromechanics-based modeling technique using stress invariants was proposed by Tran et al. [24]. Mayes et al. [25] suggested micromechanics depended on multi-continuum theory (MCT), in which the stress and strain fields are obtained using a FE analysis. Ha et al. [26] proposed a damage model based upon the micromechanics of failure (MMF) technique for the estimation of failure in composite material’s constituents using stress amplification factors. MMF is a constituent-based progressive damage model in which the data is transferred from the microscale model to the mesoscale model and vice versa. Lei et al. utilized MMF to accurately forecast the ultimate failure strength of woven configuration [27] and braided configuration [28] textile composites.
Considering the advantage of MMF in prediction accuracy, the MMF model is utilized to estimate the failure strength and stiffness of coupons fabricated from unidirectionally reinforced composite (IM7/977-3) in this work. The prediction is conducted for unnotched and open-hole specimens under tensile as well as compressive loading conditions. The coupons were categorized as unnotched tension (UNT), unnotched compression (UNC), open-hole tension (OHT), and open-hole compression (OHC). The investigated layups used in the DTDP program are [30, 60, 90, −60, −30]2S, [0, 45, 45, 90]2S, and [60, 0, −60]3S.
A micromechanics-based FE model employing representative volume element (RVE) is developed to estimate the ultimate strength of composite laminates subjected to tensile and compressive loading. The mesoscale-based structural analysis was used to give input to the RVE regarding the mesoscale stresses. Damages that occurred in different composite constituents are observed based on the failure model for each constituent. The volume-based damage homogenization technique is used to avoid damage localization throughout the total matrix zone. The predicted results using MMF are also assessed by comparison with the experimental tests data as well as the results from the other analysis methods.

2. Micromechanics Based Progressive Damage Model

2.1. Multi-Scale Analysis

The idealized model of the composite laminate is essential to predict their mechanics-based behavior subjected to static loading, prior to multi-scale analysis. The schematic of damage theory can be split up into macroscale procedure and microscale procedure. In the homogenized framework of multiscale analysis, the microscale sub-domain is a representative volume element (RVE) [29], from which micromechanical results are used to yield results from the macro-scale model. The RVE of a lamina is developed in the terms of the fiber, matrix, and interface to express microstructure as illustrated in Figure 1.
The Representative Volume Element (RVE), a micro unit cell of hexagonal fiber shape as shown in Figure 1, may be utilized to correlate the behavior of ply to the behavior of fiber and matrix constituents [30]. Stress amplification factors (SAFs) [31], denoted by M and A matrices, are used for the transformation between the ply level macro stresses ( σ ¯ ), and the fiber-based stresses (σf) and matrix-based stresses (σm) at the micro level. The transformation can be computed with the help of the following Equations (1) and (2).
σ f = M f σ ¯ + A f Δ T
σ m = M m σ ¯ + A m Δ T  
where ΔT is the difference between ambient temperature and zero temperature. Mf and Mm are the SAFs (mechanical) for the fiber and matrix, respectively, whereas Af and Am are the SAFs (thermal) for the fiber and matrix, respectively. Mf and Mm are 6 × 6 matrices, and Af and Am are 6 × 1 matrices. Equation (1) shows that the correlation between stresses at macro level and micro level is linear with a uniform fiber distribution.
Multiscale analysis is performed in three steps. First, the SAFs are computed by applying FE analysis on unit cell as shown in Figure 1. Second, Equation (1) is used to determine the micro level stresses of selected key points in the fiber and matrix using stresses in the ply. In the third step, micro stresses are then substituted in the constituent depended progressive damage model to determine the damage factors in the fiber and matrix. Later, as a result, the stiffness of each component, i.e., the fiber and matrix, was degraded, based on the calculated damages. Subsequently, the effective properties of the ply were again computed with the help of RVE.

2.2. Stress Amplification Factors

The diamond array of a fiber was selected to calculate the SAFs as shown in Figure 1. The stress state of unit cell is described by selecting 18 key points in the fiber and 19 key points in the matrix [28]. The majority of the selected points are positioned around the interface area where the maximum stress concentration takes place, and the remaining key points are situated in the center of the fiber and matrix. SAFs for the macro stress M and for the temperature increment A were calculated through direct FE analysis of a unit cell with appropriate repeated boundary conditions to ensure the correct behavior of a cell model [32]. The comprehensive structure of Equation (1) is as follows:
( σ 1 σ 2 σ 3 σ 4 σ 5 σ 6 ) = [ M 11 M 12 M 13 M 14 0 0 M 21 M 22 M 23 M 24 0 0 M 31 M 32 M 33 M 34 0 0 M 41 M 42 M 43 M 44 0 0 0 0 0 0 M 55 M 56 0 0 0 0 M 65 M 66 ] ( σ ¯ 1 σ ¯ 2 σ ¯ 3 σ ¯ 4 σ ¯ 5 σ ¯ 6 ) + ( A 1 A 2 A 3 A 4 A 5 A 6 ) Δ T
Stress amplification factors were determined by applying the uniformly distributed unidirectional unit loads one at a time at each direction to the unit cell shown in Figure 1.
For example, applying a unit load having uniform distribution at x direction of unit cell with no thermal load will result in the simplification of Equation (3). The simplified Equation (4) is given below:
( σ 1 σ 2 σ 3 σ 4 σ 5 σ 6 ) = [ M 11 M 12 M 13 M 14 0 0 M 21 M 22 M 23 M 24 0 0 M 31 M 32 M 33 M 34 0 0 M 41 M 42 M 43 M 44 0 0 0 0 0 0 M 55 M 56 0 0 0 0 M 65 M 66 ] ( 1 0 0 0 0 0 )
The solution of linear system given in Equation (4) yields Equation (5):
( σ 1 σ 2 σ 3 σ 4 σ 5 σ 6 ) = [ M 11 M 21 M 31 M 41 0 0 ]
Equation (5) shows that the calculated SAFs are the micromechanical stresses, e.g., M11 is basically the σ 1 and so on. The same method was applied to other directions to calculate all the other SAFs.

2.3. Damage Model of Constituents

A separate damage model is used for each individual constituent of the composite material, i.e., fiber and matrix, due to their different mechanical properties and behaviors. Usually, the glass fiber behaves as an isotropic material, and carbon fiber behaves as a transversely isotropic. Therefore, the behavior of both carbon and glass fibers are presumed to be linear elastic and brittle. On the other hand, the polymetric matrix behaves as an isotropic and ductile material, demonstrating non-linear behavior when it is damaged.

2.4. Fiber Failure Criteria

For fiber material, the maximum longitudinal stress damage model was used as (Equation (6)):
C f < σ f 11 < T f
The symbol σ f 11 indicates the longitudinal stress at the micro level for the fiber, and notations Tf and Cf presents the tensile strength (longitudinal) and compressive strength (longitudinal) of the fiber constituent, respectively. As soon as the produced stress in the longitudinal direction of fiber element reaches the Tf or Cf, then the fiber will fail. As a result of failure, the stiffness of fiber constituent will be significantly degraded as showed in the progressive model Figure 2.
Therefore, once the fiber failure occurs, the localized damage factor, Df, is set to 0.9 to ensure computational stability [33], and if fiber failure does not occur, then the Df is assigned value of zero. The fiber materials in this study are taken as brittle; therefore, complete damage of single element is considered catastrophic towards the total area of the fiber. As a result, the global factor for damage in the fiber was described as the maximum Df of all the calculated Df of elements (Equation (7)).
D ¯ f = M a x i m u m   ( D f ( j ) )
where D f ( j ) is the damage factor produced in jth element of fiber. Therefore, a fundamental correlation for the stresses induced in the fiber was then calculated as (Equation (8)):
σ f = ( 1 D ¯ f )   C f ε f
In Equation (8), Cf indicates the stiffness matrix of the total fiber area inside the RVE.

2.5. Matrix Failure Criteria

In the majority of cases, the matrix materials are isotropic, with different tensile and compressive strengths. Both types of strength contribute in a different way to matrix failure. Raghava et al. [34] and Ha et al. [26] presented a tailored adaptation of the von-Mises criterion for failure, using tensile strength Tm and compressive strength Cm, which is as follows (Equation (9)):
σ V M 2 C m T m + ( 1 T m 1 C m ) I 1 = 1
I1 and σVM may be computed by utilizing stress-based components with the help of the subsequent Equations (10)–(12):
I 1 = σ m 1 + σ m 2 + σ m 3
σ V M = 1 2 [ [ ( σ m 1 σ m 2 ) 2 ] + [ ( σ m 2 σ m 3 ) 2 ] + [ ( σ m 3 σ m 1 ) 2 ] + σ m 4 2 + σ m 5 2 + σ m 6 2 ]
The failure criterion mentioned in Equation (9) is the same as the condition that an equivalent stress σ e q achieves the initial Tm,
σ e q = ( β 1 )   I 1 + ( β 1 )   2   I 1 2 + 4 β σ V M 2 2 β
In Equation (12), β for the matrix is defined as the ratio between compressive strength and tensile strength. The equivalent stress σ e q , frequently termed stasis stress, is established upon the stress–strain correlation. The σ e q can be translated into the equivalent strain ε e q with the help of subsequent Equation (13).
ε e q = ( β 1 ) J 1 + ( β 1 ) 2 J 1 2 + ( 2 4 ν 1 + ν ) 2 β ε V M 2 2 β ( 1 2 ν )
In Equation (13), J1 represents first-strain invariant, ε V M represents the von Mises equivalent strain and ν is the Poisson’s ratio. J1 and ε V M can be computed with the help of strain-based components using the subsequent Equations (14) and (15) [35]:
J 1 = ε m 1 + ε m 2 + ε m 3
ε V M = 1 2 [ [ ( ε m 1 ε m 2 ) 2 ] + [ ( ε m 2 ε m 3 ) 2 ] + [ ( ε m 3 ε m 1 ) 2 ] + ε m 4 2 + ε m 5 2 + ε m 6 2 ]
The damage evolution in the matrix material was computed with the help of the equivalent strain calculated from Equation (13). Using σ e q and ε e q , a multilinear failure model for a matrix, based on stress strain, was suggested by Ha et al. [28], as shown in Figure 3. The response of the matrix, before damage, follows the linear stress–strain relation. As soon as the damage occurs, the matrix reveals hardening behavior, subsequently followed by softening behavior. The stiffness degradation occurs on the basis of the damage factor Dm in both hardening and softening behavior. The overall damage factor in the matrix Dm can be computed using the following Equation (16).
D m = 1 ( ε y ( i ) ε e q ) σ y ( i 1 ) + ( ε e q ε y ( i 1 ) ) σ y ( i ) E 0 ( ε e q ε y ( i 1 ) ) ε e q ( ε y ( i 1 ) < ε e q < ε y ( i ) )
where σ y ( i 1 ) and ε y ( i 1 ) indicate the yielding stress and strain at the initial point of the ith damage state; σ y ( i ) and ε y ( i ) indicate the yielding stress and strain at the ending point of the ith damage state; ε e q represents existing equivalent strain in the i-th damage state, and E0 indicates undamaged matrix’s stiffness. To calculate the yielding stress and strain at the start and end point of the damage state, the tensile stress–strain curve of the polymer matrix (i.e., epoxy 977-3) reported in the literature [36] is used in this work.
The stress–strain curve is divided into 10 stages. In each stage, the yield stress and corresponding yield strain are used to evaluate the matrix damage Dm. The details regarding the yield stress–strain data of each stage are described in Section 3.3. The maximum Dm was selected as 0.9 to ensure computational stability [33]. After the computation of each element’s failure state, the damage localization was prevented with the help of the volume-based damage homogenization technique [37], given by the following Equation (17).
D ¯ m = V [ D m ] p   d   V V   d   V
In Equation (17), p is a local damage’s positive weighting parameter and V is defined as the total volume of the matrix area inside the RVE. Considering the stiffness degradation of the matrix, the constitutive relation was calculated by (Equation (18)):
σ m = ( 1 D ¯ m )   C m ε m
In Equation (18), Cm indicates the stiffness matrix of the total matrix area in the RVE.

2.6. Numerical Implementation of Constituent Damage Models

Figure 4 shows the flow diagram of the algorithm, established for the methodology combining MMF and progressive damage model. The detailed methodology is as follows:
  • The overall global strain produced at the global time n, was calculated through the addition of the global strain increment at n and at the previous step n-1.
  • With regard to the composite laminate, the macro level stresses of each constituent element, i.e., the fiber and matrix, were calculated with the help of previous effective stiffness properties.
  • Subsequently, the micro level stress of each constituent was computed from the macro stresses with the help of SAFs.
  • The damage model of constituents was then employed to both fiber and matrix to calculate the damage factor in the fiber and matrix areas, symbolized by Dm and Df.
  • The total damage factor produced for matrix and fiber zones were then analyzed using their corresponding damage methods, i.e., the maximum damage for fiber and the damage homogenization using volume-based technique for the matrix.
  • The stiffness degradation of matrix and fiber were calculated using the status of the total damage factor. Subsequently, the ply level effective properties were calculated for the following time increment.
  • The numerical execution of MMF was implemented using ABAQUS, combined with the user subroutine USDFLD.

3. Methodology

An implicit scheme was employed to execute the simulation of the static loading with the help of the finite element analysis software package ABAQUS/Implicit.

3.1. Models Geometries and Boundary Conditions

The DTDP exercise program by AFRL included UNT, UNC, OHT, and OHC, and there were three layups for each coupon, resulting in twelve multi-directional configurations. The investigated layups were soft layup [30, 60, 90, −60, −30]2S, quasi-isotropic layup [0, 45, 90, −45]2S, and [60, 0, −60]3S with 0° developed in parallel with the loading direction. The model geometries along with dimensions and boundary conditions for UNT and UNC is given Figure 5 and for OHT and OHC is given in Figure 6.
To avoid unnecessary complexity to the finite element models and artificial stress concentrations, the geometrical models were simplified by exclusion of the fixtures. In the case of UNT ad UNC, only the gage section was modeled, and boundary conditions were applied to the gage section as shown in Figure 5. For all developed models, the ply thickness was taken as 0.12954 mm.

3.2. Mesh Details

The successful simulation of physical damage mechanisms occurring in composite structures as a result of finite element-based failure theory is highly reliant on the consistency of the mesh discretization [38]. Modeling and meshing laminate using single element in the direction of stack, demonstrating entire composite laminate, will no longer be able to simulate the state of interlaminar damage. Ideally, the finite element model with structured mesh type will have complete reliability to catch and identify all the failure modes; however, using this type of mesh will result in the overcapacity of computing resources as well as consumption of excessive time to complete the simulation with full load history. Therefore, a workable approach to create a finite element model is one that produces precise simulations without being excessively cumbersome.
In this study, solid element type is used with through-thickness discretization techniques, which utilize non-layered elements. In structured mesh, non-layered types of elements have a single ply per element and are given preference when models have stress concentrations causing large transverse shear and normal stress gradients. Non-layer types of elements are also capable of simulating the localized softening behavior that arises because of failure occurs in the fiber and matrix [39]. As half of the specimens investigated in the DTDP exercise by AFRL contains stress concentrations due to notches, the meshes of all the models used non-layered elements, i.e., each ply was meshed with the help of separate layer.
In finite element models, the plies were modelled using the Abaqus element type C3D8R, which is an eight-node linear brick, reduced integration, hourglass control element [40]. The reduced integration elements are efficient, specifically in the case of progressive damage simulations as the entire element is represented with the help of a single integration point [41]. The adapted mesh developed for the unnotched models is illustrated in Figure 6. The mesh discretization for the notched models as illustrated in the Figure 7, is based on “notch centric” design, characterized with the help of rounded mesh all around the region of notch that rapidly converts into a rectangular-based mesh shifting farther from the notch.
A zoomed view of the mesh is also shown in Figure 7 and Figure 8. Mesh sensitivity analysis was performed for all the cases, and we established that the final adapted mesh was to be a compromise between the output quality and time required for computation for each case. In particular, for the mesh sensitive OHT and OHC coupons, five cases were investigated with the numbers of elements of 40, 960; 67, 840; 97, 280; 179, 200; and 430, 080 elements, respectively. The mesh convergence study shows that there is no major difference in the strength values as the number of elements are increased; therefore, the case with number of elements of 97, 280, was chosen for the study, keeping in mind the time required and the output quality.

3.3. Material Properties

In this study, the matrix (977-3) is presumed as an isotropic material with two independent material properties, and the fiber (IM7) behaves as transversely isotropic with five independent material properties. Initially, the mechanical properties of IM7/977-3 ply were taken from experimental results provided by the AFRL [14]. The lamina experimental tests were all performed employing ASTM test standards 85. During the benchmark study, the initial measurement of the compressive failure strength of laminate that was delivered to the analysts did not keep into account the specimen buckling, resulting in an incorrectly measured failure load. However, after the discovery of buckling failure mode, analysts were instructed to increase the laminate compressive failure strength input property to 1680 MPa for recalibration [14].
The calibration of material properties for the carbon fiber (IM7) and matrix (977-3) was performed using the micromechanics approach, which minimizes the difference between the numerical model and experimental data. During calibration, the fiber and matrix properties were determined using the micromechanics model to estimate the composite laminate elastic properties that closely matched the experimental results. The stiffness and strength properties of the carbon-based fiber constituent [14] and matrix listed in Table 1. The tensile yield stress and yield strain data of matrix (977-3) [36] at different stages for matrix multi-linear damage model is given in Table 2. The effective properties of the ply were computed from the micro unit cell, using the micromechanics approach and are listed in Table 3.

4. Results and Discussion

The ability to precisely predict the damage onset and damage propagation have a potential role in designing and manufacturing composite structures. Furthermore, accurate progressive damage models could substantially minimize the constituents damage tolerance tests, reducing the time and cost to verify a composite structure. This section investigates the performance of the proposed micromechanics-based failure model (MMF) with experimental results [14] and other failure models [21]. First, the finite element simulations were performed for all four test specimens, i.e., UNT, UNC, OHT, and OHC, with three distinct multi-directional layups ([30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S) per coupon. The findings were then compared with the experimental results and other failure techniques to investigate the efficiency of the MMF.

4.1. Unnotched Tension

The comparison between the experimental and simulated results for static strength at failure and stiffness for different unnotched tension specimen/layup configurations is given in Table 4. The comparison of results of MMF and other analysis methods [21] with experimental results [14] for unnotched static strength and stiffness in terms of percentage error is given in the form of bar charts in Figure 9.
The results given in Table 4 and Figure 9 show that the MMF presented a good agreement for quasi-isotropic [0, 45, 90, −45]2S configuration and [60/0/−60]3S configuration with a percentage error of 3.12% and 10.45%, respectively. Whereas, for soft layup [30/60/90/−60/−30]2S, the MMF underpredicted the results with a percentage error of 18.82 %. The higher percentage error in the simulation response of this soft layup is attributed to the fact that the real failure mechanism for the soft layup as per experimental investigation is too complicated to detect using modeling techniques.
This issue was also discussed by Dalgarno et al. [17]. The investigation of experimental results and failure images [14] indicate that the laminate exhibited considerable delamination. Currently, MMF does not account for the delamination effect due to computational simplicity and execution speed. To clearly understand the overall progressive damage simulation response, Figure 10 shows a comparison of the stress–strain curves obtained by MMF and the experimental results for the UNT case for all three layups.
In addition to the unnotched tension strength, the comparison of the unnotched tension stiffness results with experimental results given in Table 4 and Figure 9 shows that the MMF exhibited excellent agreement for [30/60/90/−60/−30]2S, [0, 45, 90, −45]2S, and [60/0/−60]3S layups with percentage errors of 0.53%, 0.08%, and 0.83%, respectively.
The comparison of stiffness results in Figure 9 shows that the MMF performed better than the other failure methods. However, for the strength prediction, soft layup [30, 60, 90, −60, −30]2S, MMF, and other analysis methods, like HELIUS PFA, EST, BSAM/MIC, and EHM, have higher percentage errors. For the other two layups, i.e., [0, 45, 90, −45]2S and [60, 0, −60]3S, MMF performed reasonably well to predict the unnotched static strength.

4.2. Unnotched Compression

The results of the unnotched compression static strength and stiffness are given in Table 5. A comparison was performed between the experimental and simulated results for static strength and stiffness for three different layups, which were [30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S and [60, 0, −60]3S. For the UNC case, the comparison of MMF and other analysis methods with experimental results for the three investigated layups is given in Figure 11.
Examination of the values of static strength at failure for different configurations given in Table 5, and Figure 11 reveals a good agreement between MMF and the experimental results for [30, 60, 90, −60, −30]2S and [0, 45, 90, −45]2S with percentage errors of 11.78% and 9.29%, respectively. However, for [60, 0, −60]3S, the results of MMF are underpredicted with a percentage error of 28.37%. In addition, the comparison of simulated response by MMF and experimental response in terms of the stress–strain curves for the UNC case for all three layups is shown in Figure 12.
On the other hand, for unnotched compression stiffness, the results of MMF are well aligned with the experimental results with the percentage errors of 5.15%, 8.71%, and 6.61% for the [30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S layups, respectively.
For the stiffness comparison shown in Figure 11, the performance of MMF and other analysis methods is excellent except for the MDS-C, in which the percentage error is very high for all investigated layups. For the strength prediction in the UNC case, the comparison showed a good performance of MMF for the [30, 60, 90, −60, 30]2S and [0, 45, 90, −45]2S layups. However, for the [60, 0, −60]3S laminate sequence, most of the analysis methods, including MMF, have higher percentage errors as shown in Figure 11.

4.3. Notched Tension

The experimental and simulated results for the notched pattern coupons are given in Table 6. The results show the comparison between the experimental and simulated results for three investigated multidirectional composite layups. The comparison of MMF and other failure criterion with MMF for open hole tension static strength at failure and stiffness is illustrated in the form of bar charts in Figure 13.
From the results presented in Table 6 and Figure 13, it is concluded that the MMF model correlated better to experimental findings for the open hole tension static strength for the soft layup [30, 60, 90, 60, −30]2S and quasi-isotropic lay [0, 45, 90,−45]2S with a percentage error of 1.71% and 4.51%, respectively. However, for the [60, 0, −60]3S, the results of static strength at failure are underpredicted with a percentage error of 14.18%. To clearly recognize the overall progressive damage simulation response and accumulated simulated damage at failure, Figure 14 shows a comparison of stress–strain curves obtained by MMF and experimental results for the OHT case for all three layups along with the images of simulated fiber and matrix damage distribution at failure.
On the other hand, the results of MMF for notched tension stiffness are well aligned with the experimental results for [30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S layups with percentage errors of 1.25%, 0.62%, and 0.41%, respectively.
Like the previous two cases of UNT and UNC, MMF performed well compared to the other analysis methods for OHT stiffness. For the strength prediction case, the MMF showed good agreement for the [30, 60, 90, −60, −30]2S and [0, 45, 90, −45]2S layups. However, for the [60, 0, −60]3S layup, like the UNC case, nine out of ten analysis methods, including MMF, showed high percentage errors, as shown in Figure 13.

4.4. Notched Compression

Table 7 shows a comparison of overall open hole compression static strength and stiffness responses of simulation with the experimental results for the specimen of soft layup [30, 60, 90, −60, −30]2S, quasi-isotropic layup [0, 45, 90, −45]2S, and [60, 0, −60]3S layup. For the OHC case, Figure 15 compares percentage errors for the MMF and other analysis methods with experimental results for the three investigated composite layups.
The results given in Table 7 and Figure 15 show that the MMF exhibited good agreement with the obtained experimental results for all the experimentally investigated composite layups. For open hole compression static strength, the percentage errors for [30, 60, 90, 60, −30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S layups are 4.41%, 7.33%, and 8.38%, respectively. On the other hand, for the open hole compression stiffness, the percentages for the layups are 7.67%, 5.23%, and 6.36%, respectively. In addition, the comparison of the stress–strain curves obtained by MMF and the experimental results for the OHC case for all three layups along with the images of simulated fiber and matrix damage distribution at failure is shown in Figure 16.
The comparison given in Figure 15 showed that the MMF performed reasonably well to predict the notched compression static strength and stiffness compared to most of the other analysis methods. For OHC stiffness prediction, the GENOA, DCN, and MDS-C analysis methods had high average percentage errors. For the OHC strength prediction, there was no clear pattern for the prediction of the result, as shown in Figure 15. However, MMF performed well for the prediction of strength for all three investigated layups.

5. Conclusions

In this study, a micromechanics-based progressive damage model MMF was used to predict the ultimate strength at failure and stiffness of four different multi-directional laminated composite test specimens, i.e., UNT, UNC, OHT, and OHC. For each case, three different composite laminate layups were used: [30, 60, 90, −60, −30]2S, [0, 45, 90, −45]2S, and [60, 0, −60]3S. The results of MMF were then compared with the experimental results provided by AFRL and the results of nine different analysis methods that participated in the DTDP program.
The simulations were performed using ABAQUS, and the numerical implementation of MMF was done using the user subroutine USDFLD. The simulated results of stiffness for all 12 specimens/layups accurately aligned with the experimental results. For the ultimate strength at failure predictions, the simulated results agreed extremely well with the experimental results for the quasi-isotropic [0, 45, 90, −45]2S layup with the percentage errors within 10% of the measured failure loads.
For the [60, 0, −60]3S configuration, the findings of MMF were underpredicted with percentage errors of 28.37% and 14.18% for the UNC and OHT, respectively. For the soft layup [30, 60, 90, −60, −30]2S, the simulated results underpredicted the strength of the UNT specimen by 18.82%. For this specific laminate, the experimentally detected failure pattern observed substantial delamination during investigation. The proposed modeling technique does not take into consideration the delamination, which resulted in the higher percentage error. We concluded that MMF demonstrated the ability to correctly predict the ultimate strength at failure for a range of laminates under tensile and compressive loading.

Author Contributions

Conceptualization, M.M.; Data curation, M.U.A.; Formal analysis, A.M. and G.U.R.; Funding acquisition, M.S.M.A.K. and A.M.; Investigation, M.M. and M.U.A.; Methodology, M.M., S.J.H., A.M., and S.J.; Project administration, A.T.; Software, M.M.; Supervision, M.S.M.A.K., A.M., and A.M.; Validation, M.H.A.; Writing—original draft, M.M., M.U.A., and S.J.; Writing—review and editing, M.S.M.A.K., A.M., S.J.H., A.M., M.U.A., A.T., and M.H.A. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the Sultan Qaboos University, Muscat for funding this research under an internal research funding grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jefferson, A.J.; Arumugam, V.; Dhakal, H.N. 1—Introduction and context. In Repair of Polymer Composites; Jefferson, A.J., Arumugam, V., Dhakal, H.N., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2018; pp. 1–43. ISBN 978-0-08-102263-4. [Google Scholar]
  2. Rajak, D.K.; Pagar, D.D.; Menezes, P.L.; Linul, E. Fiber-Reinforced Polymer Composites: Manufacturing, Properties, and Applications. Polymers 2019, 11, 1667. [Google Scholar] [CrossRef] [Green Version]
  3. Rohwer, K. Predicting fiber composite damage and failure. J. Compos. Mater. 2014, 49, 2673–2683. [Google Scholar] [CrossRef] [Green Version]
  4. Talreja, R. Assessment of the fundamentals of failure theories for composite materials. Compos. Sci. Technol. 2014, 105, 190–201. [Google Scholar] [CrossRef]
  5. Zhao, Z.; Dang, H.; Xing, J.; Li, X.; Zhang, C.; Binienda, W.K.; Li, Y. Progressive Failure Simulation of Notched Tensile Specimen for Triaxially-Braided Composites. Materials 2019, 12, 833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Shah, S.Z.H.; Megat-Yusoff, P.S.M.; Karuppanan, S.; Choudhry, R.S.; Sajid, Z. Multiscale damage modelling of 3D woven composites under static and impact loads. Compos. Part A Appl. Sci. Manuf. 2021, 151, 106659. [Google Scholar] [CrossRef]
  7. Sun, C.T.; Tao, J. Prediction of failure envelopes and stress/strain behaviour of composite laminates1This article represents the authors’ contribution to a world-wide exercise to confirm the state-of-the-art for predicting failure in composites, organised by Hinton and Sod. Compos. Sci. Technol. 1998, 58, 1125–1136. [Google Scholar] [CrossRef]
  8. Shah, S.Z.H.; Karuppanan, S.; Megat-Yusoff, P.S.M.; Sajid, Z. Impact resistance and damage tolerance of fiber reinforced composites: A review. Compos. Struct. 2019, 217, 100–121. [Google Scholar] [CrossRef]
  9. Rosen, B. Mechanics of composite strengthening. In Fiber Composite Materials; American Society for Metals: Metals Park, OH, USA, 1965; pp. 37–75. [Google Scholar]
  10. Tsai, S.W.; Wu, E.M. A General Theory of Strength for Anisotropic Materials. J. Compos. Mater. 1971, 5, 58–80. [Google Scholar] [CrossRef]
  11. Hinton, M.J.; Kaddour, A.S. The background to the Second World-Wide Failure Exercise. J. Compos. Mater. 2012, 46, 2283–2294. [Google Scholar] [CrossRef]
  12. Kaddour, A.S.; Hinton, M.J.; Smith, P.A.; Li, S. The background to the third world-wide failure exercise. J. Compos. Mater. 2013, 47, 2417–2426. [Google Scholar] [CrossRef]
  13. Clay, S.B.; Engelstad, S.P. Benchmarking of composite progressive damage analysis methods: The background. J. Compos. Mater. 2016, 51, 1325–1331. [Google Scholar] [CrossRef]
  14. Clay, S.B.; Knoth, P.M. Experimental results of quasi-static testing for calibration and validation of composite progressive damage analysis methods. J. Compos. Mater. 2016, 51, 1333–1353. [Google Scholar] [CrossRef]
  15. Cui, X.; Fang, E.; Lua, J. A discrete crack network toolkit for Abaqus for damage and residual strength prediction of laminated composites. J. Compos. Mater. 2016, 51, 1355–1378. [Google Scholar] [CrossRef]
  16. Naghipour, P.; Arnold, S.M.; Pineda, E.J.; Stier, B.; Hansen, L.; Bednarcyk, B.A.; Waas, A.M. Multiscale static analysis of notched and unnotched laminates using the generalized method of cells. J. Compos. Mater. 2016, 51, 1433–1454. [Google Scholar] [CrossRef]
  17. Dalgarno, R.W.; Action, J.E.; Robbins, D.H.; Engelstad, S.P. Failure simulations of open-hole and unnotched IM7/977-3 coupons subjected to quasi-static loading using Autodesk Helius PFA. J. Compos. Mater. 2016, 51, 1421–1432. [Google Scholar] [CrossRef]
  18. Yuan, Z.; Crouch, R.; Wollschlager, J.; Shojaei, A.; Fish, J. Assessment of Altair Multiscale Designer for damage tolerant design principles (DTDP) of advanced composite aircraft structures. J. Compos. Mater. 2016, 51, 1379–1391. [Google Scholar] [CrossRef]
  19. Godines, C.; DorMohammadi, S.; Abdi, F.; Villa Montero, M.; Huang, D.; Minnetyan, L. Damage tolerant composite design principles for aircraft components under static service loading using multi-scale progressive failure analysis. J. Compos. Mater. 2016, 51, 1393–1419. [Google Scholar] [CrossRef]
  20. Bogdanor, M.J.; Oskay, C. Prediction of progressive fatigue damage and failure behavior of IM7/977-3 composites using the reduced-order multiple space-time homogenization approach. J. Compos. Mater. 2016, 51, 2101–2117. [Google Scholar] [CrossRef]
  21. Engelstad, S.P.; Clay, S.B. Comparison of composite damage growth tools for static behavior of notched composite laminates. J. Compos. Mater. 2016, 51, 1493–1524. [Google Scholar] [CrossRef]
  22. Moncada, A.M.; Chattopadhyay, A.; Bednarcyk, B.A.; Arnold, S.M. Micromechanics-based progressive failure analysis of composite laminates using different constituent failure theories. J. Reinf. Plast. Compos. 2012, 31, 1467–1487. [Google Scholar] [CrossRef]
  23. Gosse, J.; Christensen, S. Strain invariant failure criteria for polymers in composite materials. In Proceedings of the 19th AIAA Applied Aerodynamics Conference, Fluid Dynamics and Co-located Conferences, Anaheim, CA, USA, 11–14 June 2001. [Google Scholar]
  24. Tran, T.D.; Kelly, D.; Prusty, B.G.; Gosse, J.H.; Christensen, S. Micromechanical modelling for onset of distortional matrix damage of fiber reinforced composite materials. Compos. Struct. 2012, 94, 745–757. [Google Scholar] [CrossRef]
  25. Mayes, J.S.; Hansen, A.C. Composite laminate failure analysis using multicontinuum theory. Compos. Sci. Technol. 2004, 64, 379–394. [Google Scholar] [CrossRef]
  26. Ha, S.K.; Jin, K.K.; Huang, Y. Micro-Mechanics of Failure (MMF) for Continuous Fiber Reinforced Composites. J. Compos. Mater. 2008, 42, 1873–1895. [Google Scholar]
  27. Xu, L.; Huang, Y.; Zhao, C.; Ha, S. Progressive failure prediction of woven fabric composites using a multi-scale approach. Int. J. Damage Mech. 2016, 27, 97–119. [Google Scholar] [CrossRef]
  28. Xu, L.; Jin, C.; Ha, S. Ultimate strength prediction of braided textile composites using a multi-scale approach. J. Compos. Mater. 2014, 49, 477–494. [Google Scholar] [CrossRef]
  29. Li, S.; Sitnikova, E. Chapter 6—Formulation of unit cells. In Representative Volume Elements and Unit Cells; Li, S., Sitnikova, E., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2020; pp. 107–187. ISBN 978-0-08-102638-0. [Google Scholar]
  30. Li, S.; Sitnikova, E. Chapter 13—Application of unit cells to problems of finite deformation. In Representative Volume Elements and Unit Cells; Li, S., Sitnikova, E., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2020; pp. 417–438. ISBN 978-0-08-102638-0. [Google Scholar]
  31. Liu, Z.; Guan, Z.; Tan, R.; Xu, J.; Li, X. Multiscale Analysis of CFRP Laminates with MMF3 Criterion under Different Off-Axis Loading Conditions. Materials 2018, 11, 2255. [Google Scholar] [CrossRef] [Green Version]
  32. Li, S.; Sitnikova, E. Chapter 7—Periodic traction boundary conditions and the key degrees of freedom for unit cells. In Representative Volume Elements and Unit Cells; Li, S., Sitnikova, E., Eds.; Woodhead Publishing Series in Composites Science and Engineering; Woodhead Publishing: Cambridge, UK, 2020; pp. 189–221. ISBN 978-0-08-102638-0. [Google Scholar]
  33. Hashimoto, M.; Okabe, T.; Sasayama, T.; Matsutani, H.; Nishikawa, M. Prediction of tensile strength of discontinuous carbon fiber/polypropylene composite with fiber orientation distribution. Compos. Part A Appl. Sci. Manuf. 2012, 43, 1791–1799. [Google Scholar] [CrossRef]
  34. Raghava, R.; Caddell, R.M.; Yeh, G.S.Y. The macroscopic yield behaviour of polymers. J. Mater. Sci. 1973, 8, 225–232. [Google Scholar] [CrossRef] [Green Version]
  35. Jin, K.-K.; Huang, Y.; Lee, Y.-H.; Ha, S.K. Distribution of Micro Stresses and Interfacial Tractions in Unidirectional Composites. J. Compos. Mater. 2008, 42, 1825–1849. [Google Scholar] [CrossRef]
  36. DorMohammdi, S.; Godines, C.; Abdi, F.; Huang, D.; Repupilli, M.; Minnetyan, L. Damage-tolerant composite design principles for aircraft components under fatigue service loading using multi-scale progressive failure analysis. J. Compos. Mater. 2017, 51, 2181–2202. [Google Scholar] [CrossRef]
  37. Allen, D.H. Homogenization principles and their application to continuum damage mechanics. Compos. Sci. Technol. 2001, 61, 2223–2230. [Google Scholar] [CrossRef]
  38. David Müzel, S.; Bonhin, E.P.; Guimarães, N.M.; Guidi, E.S. Application of the Finite Element Method in the Analysis of Composite Materials: A Review. Polymers 2020, 12, 818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Choudhry, R.S.; Hassan, S.F.; Li, S.; Day, R. Damage in single lap joints of woven fabric reinforced polymeric composites subjected to transverse impact loading. Int. J. Impact Eng. 2015, 80, 76–93. [Google Scholar] [CrossRef] [Green Version]
  40. Systèmes, D. Abaqus Standard and Abaqus Documentation, Version 6.14; Dassault Systèmes Simulia Corp.: Johnston, RI, USA, 2014. [Google Scholar]
  41. Documentation, A. User Manual, Version 6.14.; Dassault Systèmes Simulia Corp.: Johnston, RI, USA, 2010. [Google Scholar]
Figure 1. Representative volume element (hexagonal fiber).
Figure 1. Representative volume element (hexagonal fiber).
Polymers 13 03491 g001
Figure 2. Material damage model for fiber.
Figure 2. Material damage model for fiber.
Polymers 13 03491 g002
Figure 3. Material damage model for the matrix [27].
Figure 3. Material damage model for the matrix [27].
Polymers 13 03491 g003
Figure 4. Algorithm combining the MMF and the progressive damage model.
Figure 4. Algorithm combining the MMF and the progressive damage model.
Polymers 13 03491 g004
Figure 5. Geometries and boundary conditions (a) UNT (b) UNC.
Figure 5. Geometries and boundary conditions (a) UNT (b) UNC.
Polymers 13 03491 g005
Figure 6. Geometries and boundary conditions (a) OHT (b) OHC.
Figure 6. Geometries and boundary conditions (a) OHT (b) OHC.
Polymers 13 03491 g006
Figure 7. Mesh adapted for UNT and UNC models.
Figure 7. Mesh adapted for UNT and UNC models.
Polymers 13 03491 g007
Figure 8. Mesh adapted for OHT and OHC models.
Figure 8. Mesh adapted for OHT and OHC models.
Polymers 13 03491 g008
Figure 9. Comparison of simulation with experimental results for (a) unnotched tension stiffness and (b) unnotched tension static strength.
Figure 9. Comparison of simulation with experimental results for (a) unnotched tension stiffness and (b) unnotched tension static strength.
Polymers 13 03491 g009
Figure 10. Comparison of the experimental and simulation responses for UNT static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Figure 10. Comparison of the experimental and simulation responses for UNT static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Polymers 13 03491 g010
Figure 11. Comparison of simulation with experimental results for (a) unnotched compression stiffness and (b) unnotched compression static strength.
Figure 11. Comparison of simulation with experimental results for (a) unnotched compression stiffness and (b) unnotched compression static strength.
Polymers 13 03491 g011
Figure 12. Comparison of the experimental and simulation responses for UNC static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Figure 12. Comparison of the experimental and simulation responses for UNC static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Polymers 13 03491 g012
Figure 13. Comparison of simulation with experimental results for (a) open hole tension stiffness and (b) open hole tension static strength.
Figure 13. Comparison of simulation with experimental results for (a) open hole tension stiffness and (b) open hole tension static strength.
Polymers 13 03491 g013
Figure 14. Comparison of the experimental and simulation responses along with images of simulated fiber and matrix damage distribution at failure for OHT static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Figure 14. Comparison of the experimental and simulation responses along with images of simulated fiber and matrix damage distribution at failure for OHT static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Polymers 13 03491 g014
Figure 15. Comparison of simulation with experimental results for (a) open hole compression stiffness and (b) open hole compression static strength.
Figure 15. Comparison of simulation with experimental results for (a) open hole compression stiffness and (b) open hole compression static strength.
Polymers 13 03491 g015
Figure 16. Comparison of the experimental and simulation responses along with images of simulated fiber and matrix damage distribution at failure for OHC static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Figure 16. Comparison of the experimental and simulation responses along with images of simulated fiber and matrix damage distribution at failure for OHC static strength of (a) [30/60/90/−60/−30]2S, (b) [0/45/90/−45]2S, (c) [60/0/−60]3S coupons.
Polymers 13 03491 g016
Table 1. Material properties of the carbon fiber (IM7) and matrix (977-3).
Table 1. Material properties of the carbon fiber (IM7) and matrix (977-3).
Material PropertiesValue
Carbon Fiber (IM7)
Longitudinal Tensile Modulus Ef1 (GPa)276
Transverse Tensile Modulus Ef2 (GPa)21
Longitudinal Compressive Modulus Ef1 (GPa)204
Transverse Compressive Modulus Ef2 (GPa)21
In-plane Shear Modulus Gf12 (GPa)29
Transverse Modulus Gf23 (GPa)7
Poisson’s ratio vf120.31
Poisson’s ratio vf230.28
Longitudinal Tensile strength XF (MPa)4082
Longitudinal Compressive strength XFP (MPa)2682
Fiber volume fraction0.592
Matrix (977-3)
Elastic modulus Em (GPa)3.5
Elastic Poisson’s ratio vm0.35
Tensile strength XM (MPa)81
Compressive strength XMP (MPa)131
Table 2. Tensile yield stress and yield strain data of the matrix (977-3) at different stages for the multi-linear damage model for the matrix.
Table 2. Tensile yield stress and yield strain data of the matrix (977-3) at different stages for the multi-linear damage model for the matrix.
StagesYield Stress
(MPa)
Yield Strain
(m/m)
Stage 182.160.03
Stage 2100.830.08
Stage 3113.150.16
Stage 4119.500.30
Stage 5124.730.45
Stage 6131.450.62
Stage 7135.930.76
Stage 8142.660.92
Stage 9146.761.06
Stage 10153.861.21
Table 3. Material properties of laminae (IM7/977-3).
Table 3. Material properties of laminae (IM7/977-3).
Material PropertiesValue
Longitudinal Tensile Modulus Ef1 (GPa)164
Transverse Tensile Modulus Ef2 (GPa)8.9438
Longitudinal Compressive Modulus Ef1 (GPa)121.6
Transverse modulus Ef2 (GPa)8.9236
In-plane Shear Modulus Gf12 (GPa)4.3516
Transverse Modulus Gf23 (GPa)2.9906
Poisson’s ratio vf120.324
Poisson’s ratio vf230.417
Table 4. Unnotched tension results.
Table 4. Unnotched tension results.
Unnotched Tension Static StrengthUnnotched Tension Stiffness
TeamStatic
Strength
[30, 60, 90, −60, −30]2S[0, 45, 90, –45]2S[60, 0, −60]3SElastic Modulus[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3S
Experimentσmax
(Mpa)
4738661005E (Gpa)386060
MMFσmax
(Mpa)
562893900E (Gpa)37.859.9560.5
GENOAσmax
(Mpa)
508828944E (Gpa)43.065.967.0
Helius PFAσmax
(Mpa)
679879877E (Gpa)40.261.361.6
DCNσmax
(Mpa)
506944997E (Gpa)40.561.661.9
MDS-Cσmax
(Mpa)
462890987E (Gpa)42.256.162.8
MAC/GMCσmax
(Mpa)
474897951E (Gpa)39.059.659.8
ESTσmax
(Mpa)
6038251009E (Gpa)40.561.161.0
NCYLσmax
(MPa)
4628561011E (GPa)39.760.661.5
BSAM/MICσmax
(MPa)
4328581113E (GPa)40.661.561.8
EHMσmax
(MPa)
5229111014E (GPa)39.960.461.7
Table 5. Unnotched compression results.
Table 5. Unnotched compression results.
Unnotched Compression Static StrengthUnnotched Compression Stiffness
TeamStatic Strength[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3SElastic Modulus[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3S
Experimentσmax (MPa)382603765E (GPa)334849
MMFσmax (MPa)427547548E (GPa)34.752.1852.24
GENOAσmax (MPa)439610479E (GPa)30.244.445.2
Helius PFAσmax (MPa)531619622E (GPa)38.150.950.9
DCNσmax (MPa)462560557E (GPa)36.553.851.4
MDS-Cσmax (MPa)405625836E (GPa)42.762.963.1
MAC/GMCσmax (MPa)350583649E (GPa)33.35151.2
ESTσmax (MPa)454618720E (GPa)35.252.352.3
NCYLσmax (MPa)428634632E (GPa)34.852.752.3
BSAM/MICσmax (MPa)485581551E (GPa)35.252.852.4
EHMσmax (MPa)425605602E (GPa)34.552.852.2
Table 6. Open hole tension results.
Table 6. Open hole tension results.
Open Hole Tension Static StrengthOpen Hole Tension Stiffness
TeamStatic Strength[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3SElastic Modulus[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3S
Experimentσmax (MPa)409554543E (GPa)324849
MMFσmax (MPa)402529466E (GPa)32.448.348.8
GENOAσmax (MPa)405543435E (GPa)38.656.457.2
Helius PFAσmax (MPa)400524485E (GPa)3550.650.9
DCNσmax (MPa)386458428E (GPa)35.250.951.2
MDS-Cσmax (MPa)425550558E (GPa)34.751.851
MAC/GMCσmax (MPa)378464433E (GPa)33.749.148.9
ESTσmax (MPa)409528489E (GPa)35.451.151.1
NCYLσmax (MPa)373557502E (GPa)34.550.350.9
BSAM/MICσmax (MPa)388553551E (GPa)34.148.749.3
EHMσmax (MPa)449558502E (GPa)34.650.850.8
Table 7. Open hole compression results.
Table 7. Open hole compression results.
Open Hole Compression Static StrengthOpen Hole Compression Stiffness
TeamStatic Strength[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3SElastic Modulus[30, 60, 90, −60, −30]2S[0, 45, 90, −45]2S[60, 0, −60]3S
Experimentσmax (MPa)295341358E (GPa)304444
MMFσmax (MPa)308366388E (GPa)32.346.346.8
GENOAσmax (MPa)323363380E (GPa)27.438.439.1
Helius PFAσmax (MPa)283308299E (GPa)30.543.443.6
DCNσmax (MPa)297329296E (GPa)36.152.252.5
MDS-Cσmax (MPa)271317352E (GPa)34.734.145.3
MAC/GMCσmax (MPa)330368320E (GPa)29.241.641.9
ESTσmax (MPa)296347331E (GPa)30.943.943.7
NCYLσmax (MPa)304341299E (GPa)29.841.941.9
BSAM/MICσmax (MPa)336380388E (GPa)29.641.741.9
EHMσmax (MPa)360393368E (GPa)31.845.244.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Masud, M.; Kharusi, M.S.M.A.; Ali, M.U.; Mubashar, A.; Hussain, S.J.; Tariq, A.; Rehman, G.U.; Akhtar, M.H.; Javeed, S. Prediction of the Ultimate Strength of Notched and Unnotched IM7/977-3 Laminated Composites Using a Micromechanics Approach. Polymers 2021, 13, 3491. https://doi.org/10.3390/polym13203491

AMA Style

Masud M, Kharusi MSMA, Ali MU, Mubashar A, Hussain SJ, Tariq A, Rehman GU, Akhtar MH, Javeed S. Prediction of the Ultimate Strength of Notched and Unnotched IM7/977-3 Laminated Composites Using a Micromechanics Approach. Polymers. 2021; 13(20):3491. https://doi.org/10.3390/polym13203491

Chicago/Turabian Style

Masud, Manzar, Moosa S. M. Al Kharusi, Muhammad Umair Ali, Aamir Mubashar, Shaik Javeed Hussain, Adnan Tariq, Gulfam Ul Rehman, Mahmood Hassan Akhtar, and Shama Javeed. 2021. "Prediction of the Ultimate Strength of Notched and Unnotched IM7/977-3 Laminated Composites Using a Micromechanics Approach" Polymers 13, no. 20: 3491. https://doi.org/10.3390/polym13203491

APA Style

Masud, M., Kharusi, M. S. M. A., Ali, M. U., Mubashar, A., Hussain, S. J., Tariq, A., Rehman, G. U., Akhtar, M. H., & Javeed, S. (2021). Prediction of the Ultimate Strength of Notched and Unnotched IM7/977-3 Laminated Composites Using a Micromechanics Approach. Polymers, 13(20), 3491. https://doi.org/10.3390/polym13203491

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