Modelling the nonlinear behaviour and fracture process of AS4/PEKK thermoplastic composite under shear loading

The accurate determination of non-linear shear behaviour and fracture toughness of continuous carbon-fibre/polymer composites remains a considerable challenge. These measurements are often necessary to generate material parameters for advanced computational damage models. In particular, there is a dearth of detailed shear fracture toughness characterisation for thermoplastic composites which are increasingly generating renewed interest within the aerospace and automotive sectors. In this work, carbon fibre (AS4)/ thermoplastic Polyetherketoneketone (PEKK) composite V-notched cross-ply specimens were manufactured to investigate their non-linear response under pure shear loading. Both monotonic and cyclic loading were applied to study the shear modulus degradation and progressive failure. For the first time in the reported literature, we use the essential work of fracture approach to measure the shear fracture toughness of continuous fibre reinforced composite laminates. Excellent geometric similarity in the load-displacement curves was observed for ligament-scaled specimens. The laminate fracture toughness was determined by linear regression, of the specific work of fracture values, to zero ligament thickness, and verified with computational models. The matrix intralaminar fracture toughness (ply level fracture toughness), associated with shear loading was determined by the area method. This paper also details the numerical implementation of a new three-dimensional phenomenological model for carbon fibre thermoplastic composites using the measured values, which is able to accurately represent the full non-linear mechanical response and fracture process. The constitutive model includes a new non-linear shear profile, shear modulus degradation and load reversal. It is combined with a smeared crack model for representing ply-level damage initiation and propagation. The model is shown to accurately predict the constitutive response in terms of permanent plastic strain, degraded modulus as well as load reversal. Predictions are also shown to compare favourably with the evolution of damage leading to final fracture.

Shear Strength Model for FRP-Strengthened RC Beams with Adverse FRP-Steel Interaction

RC beams shear strengthened with externally bonded fiber-reinforced polymer (FRP) U strips or side strips usually fail owing to debonding of the bonded FRP shear reinforcement. Because such debonding usually occurs in a brittle manner at relatively small shear crack widths, some of the internal steel stirrups intersected by the critical shear crack may not have reached yielding at beam shear failure. Consequently, the yield stress of internal steel stirrups in such a strengthened RC beam cannot be fully utilized. This adverse shear interaction between the internal steel shear reinforcement and the external FRP shear reinforcement may significantly reduce the benefit of the shear strengthening FRP but has not been considered explicitly by any of the shear strength models in the existing design guidelines. This paper presents a new shear strength model considering this adverse shear interaction through the introduction of a shear interaction factor. A comprehensive evaluation of the proposed model, as well as three other shear strength models, is conducted using a large test database. It is shown that the proposed shear strength model performs the best among the models compared, and the performance of the other shear strength models can be significantly improved by including the proposed shear interaction factor. Finally, a design recommendation is presented.

Modelling Static and Dynamic FRP-Concrete Bond Behaviour Using a Local Concrete Damage Model

This paper presents a study on the bond behaviour of FRP-concrete bonded joints under static and dynamic loadings, by developing a meso-scale finite element model using the K&C concrete damage model in LS-DYNA. A significant number of single shear experiments under static pull-off loading were modelled with an extensive parametric study covering key factors in the K&C model, including the crack band width, the compressive fracture energy and the shear dilatation factor. It is demonstrated that the developed model can satisfactorily simulate the static debonding behaviour, in terms of mesh objectivity, the load-carrying capacity and the local bond-slip behaviour, provided that proper consideration is given to the selection of crack band width and shear dilatation factor. A preliminary study of the effect of the dynamic loading rate on the debonding behaviour was also conducted by considering a dynamic increase factor (DIF) for the concrete strength as a function of strain rate. It is shown that a higher loading rate leads to a higher load-carrying capacity, a longer effective bond length, and a larger damaged area of concrete in the single shear loading scenario.

Improved shear strength model for FRP-strengthened RC beams accounting for adverse FRP-steel interaction

RC beams shear-strengthened with externally-bonded FRP side strips or U-strips usually fail by debonding. As such debonding occurs in a brittle manner at relatively small shear crack widths, some of the internal steel stirrups may not have reached yielding at beam shear failure. Consequently, the internal steel stirrups cannot be fully utilized. This adverse shear interaction between internal steel stirrups and external FRP strips may significantly reduce the benefit of shear-strengthening FRP but has not been considered by any of the existing FRP strengthening design guidelines. In this paper, an improved shear strength model capable of accounting for the effect of the above shear interaction is first presented, in which the unfavorable effect of shear interaction is reflected through a reduction factor (i.e. shear interaction factor). Using a large test database established in the present study, the performance of the proposed model as well as that of three other shear strength models is then assessed. This assessment shows that the proposed shear strength model performs better than the three existing models. The assessment also shows that the inclusion of the proposed shear interaction factor in the existing models can significantly improve their performance.

Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.