On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Elastic-plastic analysis of complete building structures

In this Thesis, details of a proposed method for the elastic-plastic failure load analysis of complete building structures are given. In order to handle the problem, a computer programme in Atlas Autocode is produced. The structures consist of a number of parallel shear walls and intermediate frames connected by floor slabs. The results of an experimental investigation are given to verify the theoretical results and to demonstrate various factors that may influence the behaviour of these structures. Large full scale practical structures are also analysed by the proposed method and suggestions are made for achieving design economy as well as for extending research in various aspects of this field. The existing programme for elastic-plastic analysis of large frames is modified to allow for the effect of composite action of structural members, i.e. reinforced concrete floor slabs and the supporting steel beams. This modified programme is used to analyse some framed type structures with composite action as well as those which incorporate plates and shear walls. The results obtained are studied to ascertain the influence of composite action and other factors on the load carrying capacity of both bare frames and complete building structures. The theoretical failure load presented in this thesis does not predict the overall failure load of the structure nor does it predict the partial failure load of the shear walls and slabs but it merely predicts the partial failure load of a single frame and assumes that the loss of stiffess of such a frame renders the overall structure unusable. For most structures the analysis proposed in this thesis is likely to break down prematurely due to the failure of the slab and shear wall system and this factor must be taken into account in any future work on such structures. The experimental work reported in this thesis is acknowledged to be unsatisfactory as a verification of the limited theory proposed. In particular perspex was not found to be a suitable material for testing at high loads, micro-concrete may be more suitable.

Dynamics of Granular Crystals with Elastic-Plastic Contacts

We study the behavior of granular crystals subjected to impact loading that creates plastic deformation at the contacts between constituent particles. Granular crystals are highly periodic arrangements of spherical particles, arranged into densely packed structures resembling crystals. This special class of granular materials has been shown to have unique dynamics with suggested applications in impact protection. However, previous work has focused on very low amplitude impacts where every contact point can be described using the Hertzian contact law, valid only for purely elastic deformation. In this thesis, we extend previous investigation of the dynamics of granular crystals to significantly higher impact energies more suitable for the majority of applications. Additionally, we demonstrate new properties specific to elastic-plastic granular crystals and discuss their potential applications as well. We first develop a new contact law to describe the interaction between particles for large amplitude compression of elastic-plastic spherical particles including a formulation for strain-rate dependent plasticity. We numerically and experimentally demonstrate the applicability of this contact law to a variety of materials typically used in granular crystals. We then extend our investigation to one-dimensional chains of elastic-plastic particles, including chains of alternating dissimilar materials. We show that, using the new elastic-plastic contact law, we can predict the speed at which impact waves with plastic dissipation propagate based on the material properties of the constituent particles. Finally, we experimentally and numerically investigate the dynamics of two-dimensional and three-dimensional granular crystals with elastic-plastic contacts. We first show that the predicted wave speeds for 1D granular crystals can be extended to 2D and 3D materials. We then investigate the behavior of waves propagating across oblique interfaces of dissimilar particles. We show that the character of the refracted wave can be predicted using an analog to Snell's law for elastic-plastic granular crystals and ultimately show how it can be used to design impact guiding "lenses" for mitigation applications.

Effects of 'sinking in' and 'piling up' on estimating the contact area under load in indentation

The phenomena of the 'piling up' and 'sinking-in' of surface profiles in conical indentation in elastic-plastic solids with work hardening are studied using dimensional and finite-element analysis. The degree of sinking in and piling up is shown to depend on the ratio of the initial yield strength Y to Young's modulus E and on the work-hardening exponent n. The widely used procedure proposed by Oliver and Pharr for estimating contact depth is then evaluated systematically. By comparing the contact depth obtained directly from finite-element calculations with that obtained from the initial unloading slope using the Oliver-Pharr procedure, the applicability of the procedure is discussed.

Can stress-strain relationships be obtained from indentation curves using conical and pyramidal indenters?

Applying the scaling relationships developed recently for conical indentation in elastic-plastic solids with work-hardening, we examine the question of whether stress-strain relationships of such solids can be uniquely determined by matching the calculated loading and unloading curves with that measured experimentally. We show that there can be multiple stress-strain curves for a given set of loading and unloading curves. Consequently, stress-strain relationships may not be uniquely determined from loading and unloading curves alone using a conical or pyramidal indenter.

An Energy-Based Method for Analyzing Instrumented Spherical Indentation Experiments

Using dimensional analysis and finite element calculation, we studied spherical indentation in elastic-plastic solids with work hardening. We report two previously unknown relationships between hardness, reduced modulus, indentation depth, indenter radius, and work of indentation. These relationships, together with the relationship between initial unloading stiffness and reduced modulus, provide an energy-based method for determining contact area, reduced modulus, and hardness of materials from instrumented spherical indentation measurements. This method also provides a means for calibrating the effective radius of imperfectly shaped spherical indenters. Finally, the method is applied to the analysis of instrumented spherical indentation experiments on copper, aluminum, tungsten, and fused silica.

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001 a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.

Scaling, Dimensional Analysis, and Indentation Measurements

