Gradient plasticity with the tangential-subloading surface model and the prediction of shear-band thickness of granular materials

An extended gradient elastoplastic constitutive equation is formulated, which is capable of describing the plastic strain rate due to the rate of stress inside the yield surface and the inelastic strain rate due to the stress rate component tangential to the subloading surface by incorporating the tangential-subloading surface model. Based on the extended constitutive equation, the post-localization analysis of granular materials is performed to predict the shear-band thickness. It is revealed that the shear-band thickness is almost determined by the gradient coefficient characterizing the inhomogeneity of deformation, although the stress–strain curve is strongly dependent on material properties.     

The formation of tabular compaction-band arrays: Theoretical and numerical analysis

The bifurcation analysis of compaction banding is extended to the formation of a tabular discrete compaction-band array. This analysis, taken together with the results of finite-difference simulations, shows that the bifurcation results in the formation of intermittent loading (elastic-plastic) and unloading (elastic) bands. The obtained analytical solution relates the spacing parameter χ (the ratio between the band thickness to the band-to-band distance) to all constitutive and stress-state parameters. Both this solution and numerical models reveal strong dependence of χ on the hardening modulus h: χ increases with h reduction. The band thickness in the numerical models is mesh dependent, but in terms of mesh-zone-size varies only from ∼2 to 4 depending on the constitutive parameters and independently on the mesh resolution. The thickness of the “elementary” compaction bands in real granular materials is equal to a few grain sizes. It follows that one grid zone in the numerical models corresponds approximately to one grain in the real material. The numerical models reproduce both discrete and continuous propagating compaction banding observed in the rock samples. These phenomena were shown to be dependent on the evolution of h and the dilatancy factor with deformation.     

Simulation and interpretation of diffuse mode bifurcation of elastoplastic solids

The diffuse mode bifurcation of elastoplastic solids at finite strain is investigated. The multiplicative decomposition of deformation gradient and the hyperelasto-plastic constitutive relationship are adapted to the numerical bifurcation analysis of the elastoplastic solids. First, bifurcation analyses of rectangular plane strain specimens subjected to uniaxial compression are conducted. The onset of the diffuse mode bifurcations from a homogeneous state is detected; moreover, the post-bifurcation states for these modes are traced to arrive at localization to narrow band zones, which look like shear bands. The occurrence of diffuse mode bifurcation, followed by localization, is advanced as a possible mechanism to create complex deformation and localization patterns, such as shear bands. These computational diffuse modes and localization zones are shown to be in good agreement with the associated experimental ones observed for sand specimens to ensure the validity of this mechanism. Next, the degradation of horizontal sway stiffness of a rectangular specimen due to plane strain uniaxial compression is pointed out as a cause of the bifurcation of the first antisymmetric diffuse mode, which triggers the tilting of the specimen. Last, circular and punching failures of a footing on a foundation are simulated.     

Elastoplastic microscopic bifurcation and post-bifurcation behavior of periodic cellular solids

In this study, a general framework is developed to analyze microscopic bifurcation and post-bifurcation behavior of elastoplastic, periodic cellular solids. The framework is built on the basis of a two-scale theory, called a homogenization theory, of the updated Lagrangian type. We thus derive the eigenmode problem of microscopic bifurcation and the orthogonality to be satisfied by the eigenmodes. The orthogonality allows the macroscopic increments to be independent of the eigenmodes, resulting in a simple procedure of the elastoplastic post-bifurcation analysis based on the notion of comparison solids. By use of this framework, then, bifurcation and post-bifurcation analysis are performed for cell aggregates of an elastoplastic honeycomb subject to in-plane compression. Thus, demonstrating a basic, long-wave eigenmode of microscopic bifurcation under uniaxial compression, it is shown that the eigenmode has the longitudinal component dominant to the transverse component and consequently causes microscopic buckling to localize in a cell row perpendicular to the loading axis. It is further shown that under equi-biaxial compression, the flower-like buckling mode having occurred in a macroscopically stable state changes into an asymmetric, long-wave mode due to the sextuple bifurcation in a macroscopically unstable state, leading to the localization of microscopic buckling in deltaic areas.     

Prediction of ductile fracture in tension by bifurcation,localization, and imperfection analyses

