Next Article in Journal
Experimental Study on the Thermoplastic Dripping and Flame Spread Behaviors of Energized Electrical Wire under Reduced Atmospheric Pressure
Next Article in Special Issue
An Extended Thickness-Dependent Moisture Absorption Model for Unidirectional Carbon/Epoxy Composites
Previous Article in Journal
Fabrication of an Extremely Cheap Poly(3,4-ethylenedioxythiophene) Modified Pencil Lead Electrode for Effective Hydroquinone Sensing
Previous Article in Special Issue
Adhesion Improvement between PE and PA in Multilayer Rotational Molding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linear-Nonlinear Stiffness Responses of Carbon Fiber-Reinforced Polymer Composite Materials and Structures: A Numerical Study

1
Department of Aerospace Engineering, Universiti Putra Malaysia, UPM Serdang 43400, Selangor Darul Ehsan, Malaysia
2
Institute for Nanomaterials, Advanced Technologies and Innovation (CXI), Technical University of Liberec (TUL), Studentska 2, 461 17 Liberec, Czech Republic
3
School of Mechanical Engineering, Universiti Teknologi Malaysia, Johor Bahru 81310, Malaysia
4
Department of Mechanical Engineering, Universiti Putra Malaysia, UPM Serdang 43400, Selangor Darul Ehsan, Malaysia
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(3), 344; https://doi.org/10.3390/polym13030344
Submission received: 4 January 2021 / Revised: 13 January 2021 / Accepted: 15 January 2021 / Published: 22 January 2021
(This article belongs to the Special Issue Damage Mechanics of Polymer Composites)

Abstract

:
The stiffness response or load-deformation/displacement behavior is the most important mechanical behavior that frequently being utilized for validation of the mathematical-physical models representing the mechanical behavior of solid objects in numerical method, compared to actual experimental data. This numerical study aims to investigate the linear-nonlinear stiffness behavior of carbon fiber-reinforced polymer (CFRP) composites at material and structural levels, and its dependency to the sets of individual/group elastic and damage model parameters. In this regard, a validated constitutive damage model, elastic-damage properties as reference data, and simulation process, that account for elastic, yielding, and damage evolution, are considered in the finite element model development process. The linear-nonlinear stiffness responses of four cases are examined, including a unidirectional CFRP composite laminate (material level) under tensile load, and also three multidirectional composite structures under flexural loads. The result indicated a direct dependency of the stiffness response at the material level to the elastic properties. However, the stiffness behavior of the composite structures depends both on the structural configuration, geometry, lay-ups as well as the mechanical properties of the CFRP composite. The value of maximum reaction force and displacement of the composite structures, as well as the nonlinear response of the structures are highly dependent not only to the mechanical properties, but also to the geometry and the configuration of the structures.

Graphical Abstract

1. Introduction

As a widely used material in advanced industries—such as aerospace, automotive, etc.—fiber-reinforced polymer (FRP) composites have been the subject of many studies [1,2]. Most of these studies have employed a combination of experimental and numerical methods to predict the material or structural behaviors of composites under different loading conditions with respect to their industrial applications [3,4]. Finite element (FE) simulation is the mostly used numerical method because of the availability of well-stablished models and the accuracy of the results [5]. In real-time simulation of structures, the necessary requirement including the precise material properties and constitutive model, as well as correct model configuration, boundary, and loading conditions [1,6,7,8,9]. On the other hand, every FE model and simulation process should be validated through a logically accepted comparison of the FE results with the actual data obtained through experiment or a valid theoretical model [2,4,10,11,12].
In many case studies—including material or structures at different micro-macro scales—the validation of FE models has been frequently done by comparing the stiffness curve of solid objects that obtained through FE simulation and compared with their actual behavior measured in the experiments or numerical approach [1,11,13]. This curve is called either stiffness response or load-deformation/displacement, that generally comprise of two parts, that start with an initial linear response to a maximum reaction force at specific displacement, and continued with a nonlinear part to a stage where the load drop is seen [10,14,15]. Some researchers employed the load-displacement curve to validate the numerical simulation of the material behavior [9,16], while others validate the structural response [1,17,18,19]. In this regard, different modes of loading such as tensile, compression, and bending has been investigated [1,17,18,19].
Numerical modelling of the elastic-damage behavior of composite materials and structure is of importance, because the occurrence of invisible damage inside the composite materials and structures such as delamination, fiber breakage, and matrix cracking could cause catastrophic failure of the whole structure [20,21]. Several failure criteria have been derived for damage and failure of FRP composite structures including Chang-Lessard [22] and Greszczuk [23] fiber bucking models, Puck [24] fiber breakage models, Lee [25] fiber failure in both tension and compression, Hashin-Rotem [26] and Shahid-Chang [27] matrix cracking model, Maimí et al. [28] matrix crushing model, etc. In continuum damage mechanics [29,30], the necessity of damage evolution law to be coupled with a suitable failure criteria, in order to estimate the degradation of mechanical properties to fracture is highlighted. Xue and Kirane 2020 [31] and R. Koloor et al., 2018 [1] have developed damage evolution law that enables the prediction of damage initiation and propagation in FRP composites, that are very suitable for estimation of the stiffness response of composite materials and structures.
Although the theory of composite materials and damage models are used to obtain the stiffness behavior of the structure; however, no studies have been implemented to investigate and describe the stiffness response of the composites in detail and its dependency to the elastic-damage properties at material and structural levels. In this regard, the suitability and applicability of a composite structures to be used for validation of a new damage model or characterization of the mechanical properties through numerical approach, is not determined. On the other hand, the effectiveness of each properties in determination of the mechanical behavior at structural level, has not been studied. Therefore, this study aims to investigate the linear-nonlinear stiffness responses of a carbon fiber-reinforced polymer (CFRP) composites at two levels of material and structure, and describe its dependency to elastic and damage model parameters. For this purpose, four FE models of different composite cases (material and structural levels) are selected for simulation at full damage and failure states, in which the effect of sets of individual and group of elastic and damage parameters are examined. In total, 48 FE models representing the combination of different specimens and properties are run and the results are compared systematically, to provide a comprehensive analysis of stiffness response dependency to the mechanical properties. The results provide important insight into the analysis of composite structures when used for examination of design cases, validation of new theoretical models, and characterization of the mechanical properties.

2. Damage Model of FRP Composite Material

