Deformation field evolution in indentation of a porous brittle solid

An experimental study of plane strain wedge indentation of a model porous brittle solid has been made to understand the effect of indentation parameters on the evolution of the deformation field and the accompanying volume change. A series of high-speed, high-resolution images of the indentation region and simultaneous measurements of load response were captured through the progression of the indentation process. Particle image velocimetry analysis of the images facilitated in situ measurement of the evolution of the resulting plastic zone in terms of incremental material displacement (velocity), strain rate, strain and volume change (e.g., local pore collapse). These measurements revealed initiation and propagation of flow localizations and fractures, as well as enabled estimate of volume changes occurring in the deformation zone. The results were directly compared with theoretical estimates of indentation pressure and deformation zone geometry and were used to validate a modified cavity expansion solution that incorporates effects of volume changes in the plastic zone. (C) 2015 Elsevier Ltd. All rights reserved.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

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.

The investigation on dynamic fracture behavour of materials under compressive-shear combined stress waves

In this paper, an improved plate impact experimental technique is presented for studying dynamic fracture mechanism of materials, under the conditions that the impacting loading is provided by a single pulse and the loading time is in the sub-microsecond range. The impacting tests are carried out on the pressure-shear gas gun. The loading rate achieved is dK/dt similar to 10(8) MPa m(1/2) s(-1). With the elimination of influence of the specimen boundary, the plane strain state of a semi-infinite crack in an infinite elastic plate is used to simulate the deformation fields of crack tip. The single pulses are obtained by using the "momentum trap" technique. Therefore, the one-time actions of the single pulse are achieved by eradicating the stress waves reflected from the specimen boundary or diffracted from the crack surfaces. In the current study, some important phenomena have been observed. The special loading of the single pulse can bring about material damage around crack tip, and affect the material behavior, such as kinking and branching of the crack propagation. Failure mode transitions from mode I to mode II crack are observed under asymmetrical impact conditions. The mechanisms of the dynamic crack propagation are consistent with the damage failure model.

Elastodynamic erosion of thermal barrier coatings

The finite element method is used to analyze the elastodynamic response of a columnar thermal barrier coating due to normal impact and oblique impact by an erosive particle. An assessment is made of the erosion by crack growth from preexisting flaws at the edge of each column: it is demonstrated that particle impacts can be sufficiently severe to give rise to columnar cracking. First, the transient stress state induced by the normal impact of a circular cylinder or a sphere is calculated in order to assess whether a 2D calculation adequately captures the more realistic 3D behavior. It is found that the transient stress states for the plane strain and axisymmetric models are similar. The sensitivity of response to particle diameter and to impact velocity is determined for both the cylinder and the sphere. Second, the transient stress state is explored for 2D oblique impact by a circular cylindrical particle and by an angular cylindrical particle. The sensitivity of transient tensile stress within the columns to particle shape (circular and angular), impact angle, impact location, orientation of the angular particle, and to the level of friction is explored in turn. The paper concludes with an evaluation of the effect of inclining the thermal barrier coating columns upon their erosion resistance. © 2011 The American Ceramic Society.

On The Interaction Of Channeling/Segmentation Cracks Of Coating: Study On The Critical Spacing

Channeling/segmentation cracks may arise in the coating subjected to in-plane tensile stress. The interaction between these multiple cracks, say the effect of the spacing between two adjacent cracks oil the behaviors of channels themselves and the interface around the interface corners, attracts wide interest. However, if the spacing is greater than a specific magniture,, namely the Critical Spacing (CS), there should be no interaction between such channeling/segmentation cracks. In this study, file mechanism of the effect of the crack spacing oil the interfacial stress around the interface corner will be Interpreted firstly. Then the existence of the CS will be verified and the relationship between the CS and the so-called stress transfer length Ill coating will be established for plane strain condition. Finally, the dependence of the stress transfer length, simultaneously of the CS, on the sensitive parameters will be investigated with finite element method and expressed with a simple empirical formula. (C) 2007 Elsevier Ltd. All rights reserved.

