Atomistic simulations of crack nucleation and intergranular fracture in bulk nanocrystalline nickel

We report large scale molecular dynamics simulations of dynamic cyclic uniaxial tensile deformation of pure, fully dense nanocrystalline Ni, to reveal the crack initiation, and consequently intergranular fracture is the result of coalescence of nanovoids by breaking atomic bonds at grain boundaries and triple junctions. The results indicate that the brittle fracture behavior accounts for the transition from plastic deformation governed by dislocation to one that is grain-boundary dominant when the grain size reduces to the nanoscale. The grain-boundary mediated plasticity is also manifested by the new grain formation and growth induced by stress-assisted grain-boundary diffusion observed in this work. This work illustrates that grain-boundary decohesion is one of the fundamental deformation mechanisms in nanocrystalline Ni.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

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.

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.

Fatigue-crack closure measurements on 2024-t3 sheet specimens

A computer-controlled procedure has been developed for automatic measurement of the crack opening stress S-op during fatigue tests. A crack opening displacement gauge (GOD meter) is used to obtain digital data on the load versus COD curves. Three methods for deriving S-op from the data sets are compared: (1) a slope method, (2) a tangent lines intersecting method, and (3) a tangent point method. The effect of the position of the COD meter with respect to the crack tip on S-op is studied in tests of 2024-T3 specimens. Results of crack growth and S-op are presented for CA loading with an overload, and with an overload followed by an underload.

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.