In modelling aspect, each lamina is treated as a homogeneous orthotropic layer. Classical lamina theory is used to model composite lamina [32]. To describe the material properties and constitutive model of the unidirectional (UD) lamina, two coordinate systems are required, the material (1-, 2-, 3-axis) and global (x-, y-, z-axis) coordinates. In the material coordinate system, axis 1 and 2 are in the plan of lamina, where axis 1 is along the fibers and axis 2 is normal to the fibers. Therefore, the stress in an equivalent thin lamina (i.e., material coordinate), which is under plane stress condition, can be written in terms of the global strain as [17]
{ σ 11 σ 22 τ 12 } { σ ^ } = [ cos 2 θ sin 2 θ 2 sin θ cos θ sin 2 θ cos 2 θ 2 sin θ cos θ sin θ cos θ sin θ cos θ cos 2 θ sin 2 θ ] [ T ] × [ Q ¯ 11 Q ¯ 12 Q ¯ 16 Q ¯ 12 Q ¯ 22 Q ¯ 26 Q ¯ 16 Q ¯ 26 Q ¯ 66 ] [ Q ¯ i j ] × { ε x ε y γ x y } { ε }
where { σ } is the material (local) stress tensor, [ T ] is the axis rotation matrix, [ Q ¯ i j ] is the stiffness matrix which determines the mechanical properties of the lamina, { ε } is the strain tensor in the global coordinate, and θ is the lamina angle. In this equation, mainly the elastic properties that are obtained through standard tests [33], are used to compute the linear behavior of the composite laminate [17,34].
There are several constitutive damage models for defining the FRP composites failure [3,4,35,36,37,38]. A constitutive damage model of lamina [1] is applied for the simulation of elastic-damage behavior of the FRP composite materials and structures. In this model, each damage mode is defined by a bilinear curve that represents the elastic behavior and the softening process in each mode of loading (refer to Section 3).

2.1. Damage Initiation

The initiation of damage is predicted by the Hashin’s failure model in laminas at meso-scale [35]. The stress-based Hashin’s criteria are described for fiber and matrix phases in the lamina under plane stress condition, as follows:
For the fiber fracture, buckling or kinking:
( σ ^ 11 X T ) 2 + ( τ ^ 12 S L ) 2 = d f t ;     if   σ ^ 11 0   ( Tension )
( σ ^ 11 X C ) 2 = d f c ;     if   σ ^ 11 < 0   ( Compression )
For the matrix cracking and crushing:
( σ ^ 22 Y T ) 2 + ( τ ^ 12 S L ) 2 = d m t ;     if   σ ^ 22 0 ( Tension )
( σ ^ 22 2 S T ) 2 + [ ( Y C 2 S T ) 2 1 ] ( σ ^ 22 Y C ) + ( τ ^ 12 S L ) 2 = d m c ; if   σ ^ 22 < 0   ( Compression )
In these equations, σ ^ i j are the effective stresses in the meso-scale lamina, XT, YT, XC, YC, SL, and ST are the strength properties, and the parameters d f t , d f c and d m t , d m c shows the internal damage variables, in which the index represents the fiber (f) and matrix (m) phases, and the power expresses tensile (t) and compression (c) loading modes. The strength properties are normally extracted through standard test method such as tension, compression, and shear experiment on 0o and 90o standard samples [33].

2.2. Post-Damage Initiation

