Discrete element modelling of unidirectional fibre-reinforced polymers under transverse tension

The mechanical behaviour of unidirectional fibre-reinforced polymer composites subjected to transverse tension was studied using a two dimensional discrete element method. The Representative Volume Element (RVE) of the composite was idealised as a polymer matrix reinforced with randomly distributed parallel fibres. The matrix and fibres were constructed using disc particles bonded together using parallel bonds, while the fibre/matrix interfaces were represented by a displacement-softening model. The prevailing damage mechanisms observed from the model were interfacial debonding and matrix plastic deformation. Numerical simulations have shown that the magnitude of stress is significantly higher at the interfaces, especially in the areas with high fibre densities. Interface fracture energy, stiffness and strength all played important roles in the overall mechanical performance of the composite. It was also observed that tension cracks normally began with interfacial debonding. The merge of the interfacial and matrix micro-cracks resulted in the final catastrophic fracture.     

Incremental damage theory of particulate-reinforced composites with a ductile interphase

This paper deals with a new micromechanics model of particulate-reinforced composites (PRCs) describing the evolution of debonding damage, matrix plasticity and particle size effect on the deformation. A ductile interphase was considered in the frame of incremental damage theory to analyze the dependence of elastic–plastic–damage behavior on particle size. Progressive debonding damage was controlled by a critical energy criterion for particle–matrix interfacial separation. The equivalent stresses of the matrix and interphase were determined by field fluctuation method. The influences of progressive debonding damage, particle size and interphase properties on the overall stress–strain response of PRC were explained simultaneously. Due to the existence of a ductile interphase, stress transfer and plastic initiation in PRC become very complicated, and thus a unit-cell (UC) based FEM was used to simulate their evolutions and demonstrate the role of the interphase. Finally, particle size effect on the mechanical behaviors of composites was interpreted.     

Mechanical properties of epoxy composites filled with silane-functionalized graphene oxide

In this work, the effects of as-produced GO and silane functionalized GO (silane-f-GO) loading and silane functionalization on the mechanical properties of epoxy composites are investigated and compared. Such silane functionalization containing epoxy ended-groups is found to effectively improve the compatibility between the silane-f-GO and the epoxy matrix. Increased storage modulus, glass transition temperature, thermal stability, tensile and flexural properties and fracture toughness of epoxy composites filled with the silane-f-GO sheets are observed compared with those of the neat epoxy and GO/epoxy composites. These findings confirm the improved dispersion and interfacial interaction in the composites arising from covalent bonds between the silane-f-GO and the epoxy matrix. Moreover, several possible fracture mechanisms, i.e. crack pinning/deflection, crack bridging, and matrix plastic deformation initiated by the debonding/delamination of GO sheets, were identified and evaluated.     

连续石墨纤维增强铝基复合材料准静态拉伸损伤演化与断裂力学行为

针对连续石墨纤维增强铝基(CF/Al)复合材料,采用三种纤维排布方式的代表体积单元(RVE)建立了其细观力学有限元模型,采用准静态拉伸试验与数值模拟结合的方法,研究了其在轴向拉伸载荷下的渐进损伤与断裂力学行为。结果表明,采用基体合金和纤维原位力学性能建立的细观力学有限元模型,对轴向拉伸弹性模量和极限强度的计算结果与实验结果吻合良好,而断裂应变计算值较实验结果偏低。轴向拉伸变形中首先出现界面和基体合金损伤现象,随应变增加界面发生失效并诱发基体合金的局部失效,最后复合材料因纤维发生失效而破坏,从而出现界面脱粘后纤维拔出与基体合金撕裂共存的微观形貌。细观力学有限元分析结果表明,在复合材料制备后纤维性能衰减而强度较低条件下,改变界面强度和刚度对复合材料轴向拉伸弹塑性力学行为的影响较小,复合材料中纤维强度水平是决定该复合材料轴向拉伸力学性能的主要因素。     

Mechanical properties and toughening mechanisms of polypropylene/barium sulfate composites

The effect of interfacial interaction on the mechanical performance of a group of polypropylene (PP)/barium sulfate (BaSO₄) composites were studied. It was found that PP can be toughened with specially treated BaSO₄ particles. The interfacial modification contributes to the toughening in two aspects. The first is to provide a proper interfacial adhesion and control the interfacial debonding occurs at well-timed stages. This ensures the inorganic particles transfer the stress and stabilizes the cracks at the initial stage of the deformation, and satisfy the stress conditions for plastic deformation of matrix ligaments subsequently via debonding. The second is that the modified interface between PP matrix and filler particles increases the nucleating ability of the fillers and retards the motion of the PP chains. This leads to the formation of PP crystals with less perfection and smaller size in the matrix and promotes plastic deformation of the matrix after the debonding occurs.     