Non-Slipping Jkr Model For Transversely Isotropic Materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers. (c) 2007 Elsevier Ltd. All rights reserved.

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.

Transient effect on ideally plastic stationary crack-tip fields

This paper analyses the transient effect on ideally plastic stationary crack-tip fields under mode I plane strain conditions, when the inertial forces are not negligible. It is shown that the governing equation for such a problem can be expressed in formal simplicity when referred to a system of moving curvilinear coordinates, which is a generalization of the system defined by the slip-line field in quasi-static plasticity. A perturbation method of solving the equations is described and illustrated by application to problems of ideally plastic stationary crack-tip fields when the inertia forces are not negligible.

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.

Transient effect on ideally plastic stationary crack-tip fields

This paper analyses the transient effect on ideally plastic stationary crack tip fields under mode I plane strain conditions, when the inertial forces are not negligible. It is shown that the governing equation for such a problem can be expressed in formal simplicity when referred to a system of moving curvilinear coordinates, which is a generalization of the system defined by the slip-line field in quasi-static plasticity. A perturbation method of solving the equations is described and illustrated by application to problems of ideally plastic stationary crack tip fields when the inertial forces are not negligible.

Higher-order analysis of crack-tip fields in power-law hardening materials

In this paper, we present an exact higher-order asymptotic analysis on the near-crack-tip fields in elastic-plastic materials under plane strain, Mode I. A four- or five-term asymptotic series of the solutions is derived. It is found that when 1.6 < n less-than-or-equal-to 2.8 (here, n is the hardening exponent), the elastic effect enters the third-order stress field; but when 2.8< n less-than-or-equal-to 3.7 this effect turns to enter the fourth-order field, with the fifth-order field independent. Moreover, if n>3.7, the elasticity only affects the fields whose order is higher than 4. In this case, the fourth-order field remains independent. Our investigation also shows that as long as n is larger than 1.6, the third-order field is always not independent, whose amplitude coefficient K3 depends either on K1 or on both K1 and K2 (K1 and K2 arc the amplitude coefficients of the first- and second-order fields, respectively). Firmly, good agreement is found between our results and O'Dowd and Shih's numerical ones[8] by comparison.

Higher-order analysis of crack tip fields in elastic power-law hardening materials

A HIGHER-ORDER asymptotic analysis of a stationary crack in an elastic power-law hardening material has been carried out for plane strain, Mode 1. The extent to which elasticity affects the near-tip fields is determined by the strain hardening exponent n. Five terms in the asymptotic series for the stresses have been derived for n = 3. However, only three amplitudes can be independently prescribed. These are K1, K2 and K5 corresponding to amplitudes of the first-, second- and fifth-order terms. Four terms in the asymptotic series have been obtained for n = 5, 7 and 10; in these cases, the independent amplitudes are K1, K2 and K4. It is found that appropriate choices of K2 and K4 can reproduce near-tip fields representative of a broad range of crack tip constraints in moderate and low hardening materials. Indeed, fields characterized by distinctly different stress triaxiality levels (established by finite element analysis) have been matched by the asymptotic series. The zone of dominance of the asymptotic series extends over distances of about 10 crack openings ahead of the crack tip encompassing length scales that are microstructurally significant. Furthermore, the higher-order terms collectively describe a spatially uniform hydrostatic stress field (of adjustable magnitude) ahead of the crack. Our results lend support to a suggestion that J and a measure of near-tip stress triaxiality can describe the full range of near-tip states.

Singular behaviour near the tip of a sharp V-notch in a power law hardening material

This paper presents an asymptotic analysis of the near-tip stress and strain fields of a sharp V-notch in a power law hardening material. First, the asymptotic solutions of the HRR type are obtained for the plane stress problem under symmetric loading. It is found that the angular distribution function of the radial stress sigma(r) presents rapid variation with the polar angle if the notch angle beta is smaller than a critical notch angle; otherwise, there is no such phenomena. Secondly, the asymptotic solutions are developed for antisymmetric loading in the cases of plane strain and plane stress. The accurate calculation results and the detailed comparisons are given as well. All results show that the singular exponent s is changeable for various combinations of loading condition and plane problem.