We provide an overview of the basic concepts of scaling and dimensional analysis, followed by a review of some of the recent work on applying these concepts to modeling instrumented indentation measurements. Specifically, we examine conical and pyramidal indentation in elastic-plastic solids with power-law work-hardening, in power-law creep solids, and in linear viscoelastic materials. We show that the scaling approach to indentation modeling provides new insights into several basic questions in instrumented indentation, including, what information is contained in the indentation load-displacement curves? How does hardness depend on the mechanical properties and indenter geometry? What are the factors determining piling-up and sinking-in of surface profiles around indents? Can stress-strain relationships be obtained from indentation load-displacement curves? How to measure time dependent mechanical properties from indentation? How to detect or confirm indentation size effects? The scaling approach also helps organize knowledge and provides a framework for bridging micro- and macroscales. We hope that this review will accomplish two purposes: (1) introducing the basic concepts of scaling and dimensional analysis to materials scientists and engineers, and (2) providing a better understanding of instrumented indentation measurements.

A constitutive model for cavitation and cavity growth in rubber-like materials under arbitrary tri-axial loading

In this paper, we attempted to construct a constitutive model to deal with the phenomenon of cavitation and cavity growth in a rubber-like material subjected to an arbitrary tri-axial loading. To this end, we considered a spherical elementary representative volume in a general Rivlin's incompressible material containing a central spherical cavity. The kinematics proposed by [Hou, H.S., Abeyaratne, R., 1992. Cavitation in elastic and elastic-plastic solids. J. Mech. Phys. Solids 40, 571-722] was adopted in order to construct an approximate but optimal field. In order to establish a suitable constitutive law for this class of materials, we utilized the homogenisation technique that permits us to calculate the average strain energy density of the volume. The cavity growth was considered through a physically realistic failure criterion. Combination of the constitutive law and the failure criterion enables us to describe correctly the global behaviour and the damage evolution of the material under tri-axial loading. It was shown that the present models can efficiently reproduce different stress states, varying from uniaxial to tri-axial tensions, observed in experimentations. Comparison between predicted results and experimental data proves that the proposed model is accurate and physically reasonable. Another advantage is that the proposed model does not need special identification work, the initial Rivlin's law for the corresponding incompressible material is sufficient to form the new law for the compressible material resulted from cavitation procedure. (C) 2007 Elsevier Ltd. All rights reserved.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Ultrahigh compression of water using intense heavy ion beams: laboratory planetary physics

Intense heavy ion beams offer a unique tool for generating samples of high energy density matter with extreme conditions of density and pressure that are believed to exist in the interiors of giant planets. An international accelerator facility named FAIR (Facility for Antiprotons and Ion Research) is being constructed at Darmstadt, which will be completed around the year 2015. It is expected that this accelerator facility will deliver a bunched uranium beam with an intensity of 5x10(11) ions per spill with a bunch length of 50-100 ns. An experiment named LAPLAS (Laboratory Planetary Sciences) has been proposed to achieve a low-entropy compression of a sample material like hydrogen or water (which are believed to be abundant in giant planets) that is imploded in a multi-layered target by the ion beam. Detailed numerical simulations have shown that using parameters of the heavy ion beam that will be available at FAIR, one can generate physical conditions that have been predicted to exist in the interior of giant planets. In the present paper, we report simulations of compression of water that show that one can generate a plasma phase as well as a superionic phase of water in the LAPLAS experiments.

Saturated duration for dynamic structural plastic response and fundamental periods of elastic vibration

The relationship is determined between saturated duration of rectangular pressure pulses applied to rigid, perfectly plastic structures and their fundamental periods of elastic vibration. It is shown that the ratio between the saturated duration and the fundamental period of elastic vibration of a structure is dependent upon two factors: the first one is the slenderness or thinness ratio of the structure; and the second one is the square root of ratio between the Young's elastic modulus and the yield stress of the structural material. Dimensional analysis shows that the aforementioned ratio is one of the basic similarity parameters for elastic-plastic modeling under dynamic loading.

Size Effect and Geometrical Effect of Solids in Micro-Indentation Test

Micro-indentation tests at scales of the order of sub-micron show that the measured hardness increases strongly with decreasing indent depth or indent size, which is frequently referred to as the size effect. At the same time, at micron or sub-micron scale, another effect, which is referred to as the geometrical size effects such as crystal grain size effect, thin film thickness effect, etc., also influences the measured material hardness. However, the trends are at odds with the size-independence implied by the conventional elastic-plastic theory. In the present research, the strain gradient plasticity theory (Fleck and Hutchinson) is used to model the composition effects (size effect and geometrical effect) for polycrystal material and metal thin film/ceramic substrate systems when materials undergo micro-indenting. The phenomena of the "pile-up" and "sink-in" appeared in the indentation test for the polycrystal materials are also discussed. Meanwhile, the micro-indentation experiments for the polycrystal Al and for the Ti/Si_3N_4 thin film/substrate system are carried out. By comparing the theoretical predictions with experimental measurements, the values and the variation trends of the micro-scale parameter included in the strain gradient plasticity theory are predicted.

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.