Fatigue damage behaviors of carbon fiber-reinforced epoxy composites containing nanoclay

The effects of nanoclay inclusion on cyclic fatigue behavior and residual properties of carbon fiber-reinforced composites (CFRPs) after fatigue have been studied. The tension–tension cyclic fatigue tests are conducted at various load levels to establish the S-N curve. The residual strength and modulus are measured at different stages of fatigue cycles. The scanning electron microscopy (SEM) and scanning acoustic microscopy (SAM) are employed to characterize the underlying fatigue damage mechanisms and progressive damage growth. The incorporation of nanoclay into CFRP composites not only improves the mechanical properties of the composite in static loading, but also the fatigue life for a given cyclic load level and the residual mechanical properties after a given period of cyclic fatigue. The corresponding fatigue damage area is significantly reduced due to nanoclay. Nanoclay serves to suppress and delay delamination damage growth and eventual failure by improving the fiber/matrix interfacial bond and through the formation of nanoclay-induced dimples.     

Fracture mechanisms of epoxy-based ternary composites filled with rigid-soft particles

Epoxy composites filled with different amounts of aggregate-free silica nanoparticles and phase-separated submicron rubber particles were fabricated to study the synergistic effect of multi-phase particles on mechanical properties of the composites. Compared with binary composites with single-phase particles, the ternary composites with both rigid and soft particles offer a good balance in stiffness, strength and fracture toughness, showing capacities in tailoring the mechanical properties of modified epoxy resins. It was observed that debonding of silica nanoparticles from matrix in the ternary composites was less pronounced than that in the binary composites. Moreover, the rubber particles became smaller and their shape tends to be irregular, affected by the presence of rigid silica nanoparticles. The toughening mechanisms in the epoxy composites were evaluated, and the enlarged plastic deformation around the crack tip, induced by the combination of rigid and soft particles, seems to be a dominant factor in enhancing fracture toughness of the ternary composites.     

Flexural and compression response of woven E-glass/polyester–CNF nanophased composites

A significant improvement in fiber reinforced polymeric composite (FRPC) materials can be obtained by incorporating a very small amount of nanofiller in the matrix material. In this work, an ultrasonic liquid processor was used to infuse carbon nanofiber (CNF) into the polyester matrix which was then mixed with catalyst using a mechanical agitator. Both conventional and CNF-filled glass-fiber reinforced polyester composites (GRPC) were fabricated using the vacuum assisted resin transfer molding (VARTM) process. Excellent dispersion of CNFs into the polyester resin was observed from the scanning electron microscopy (SEM) micrographs. Flexural and quasi-static tests were performed for investigating the mechanical responses. Fracture surface was examined using optical microscopy (OM) and SEM. Flexure tests performed on the conventional GRPC, 0.1–0.4 wt.% CNF-filled GRPC showed up to 49% and 31% increase in the flexural strength and modulus, respectively, compared to the conventional one with increasing loading of CNFs up to 0.2 wt.%. Similar trend was seen in quasi-static compression properties. SEM evaluation revealed relatively less damage in the tested fracture surfaces of the nanophased composites in terms of matrix failure, fiber breakage, matrix–fiber debonding, and delamination, compared to the conventional one. This might be the result of better interfacial interaction between matrix and fibers, due to the presence of CNFs.     

含界面相三维全五向编织复合材料拉伸行为模拟及失效机制研究

基于三维全五向（Q5D）编织复合材料的细观结构模型,通过引入界面相单元,建立了含界面相Q5D编织复合材料单轴拉伸损伤失效分析模型。应用Python语言实现对ABAQUS的二次开发,将Linde等提出的失效准则和Von-Mises应力准则分别用于纱线和基体的渐进损伤判断,并确定材料的整体失效模式;对于界面相,采用Quads准则进行损伤判断。利用周期性位移边界条件,对含界面相Q5D编织复合材料的纵向拉伸应力-应变行为进行了渐进损伤数值模拟,详细讨论了在纵向拉伸载荷作用下材料的细观损伤起始、扩展和最终失效的演化过程,分析了材料的细观损伤失效机制,预测了材料的极限破坏强度,并研究了界面相性能对材料整体力学行为的影响规律。研究结果表明,数值模拟结果与实验值吻合较好,验证了渐进损伤模型的有效性,为该类材料的力学分析和优化设计奠定了基础。     

Measurement of resistance curves in the longitudinal failure of composites using digital image correlation

