Discrete dislocation plasticity analysis of crack-tip fields in polycrystalline materials

Small scale yielding around a mode I crack is analysed using polycrystalline discrete dislocation plasticity. Plane strain analyses are carried out with the dislocations all of edge character and modelled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, nucleation, interaction with obstacles and annihilation are incorporated through a set of constitutive rules. Grain boundaries are modelled as impenetrable to dislocations. The polycrystalline material is taken to consist of two types of square grains, one of which has a bcc-like orientation and the other an fcc-like orientation. For both orientations there are three active slip systems. Alternating rows, alternating columns and a checker-board-like arrangement of the grains is used to construct the polycrystalline materials. Consistent with the increasing yield strength of the polycrystalline material with decreasing grain size, the calculations predict a decrease in both the plastic zone size and the crack-tip opening displacement for a given applied mode I stress intensity factor. Furthermore, slip-band and kink-band formation is inhibited by all grain arrangements and, with decreasing grain size, the stress and strain distributions more closely resemble the HRR fields with the crack-tip opening approximately inversely proportional to the yield strength of the polycrystalline materials. The calculations predict a reduction in fracture toughness with decreasing grain size associated with the grain boundaries acting as effective barriers to dislocation motion.

Discrete dislocation plasticity analysis of size effects in single crystals

The effect of size and slip system configuration on the tensile stress-strain response of micron-sized planar crystals as obtained from discrete dislocation plasticity simulations is presented. The crystals are oriented for either single or symmetric double slip. With the rotation of the tensile axis unconstrained, there is a strong size dependence, with the flow strength increasing with decreasing specimen size. Below a certain specimen size, the flow strength of the crystals is set by the nucleation strength of the initially present Frank-Read sources. The main features of the size dependence are the same for both the single and symmetric double slip configurations.

Cortical thickness mapping to identify focal osteoporosis in patients with hip fracture.

BACKGROUND: Individuals with osteoporosis are predisposed to hip fracture during trips, stumbles or falls, but half of all hip fractures occur in those without generalised osteoporosis. By analysing ordinary clinical CT scans using a novel cortical thickness mapping technique, we discovered patches of markedly thinner bone at fracture-prone regions in the femurs of women with acute hip fracture compared with controls. METHODS: We analysed CT scans from 75 female volunteers with acute fracture and 75 age- and sex-matched controls. We classified the fracture location as femoral neck or trochanteric before creating bone thickness maps of the outer 'cortical' shell of the intact contra-lateral hip. After registration of each bone to an average femur shape and statistical parametric mapping, we were able to visualise and quantify statistically significant foci of thinner cortical bone associated with each fracture type, assuming good symmetry of bone structure between the intact and fractured hip. The technique allowed us to pinpoint systematic differences and display the results on a 3D average femur shape model. FINDINGS: The cortex was generally thinner in femoral neck fracture cases than controls. More striking were several discrete patches of statistically significant thinner bone of up to 30%, which coincided with common sites of fracture initiation (femoral neck or trochanteric). INTERPRETATION: Femoral neck fracture patients had a thumbnail-sized patch of focal osteoporosis at the upper head-neck junction. This region coincided with a weak part of the femur, prone to both spontaneous 'tensile' fractures of the femoral neck, and as a site of crack initiation when falling sideways. Current hip fracture prevention strategies are based on case finding: they involve clinical risk factor estimation to determine the need for single-plane bone density measurement within a standard region of interest (ROI) of the femoral neck. The precise sites of focal osteoporosis that we have identified are overlooked by current 2D bone densitometry methods.

Tensile response of passivated films with climb-assisted dislocation glide