In this investigation, it is shown that the onset of ductile fracture in tension can be interpreted as the result of a supercritical bifurcation of homogeneous deformation and that this fact can be applied to predict ductile fracture initiation of materials with general imperfections or flaws. We focus on one dimensional quasi-static simple tension for rate-independent isotropic plastic materials. For deformation beyond the bifurcation point, multiple equilibrium paths appear. The homogeneous deformation, as one of the equilibrium paths, loses stability while the inhomogeneous paths are stable, thus indicating the occurrence of strain localization. This investigation also provides a physical example for the application of the Lambert W function in material localization analyses. Material instability is treated as the instability of a static system with dynamic perturbation. We also address the presence of microvoids in a power law plastic material as an unfolding of the supercritical pitchfork bifurcation. The imperfect system, idealized as spherical voids within the plastic matrix, is analyzed using the familiar Gurson model which is based on the presumption of a randomly voided material and characterized by the volume fraction of voids. If, in addition, the sizes of the microvoids are known, this then provides a length scale for the imperfection zone. In this manner, relevance to the sample size effects of strain-to-failure for ductile fracture initiation is addressed by considering separate zones with variations in void volume fractions. Fracture initiation predictions are presented and compare very well to existing experimental results.     

A constitutive model for plastically anisotropic solids with non-spherical voids

Plastic constitutive relations are derived for a class of anisotropic porous materials consisting of coaxial spheroidal voids, arbitrarily oriented relative to the embedding orthotropic matrix. The derivations are based on nonlinear homogenization, limit analysis and micromechanics. A variational principle is formulated for the yield criterion of the effective medium and specialized to a spheroidal representative volume element containing a confocal spheroidal void and subjected to uniform boundary deformation. To obtain closed form equations for the effective yield locus, approximations are introduced in the limit-analysis based on a restricted set of admissible microscopic velocity fields. Evolution laws are also derived for the microstructure, defined in terms of void volume fraction, aspect ratio and orientation, using material incompressibility and Eshelby-like concentration tensors. The new yield criterion is an extension of the well known isotropic Gurson model. It also extends previous analyses of uncoupled effects of void shape and material anisotropy on the effective plastic behavior of solids containing voids. Preliminary comparisons with finite element calculations of voided cells show that the model captures non-trivial effects of anisotropy heretofore not picked up by void growth models.     

Microscopic symmetric bifurcation condition of cellular solids based on a homogenization theory of finite deformation

In this paper, we establish a homogenization framework to analyze the microscopic symmetric bifurcation buckling of cellular solids subjected to macroscopically uniform compression. To this end, describing the principle of virtual work for infinite periodic materials in the updated Lagrangian form, we build a homogenization theory of finite deformation, which satisfies the principle of material objectivity. Then, we state a postulate that at the onset of microscopic symmetric bifurcation, microscopic velocity becomes spontaneous, yet changing the sign of such spontaneous velocity has no influence on the variation in macroscopic states. By applying this postulate to the homogenization theory, we derive the conditions to be satisfied at the onset of microscopic symmetric bifurcation. The resulting conditions are verified by analyzing numerically the in-plane biaxial buckling of an elastic hexagonal honeycomb. It is thus shown that three kinds of experimentally observed buckling modes of honeycombs i.e., uniaxial, biaxial and flower-like modes, are attained and classified as microscopic symmetric bifurcation. It is also shown that the multiplicity of bifurcation gives rise to the complex cell-patterns in the biaxial and flower-like modes.     

Instabilities in bimetallic layers

The analysis of the thickness fluctuations in the extrusion process of bimetallic tubes has motivated the present work. For an orthotropic, incrementally-linear solid, the possibility of non-uniform solutions near the bimaterial interface is investigated by assuming an initial perturbation along the interface. The bifurcation equation for the problem is established and firstly solved numerically to obtain a critical strain in terms of a corresponding wavenumber. Then the bifurcation equation is solved for strains above the critical strain in order to analyse the growth of the perturbation, by means of an instability parameter introduced in the analysis. A set of values for several selected parameters of the process and three constitutive equations are considered in the computations. Their influence in the stability of the process is discussed in detail.     

Micromechanics of coalescence in ductile fracture