This paper presents a new methodology to measure the crack resistance curves associated with fiber-dominated failure modes in polymer–matrix composites. The crack resistance curves not only characterize the fracture toughness of the material, but are also the basis for the identification of the parameters of the softening laws used in the numerical simulation of fracture in composite materials. The proposed method is based on the identification of the crack tip location using Digital Image Correlation and the calculation of the J-integral directly from the test data using a simple expression derived for cross-ply composite laminates. It is shown that the results obtained using the proposed methodology yield crack resistance curves similar to those obtained using Finite Element based methods for compact tension carbon–epoxy specimens. However, it is also shown that, while the Digital Image Correlation based technique mitigates the problems resulting from Finite Element based data reduction schemes applied to compact compression tests, the delamination that accompanies the propagation of a kink-band renders compact compression test specimens unsuitable to measure resistance curves associated with fiber kinking.     

Microstructure effects on transverse cracking in composite laminae by DEM

A particle discrete element method (DEM) was employed to simulate transverse cracking in laminated fiber reinforced composites. The microstructure of the laminates was modeled by a DEM model using different mechanical constitutive laws and materials parameters for different constituents, i.e. fiber, matrix and fiber/matrix interface. Rectangular, hexagonal and random fiber distributions were simulated to study the effect of fiber distribution on the transverse cracking. The initiation and dynamic propagation of transverse cracking and interfacial debonding were all captured by the DEM simulation, which showed similar patterns to those observed from experiments. The effect of fiber volume fraction was also studied for laminae with randomly distributed fibers. It was found that the distribution and volume fraction of fibers affected not only the transverse cracking path, but also the behavior of matrix plastic deformation and fiber/matrix interface yielding in the material.     

Failure analysis and the optimal toughness design of carbon nanotube-reinforced composites

The combined analysis of the fracture toughness enhancement of carbon nanotube (CNT)-reinforced composites is herein carried out on the basis of atomistic simulation, shear-lag theory and facture mechanics. It is found that neither longer reinforced CNTs nor stronger CNT/matrix interfaces can definitely lead to the better fracture toughness of these composites. In contrast, the optimal interfacial chemical bond density and the optimal CNT length are those making the failure mode just in the transition from CNT pull-out to CNT break. To verify our theory, an atomic/continuum finite element method (FEM) is applied to investigate the fracture behavior of CNT-reinforced composites with different interfacial chemical bond densities. Our analysis shows that the optimal interfacial chemical bond density for (6,6) CNTs is about 5–10% and that increasing the CNT length beyond 100 nm does not further improve fracture toughness, but can easily lead to the self-folding and clustering of the CNTs. The proposed theoretical model is also applicable to short fiber-reinforced composites.     

2D and 3D imaging of fatigue failure mechanisms of 3D woven composites

A detailed investigation of the failure mechanisms for angle-interlocked (AI) and modified layer-to-layer (MLL) three dimensional (3D) woven composites under tension–tension (T–T) fatigue loading has been conducted using surface optical microscopy, cross-sectional SEM imaging, and non-destructive X-ray computed tomography (CT). X-ray microCT has revealed how cracks including surface matrix cracks, transverse matrix cracks, fibre/matrix interfacial debonding or delamination develop, and has delineated the complex 3D morphology of these cracks in relation to fibre architecture. For both weaves examined, transverse cracks soon become uniformly distributed in the weft yarns. A higher crack density was found in the AI composite than the MLL composite. Transverse cracking initiates in the fibre rich regions of weft yarns rather than the resin rich regions. Delaminations in the failed MLL specimen were more extensive than the AI specimen. It is suggested that for the MLL composite that debonding between the binder yarns and surrounding material is the predominant damage mechanism.     

Three-phase incremental damage theory of particulate-reinforced composites with a brittle interphase

This paper deals with a constitutive model of particulate-reinforced composites (PRCs) which can describe the evolution of debonding damage, matrix plasticity and particle size effects on deformation and damage. An incremental damage model of PRC based on Mori–Tanaka’s mean field concept has been extended to three-phase composites for interpreting particle size effect. The interphase was perfectly incorporated into the present micromechanics model as a third phase with the help of double-inclusion model. Progressive damage was controlled by a critical energy criterion for particle–matrix interfacial separation. Based on the developed model, influences of progressive debonding damage, particle size and interphase properties on the overall stress–strain response of PRC were discussed. Finally, particle size effect on the mechanical behaviors of composites was clearly interpreted from the role of the interphase, which is different from all the existing researches.     

A numerical approach for predicting the failure locus of fiber reinforced composites under combined transverse compression and axial tension

