Deformation and failure of a film/substrate system subjected to spherical indentation: Part 1. Experimental validation of stresses and strains derived using Hankel transform technique in an elastic film/substrate system

Our concern here is to rationalize experimental observations of failure modes brought about by indentation of hard thin ceramic films deposited on metallic substrates. By undertaking this exercise, we would like to evolve an analytical framework that can be used for designs of coatings. In Part I of the paper we develop an algorithm and test it for a model system. Using this analytical framework we address the issue of failure of columnar TiN films in Part II [J. Mater. Res. 21, 783 (2006)] of the paper. In this part, we used a previously derived Hankel transform procedure to derive stress and strain in a birefringent polymer film glued to a strong substrate and subjected to spherical indentation. We measure surface radial strains using strain gauges and bulk film stresses using photo elastic technique (stress freezing). For a boundary condition based on Hertzian traction with no film interface constraint and assuming the substrate constraint to be a function of the imposed strain, the theory describes the stress distributions well. The variation in peak stresses also demonstrates the usefulness of depositing even a soft film to protect an underlying substrate.

On the hardness and elastic modulus of bulk metallic glass matrix composites

The influences of the amorphous matrix and crystalline dendrite phases on the hardness and elastic moduli of Zr/Ti-based bulk metallic glass matrix composites have been assessed. While the moduli of the composites correspond to those predicted by the rule of mixtures, the hardness of the composites is similar to that of the matrix, suggesting that the plastic flow in the composites under constrained conditions such as indentation is controlled by the flow resistance of the contiguous matrix. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

On the two-dimensional thermoelectric power in quantum wells of non-parabolic materials under magnetic quantization

We present a simplified theoretical formulation of the thermoelectric power (TP) under magnetic quantization in quantum wells (QWs) of nonlinear optical materials on the basis of a newly formulated magneto-dispersion law. We consider the anisotropies in the effective electron masses and the spin-orbit constants within the framework of k.p formalism by incorporating the influence of the crystal field splitting. The corresponding results for III-V materials form a special case of our generalized analysis under certain limiting conditions. The TP in QWs of Bismuth, II-VI, IV-VI and stressed materials has been studied by formulating appropriate electron magneto-dispersion laws. We also address the fact that the TP exhibits composite oscillations with a varying quantizing magnetic field in QWs of n-Cd3As2, n-CdGeAs2, n-InSb, p-CdS, stressed InSb, PbTe and Bismuth. This reflects the combined signatures of magnetic and spatial quantizations of the carriers in such structures. The TP also decreases with increasing electron statistics and under the condition of non-degeneracy, all the results as derived in this paper get transformed into the well-known classical equation of TP and thus confirming the compatibility test. We have also suggested an experimental method of determining the elastic constants in such systems with arbitrary carrier energy spectra from the known value of the TP. (C) 2010 Elsevier Ltd. All rights reserved.

A Bilinear Constitutive Model for Isotropic Bimodulus Materials

Proper formulation of stress-strain relations, particularly in tension-compression situations for isotropic biomodulus materials, is an unresolved problem. Ambartsumyan's model [8] and Jones' weighted compliance matrix model [9] do not satisfy the principle of coordinate invariance. Shapiro's first stress invariant model [10] is too simple a model to describe the behavior of real materials. In fact, Rigbi [13] has raised a question about the compatibility of bimodularity with isotropy in a solid. Medri [2] has opined that linear principal strain-principal stress relations are fictitious, and warned that the bilinear approximation of uniaxial stress-strain behavior leads to ill-working bimodulus material model under combined loading. In the present work, a general bilinear constitutive model has been presented and described in biaxial principal stress plane with zonewise linear principal strain-principal stress relations. Elastic coefficients in the model are characterized based on the signs of (i) principal stresses, (ii) principal strains, and (iii) on the value of strain energy component ratio ER greater than or less than unity. The last criterion is used in tension-compression and compression-tension situations to account for different shear moduli in pure shear stress and pure shear strain states as well as unequal cross compliances.

On numerical implementation of an isotropic elastic-viscoplastic constitutive model for bulk metallic glasses

An efficient algorithm within the finite deformation framework is developed for finite element implementation of a recently proposed isotropic, Mohr-Coulomb type material model, which captures the elastic-viscoplastic, pressure sensitive and plastically dilatant response of bulk metallic glasses. The constitutive equations are first reformulated and implemented using an implicit numerical integration procedure based on the backward Euler method. The resulting system of nonlinear algebraic equations is solved by the Newton-Raphson procedure. This is achieved by developing the principal space return mapping technique for the present model which involves simultaneous shearing and dilatation on multiple potential slip systems. The complete stress update algorithm is presented and the expressions for viscoplastic consistent tangent moduli are derived. The stress update scheme and the viscoplastic consistent tangent are implemented in the commercial finite element code ABAQUS/Standard. The accuracy and performance of the numerical implementation are verified by considering several benchmark examples, which includes a simulation of multiple shear bands in a 3D prismatic bar under uniaxial compression.

Interplay of spin and lattice degrees of freedom in the frustrated antiferromagnet CdCr2O4: High-field and temperature-induced anomalies of the elastic constants