Significant progress has been recently made in modelling the onset of void coalescence by internal necking in ductile materials. The aim of this paper is to develop a micro-mechanical framework for the whole coalescence regime, suitable for finite-element implementation. The model is defined by a set of constitutive equations including a closed form of the yield surface along with appropriate evolution laws for void shape and ligament size. Normality is still obeyed during coalescence. The derivation of the evolution laws is carefully guided by coalescence phenomenology inferred from micromechanical unit-cell calculations. The major implication of the model is that the stress carrying capacity of the elementary volume vanishes as a natural outcome of ligament size reduction. Moreover, the drop in the macroscopic stress accompanying coalescence can be quantified for many initial microstructures provided that the microstructure state is known at incipient coalescence. The second part of the paper addresses a more practical issue, that is the prediction of the acceleration rate δ in the Tvergaard-Needleman phenomenological approach to coalescence. For that purpose, a Gurson-like model including void shape effects is used. Results are presented and discussed in the limiting case of a non-hardening material for different initial microstructures and various stress states. Predicted values of δ are extremely sensitive to stress triaxiality and initial spacing ratio. The effect of initial porosity is significant at low triaxiality whereas the effect of initial void shape is emphasized at high triaxiality.     

Perturbation growth and localization in fluid-saturated inelastic porous media under quasi-static loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under quasi-static loadings using a perturbation approach and focussing in particular on the two limiting cases of the onset of growth and of blowing up of perturbations.For associative flow rules for the skeleton, both onset of growth and blowing up of perturbations depend only on the underlying drained properties. Unbounded growth is obtained when the condition of localization for the underlying drained deformation (singularity of the drained acoustic tensor) is approached or just passed. Onset of growth has always a divergence growth character and critical conditions are always associated to the shortwavelength regime leading to the fact that the failure mode is expected to be a localized one.For non-associative behaviour of the skeleton we show in contrast that the onset of growth and unbounded growth may be defined either by the drained or undrained properties. One or the other depends on the details of the constitutive behaviour but also on the type of loadings. In particular, unbounded growth occurs when either the condition of localization under drained or undrained conditions is first passed. Transition from decaying to growing behaviour may have a divergence character or flutter-type character. Here the critical conditions are associated either to the shortwavelength or to the longwavelength regimes and therefore the failure mode may be localized or diffuse.The hierarchy between criticality of drained and undrained properties is analysed for a general class of constitutive equations and the results are fully and explicitly illustrated for saturated porous media with skeleton obeying Drucker-Prager like constitutive model.     

Analysis of material instabilities in inelastic solids by incremental energy minimization and relaxation methods: evolving deformation microstructures in finite plasticity

We propose an approach to the definition and analysis of material instabilities in rate-independent standard dissipative solids at finite strains based on finite-step-sized incremental energy minimization principles. The point of departure is a recently developed constitutive minimization principle for standard dissipative materials that optimizes a generalized incremental work function with respect to the internal variables. In an incremental setting at finite time steps this variational problem defines a quasi-hyperelastic stress potential. The existence of this potential allows to be recast a typical incremental boundary-value problem of quasi-static inelasticity into a principle of minimum incremental energy for standard dissipative solids. Mathematical existence theorems for sufficiently regular minimizers then induce a definition of the material stability of the inelastic material response in terms of the sequentially weakly lower semicontinuity of the incremental variational functional. As a consequence, the incremental material stability of standard dissipative solids may be defined in terms of the quasi-convexity or the rank-one convexity of the incremental stress potential. This global definition includes the classical local Hadamard condition but is more general. Furthermore, the variational setting opens up the possibility to analyze the post-critical development of deformation microstructures in non-stable inelastic materials based on energy relaxation methods. We outline minimization principles of quasi- and rank-one convexifications of incremental non-convex stress potentials for standard dissipative solids. The general concepts are applied to the analysis of evolving deformation microstructures in single-slip plasticity. For this canonical model problem, we outline details of the constitutive variational formulation and develop numerical and semi-analytical solution methods for a first-level rank-one convexification. A set of representative numerical investigations analyze the development of deformation microstructures in the form of rank-one laminates in single slip plasticity for homogeneous macro-deformation modes as well as inhomogeneous macroscopic boundary-value problems. The well-posedness of the relaxed variational formulation is indicated by an independence of typical finite element solutions on the mesh-size.     

A fundamental model of cyclic instabilities in thermal barrier systems

