Numerical simulations of fracture initiation in ductile solids under mode I, dynamic loading

In this paper, an overview of some recent computational studies by the authors on ductile crack initiation under mode I, dynamic loading is presented. In these studies, a large deformation finite element procedure is employed along with the viscoplastic version of the Gurson constitutive model that accounts for the micro-mechanical processes of void nucleation, growth and coalescence. A three-point bend fracture specimen subjected to impact, and a single edge notched specimen loaded by a tensile stress pulse are analysed. Several loading rates are simulated by varying the impact speed or the rise time and magnitude of the stress pulse. A simple model involving a semi-circular notch with a pre-nucleated circular hole situated ahead of it is considered. The growth of the hole and its interaction with the notch tip, which leads to plastic strain and porosity localization in the ligament connecting them, is simulated. The role of strain-rate dependence on ductile crack initiation at high loading rates, and the specimen geometry effect on the variation of dynamic fracture toughness with loading rate are Investigated.

Dislocation mechanisms for plastic flow of nickel in the temperature range 4.2 – 1200 K

The temperature ranges of thermal and athermal deformation behaviour of nickel are identified by employing the temperature-dependence of flow-stress and strain-rate cycling data. The results are used to present a unified view of dislocation mechanisms of glide encompassing the two thermally activated and the intermediate athermal regimes of plastic flow.In the low-temperature thermally activated region (<250 K) the strain rate is found to be controlled by the repulsive intersection of glide and forest dislocations, in accordance with current ideas. The athermal stress in this region can be attributed mainly to the presence of strong attractive junctions which are overcome by means of Orowan bowing, a small contribution also coming from the elastic interactions between dislocations. The values of activation area and activation energy obtained in the high-temperature region (> 750 K) negate the operation of a diffusion-controlled mechanism. Instead, the data support a thermal activation model involving unzipping of the attractive junctions. The internal (long-range) stress contribution here results solely from the elastic interactions between dislocations. This view concerning the high-temperature plastic flow is further supported by the observation that the Cottrell–Stokes law is obeyed over large strains in the range 750–1200 K.

Energy loss due to adhesion in longitudinal impact of elastic cylinders

Adhesive interaction between impacting bodies can cause energy loss, even in an otherwise elastic impact. Adhesion force induces tensile stress in the bodies, which modifies the stress wave profile and influences the restitution behavior. We investigate this effect by developing a finite element framework, which incorporates a Lennard-Jones-type potential for modeling the adhesive interaction between volume elements. With this framework, the classical problems in contact mechanics can be revisited without the restrictive surface-force approximation. In this paper, we study the longitudinal impact of an elastic cylinder on a rigid half-space with adhesion. In the absence of adhesion, this problem reduces to the impact between two identical cylinders in which there is no energy loss. Adhesion causes a fraction of energy in the stress waves to remain in the cylinder as residual stress waves. This apparent loss in kinetic energy is shown to be a unique function of maximum tensile strain energy. We have developed a 1-D model in terms of interaction force parameters, velocity and material properties to estimate the tensile stain energy. We show that this model can be used to predict practically important phenomena like capture wherein the impacting bodies stick together. (C) 2013 Elsevier Masson SAS. All rights reserved.

General framework for acoustic emission during plastic deformation

