On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

A review of dynamic fracture mechanics in China

This paper summarizes the recent development of dynamic fracture in China. The review covers analytical and numerical results on elastodynamic crack fields in 3D and layered media; experimental and theoretical research on dynamic mechanical properties of rocks and advanced materials; transient effects on ideally plastic crack-tip fields when the inertia forces are not negligible.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

Composite materials dynamic fracture studies by generalized Shmuely difference algorithm

The generalized Shmuely Difference Algorithm (GSDA) is presented here to analyze the dynamic fracture performance of orthogonal-anisotropic composite materials, such as glass fibre reinforced phenolplast. The difference recurrence Formulae and boundary condition difference extrapolation formulae are derived and programmed. The dynamic stress intensity factors (DSIF) of the isotropic and anisotropic centrally cracked plates are computed respectively using GSDA and compared with that published previously. GSDA is proved effective and reliable. Copyright (C) 1996 Elsevier Science Ltd.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Fracture criterion based on the higher-order asymptotic fields

A complete development for the higher-order asymptotic solutions of the crack tip fields and finite element calculations for mode I loading of hardening materials in plane strain are performed. The results show that in the higher-order asymptotic solution (to the twentieth order), only three coefficients are independent. These coefficients are determined by matching with the finite element solutions carried out in the present paper (our attention is focused on the first five terms of the higher-order asymptotic solution). We obtain an analytic characterization of crack tip fields, which conform very well to the finite element solutions over wide range. A modified two parameter criterion based on the asymptotic solution of five terms is presented. The upper bound and lower bound fracture toughness curves predicted by modified two parameter criterion are given. These two curves agree with most of the experimental data and fully capture the proper trend.

An analysis of the effects of symmetrical tilt grain boundaries on the plastic deformation

A numerical investigation on the simple polycrystals containing three symmetrical tilt grain boundaries (GBs) is carried out within the framework of crystal plasticity which precisely considers the finite deformation and finite lattice rotation as well as elastic anisotropy. The calculated results show that the slip geometry and the redistribution of stresses arising from the anisotropy and boundary constraint play an important role in the plastic deformation in the simple polycrystals. The stress level along GB is sensitive to the load level and misorientation, and the stresses along QB are distributed nonuniformly. The GB may exhibit a softening or strengthening feature, which depends on the misorientation angle. The localized deformation bands usually develop accompanying the GB plastic deformation, the impingement of the localized band on the GB may result in another localized deformation band. The yield stresses with different misorientation angles are favorably compared with the experimental results.

The interlaminar fracture-behavior and toughening mechanisms of new carbon-fiber-reinforced bismaleimide composites

The interlaminar fracture behaviour of carbon fibre-reinforced bismaleimide (BMI) composites prepared by using a new modified BMI matrix has been investigated by various methods. Laminates of three typical stacking sequences were evaluated. Double cantilever beam, end-notch flexure and edge-delamination tension tests were conducted under conventional conditions and in a scanning electron microscope. The strain energy release rates in Mode I and Mode III G(lc) and G(llc), as well as the total strain energy release rate, G(mc), have been determined and found to be higher than those for laminates with an epoxy matrix. Dynamic delamination propagation was also studied. The toughening mechanisms are discussed.

Fatigue fracture-behavior of carbon-fiber-reinforced modified bismaleimide composites

The fracture toughness and fatigue fracture behaviour of carbon-fiber-reinforced modified bismaleimide (BMI) composites have been studied. These composites were found to have higher fracture toughnes, better damage tolerance and longer fatigue life than carbon-fiber composites with epoxy matrices. Delamination is the major mode of failure in fatigue and it is controlled by the properties of the matrix and interface. The improved performance is dire to the presence of thermoplastic particles in the modified BMI matrix which gives rise to enhanced fiber/matrix adhesion and more extensive plastic deformation. The fatigue behaviour also depends on the stacking sequence, with the multidirectional [45/90/-45/0] fiber-reinforced modified BMI composite having a lower crack propagation rate and longer fatigue life than the unidirectional laminate. This arises because of the constraint on the damage processes due to the different fiber orientation in the plies.

Statistical description of pattern evolution in damage-fracture

A model of dynamical process and stochastic jump has been put forward to study the pattern evolution in damage-fracture. According to the final states of evolution processes, the evolution modes can be classified as globally stable modes (GS modes) and evolution induced catastrophic modes (ElC modes); the latter are responsible for fracture. A statistical description is introduced to clarify the pattern evolution in this paper. It is indicated that the appearance of fracture in disordered materials should be depicted by probability distribution function.