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.

Anisotropic, gradient and metal-like mechanical behaviors of teeth and their implications on tooth functions

The anisotropy and gradient of the elastic modulus and the hardness of teeth were investigated by means of instrumented indentation method. Such properties are attributed to the unique microstructures of teeth based on scanning electron microscopic analysis. By comparing the relationship between the ratio of hardness to the reduced elastic modulus and the ratio of elastic unloading work to the total work of teeth in course of indentation to those of other materials, we found that the material behaviors of teeth display metal-like characteristics rather than ceramics as considered traditionally. These material behaviors and relevant functions are discussed briefly.

The size-dependent bending elastic properties of nanobeams with surface effects

A theoretical model is presented to investigate the size-dependent bending elastic properties of a nanobeam with the influence of the surface relaxation and the surface tension taken into consideration. The surface layer and its thickness of a nanostructure are defined unambiguously. A three-dimensional (3D) crystal model for a nanofilm with n layers of relaxed atoms is investigated. The four nonzero elastic constants of the nanofilm are derived, and then the Young's modulus for simple tension is obtained. Using the relation of energy equilibrium, the size-dependent effective elastic modulus and effective flexural rigidity of a nanobeam with two kinds of cross sections are derived, and their dependence on the surface relaxation and the surface tension is analysed.

Validation of component assembly model and extension to plasticity

Material potential energy is well approximated by '' pair-functional '' potentials. During calculating potential energy, the orientational and volumetric components have been derived from pair potentials and embedding energy, respectively. Slip results in plastic deformation, and slip component has been proposed accordingly. Material is treated as a component assembly, and its elastic, plastic and damage properties are reflected by different components respectively. Material constitutive relations are formed by means of assembling these three kinds of components. Anisotropy has been incorporated intrinsically via the concept of component. Theoretical and numerical results indicate that this method has the capacity of reproducing some results satisfactorily, with the advantages of physical explicitness, etc. (c) 2007 Elsevier Ltd. All rights reserved.

Simple phenomenological determination of contact stiffness and elastic modulus of Ce-based bulk metallic glasses through nanoindentation

The viscoelastic deformation of Ce-based bulk metallic glasses (BMGs) with low glass transition temperature is investigated at room temperature. Contact stiffness and elastic modulus of Ce-based BMGs cannot be derived using the conventional Oliver-Pharr method [W. C. Oliver and G. M. Pharr, J. Mater. Res. 7, 1564 (1992)]. The present work shows that the time dependent displacement of unloading segments can be described well by a generalized Kelvin model. Thus, a modified Oliver-Pharr method is proposed to evaluate the contact stiffness and elastic modulus, which does, in fact, reproduce the values obtained via uniaxial compression tests. (c) 2007 American Institute of Physics.

On the initial unloading slope in indentation of elastic-plastic solids by an indenter with an axisymmetric smooth profile

A simple relationship between the initial unloading slope, the contact area, and the elastic modulus is derived for indentation in elastic-plastic solids by an indenter with an arbitrary axisymmetric smooth profile. Although the same expression was known to hold for elastic solids, the new derivation shows that it is also true for elastic-plastic solids with or without work hardening and residual stress. These results should provide a sound basis for the use of the relationship for mechanical property determination using indentation techniques. (C) 1997 American Institute of Physics.

Dynamic analysis of arrest of buckle propagation on a beam on a nonlinear elastic foundation by FEM

Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.

Numerical analyses for the material bifurcation in plane sheet under tension

The various patterns (shear banding, surface wrinkling and necking) of material bifurcation in plane sheet under tension are investigated in this paper by means of a numerical method. It is found that numerical analysis can provide better ground for searching for the lowest critical loads. The inhomogeneity caused by void damage and the nonuniformity in the stress distribution across sheet thickness are proved to have detrimental effects on the material bifurcation. Nevertheless, material stability can be promoted by any means of depressing void damage or alleviating stress, even locally across the thickness. Besides, the peculiar behaviour of material bifurcation under slight biaxiality state is demonstrated. Copyright (C) 1996 Elsevier Science Ltd

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals .1. Theoretical developments

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals is proposed. The anisotropic hardening of single crystals, at first, is discussed with the viewpoint of yield surface and a new formulation of it is proposed. Then, a model for the anisotropic hardening of polycrystals is suggested by increasing the number of slip systems and incorporating the interaction of all slip systems. The interaction of grains through grain boundaries is shown to be similar to, and incorporated into, the interaction of slip systems in grains. The numerical predictions and their comparisons with experiments will follow in Part II of this paper.

A new expression of hardening coefficients for fcc-crystal and calibration of the material constants

In order to describe the effect of latent hardening on the macro-plastic behavior of foc-crystal, a new expression for hardening coefficient is proposed in which there are 12 material constants, each having clear physical meaning. And a method of material constant calibration is suggested and used to determine the material constants of copper and aluminum crystal. The simulated load-elongation curves along various crystallographic orientations are comparable with the experimental ones.

Dislocation nucleation and emission from crack-tip

The problems of dislocation nucleation and emission from a crack tip are analysed based on Peierls model. The concept adopted here is essentially the same as that proposed by Rice. A slight modification is introduced here to identify the pure linear elastic response of material. A set of new governing equations is developed, which is different from that used by Beltz and Rice. The stress field and the dislocation density field can be expressed as the first and second Chebyshev polynomial series respectively. Then the opening and slip displacements can be expanded as the trigonometric series. The Newton-Raphson Method is used to solve a set of nonlinear algebraic equations. The new governing equations allow us to extend the analyses to the case of dislocation emission. The calculation results for pure shearing, pure tension and combined tension and shear loading are given in detail.

Unloading characteristics of anti-plane shear crack in softening materials

The local characteristics of the anti-plane shear stress and strain field are determined for a material where the stress increases linearly with strain up to a limit and then softens nonlinearly. Two unloading models are considered such that the unloading path always returns to the origin while the other assumes the unloading modulus to be that of the initial shear modulus. As the applied shear increases, an unloading zone is found to prevail between a zone in which the material softens and another zone in which the material is linear-elastic even though the crack does not propagate. The divisions of these zones are displayed graphically.

The inhomogeneity of plastic-deformation in ductile single-crystals

A discrete slip model which characterizes the inhomogeneity of material properties in ductile single crystals is proposed in this paper. Based on this model rate-dependent finite element investigations are carried out which consider the finite deformation, finite rotation, latent hardening effect and elastic anisotropy. The calculation clearly exhibits the process from microscopic inhomogeneous and localized deformation to necking and the formation of LSBS and reveals several important features of shear localization. For example, the inhomogeneous deformation is influenced by the imperfections and initial non-uniformities of material properties. The inhomogeneous deformation may either induce necking which results in the lattice rotation and leads to geometrical softening, which in turn promotes the formation of CSBS, or induces heavily localized deformation. The microscopic localized deformation eventually develops into the LSBS and results in a failure. These results are in close agreement with experiment. Our calculations also find that the slip lines on the specimen's surface at necking become curved and also find that if the necking occurs before the formation of LSBS, this band must be misoriented from the operative slip systems. In this case, the formation of LSBS must involve non-crystallographic effects. These can also be indirectly confirmed by experiment. All these suggest that our present discrete slip model offers a correct description of the inhomogeneous deformation characterization in ductile crystals.

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.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.