Despite the long history, so far there is no general theoretical framework for calculating the acoustic emission spectrum accompanying any plastic deformation. We set up a discrete wave equation with plastic strain rate as a source term and include the Rayleigh-dissipation function to represent dissipation accompanying acoustic emission. We devise a method of bridging the widely separated time scales of plastic deformation and elastic degrees of freedom. While this equation is applicable to any type of plastic deformation, it should be supplemented by evolution equations for the dislocation microstructure for calculating the plastic strain rate. The efficacy of the framework is illustrated by considering three distinct cases of plastic deformation. The first one is the acoustic emission during a typical continuous yield exhibiting a smooth stress-strain curve. We first construct an appropriate set of evolution equations for two types of dislocation densities and then show that the shape of the model stress-strain curve and accompanying acoustic emission spectrum match very well with experimental results. The second and the third are the more complex cases of the Portevin-Le Chatelier bands and the Luders band. These two cases are dealt with in the context of the Ananthakrishna model since the model predicts the three types of the Portevin-Le Chatelier bands and also Luders-like bands. Our results show that for the type-C bands where the serration amplitude is large, the acoustic emission spectrum consists of well-separated bursts of acoustic emission. At higher strain rates of hopping type-B bands, the burst-type acoustic emission spectrum tends to overlap, forming a nearly continuous background with some sharp acoustic emission bursts. The latter can be identified with the nucleation of new bands. The acoustic emission spectrum associated with the continuously propagating type-A band is continuous. These predictions are consistent with experimental results. More importantly, our study shows that the low-amplitude continuous acoustic emission spectrum seen in both the type-B and type-A band regimes is directly correlated to small-amplitude serrations induced by propagating bands. The acoustic emission spectrum of the Luders-like band matches with recent experiments as well. In all of these cases, acoustic emission signals are burstlike, reflecting the intermittent character of dislocation-mediated plastic flow.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Dislocation theory of the fracture criterion for anisotropic solids

A dislocation theory of fracture criterion for the mixed dislocation emission and cleavage process in an anisotropic solid is developed in this paper. The complicated cases involving mixed-mode loading are considered here. The explicit formula for dislocations interaction with a semi-infinite crack is obtained. The governing equation for the critical condition of crack cleavage in an anisotropic solid after a number dislocation emissions is established. The effects of elastic anisotropy, crack geometry and load phase angle on the critical energy release rate and the total number of the emitted dislocations at the onset of cleavage are analysed in detail. The analyses revealed that the critical energy release rates can increase to one or two magnitudes larger than the surface energy because of the dislocation emission. It is also found elastic anisotropy and crystal orientation have significant effects on the critical energy release rates. The anisotropic values can be several times the isotropic value in one crack orientation. The values may be as much as 40% less than the isotropic value in another crack orientation and another anisotropy parameter. Then the theory is applied to a fee single crystal. An edge dislocation can emit from the crack tip along the most highly shear stressed slip plane. Crack cleavage can occur along the most highly stressed slip plane after a number of dislocation emissions. Calculation is carried out step by step. Each step we should judge by which slip system is the most highly shear stressed slip system and which slip system has the largest energy release rate. The calculation clearly shows that the crack orientation and the load phase angle have significant effects on the crystal brittle-ductile behaviours.

An elasto-plastic constitutive relation of whisker-reinforced metal-matrix composite

A material model for whisker-reinforced metal-matrix composites is constructed that consists of three kinds of essential elements: elastic medium, equivalent slip system, and fiber-bundle. The heterogeneity of material constituents in position is averaged, while the orientation distribution of whiskers and slip systems is considered in the structure of the material model. Crystal and interface sliding criteria are addressed. Based on the stress-strain response of the model material, an elasto-plastic constitutive relation is derived to discuss the initial and deformation induced anisotropy as well as other fundamental features. Predictions of the present theory for unidirectional-fiber-reinforced aluminum matrix composites are favorably compared with FEM results.

Mode I crack growth resistance of metallic foams

A Dugdale-type cohesive zone model is used to predict the mode I crack growth resistance (R-curve) of metallic foams, with the fracture process characterized by an idealized traction-separation law that relates the crack surface traction to crack opening displacement. A quadratic yield function, involving the von Mises effective stress and mean stress, is used to account for the plastic compressibility of metallic foams. Finite element calculations are performed for the crack growth resistance under small scale yielding and small scale bridging in plane strain, with K-field boundary conditions. The following effects upon the fracture process are quantified: material hardening, bridging strength, T-stress (the non-singular stress acting parallel to the crack plane), and the shape of yield surface. To study the failure behaviour and notch sensitivity of metallic foams in the presence of large scale yielding, a study is made for panels embedded with either a centre-crack or an open hole and subjected to tensile stressing. For the centre-cracked panel, a transition crack size is predicted for which the fracture response switches from net section yielding to elastic-brittle fracture. Likewise, for a panel containing a centre-hole, a transition hole diameter exists for which the fracture response switches from net section yielding to a local maximum stress criterion at the edge of the hole.