Cyclic morphological instabilities in the thermally grown oxide (TGO) represent a source of failure in some thermal barrier systems. Observations and simulations have indicated that several factors interact to cause these instabilities to propagate: (i) thermal cycling; (ii) thermal expansion misfit; (iii) oxidation strain; (iv) yielding in the TGO and the bond coat; and (v) initial geometric imperfections. This study explores a fundamental understanding of the propagation phenomenon by devising a spherically symmetric model that can be solved analytically. The applicability of this model is addressed through comparison with simulations conducted for representative geometric imperfections and by analogy with the elastic/plastic indentation of a half space. Finite element analysis is used to confirm and extend the model. The analysis identifies the dependencies of the instability on the thermo-mechanical properties of the system. The crucial role of the in-plane growth strain is substantiated, as well as the requirement for bond coat yielding. It is demonstrated that yielding of the TGO is essential and is, in fact, the phenomenon that differentiates between cyclic and isothermal responses.     

Bounds and estimates for the effect of strain gradients upon the effective plastic properties of an isotropic two-phase composite

Predictions are made for the size effect on strength of a random, isotropic two-phase composite. Each phase is treated as an isotropic, elastic-plastic solid, with a response described by a modified deformation theory version of the Fleck-Hutchinson strain gradient plasticity formulation (Fleck and Hutchinson, J. Mech. Phys. Solids 49 (2001) 2245). The essential feature of the new theory is that the plastic strain tensor is treated as a primary unknown on the same footing as the displacement. Minimum principles for the energy and for the complementary energy are stated for a composite, and these lead directly to elementary bounds analogous to those of Reuss and Voigt. For the case of a linear hardening solid, Hashin-Shtrikman bounds and self-consistent estimates are derived. A non-linear variational principle is constructed by generalising that of Ponte Castañeda (J. Mech. Phys. Solids 40 (1992) 1757). The minimum principle is used to derive an upper bound, a lower estimate and a self-consistent estimate for the overall plastic response of a statistically homogeneous and isotropic strain gradient composite. Sample numerical calculations are performed to explore the dependence of the macroscopic uniaxial response upon the size scale of the microstructure, and upon the relative volume fraction of the two phases.     

Experimental investigation and modeling of non-monotonic creep behavior in polymers

The deformation behavior of two unfilled engineering thermoplastics, ultra high molecular weight polyethylene (UHMWPE) and polycarbonate (PC), has been investigated in creep test conditions. It has been found that a loading history (prior to the creep test) comprising of loading to a maximum stress or strain value followed by partial unloading to arrive at the target stress value can greatly modify the strain-time behavior. Under such a test protocol, while the expected increase in strain during creep (constant tensile load) is observed, at relatively low creep stresses specimens have also demonstrated a monotonic decrease in strain. In an intermediate stress range, specimens have demonstrated time dependent behavior comprising of a transition from decreasing to increasing strain during creep in tension. This paper presents experimental results to delineate these findings and explore the effect of prior strain rate on the qualitative and quantitative changes in the output (strain-time) behavior. Furthermore, modification of the viscoplasticity theory based on overstress (VBO) model into a double element configuration is introduced. These changes confer upon the model the ability to yield non-monotonic behavior in creep, and supporting simulation results have been included. These changes, therefore, allow the model to simulate strain rate sensitivity, creep, relaxation, and recovery behavior, but more importantly address the issue of non-monotonic changes in creep and relaxation when a loading history involves some degree of unloading.     

Numerical study of sliding wear caused by a loaded pin on a rotating disc

A computational approach is proposed to predict the sliding wear caused by a loaded spherical pin contacting a rotating disc, a condition typical of the so-called pin-on-disc test widely used in tribological studies. The proposed framework relies on the understanding that, when the pin contacts and slides on the disc, a predominantly plane strain region exists at the centre of the disc wear track. The wear rate in this plane strain region can therefore be determined from a two dimensional idealisation of the contact problem, reducing the need for computationally expensive three dimensional contact analyses. Periodic unit cell techniques are used in conjunction with a ratchetting-based failure criterion to predict the wear rate in the central plane strain region. The overall three dimensional wear rate of the disc is then determined by scaling the plane strain wear rate with a conversion factor related to the predicted shape of the wear track. The approach is used to predict pin-on-disc test data from an Al-Si coating using a tungsten carbide pin. The predicted results are found to be consistent with measured data.     

Kink band instability in layered structures

A recent two-dimensional prototype model for the initiation of kink banding in compressed layered structures is extended to embrace the two propagation mechanisms of band broadening and band progression. As well as interlayer friction, overburden pressure and layer bending energy, the characteristics of transverse layer compressibility and foundation stiffness are now included. Experiments on constrained layers of paper show good agreement with the predictions of angle of orientation, kink band width and post-kink load-deflection response obtained from the model.