The tensile response of single crystal films passivated on two sides is analysed using climb enabled discrete dislocation plasticity. Plastic deformation is modelled through the motion of edge dislocations in an elastic solid with a lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation incorporated through a set of constitutive rules. The dislocation motion in the films is by glide-only or by climb-assisted glide whereas in the surface passivation layers dislocation motion occurs by glide-only and penalized by a friction stress. For realistic values of the friction stress, the size dependence of the flow strength of the oxidised films was mainly a geometrical effect resulting from the fact that the ratio of the oxide layer thickness to film thickness increases with decreasing film thickness. However, if the passivation layer was modelled as impenetrable, i.e. an infinite friction stress, the plastic hardening rate of the films increases with decreasing film thickness even for geometrically self-similar specimens. This size dependence is an intrinsic material size effect that occurs because the dislocation pile-up lengths become on the order of the film thickness. Counter-intuitively, the films have a higher flow strength when dislocation motion is driven by climb-assisted glide compared to the case when dislocation motion is glide-only. This occurs because dislocation climb breaks up the dislocation pile-ups that aid dislocations to penetrate the passivation layers. The results also show that the Bauschinger effect in passivated thin films is stronger when dislocation motion is climb-assisted compared to films wherein dislocation motion is by glide-only. © 2012 Elsevier Ltd.

Modeling fatigue crack growth in crystalline solids with discrete dislocation plasticity

Analyses of crack growth under cyclic loading conditions are discussed where plastic flow arises from the motion of large numbers of discrete dislocations and the fracture properties are embedded in a cohesive surface constitutive relation. The formulation is the same as used to analyse crack growth under monotonic loading conditions, differing only in the remote loading being a cyclic function of time. Fatigue, i.e. crack growth in cyclic loading at a driving force for which the crack would have arrested under monotonic loading, emerges in the simulations as a consequence of the evolution of internal stresses associated with the irreversibility of the dislocation motion. A fatigue threshold, Paris law behaviour, striations, the accelerated growth of short cracks and the scaling with material properties are outcomes of the calculations. Results for single crystals and polycrystals will be discussed.

Plane-strain discrete dislocation plasticity with climb-assisted glide motion of dislocations

A small-strain two-dimensional discrete dislocation plasticity (DDP) framework is developed wherein dislocation motion is caused by climb-assisted glide. The climb motion of the dislocations is assumed to be governed by a drag-type relation similar to the glide-only motion of dislocations: such a relation is valid when vacancy kinetics is either diffusion limited or sink limited. The DDP framework is employed to predict the effect of dislocation climb on the uniaxial tensile and pure bending response of single crystals. Under uniaxial tensile loading conditions, the ability of dislocations to bypass obstacles by climb results in a reduced dislocation density over a wide range of specimen sizes in the climb-assisted glide case compared to when dislocation motion is only by glide. A consequence is that, at least in a single slip situation, size effects due to dislocation starvation are reduced. By contrast, under pure bending loading conditions, the dislocation density is unaffected by dislocation climb as geometrically necessary dislocations (GNDs) dominate. However, climb enables the dislocations to arrange themselves into lower energy configurations which significantly reduces the predicted bending size effect as well as the amount of reverse plasticity observed during unloading. The results indicate that the intrinsic plasticity material length scale associated with GNDs is strongly affected by thermally activated processes and will be a function of temperature. © 2013 IOP Publishing Ltd.

Climb-enabled discrete dislocation plasticity

A small strain two-dimensional discrete dislocation plasticity framework coupled to vacancy diffusion is developed wherein the motion of edge dislocations is by a combination of glide and climb. The dislocations are modelled as line defects in a linear elastic medium and the mechanical boundary value problem is solved by the superposition of the infinite medium elastic fields of the dislocations and a complimentary non-singular solution that enforces the boundary conditions. Similarly, the climbing dislocations are modelled as line sources/sinks of vacancies and the vacancy diffusion boundary value problem is also solved by a superposition of the fields of the line sources/sinks in an infinite medium and a complementary non-singular solution that enforces the boundary conditions. The vacancy concentration field along with the stress field provides the climb rate of the dislocations. Other short-range interactions of the dislocations are incorporated via a set of constitutive rules. We first employ this formulation to investigate the climb of a single edge dislocation in an infinite medium and illustrate the existence of diffusion-limited and sink-limited climb regimes. Next, results are presented for the pure bending and uniaxial tension of single crystals oriented for single slip. These calculations show that plasticity size effects are reduced when dislocation climb is permitted. Finally, we contrast predictions of this coupled framework with an ad hoc model in which dislocation climb is modelled by a drag-type relation based on a quasi steady-state solution. © 2013 Elsevier Ltd. All rights reserved.