Non-slipping Adhesive Contact between Mismatched Elastic Cylinders

要: We have recently proposed a generalized JKR model for non-slipping adhesive contact between two elastic spheres subjected to a pair of pulling forces and a mismatch strain (Chen, S., Gao, H., 2006c. Non-slipping adhesive contact between mismatched elastic spheres: a model of adhesion mediated deformation sensor. J. Mech. Phys. Solids 54, 1548-1567). Here we extend this model to adhesion between two mismatched elastic cylinders. The attention is focused on how the mismatch strain affects the contact area and the pull-off force. It is found that there exists a critical mismatch strain at which the contact spontaneously dissociates. The analysis suggests possible mechanisms by which mechanical deformation can affect binding between cells and molecules in biology.

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic–plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a “hold” period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and “hold-at-the-maximum-indenter-displacement” for determining the instantaneous modulus using spherical indenters.

Experimental Verification And Theoretical Analysis Of The Relationships Between Hardness, Elastic Modulus, And The Work Of Indentation

The relationship between hardness (H), reduced modulus (E-r), unloading work (W-u), and total work (W-t) of indentation is examined in detail experimentally and theoretically. Experimental study verifies the approximate linear relationship. Theoretical analysis confirms it. Furthermore, the solutions to the conical indentation in elastic-perfectly plastic solid, including elastic work (W-e), H, W-t, and W-u are obtained using Johnson's expanding cavity model and Lame solution. Consequently, it is found that the W-e should be distinguished from W-u, rather than their equivalence as suggested in ISO14577, and (H/E-r)/(W-u/W-t) depends mainly on the conical angle, which are also verified with numerical simulations. (C) 2008 American Institute of Physics.

Numerical verification for instrumented spherical indentation techniques in determining the plastic properties of materials

 Instrumented indentation tests have been widely adopted for elastic modulus determination. Recently, a number of indentation-based methods for plastic properties characterization have been proposed, and rigorous verification is absolutely necessary for their wide application. In view of the advantages of spherical indentation compared with conical indentation in determining plastic proper-ties, this study mainly concerns verification of spherical indentation methods. Five convenient and simple models were selected for this purpose, and numerical experiments for a wide range of materials are carried out to identify their accuracy and sensitivity characteristics. The verification results show that four of these five methods can give relatively accurate and stable results within a certain material domain, which is defined as their validity range and has been summarized for each method.

Bio-inspired mechanics of reversible adhesion: Orientation-dependent adhesion strength for non-slipping adhesive contact with transversely isotropic elastic materials

Geckos and many insects have evolved elastically anisotropic adhesive tissues with hierarchical structures that allow these animals not only to adhere robustly to rough surfaces but also to detach easily upon movement. In order to improve Our understanding of the role of elastic anisotropy in reversible adhesion, here we extend the classical JKR model of adhesive contact mechanics to anisotropic materials. In particular, we consider the plane strain problem of a rigid cylinder in non-slipping adhesive contact with a transversely isotropic elastic half space with the axis of symmetry oriented at an angle inclined to the surface. The cylinder is then subjected to an arbitrarily oriented pulling force. The critical force and contact width at pull-off are calculated as a function of the pulling angle. The analysis shows that elastic anisotropy leads to an orientation-dependent adhesion strength which can vary strongly with the direction of pulling. This study may suggest possible mechanisms by which reversible adhesion devices can be designed for engineering applications. (C) 2006 Elsevier Ltd. All rights reserved.

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.

On the velocity of buckle propagation in a beam on a nonlinear elastic-foundation

The paper revisits a simple beam model used by Chater et al. (1983, Proc. IUTAM Symp. Collapse, Cambridge University Press) to examine the dynamics of propagating buckles on it. It was found that, if a buckle is initiated at a constant pressure higher than the propagation pressure of the model (P-p), the buckle accelerates and gradually reaches a constant velocity which depends upon the pressure, while if it is initiated at P-p, the buckle propagates at a velocity which depends upon the initial imperfection. The causes for the difference are also investigated.