Temperature and magnetic field studies of the elastic constants of the chromium spinel CdCr2O4 show pronounced anomalies related to strong spin-phonon coupling in this frustrated antiferromagnet. A detailed comparison of the longitudinal acoustic mode propagating along the 111] direction with a theory based on an exchange-striction mechanism leads to an estimate of the strength of the magnetoelastic interaction. The derived spin-phonon coupling constant is in good agreement with previous determinations based on infrared absorption. Further insight is gained from intermediate and high magnetic field experiments in the field regime of the magnetization plateau. The role of the antisymmetric Dzyaloshinskii-Moriya interaction is discussed.

Symmetry-breaking transitions in equilibrium shapes of coherent precipitates: Effect of elastic anisotropy and inhomogeneity

We examine the symmetry-breaking transitions in equilibrium shapes of coherent precipitates in two-dimensional (2-D) systems under a plane-strain condition with the principal misfit strain components epsilon(xx)*. and epsilon(yy)*. For systems with cubic elastic moduli, we first show all the shape transitions associated with different values of t = epsilon(yy)*/epsilon(xx)*. We also characterize each of these transitions, by studying its dependence on elastic anisotropy and inhomogeneity. For systems with dilatational misfit (t = 1) and those with pure shear misfit (t = -1), the transition is from an equiaxed shape to an elongated shape, resulting in a break in rotational symmetry. For systems with nondilatational misfit (-1 < t < 1; t not equal 0), the transition involves a break in mirror symmetries normal to the x- and y-axes. The transition is continuous in all cases, except when 0 < t < 1. For systems which allow an invariant line (-1 less than or equal to t < 0), the critical size increases with an increase in the particle stiffness. However, for systems which do not allow an invariant line (0 < t less than or equal to 1), the critical size first decreases, reaches a minimum, and then starts increasing with increasing particle stiffness; moreover, the transition is also forbidden when the particle stiffness is greater than a critical value.

A new beam finite element for the analysis of functionally graded materials

A new beam element is developed to study the thermoelastic behavior of functionally graded beam structures. The element is based on the first-order shear deformation theory and it accounts for varying elastic and thermal properties along its thickness. The exact solution of static part of the governing differential equations is used to construct interpolating polynomials for the element formulation. Consequently, the stiffness matrix has super-convergent property and the element is free of shear locking. Both exponential and power-law variations of material property distribution are used to examine different stress variations. Static, free vibration and wave propagation problems are considered to highlight the behavioral difference of functionally graded material beam with pure metal or pure ceramic beams. (C) 2003 Elsevier Science Ltd. All rights reserved.

Fully Comprehensive Geometrically Non-Linear Dynamic Analysis of Multi-Body Beam Systems with Elastic Couplings

This paper is concerned with the dynamic analysis of flexible,non-linear multi-body beam systems. The focus is on problems where the strains within each elastic body (beam) remain small. Based on geometrically non-linear elasticity theory, the non-linear 3-D beam problem splits into either a linear or non-linear 2-D analysis of the beam cross-section and a non-linear 1-D analysis along the beam reference line. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction,results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis,the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here.The analysis methodology can be viewed as a 3-step procedure. First, the sectional properties of beams made of composite materials are determined either based on an asymptotic procedure that involves a 2-D finite element nonlinear analysis of the beam cross-section to capture trapeze effect or using strip-like beam analysis, starting from Classical Laminated Shell Theory (CLST). Second, the dynamic response of non-linear, flexible multi-body beam systems is simulated within the framework of energy-preserving and energy-decaying time integration schemes that provide unconditional stability for non-linear beam systems. Finally,local 3-D responses in the beams are recovered, based on the 1-D responses predicted in the second step. Numerical examples are presented and results from this analysis are compared with those available in the literature.

SCALE EFFECTS ON ULTRASONIC WAVE DISPERSION CHARACTERISTICS OF MONOLAYER GRAPHENE EMBEDDED IN AN ELASTIC MEDIUM

Ultrasonic wave propagation in a graphene sheet, which is embedded in an elastic medium, is studied using nonlocal elasticity theory incorporating small-scale effects. The graphene sheet is modeled as an one-atom thick isotropic plate and the elastic medium/substrate is modeled as distributed springs. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. After that, an ultrasonic type of wave propagation model is also derived. The explicit expressions for the cut-off frequencies are also obtained as functions of the nonlocal scaling parameter and the y-directional wavenumber. Local elasticity shows that the wave will propagate even at higher frequencies. But nonlocal elasticity predicts that the waves can propagate only up to certain frequencies (called escape frequencies), after which the wave velocity becomes zero. The results also show that the escape frequencies are purely a function of the nonlocal scaling parameter. The effect of the elastic medium is captured in the wave dispersion analysis and this analysis is explained with respect to both local and nonlocal elasticity. The simulations show that the elastic medium affects only the flexural wave mode in the graphene sheet. The presence of the elastic matrix increases the band gap of the flexural mode. The present results can provide useful guidance for the design of next-generation nanodevices in which graphene-based composites act as a major element.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Odd-even effect in the elastic modulii of alpha,omega-Alkanedicarboxylic acids

Nanoindentation studies on alpha,omega-alkanedicarboxylic acids reveal that the elastic modulus, E, shows an odd-even alternation in exactly the same manner as the melting temperature, T-m. The results are consistent with the hypothesis that the strained molecular conformations in the odd diacids are the reasons for these alternations in T-m. The same packing features that lower T-m in the odd acids lead to easy accommodation of the deformation during nanoindentation and hence their low E.