Composite structures are made of multidirectional (MD) fiber-reinforced polymer (FRP) composite laminates, which fail due to multiple damages in matrix, interface, and fiber constituents at different scales. The yield point of a unidirectional FRP composite is assumed as the lamina strength limit representing the damage initiation phenomena, while yielding of MD composites in structural applications are not quantified due to the complexity of the sequence of damage evolutions in different laminas dependent on their angle and specification. This paper proposes a new method to identify the yield point of MD composite structures based on the evolution of the damage dissipation energy (DDE). Such a characteristic evolution curve is computed using a validated finite element model with a mesoscale damage-based constitutive model that accounts for different matrix and fiber failure modes in angle lamina. The yield point of composite structures is identified to correspond to a 5% increase in the initial slope of the DDE evolution curve. The yield points of three antisymmetric MD FRP composite structures under flexural loading conditions are established based on Hashin unidirectional (UD) criteria and the energy-based criterion. It is shown that the new energy concept provides a significantly larger safe limit of yield for MD composite structures compared to UD criteria, in which the accumulation of energy dissipated due to all damage modes is less than 5% of the fracture energy required for the structural rupture.

It is necessary to consider the influence of moisture damage on the interlaminar fracture toughness for composite structures that are used for outdoor applications. However, the studies on the progressive variation of the fracture toughness as a function of moisture content M (%) is rather limited. In this regard, this study focuses on the characterization of mode II delamination of carbon/epoxy composites conditioned at 70 °C/85% relative humidity (RH). End-notched flexure test is conducted for specimens aged at various moisture absorption levels. Experimental results reveal that mode II fracture toughness degrades with the moisture content, with a maximum of 23% decrement. A residual property model is used to predict the variation of the fracture toughness with the moisture content. Through numerical simulations, it is found that the approaches used to estimate the lamina and cohesive properties are suitable to obtain reliable simulation results. In addition, the damage initiation is noticed during the early loading stage; however, the complete damage is only observed when the numerical peak load is achieved. Results from the present research could serve as guidelines to predict the residual properties and simulate the mode II delamination behavior under moisture attack.

The honeycomb (HC) core of sandwich structures undergoes flexural loading and carries the normal compression and shear. The mechanical properties and deformation response of the core need to be established for the design requirements. In this respect, this article describes the development of the smallest possible representative cell (RC) models for quantifying the deformation and failure process of the Nomex polymer-based hexagonal HC core structure under the out-of-plane quasi-static loadings. While the hexagonal single and multi-cell models are suitable for the tension and compression, a six-cell model is the simplest RC model developed for shear in the transverse and ribbon direction. Hashin’s matrix and fiber damage equations are employed in simulating the failure process of the orthotropic cell walls, using the finite element (FE) analysis. The FE-calculated load–displacement curves are validated with the comparable measured responses throughout the loading to failure. The location of the fracture plane of the critical cell wall in the out-of-plane tension case is well predicted. The wrinkling of the cell walls, leading to the structural buckling of the HC core specimen in the compression test, compares well with the observed failure mechanisms. In addition, the observed localized buckling of the cell wall by the induced compressive stress during the out-of-plane shear in both the transverse and ribbon direction is explained. The mesoscale RC models of the polymer hexagonal HC core structure have adequately demonstrated the ability to predict the mechanics of deformation and the mechanisms of failure.

The mechanical behavior of graphene/polymer interfaces in the graphene-reinforced epoxy nanocomposite is one of the factors that dictates the deformation and damage response of the nanocomposites. In this study, hybrid molecular dynamic (MD) and finite element (FE) simulations of a graphene/polymer nanocomposite are developed to characterize the elastic-damage behavior of graphene/polymer interfaces under a tensile separation condition. The MD results show that the graphene/epoxy interface behaves in the form of elastic-softening exponential regressive law. The FE results verify the adequacy of the cohesive zone model in accurate prediction of the interface damage behavior. The graphene/epoxy cohesive interface is characterized by normal stiffness, tensile strength, and fracture energy of 5 × 10⁻⁸ (aPa·nm⁻¹), 9.75 × 10⁻¹⁰ (nm), 2.1 × 10⁻¹⁰ (N·nm⁻¹) respectively, that is followed by an exponential regressive law with the exponent, α = 7.74. It is shown that the commonly assumed bilinear softening law of the cohesive interface could lead up to 55% error in the predicted separation of the interface.

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.

The nanoindentation test is frequently used as an alternate method to obtain the mechanical properties of ductile materials. However, due to the lack of information about the effects of the sample and indenter physical configurations, the accuracy of the extracted material properties in nanoindentation tests requires further evaluation that has been considered in this study. In this respect, a demonstrator ductile material, aluminum 1100, was tested using the Triboscope nanoindenter system with the Berkovich indenter. A 3D finite element simulation of the nanoindentation test was developed and validated through exact prediction of the structural response with measured data. The validated model was then employed to examine the effects of various test configurations on the load–displacement response of the sample material. These parameters were the different indenter edge-tip radii, different indentation depths, different sample tilts, and different friction conditions between the indenter and the material surface. Within the range of the indenter edge-tip radii examined, the average elastic modulus and hardness were 78.34 ± 14.58 and 1.6 ± 0.24 GPa, respectively. The different indentation depths resulted in average values for the elastic modulus and hardness of 77.03 ± 6.54 and 1.58 ± 0.17 GPa, respectively. The uneven surface morphology, as described by the inclination of the local indentation plane, indicated an exponential increase in the extracted values of elastic modulus and hardness, ranging from 71.83 and 1.47 GPa (for the reference case, θ = 0°) to 243.39 and 5.05 GPa at θ = 12°. The mechanical properties that were obtained through nanoindentation on the surface with 6° tilt or higher were outside the range for aluminum properties. The effect of friction on the resulting mechanical response and the properties of the material was negligible.