The orthotropic intrinsic of the FRP Composite lamina causes the occurrence of mixed-mode damage and multiple failures [17,34]. Once damage initiated, this effect could be applied by updating the elastic stress tensor shown in Equation (1), by multiplying by a damage variable matrix as
σ ^ i j =   { σ i j ,    b e f o r e   d a m a g e   i n i t i a t i o n D σ i j ,    I f   a n y   d a m a g e   i n i t i a t e d
where σ ^ i j is the effective stress and D is the damage variable matrix, which obtained from the strain equivalence hypothesis, as [7,29,34,37,39]
D = [ 1 / ( 1 d f ) 0 0 0 1 / ( 1 d m ) 0 0 0 1 / ( 1 d s ) ]
in which d f ,   d m and d s derive from the internal damage variables in the lamina (i.e., Equations (2a)–(2d))
d m = { d m t   if   σ ^ 22 0 , d m c   if   σ ^ 22 < 0 ,    d f = { d f t   if   σ ^ 11 0 , d f c   if   σ ^ 11 < 0 , d s = 1 ( 1 d m t ) ( 1 d m c ) ( 1 d f t ) ( 1 d f c )
The effective stress parameters are applied in Hashin model to predict the initiation of damage at different levels.

2.3. Damage Propagation

The propagation of damage in the FRP composite lamina is modeled by applying energy-based criteria to define the softening behavior to the final failure [34]. Based on these criteria, during the damage evolution, the relation between the equivalent stress and displacement in each failure mode is expressed by [34]:
Fiber tension ( σ ^ 11 0 ) :
σ e q . = ( ( σ 11 o ε 11 o + τ 12 o ε 12 o ) 2 ( ( L c ( σ 11 o ε 11 o + τ 12 o ε 12 o ) ) 2 G C X T ) × ( ε 11 o 2 + ε 12 o 2 ) ) × ( δ e q . 2 G C X T   ε 11 o 2 + ε 12 o 2   σ 11 o ε 11 o + τ 12 o ε 12 o )
Fiber compression ( σ ^ 11 < 0 ) :
σ e q . = ( σ 11 o 2 ( L c ε 11 o σ 11 o 2 G C X C ) ) × ( δ e q . 2 G C X C   σ 11 o )
Matrix tension ( σ ^ 22 0 ) :
σ e q . = ( ( σ 22 o ε 22 o + τ 12 o ε 12 o ) 2 ( ( L c ( σ 22 o ε 22 o + τ 12 o ε 12 o ) ) 2 G C Y T ) × ( ε 22 o 2 + ε 12 o 2 ) ) × ( δ e q . 2 G C Y T   ε 22 o 2 + ε 12 o 2 σ 22 o ε 22 o + τ 12 o ε 12 o )
Matrix compression ( σ ^ 22 < 0 ) :
σ e q . = ( ( σ 22 o ε 22 o + τ 12 o ε 12 o ) 2 ( L c ( σ 22 o ε 22 o + τ 12 o ε 12 o ) 2 G C Y C   ) ( ε 22 o 2 + ε 12 o 2 ) ) × ( δ e q . 2 G C Y C   ε 22 o 2 + ε 12 o 2   σ 22 o ε 22 o + τ 12 o ε 12 o )
where L c is the length of a first-order element on the lamina, and this form of brackets, <•>, represents the Macaulay bracket. In these equations G C X T , G C X C , G C Y T , and G C Y C are the values of fracture energy in fiber and transvers-to-fiber directions for each mode of loadings, in which has been obtained through standard [33] and innovative [40] experimentations. The parameter σ i j o , τ i j o , and ε i j o are the effective stresses, which computed at the damage initiation using the following damage evolution parameter [34,41]
d = δ eq f ( δ eq δ eq 0 ) δ eq ( δ eq f δ eq 0 ) δ eq δ eq 0
where δ eq 0 represents the equivalent displacement at the onset of damage (i.e., d = 0), and δ eq f is the equivalent displacement at failure of a material point (i.e., d = 1). Further explanation about these criteria for each level of elastic-damage behavior of the FRP composite, i.e., damage initiation, post damage initiation, and damage propagation, is found in [1,17,34].

3. Mechanical Properties and Damage Model Parameters

The mathematical formulation of the constitutive damage model describing the elastic to damage and fracture of FRP composites, is described in Section 2. The schematic bilinear curve representing the constitutive model is shown in Figure 1a. This constitutive curve includes a linear section (line OA) to model the elastic behavior, which depends to nine properties of elastic moduli, shear moduli, and Poisson’s ratio (Equation (1)) [2,32]. The maximum elastic behavior is concluded in the yielding point ( δ o , σ o ) which represents the onset of damage and depends to the six strength properties of matrix and fiber in tension and compression loading conditions (Equations (2a)–(2d)) [35,42,43]. Then, a second linear line (AB) representing the damage evolution process to fracture ( δ f ) that is calculated based on fracture energy value in each different failure modes (Equations (4a)–(4d)) [1,34,44,45]. These mechanical properties are essential to compute a full elastic-to-failure process at material point, that are obtained through standard test processes [33]. In the past three decades, many studies (e.g., [46,47,48,49,50,51,52,53]) have used these mechanical parameters to investigate the behavior of composite parts, in which a range of different value with 50–200% tolerance have been reported. Here, a schematic image of the average values of those works in the form of different set of parameters are present in Figure 1b. In this figure, the central data shown using green-line is the reference data, which are validated and presented by the same author in the previous works [1,17,44,45], as listed in Table 1. It is worth mentioning that Figure 1b provides an general view of the scattered data (elastic and damage parameters), and highlights the range of the considered values in the view of a representation of the constitutive damage model. It should be noted that the authors are aware of the fact that, these values may change for different carbon fibers or thermoset matrix constituents that are used in the manufacturing of CFRP composite structures. However, the objective is to show the variations of mechanical parameters in view of its effect in the constitutive model.
As mentioned, the main mechanical properties and damage model parameters that are used as the ref data (green-color curve, Figure 1) for the simulation of CFRP composite materials and structures, are extracted from the published work of the same authors [1,17,44,45], which are rigorously examined and validated through mechanical analysis of different structures. These properties are listed in Table 1. To define a set of elastic and damage properties to build a systematic study and investigate the effect of mechanical properties on the stiffness response of the composite materials and structures, two different ranges of properties (along the reference properties) one 50% of the ref data (one-half ref) and second 1.5 times of the ref properties (three-half ref) are selected. It is important to make a note that, mechanical properties similar to the range of one-half ref and three-half ref properties, have been reported as the properties of CFRP composite in the previous studies. These sets of different properties projected from different values reported in the previous studies (e.g., [46,47,48,49,50,51,52,53] as shown in Figure 1b), are used for the parametric study on the stiffness responses of composite materials and structures (Figure 1), as listed in Table 1. The detail information about the property combinations is provided in Figure 2, using both tables and schematic view of the bilinear damage model curves. In all cases, a full damage model comprising elastic, strength, and energy properties are used. The CFRP composite properties used in this investigation are considered from a real structure that has been manufactured using pre-impregnated laminas (prepreg) with high modulus carbon fibers with epoxy resin (M40J fibers and NCHM 6376 resin, Structil France) [1,17], in order to make a realistic numerical study.
In the first set of properties (Figure 2a), the variation of the elastic properties on the stiffness response of the composite material and structures, is examined while the strength and energy properties selected from the Ref data. In the second set of properties (Figure 2b), the variation of the strength properties is examined on the stiffness response, while the elastic and energy properties are assumed similar to the ref data. In the third set of properties (Figure 2c), the variation of the energy properties is examined, while the elastic and strength properties are defined based on the ref data. Finally, in the fourth set of properties (Figure 2d), the variation of the strength and energy properties are examined on the stiffness response of the composite cases, while the elastic properties are assumed similar to the ref data.

4. Finite Element Model and Simulation Process

In this study, the stiffness response of CFRP composite material and structures are investigated with respect to the variation of the elastic-damage properties. Composites are normally tested at different levels of material, structure, or super-structure. A list of specimens at different material and structural levels that are considered in the computing activities of the present study, are provided in Table 2, and their configuration and geometry are shown in Figure 3.
The material level tests such as UD 0° or 90° FRP composite under tensile or compressive loadings [33], are normally used for the extraction of mechanical properties, which are not a good target for this study. However, an angle laminate made of several UD lamina with an arbitrary angle which is the same in all layers, is a good case for the simulation of composite at material level. Because, the response of the material is reflected by a combination of normal and shear stresses and deformations. In this regard, a beam made of five layers of 45° UD lamina is considered for data analysis at materials level, as shown in Figure 3a. The stiffness response at structural level is investigated by using a MD composite plate under bending, as shown in Figure 3b. Two cases are considered with different fiber arrangements and thicknesses, to highlight the structural damage growth dependency to the mechanical properties. A second stage of the investigation is implemented by modelling of a super-structure under flexural loading condition. In this regard, a hat-stiffened part of composite structures that are used in aerospace application, is modelled under four-point bending, as shown in Figure 3c. The flexural load applied from the bottom part of the structure is to apply the load in a more convenient way, because in the real application, the loading will not appear on the hat section.
The four different sets of properties shown in Figure 2, are used in four models of the composite at the material and structural levels described in Figure 3. As the result of that, 48 FE models are developed, which are listed and coded in Table 3. The case codes are used to name the results of stiffness response for brevity and referment to each case.

4.1. Finite Element Simulation

In this study, CFRP laminated composite materials and structures are simulated using a multilayer construction FE model. In this model each lamina is modelled as a homogeneous orthotropic layer, that created using 8-node three-dimensional continuum shell elements in FE environment. All laminas are stacked with their specific orientation on each other, in which a weak physical interface is considered between the adjacent layers, as shown in Figure 4. This FE model that has been validated in the previous study [1] for the CFRP laminated composite with the reference constitutive model, is applied to the simulation of the material and the structures described in Figure 3. The multilayer structure is more suitable for the simulation of the progressive damage behavior of the CFRP composite laminate manufactures using the autoclave or prepreg method, in which the prepreg laminas are chemically attached together in the curing procedure. The detailed descriptions about the model creation, simulation process, mesh convergence study, and the effect of the intralayer interface can be found in [17,34].
The geometry of the materials and structural models were created based on the dimensions of the specimens listed in Table 2. As an example, the 3D geometry, load, and boundary condition of the composite structure (Profile, Table 2) with the mesh configuration are demonstrated in Figure 5. The specimens were meshed using SC8R elements [37], in which the average element size is 0.2 mm in the damage areas. For the loadings tools and supports a frictionless contact with the specimen is considered in which rigid body element with R3D4 types are applied [37]. In this simulation the interfaces between the plies are perfectly bonded, in which zero relative motion is considered. To eliminate the effect of mesh size on the results of the simulations, mesh convergence study is done in two levels; firstly, the optimum element size is investigated for the pure elastic simulation, and then the final element size is found for the elastic-damage model, in which the detail is provided elsewhere [17]. The mesh configurations of the other cases are created with the same procedure according to the geometry shown in Figure 3, in which example of them could be find previous published work [1,45].

4.2. Validation of the Finite Element Model

Numerical method has become one of the important tools for design and analysis of composite structures. An accurate numerical analysis of FRP composite materials and structures, very much depend on the credibility of the input-data, material model, and simulation processes. The material properties used in this study, has been measured experimentally and examined thorough several cases of multidirectional CFRP composite structure under tensile, bending and flexural loading, as reported in the published literature [33,34]. The material model, and FE simulation process are also used in model creation of several CFRP composite cases which have been validated experimentally [1,8,17,44,45,54,55]. In this study, the same material properties, and simulation processes have been utilized to create three distinct composite cases at different levels, to highlight the dependency of the mechanical properties to the global stiffness response.

5. Results and Discussion

The FE results of the material and structural stiffness in term of global load–displacement curve are presented quantitatively and discussed in this section. Results are divided into two main parts, firstly the influences of elastic properties are shown, and secondly the effect of damage parameters will be discussed.

5.1. Stiffness Response with Respect to the Variation of Elastic Properties

The load-displacement/deflection responses of the CFRP composite under tensile and flexural loads are presented in Figure 6. The legend provided in each picture has been given in Table 3. The reference curve that has been validated by the experiment is shown using green color, while the results obtained from one-half and third-half elastic properties are shown using blue and red color curves, respectively. The effect of elastic properties appears significant in changing the stiffness response of the structures especially in the linear region. Since the composite laminate is selected with 45° laminas angle specification, both elastic and shear moduli affect the slope of the global load-displacement response. However, the values of shear moduli (1.2–5.1 MPa) are almost 3–30 time less than elastic moduli values (3.6–158 MPa), as given in Table 1, which have led to the disproportional slope in one-half and third-half responses compare to the ref data.
Considering that in this section, the focus is only on the variation of elastic properties while the damage properties are kept the same in all cases, it is worth to highlight that while the maximum reaction force of the specimens are computed similarly (≈700 N) at the material level (Figure 6a), it has been changed significantly at the structural level (Figure 6b–d). In addition, it should be noted that, the deformation at maximum load is also affected to be larger for thick plates (Figure 6c) and profile cases (Figure 6d). Such effects were seen as the response to geometry changes in the structure cases that is intensified for the structures with larger geometry.

5.2. Stiffness Response with Respect to Variation of Damage Parameters

In this section, the stiffness responses of CFRP composite at material and structural levels are investigated with respect to the variation of the damage properties (strength and fracture energy), while the elastic properties are kept similar to the ref data. The FE result of each specimen is provided separately and discussed. The results of each specimen are shown using three graphs that illustrate first the effect of strength variation, second the effect of fracture energy variations, and third the effect of both strength and energy variations, that provides a systematic understanding about the effect of the damage parameters in detail. The reference response is shown using green color curve, while one-half damage properties and three-half properties are shown using blue and red color curves, respectively. The code of cases provided in the legend, has been given in Table 3. The values of horizontal and vertical axes are unified for all cases, to show a consistent comparison between the results.

5.2.1. CFRP Composite Beam under Tension

The load-displacement responses of the CFRP composite beam under tensile load are presented in Figure 7. Although the results seem to be simply understandable, however it will help to comprehend the variation of stiffness response in complex condition at structural level shown in Figure 8, Figure 9 and Figure 10.
Figure 7a indicates a direct effect on the maximum reaction force of the composite beam based-on the variation of strength property. Considering a bigger or smaller strength value, causes the bilinear curve to have a shorter or longer softening curve (Figure 2b), respectively, and therefore the failure point ( δ f ) will be computed in a smaller or bigger range, respectively. The result of such assumption has affected the nonlinear part of the stiffness responses to gradually (B2-S-1) or sharply (Be-S-2) be reduced (Figure 7a).
Figure 7b shows the effect of fracture energy variation on the nonlinear part of the stiffness curve, to sharply degrade as a result of smaller energy value, or smoothly degrade as a result of a bigger energy value. Therefore, the composite beam with the assumption of bigger fracture energy fails later due to the larger δ f value (Figure 2c).
Figure 7c illustrates the effect of both strength and energy values on the stiffness response of the composite beam. Considering the simultaneous variation of these two parameters, a similar value of displacement at failure ( δ f ) is calculated, as shown in Figure 2d. A proportional variation of these parameters with respect to ref data resulted in a smooth degradation of the stiffness response in the nonlinear section.

5.2.2. Thin Composite Plate under Three-Point Bending

The load-deflection responses of a thin CFRP composite plate under bending load are presented based on different damage properties, and shown in Figure 8. The results are interpreted with respect to the shapes of the bilinear models that varied based on different property assumptions as shown in Figure 2. It should be noted that, in thin composite plate, the specimen shows a better deformation flexibility before failure. The results are shown for the maximum deflection of 40 mm, where a negative slope in the stiffness response is computed.
Figure 8a shows the effect of variation on the strength property, which is obvious in the P1-S-1 with 50% smaller strength value compare to ref data. However, a larger value of strength (P1-S-2) has not significantly changed the global stiffness response, which is due to the less capacity of the composite structure (thin plate) in generation of stress value to the strength level. It is notable that, the stress distribution is higher at the middle section of the plate where normally shows a bigger deflection; however, this value changes from a location to another. Therefore, in some locations the stress may exceed the strength, however, those locations should be sufficiently big to change the global stiffness response of the structure. This fact is more effective when using smaller strength value, that increased the area of the regions where stress exceeds the strength of the structure, as highlighted in the result of P1-S-1 case.
Figure 8b indicates the effect of fracture energy variation on the global stiffness response of the thin composite plate. Results indicated no changes in the stiffness response while the fracture energy value is significantly changed. From the damage mechanics point of view, the effect of fracture energy will appear if the damage initiation criteria are satisfied (Equations (2a)–(2d)), then the damage propagation is computed in the simulation process. In addition, the region undergone damage propagation should be large enough to affect the global stiffness response of the structures, which is not applicable in the considered thin plate. Therefore, the reason for the degradation of the stiffness response to a negative value, is due to the structural behavior of the thin plate and its flexibility. As the thin plate is loaded under three-point bending, a bigger area of the plate moved in between the support-span through the plate/roller sliding process (as the deflection getting larger ≈ 30–40 mm). Such a movement increases the effective length of the plate against bending, and results in lowering the global stiffness response. It is worth to note that, this plate is not a good case for research when the objective is the constitutive modeling or damage characterization of the composite materials, as the structure is not capable to reflect the structural response based on the material property variation.
Figure 8c demonstrates the effect of both strength and energy values on the global stiffness of the thin composite plate. As explained before, the single variation of strength and fracture energy have not been able to change the stiffness response, which is also seen when both parameters are changed. The thin composite plate is introduced as a good example, that should not be used for the damage modeling and characterization of composite materials.

5.2.3. Thick Composite Plate under Three-Point Bending

The load-deflection responses of a thick CFRP composite plate under three-point bending load based on various damage properties, are presented in Figure 9. The results are interpreted with respect to the shapes of the bilinear model that varied based on different property assumptions as shown in Figure 2. It should be noted that, the reaction force obtained from the thick composite plate is much higher than thin plate (Section 5.2.2). Therefore, it is expected to capture very much different structural responses compare to the previous case, which helps to highlight the objective of this study. It is worth mentioning that a notable higher global stiffness is computed at 4–7 mm deflection in all cases, which is the result of saddle deformation of the composite plate during the bending process [1,17]. A small magnitude of this effect is also seen in the thin plate case (Figure 8, at 11 mm deflection). This effect is normally bigger for specimen with larger width.
Figure 9a indicates a significant influence of the strength variation, which not only changes the maximum reaction force of the structure, but also notably changes the global stiffness (slope of the curve) of the structure (although the elastic properties of all three cases are kept similar to ref data). The deviation of the stiffness value from the ref data is visible at around 500 N force in P2-S-1 case, and around 1000 N for P2-S-2 case. Such deviation is smaller for P2-S-1 case with lower strength value, because that is the point where the stress exceeds from the strength in the middle of the structure which causes the saddle deformation. It should be noted that, the stress exceeding the strength is a sign of damage propagation, that could also have effect on the reduction of the structural stiffness.
Figure 9b shows the effect of fracture energy variation on the global stiffness response of the thick composite plate. It is very interesting to see that, although the strength values are kept similar to ref data in all cases, but the maximum reaction force of the structures has changed for +10% in P2-E-1 case, and −6% in P2-E-2 case compared to ref data. This effect is due to a steeper softening curve and smaller displacement at failure ( δ f ) for P2-E-1 case, as shown in Figure 2c, which result in the propagation of damage in a bigger area.
Figure 9c illustrates the effect of both strength and energy values on the global stiffness of the thick composite plate. The influence of single damage parameters variation (Figure 9a,b), is visible with slightly higher magnitude. The maximum reaction force and displacement of the P2-SE-2 case are notably higher than P2-S-2. Oppositely, the maximum reaction force and displacement of the P2-SE-1 case are notably lower than P2-S-1.
It is worth to mention that, similar example to the thick composite plate in this study, could be a good case for researches related to a new damage modeling and characterization of FRP laminated composites. Using the thick composite plate, the existence of any inaccuracy in the determination of material properties or damage model could be highlighted, as the property changes result in the variation of the linear and nonlinear response of the composite structure.

5.2.4. Composite Profile under Flexural Loading

The load–deflection responses of hat-stiffened composite structure under the four-point bending load, based on various damage properties, are presented in Figure 10. The results are construed with respect to various shapes of the bilinear models as shown in Figure 2. Figure 10a demonstrates the effect of the strength variation on the global stiffness response, maximum reaction force, and displacement of the profile structure. Similar to the case of thick composite plate (Figure 9), a notable change in the global stiffness (slope of the linear part of the load–displacement curve) is computed at 8–9 mm deflection in all cases compare to ref data. The deviation of the stiffness value from the ref data is visible at around 400 N force in Pr-S-1 case, which is at around 800 N for Pr-S-2 case. Similar to the case of thick composite plate, the deviation is smaller for Pr-S-1 case as the strength is lower, which causes early stress exceeds of the strength in the mid-section of the structure, that results in lower capacity of the structure to resist the load. It is worth to mention that, the nonlinear portion of the Pr-S-2 is very short (displacement 54–57.5 mm), which is due to the steeper softening curve and smaller value of displacement at failure ( δ f ), as described in Figure 2b.
Figure 10b indicates the effect of fracture energy variation on the global stiffness response of the composite structure. It is important to highlight that, although the strength values are kept similar in all cases, however, the maximum reaction force of the structure is computed 17% lower in Pr-E-1 case compared to ref data. On the other hand, the nonlinear part of the stiffness response is shown around 300 N higher reaction force compared to ref data, before the final failure of the structure. As mentioned before, this effect is due to a steeper softening curve (and smaller displacement at failure ( δ f )) of the ref data compared to Pr-E-2 case, as shown in Figure 2c.
Figure 10c illustrates the effect of both strength and energy values on the global stiffness of the composite structure. The effect of both damage parameters that were individually illustrated in Figure 10a,b, is visible with slightly higher magnitude. Similarly, the maximum reaction force, as well as the displacement range in the nonlinear section (54–77 mm) of the Pr-SE-2 case is very much bigger than Pr-S-2. Notably, the maximum reaction force of the Pr-SE-1 case is lower than Pr-S-1.
The investigation on this structure provides important data on the dependency of the stiffness response of the structure to the damage model and parameters. This is important as the design of composite structures inevitably depends on the damage investigation and the determination of the inelastic behavior of the structure with respect to the stiffness response.

6. Concluding Remarks

This study focused on the dependency of the stiffness response or load-deformation/displacement behavior of CFRP composite materials and structures, to the elastic and damage properties. This is because the stiffness behavior is the most important response that frequently has been utilized for the validation of the theoretical constitutive models, damage theories and mechanical characterization of composites, as well as for the design of composite structures under severe loading. In this regard, the effects of sets of individual/group of mechanical properties including the elastic and shear moduli, strengths, and fracture energies, are studied on the stiffness response of CFRP composite at material and structural levels. These properties are responsible to model different intralaminar failure modes of matrix cracking/crushing, lamina shear failure, and fiber breakage/buckling in composite lamina at meso scale. Four different specimens including a [45]5 UD laminated composite beam under tension in material level, two thin and thick CFRP composite plates under three-point bending in structural level, and an aerospace case of hat-stiffened composite structure under flexural loading, were used in a FE environment to apply the elastic-damage model and compute the stiffness response based on different mechanical properties. A set of validated mechanical properties, model and simulation process were used to obtain a reference data, that was utilized to examine different possible scenarios of material and structural behaviors.
The result indicated a direct dependency of the stiffness response of the composite at material level to the mechanical properties. However, the stiffness behavior of the composite structures depends both on the structural configuration, geometry, lay-ups as well as the mechanical properties of the CFRP composite. The elastic properties directly affect the initial slope and linear response of the composite at material level. However, the composite structures showed two slopes in the linear part of the stiffness curve, in which the initial slope (e.g., Figure 9c, within the deflection values of 0–5 mm) depends to the elastic properties, while the secondary slope (e.g., Figure 9c, within the deflection value of 5 mm to the deflection at maximum load) depends to the geometry of the structure and elastic-damage properties. In this regard, the value of maximum reaction force and displacement of the composite structures are highly dependent not only to the mechanical properties, but also to the geometry and the configuration of the structure.
Since the stiffness curve is frequently used in the design and analysis of advanced structures, structure cases similar to the case in this studies are introduced as good examples, if the objectives of a research is the investigation of a new theoretical modeling or characterization of FRP composites to highlight the accuracy of material model and parameters (e.g., thick composite plate, Section 5.2.3), or if the objective is the design of composite part and defining the tolerance of the structure by using a mean to the mechanical performance (e.g., hat-stiffened composite structure, Section 5.2.4).

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

This research was supported by Universiti Putra Malaysia under Putra Grant no. GP/2018/9635100. Also, it was supported by the Ministry of Education, Youth, and Sports of the Czech Republic and the European Union (European Structural and Investment Funds Operational Program Research, Development, and Education) in the framework of the project “Modular platform for autonomous chassis of specialized electric vehicles for freight and equipment transportation”, reg. no. CZ.02.1.01/0.0/0.0/16_025/0007293. In addition, the research was supported by the Research University Grant from Universiti Teknologi Malaysia (grant no. HIR—Q.J130000.2424.04G40), in which the authors acknowledge the Ministry of Higher Education (MOHE) and Universiti Teknologi Malaysia for the approved fund which makes this important research viable and effective.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koloor, S.S.R.; Khosravani, M.R.; Hamzah, R.; Tamin, M. FE model-based construction and progressive damage processes of FRP composite laminates with different manufacturing processes. Int. J. Mech. Sci. 2018, 141, 223–235. [Google Scholar] [CrossRef]
  2. Rangappa, S.M.; Parameswaranpillai, J.; Siengchin, S.; Kroll, L. Lightweight Polymer Composite Structure: Design and Manufacturing Techniques; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  3. Kaddour, A.; Hinton, M.; Smith, P.; Li, S. The background to the third world-wide failure exercise. J. Compos. Mater. 2013, 47, 2417–2426. [Google Scholar] [CrossRef]
  4. Hinton, M.J.; Kaddour, A.S.; Soden, P.D. Failure Criteria in Fibre Reinforced Polymer Composites: The World-Wide Failure Exercise; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  5. Pidaparti, R.M. Engineering Finite Element Analysis; Morgan & Claypool Publishers: San Rafael, CA, USA, 2017. [Google Scholar]
  6. Schuecker, C.; Pettermann, H.E. A continuum damage model for fiber reinforced laminates based on ply failure mechanisms. Compos. Struct. 2006, 76, 162–173. [Google Scholar] [CrossRef]
  7. Hallett, S.R.; Jiang, W.-G.; Khan, B.; Wisnom, M.R. Modelling the interaction between matrix cracks and delamination damage in scaled quasi-isotropic specimens. Compos. Sci. Technol. 2008, 68, 80–89. [Google Scholar] [CrossRef]
  8. Koloor, S.S.R.; Tamin, M. Mode-II interlaminar fracture and crack-jump phenomenon in CFRP composite laminate materials. Compos. Struct. 2018, 204, 594–606. [Google Scholar] [CrossRef]
  9. Lopes, B.; Arruda, M.R.T.; Almeida-Fernandes, L.; Castro, L.; Silvestre, N.; Correia, J.R. Assessment of mesh dependency in the numerical simulation of compact tension tests for orthotropic materials. Compos. Part C Open Access 2020, 1, 100006. [Google Scholar] [CrossRef]
  10. Ng, T.P.; Koloor, R.S.S.; Djuansjah, J.R.P.; Abdul, K.M.R. Assessment of compressive failure process of cortical bone materials using damage-based model. J. Mech. Behav. Biomed. Mater. 2017, 66, 1–11. [Google Scholar] [CrossRef]
  11. Karimzadeh, A.; Koloor, S.S.R.; Ayatollahi, M.R.; Bushroa, A.R.; Yahya, M.Y. Assessment of nano-indentation method in mechanical characterization of heterogeneous nanocomposite materials using experimental and computational approaches. Sci. Rep. 2019, 9, 1–14. [Google Scholar] [CrossRef] [Green Version]
  12. Rahimian, K.S.S.; Karimzadeh, A.; Tamin, M.N.; Abd, S.M.H. Effects of sample and indenter configurations of nanoindentation experiment on the mechanical behavior and properties of ductile materials. Metals 2018, 8, 421. [Google Scholar] [CrossRef] [Green Version]
  13. Koloor, S.S.R.; Rahimian-Koloor, S.M.; Karimzadeh, A.; Hamdi, M.; Petrů, M.; Tamin, M. Nano-level damage characterization of graphene/polymer cohesive interface under tensile separation. Polymers 2019, 11, 1435. [Google Scholar] [CrossRef] [Green Version]
  14. Bassiri, N.A.; Farokhi, N.A.; Xin, L.; Ayob, A.; Yahya, M.Y.; R. Koloor, S.S.; Petrů, M.; Abu, H.S. Dynamic response of aluminium sheet 2024-T3 subjected to close-range shock wave: Experimental and numerical studies. J. Mater. Res. Technol. 2021, 10, 349–362. [Google Scholar] [CrossRef]
  15. Khan, M.S.; Abdul-Latif, A.; Koloor, S.S.R.; Petrů, M.; Tamin, M.N. Representative Cell Analysis for Damage-Based Failure Model of Polymer Hexagonal Honeycomb Structure under the Out-of-Plane Loadings. Polymers 2020, 13, 52. [Google Scholar] [CrossRef] [PubMed]
  16. Congro, M.; Roehl, D.; Mejia, C. Mesoscale computational modeling of the mechanical behavior of cement composite materials. Compos. Struct. 2021, 257, 113137. [Google Scholar] [CrossRef]
  17. Koloor, S.S.R.; Ayatollahi, M.R.; Tamin, M.N. Elastic-damage deformation response of fiber-reinforced polymer composite laminates with lamina interfaces. J. Reinf. Plast. Compos. 2017, 36, 832–849. [Google Scholar] [CrossRef]
  18. Kolanu, N.R.; Raju, G.; Ramji, M. A unified numerical approach for the simulation of intra and inter laminar damage evolution in stiffened CFRP panels under compression. Compos. Part B Eng. 2020, 190, 107931. [Google Scholar] [CrossRef]
  19. Lapczyk, I.; Hurtado, J.A. Progressive damage modeling in fiber-reinforced materials. Compos. Part A Appl. Sci. Manuf. 2007, 38, 2333–2341. [Google Scholar] [CrossRef]
  20. Eyer, G.; Montagnier, O.; Hochard, C.; Charles, J.P. Effect of matrix damage on compressive strength in the fiber direction for laminated composites. Compos. Part A Appl. Sci. Manuf. 2017, 94, 86–92. [Google Scholar] [CrossRef]
  21. Higuchi, R.; Okabe, T.; Nagashima, T. Numerical simulation of progressive damage and failure in composite laminates using XFEM/CZM coupled approach. Compos. Part A Appl. Sci. Manuf. 2017, 95, 197–207. [Google Scholar] [CrossRef]
  22. Chang, F.-K.; Lessard, L.B. Damage tolerance of laminated composites containing an open hole and subjected to compressive loadings: Part I—Analysis. J. Compos. Mater. 1991, 25, 2–43. [Google Scholar] [CrossRef]
  23. Greszczuk, L. Microbuckling of lamina-reinforced composites. In Proceeding of the Composite Materials: Testing and Design (Third Conference), Williamsburg, VA, USA, 21–22 March 1973; ASTM International: West Conshohocken, PA, USA, 1974. [Google Scholar]
  24. Puck, A.; Schürmann, H. Failure analysis of FRP laminates by means of physically based phenomenological models. Compos. Sci. Technol. 2002, 62, 1633–1662. [Google Scholar] [CrossRef]
  25. Lee, J.D. Three dimensional finite element analysis of damage accumulation in composite laminate. Comput. Struct. 1982, 15, 335–350. [Google Scholar] [CrossRef]
  26. Hashin, Z.; Rotem, A. A fatigue failure criterion for fiber reinforced materials. J. Compos. Mater. 1973, 7, 448–464. [Google Scholar] [CrossRef] [Green Version]
  27. Shahid, I.; Chang, F.-K. An accumulative damage model for tensile and shear failures of laminated composite plates. J. Compos. Mater. 1995, 29, 926–981. [Google Scholar] [CrossRef]
  28. Maimí, P.; Camanho, P.P.; Mayugo, J.; Dávila, C. A continuum damage model for composite laminates: Part I–Constitutive model. Mech. Mater. 2007, 39, 897–908. [Google Scholar] [CrossRef]
  29. Murakami, S. Continuum Damage Mechanics: A Continuum Mechanics Approach to the Analysis of Damage and Fracture; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 185. [Google Scholar]
  30. Kachanov, L. Introduction to Continuum Damage Mechanics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; Volume 10. [Google Scholar]
  31. Xue, J.; Kirane, K. Effect of the Shape of the Softening Damage Law on the Predicted Fracturing Behavior of Composites. In Proceedings of the American Society for Composites—Thirty-Fifth Technical Conference, Hoboken, NJ, USA, 14–17 September 2020. [Google Scholar]
  32. Kaw, A.K. Mechanics of Composite Materials; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  33. ASTM. D4762–08, Standard Guide for Testing Polymer Matrix Composite Materials. 2008. Available online: https://www.astm.org/Standards/D4762.htm (accessed on 10 December 2020).
  34. Koloor, R.S.S. Simulation Methodology for Fracture Processes of Composite Laminates Using Damage-Based Models. Ph.D. Degree, Universiti Teknologi Malaysia, Johor, Malaysia, 2016. [Google Scholar]
  35. Hashin, Z. Failure criteria for unidirectional fiber composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  36. Orifici, A.; Herszberg, I.; Thomson, R. Review of methodologies for composite material modelling incorporating failure. Compos. Struct. 2008, 86, 194–210. [Google Scholar] [CrossRef]
  37. ABAQUS-6.9EF. Theory Manual. 2010. Available online: http://130.149.89.49:2080/v6.9ef/index.html (accessed on 21 January 2021).
  38. Hinton, M.; Kaddour, A. The background to the Second World-Wide Failure Exercise. J. Compos. Mater. 2012, 46, 2283–2294. [Google Scholar] [CrossRef]
  39. Chaboche, J. The Concept of Effective Stress Applied Elasticity and to Viscoplasticity in Presence of Anisotropic Damage; Royal Aircraft Establishment: Farnborough, UK, 1979; pp. 737–760. [Google Scholar]
  40. Pinho, S.T. Modelling Failure of Laminated Composites Using Physically-Based Failure Models. Ph.D. Thesis, Imperial College, London, UK, 2005. [Google Scholar]
  41. Camanho, P.P.; Dávila, C.G. Mixed-mode decohesion finite elements for the simulation of delamination in composite materials. NASA Tech. Pap. 2002, 211737, 33. [Google Scholar]
  42. Matzenmiller, A.; Lubliner, J.; Taylor, R. A constitutive model for anisotropic damage in fiber-composites. Mech. Mater. 1995, 20, 125–152. [Google Scholar] [CrossRef]
  43. Kaddour, A.; Hinton, M.; Smith, P.; Li, S. A comparison between the predictive capability of matrix cracking, damage and failure criteria for fibre reinforced composite laminates: Part A of the third world-wide failure exercise. J. Compos. Mater. 2013, 47, 2749–2779. [Google Scholar] [CrossRef]
  44. Koloor, R.S.S.; Tamin, M. Effects of lamina damages on flexural stiffness of CFRP composites. In Proceedings of the 8th Asian-Australasian Conference on Composite Materials 2012, ACCM 2012—Composites: Enabling Tomorrow’s Industry Today, Kuala Lumpur, Malaysia, 6–8 November 2012; Asian Australasian Association for Composite Materials (AACM): Red Hook, NY, USA, 2013; Volume 1, pp. 237–243. [Google Scholar]
  45. Koloor, R.S.S.; Karimzadeh, A.; Yidris, N.; Petrů, M.; Ayatollahi, M.R.; Tamin, M.N. An energy-based concept for yielding of multidirectional FRP composite structures using a mesoscale lamina damage model. Polymers 2020, 12, 157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Varvani-Farahani, A.; Haftchenari, H.; Panbechi, M. An energy-based fatigue damage parameter for off-axis unidirectional FRP composites. Compos. Struct. 2007, 79, 381–389. [Google Scholar] [CrossRef]
  47. Wang, Y.; Wong, P.; Kodur, V. An experimental study of the mechanical properties of fibre reinforced polymer (FRP) and steel reinforcing bars at elevated temperatures. Compos. Struct. 2007, 80, 131–140. [Google Scholar] [CrossRef]
  48. Vermeulen, B.; Van Tooren, M. Design case study for a comparative performance analysis of aerospace materials. Mater. Des. 2006, 27, 10–20. [Google Scholar] [CrossRef]
  49. Camanho, P.P.; Fink, A.; Obst, A.; Pimenta, S. Hybrid titanium–CFRP laminates for high-performance bolted joints. Compos. Part A Appl. Sci. Manuf. 2009, 40, 1826–1837. [Google Scholar] [CrossRef] [Green Version]
  50. Oya, N.; Hamada, H. Mechanical properties and failure mechanisms of carbon fibre reinforced thermoplastic laminates. Compos. Part A Appl. Sci. Manuf. 1997, 28, 823–832. [Google Scholar] [CrossRef]
  51. Kumar, S.B.; Sridhar, I.; Sivashanker, S.; Osiyemi, S.; Bag, A. Tensile failure of adhesively bonded CFRP composite scarf joints. Mater. Sci. Eng. B 2006, 132, 113–120. [Google Scholar] [CrossRef]
  52. Wang, S.-X.; Wu, L.-Z.; Ma, L. Low-velocity impact and residual tensile strength analysis to carbon fiber composite laminates. Mater. Des. 2010, 31, 118–125. [Google Scholar] [CrossRef]
  53. Fan, J.; Guan, Z.; Cantwell, W. Numerical modelling of perforation failure in fibre metal laminates subjected to low velocity impact loading. Compos. Struct. 2011, 93, 2430–2436. [Google Scholar] [CrossRef]
  54. Koloor, R.S.S.; Abdul-Latif, A.; Tamin, M.N. Mechanics of Composite Delamination under Flexural Loading. Key Eng. Mater. 2011, 462, 726–731. [Google Scholar] [CrossRef]
  55. Koloor, R.S.S.; Redzuan, N.; Tamin, M.N. Shear-dominated Interlaminar Fracture Process in CFRP Composite Laminates. In Proceedings of the 9th International Conference on Fracture & Strength of Solids, Jeju Island, Korea, 9–13 June 2013. [Google Scholar]
Figure 1. Bilinear damage model of FRP composites (a) as described by the theory, and (b) the variation of this model based on different mechanical and damage properties reported in the previous studies.
Figure 1. Bilinear damage model of FRP composites (a) as described by the theory, and (b) the variation of this model based on different mechanical and damage properties reported in the previous studies.
Polymers 13 00344 g001
Figure 2. Bilinear damage model of FRP composites based on different (a) elastic, (b) strength, (c) fracture energy, and (d) strength and energy properties.
Figure 2. Bilinear damage model of FRP composites based on different (a) elastic, (b) strength, (c) fracture energy, and (d) strength and energy properties.
Polymers 13 00344 g002aPolymers 13 00344 g002b
Figure 3. Specimens configuration at different levels of material level (a), structural level (b), and super-structural level (c) as an example of hat-stiffened composite section in the Fuselage section of Boeing 787.
Figure 3. Specimens configuration at different levels of material level (a), structural level (b), and super-structural level (c) as an example of hat-stiffened composite section in the Fuselage section of Boeing 787.
Polymers 13 00344 g003
Figure 4. FE model of multi-layer construction of laminated composite.
Figure 4. FE model of multi-layer construction of laminated composite.
Polymers 13 00344 g004
Figure 5. Mesh configuration of the 3D geometry and boundary condition of the sample.
Figure 5. Mesh configuration of the 3D geometry and boundary condition of the sample.
Polymers 13 00344 g005
Figure 6. Stiffness response of CFRP composites at (a) material, (b,c) structure, and (d) super-structure levels, to the variation of the elastic properties.
Figure 6. Stiffness response of CFRP composites at (a) material, (b,c) structure, and (d) super-structure levels, to the variation of the elastic properties.
Polymers 13 00344 g006
Figure 7. Stiffness response of CFRP composites beam (material level) under tensile load. The load–displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Figure 7. Stiffness response of CFRP composites beam (material level) under tensile load. The load–displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Polymers 13 00344 g007aPolymers 13 00344 g007b
Figure 8. Stiffness response of thin CFRP composites plate (structural level) under three-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Figure 8. Stiffness response of thin CFRP composites plate (structural level) under three-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Polymers 13 00344 g008
Figure 9. Stiffness response of thick CFRP composites plate (structural level) under three-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Figure 9. Stiffness response of thick CFRP composites plate (structural level) under three-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Polymers 13 00344 g009
Figure 10. Stiffness response of hat-stiffened structure made of CFRP composites (super-structural level) under four-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Figure 10. Stiffness response of hat-stiffened structure made of CFRP composites (super-structural level) under four-point bending load. The load-displacement response based on (a) strength, (b) energy, and (c) both strength and energy properties variations.
Polymers 13 00344 g010
Table 1. Elastic and damage properties of UD CFRP composite lamina.
Table 1. Elastic and damage properties of UD CFRP composite lamina.
Elastic PropertiesDamage Model Parameters
Sets Polymers 13 00344 i0011/2RefRef3/2Ref 1/2RefRef3/2Ref
E11, GPa52.75105.5158.25Longitudinal tensile strength, MPaXT67013402010
E22, GPa3.67.210.8Longitudinal compressive strength, MPaXC59611921788
E33, GPa3.67.210.8Transverse tensile strength, MPaYT9.819.629.4
G12, GPa1.73.45.1Transverse compressive strength, MPaYC46.1592.3138.45
G13, GPa1.73.45.1Longitudinal shear strength, MPaSL25.55176.5
G23, GPa1.262.523.78Transverse shear strength, MPaST11.52334.5
ν120.340.340.34Longitudinal tensile fracture energy, N/mmGXT24.248.472.6
ν130.340.340.34Longitudinal compressive fracture energy, N/mmGXC30.1560.390.45
ν230.3780.3780.378Transverse tensile fracture energy, N/mmGYT2.254.56.75
Transverse compressive fracture energy, N/mmGYC4.258.512.75
Table 2. List of FE models of CFRP composite materials and structures.
Table 2. List of FE models of CFRP composite materials and structures.
No.LevelsComposite CaseSpecificationLoadDimension
1MaterialBeam[45]5Tensile200 × 20 × 1 mm3
2StructureThinPlate 1[60/45/90/−45/30]Three-point bending150 × 60 × 1 mm3
ThickPlate 2[453/903/03/−453/453]Three-point bending150 × 60 × 3 mm3
3Super-structureProfile Hat structure: [45/−45/90/45/90]
Plate structure: [45/−45/45/−45/45]
Four-point bendingRefer to Figure 3
Table 3. Different sets of elastic-damage material properties.
Table 3. Different sets of elastic-damage material properties.
No.PropertiesCasesProperty CodeFE ModelCase Code
1Elastic Properties
(Figure 2a)
1Ep-1 BeamBe-Ep-1
Plate1P1-Ep-1
Plate2P2-Ep-1
ProfilePr-Ep-1
RefEp-RefBeamBe-Ep-Ref *
Plate1P1- Ep-Ref **
Plate2P2- Ep-Ref ***
ProfilePr-Ep-Ref ****
2Ep-2BeamBe-Ep-2
Plate1P1-Ep-2
Plate2P2-Ep-2
ProfilePr-Ep-2
2Strength values
(Figure 2b)
1S-1BeamBe-S-1
Plate1P1-S-1
Plate2P2-S-1
ProfilePr-S-1
RefS-RefBeamBe-S-Ref *
Plate1P1-S-Ref **
Plate2P2-S-Ref ***
ProfilePr-S-Ref ****
2S-2BeamBe-S-2
Plate1P1-S-2
Plate2P2-S-2
ProfilePr-S-2
3Energy values
(Figure 2c)
1E-1BeamBe-E-1
Plate1P1-E-1
Plate2P2-E-1
ProfilePr-E-1
RefE-RefBeamBe-E-Ref *
Plate1P1-E-Ref **
Plate2P2-E-Ref ***
ProfilePr-E-Ref ****
2E-2BeamBe-E-2
Plate1P1-E-2
Plate2P2-E-2
ProfilePr-E-2
4Strength and Energy values
(Figure 2d)
1SE-1BeamBe-SE-1
Plate1P1-SE-1
Plate2P2-SE-1
ProfilePr-SE-1
RefSE-RefBeamBe-SE-Ref *
Plate1P1-SE-Ref **
Plate2P2-SE-Ref ***
ProfilePr-SE-Ref ****
2SE-2BeamBe-SE-2
Plate1P1-SE-2
Plate2P2-SE-2
ProfilePr-SE-2
Note; Each set of data presented by *, **, ***, and **** are same. Detail info of the FE models and geometry are provided in Table 2, and Figure 3.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

R. Koloor, S.S.; Karimzadeh, A.; Abdullah, M.R.; Petrů, M.; Yidris, N.; Sapuan, S.M.; Tamin, M.N. Linear-Nonlinear Stiffness Responses of Carbon Fiber-Reinforced Polymer Composite Materials and Structures: A Numerical Study. Polymers 2021, 13, 344. https://doi.org/10.3390/polym13030344

AMA Style

R. Koloor SS, Karimzadeh A, Abdullah MR, Petrů M, Yidris N, Sapuan SM, Tamin MN. Linear-Nonlinear Stiffness Responses of Carbon Fiber-Reinforced Polymer Composite Materials and Structures: A Numerical Study. Polymers. 2021; 13(3):344. https://doi.org/10.3390/polym13030344

Chicago/Turabian Style

R. Koloor, S. S., A. Karimzadeh, M. R. Abdullah, M. Petrů, N. Yidris, S. M. Sapuan, and M. N. Tamin. 2021. "Linear-Nonlinear Stiffness Responses of Carbon Fiber-Reinforced Polymer Composite Materials and Structures: A Numerical Study" Polymers 13, no. 3: 344. https://doi.org/10.3390/polym13030344

APA Style

R. Koloor, S. S., Karimzadeh, A., Abdullah, M. R., Petrů, M., Yidris, N., Sapuan, S. M., & Tamin, M. N. (2021). Linear-Nonlinear Stiffness Responses of Carbon Fiber-Reinforced Polymer Composite Materials and Structures: A Numerical Study. Polymers, 13(3), 344. https://doi.org/10.3390/polym13030344

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