A computational model based on the finite element method is presented for the estimation of strength of a fiber-reinforced lamina subjected to a combination of the transverse compression and axial tension. A complex damage mechanism including fiber breakage, fiber/matrix debonding and matrix plastic deformation is reproduced in the proposed model by using appropriate constitutive equations. The numerical simulation of mechanical response of the unidirectional lamina under biaxial loading is used to obtained the failure locus. Subsequently, the model is verified against an analytical solution and experimental data. It was found that the numerical calculations agree better with experimental results than analytical predictions.     

Influence of thermal distortion to compression failure of polymer matrix composite panels in fire

This paper discusses the influence of thermal distortion to compression failure load of polymer matrix composite panels in fire. The thermal gradient along the thickness of the panel, caused by low thermal conductivity, induces transverse deformation field even without external mechanical loads, known as thermal distortion. The deformation field of thermal distortion obtained from combined thermal–mechanical analysis is presented in terms of non-dimensional parameters that represent thermal and mechanical properties, including those in phase transition at high temperature range, and dimensional lengths. A failure model is proposed to assess the transition from Euler buckling to plastic compression failure under the influence of thermal distortion. For large slenderness ratio, the panel is failed by Euler buckling. For small slenderness ratio, the panel is failed by plastic failure mechanism, where the analysis quantifies the influence of thermal distortion to the reduction of load-bearing capacity. The compression failure mechanism under the influence of thermal distortion is discussed for developing quantitative approach to design polymer matrix composite structures in fire.     

Dynamic compressive behavior of unidirectional IM7/5250-4 laminate after thermal oxidation

Fiber reinforced high temperature polymer matrix composites are currently gaining wide usage in aircraft structures, especially in airframe and engine inlet casing. The failure of composites in worst-case operational conditions mandates the extensive investigation of the mechanical behavior, and the durability in long-term performance and service life under thermal oxidation. In this work, unidirectional IM7 carbon fiber reinforced high-temperature BMI resin composite (IM7/5250-4) were isothermally aged in air for 2 months at 195 °C and 245 °C, respectively. The dynamic behavior of thermally aged composites was investigated on a split Hopkinson pressure bar (SHPB) in three principal directions. The results indicate that thermal oxidation leads to significant reduction in both stiffness and strength of the composites. Optical micrographs of fracture surface and failure pattern of composite after SHPB impact reveals oxidation induced debonding along the fiber–matrix interface due to oxygen diffusion under long-term exposure to elevated temperatures.     

Progressive damage modeling in fiber-reinforced materials

This paper presents an anisotropic damage model suitable for predicting failure and post-failure behavior in fiber-reinforced materials. In the model the plane stress formulation is used and the response of the undamaged material is assumed to be linearly elastic. The model is intended to predict behavior of elastic-brittle materials that show no significant plastic deformation before failure. Four different failure modes – fiber tension, fiber compression, matrix tension, and matrix compression – are considered and modeled separately. The onset of damage is predicted using Hashin’s initiation criteria [Hashin Z, Rotem A. A fatigue failure criterion for fiber-reinforced materials. J Compos Mater 1973;7:448; Hashin Z. Failure criteria for unidirectional fiber composites. J Appl Mech 1980;47:329–34] and the progression of damage is controlled by a new damage evolution law, which is easy to implement in a finite element code. The evolution law is based on fracture energy dissipation during the damage process and the increase in damage is controlled by equivalent displacements. The issues related to numerical implementation, such as mesh sensitivity and convergence in the softening regime, are also addressed.     

Split Hopkinson pressure bar testing of 3D multi-axial warp knitted carbon/epoxy composites

High-strain-rate compression experiments were performed on 3D MWK carbon/epoxy composites with different fiber architectures at room and elevated temperature using an SHPB apparatus. Macro-fracture and SEM micrographs were examined to understand the failure mechanism. The results show the dynamic properties increase with the strain rate and show a high-strain-rate sensitivity. Meanwhile, composites with [0°/0°/0°/0°] have higher properties. Moreover, composites show temperature sensitivity and the properties decrease significantly, especially for composites with [0°/90°/+45°/−45°]. The results also indicate composites take on more serious damage and failure with the strain rate. The failure of composites with [0°/0°/0°/0°] behaves in multiple delaminating, overall expansion and 0° fibers tearing. While that of composites with [0°/90°/+45°/−45°] is mainly interlaminar delaminating, local fibers tearing and fracture on different fiber layers. In addition, with increasing the temperature, the composite shows less fracture and becomes more plastic. The damage of matrix yielding, interface debonding and twisting of fibers increase significantly.