Novel Mechanical Behavior Of Zno Nanorods

A novel stress-strain relation with two stages of linear elastic deformation is observed in [0 0 0 1]-oriented ZnO nanorods under uniaxial tensile loading. This phenomenon results from a phase transformation from wurtzite (WZ, P6(3)mc space group) to a body-centered tetragonal structure with four-atom rings (denoted as BCT-4) belonging to the P4(2)/mnm space group. The analysis here focuses on the effects of nanorod size and temperature on the phase transformation and the associated mechanical behavior. It is found that as size is increased from 19.5 to 45.5 angstrom, the critical stress for nucleation of the transformation decreases by 25% from 21.90 to 16.50 GPa and the elastic moduli of the WZ- and BCT-4-structured nanorods decrease by 24% (from 299.49 to 227.51 GPa) and 38% (from 269.29 to 166.86 GPa), respectively. A significant temperature effect is also observed, with the critical stress for transformation initiation decreasing 87.8% from 17.89 to 2.19 GPa as temperature increases from 300 to 1500 K. (c) 2007 Elsevier B.V. All rights reserved.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

Micromechanics Modeling Of Strength For Nanocrystalline Copper

We present a model in this paper for predicting the inverse Hall-Petch phenomenon in nanocrystalline (NC) materials which are assumed to consist of two phases: grain phase of spherical or spheroidal shapes and grain boundary phase. The deformation of the grain phase has an elasto-viscoplastic behavior, which includes dislocation glide mechanism, Coble creep and Nabarro-Herring creep. However the deformation of grain boundary phase is assumed to be the mechanism of grain boundary diffusion. A Hill self-consistent method is used to describe the behavior of nanocrystalline pure copper subjected to uniaxial tension. Finally, the effects of grain size and its distribution, grain shape and strain rate on the yield strength and stress-strain curve of the pure copper are investigated. The obtained results are compared with relevant experimental data in the literature.

The Characterization Of Creep And Time-Dependent Properties Of Bulk Metallic Glasses Using Nanoindentation

Viscoelastic deformation and creep behavior of La- and Ce-based bulk metallic glasses (BMGs) with low glass transition temperature are investigated through nanoindentation at room temperature. Creep compliance and retardation spectra are derived to study the creep mechanism. The time-dependent displacement can be well described by a generalized Kelvin model. A modification is proposed to determine the elastic modulus from the generalized Kelvin model. The results are in excellent agreement with the elastic modulus determined by uniaxial compression tests. (c) 2007 Published by Elsevier B.V.

Deformation Behavior Of Fe-Based Bulk Metallic Glass During Nanoindentation

Fe-based bulk metallic glasses (BMGs) normally exhibit super high strength but significant brittleness at ambient temperature. Therefore, it is difficult to investigate the plastic deformation behavior and mechanism in these alloys through conventional tensile and compressive tests due to lack of distinct macroscopic plastic strain. In this work, the deformation behavior of Fe52Cr15Mo9Er3C15B6 BMG was investigated through instrumented nanoindentation and uniaxial compressive tests. The results show that serrated flow, the typical plastic deformation feature of BMGs, could not be found in as-cast and partially crystallized samples during nanoindentation. In addition, the deformation behavior and mechanical properties of the alloy are insensitive to the applied loading rate. The mechanism for the appearance of the peculiar deformation behavior in the Fe-based BMG is discussed in terms of the temporal and spatial characteristics of shear banding during nanoindentation.

Molecular Dynamics Study Of The Mechanical Behavior Of Nickel Nanowire: Strain Rate Effects

We present the analysis of uniaxial deformation of nickel nanowires using molecular dynamics simulations, and address the strain rate effects on mechanical responses and deformation behavior. The applied strain rate is ranging from 1 x 10(8) s(-1) to 1.4 x 10(11) s(-1). The results show that two critical strain rates, i.e., 5 x 10(9) s(-1) and 8 x 10(10) s(-1), are observed to play a pivotal role in switching between plastic deformation modes. At strain rate below 5 x 10(9) s(-1), Ni nanowire maintains its crystalline structure with neck occurring at the end of loading, and the plastic deformation is characterized by {111} slippages associated with Shockley partial dislocations and rearrangements of atoms close to necking region. At strain rate above 8x10(10) s(-1), Ni nanowire transforms from a fcc crystal into a completely amorphous state once beyond the yield point, and hereafter it deforms uniformly without obvious necking until the end of simulation. For strain rate between 5 x 10(9) s(-1) and 8 x 10(10) s(-1), only part of the nanowire exhibits amorphous state after yielding while the other part remains crystalline state. Both the {111} slippages in ordered region and homogenous deformation in amorphous region contribute to the plastic deformation. (C) 2007 Published by Elsevier B.V.

Effect Of Load Triaxiality On Polymorphic Transitions In Zinc Oxide

Molecular dynamics (MD) simulations and first-principles calculations are carried out to analyze the stability of both newly discovered and previously known phases of ZnO under loading of various triaxialities. The analysis focuses on a graphite-like phase (FIX) and a body-centered-tetragonal phase (BCT-4) that were observed recently in [0 1 (1) over bar 0]- and [0 0 0 1]-oriented nanowires respectively under uniaxial tensile loading as well as the natural state of wurtzite (WZ) and the rocksalt (RS) phase which exists under hydrostatic pressure loading. Equilibrium critical stresses for the transformations are obtained. The WZ -> HX transformation is found to be energetically favorable above a critical tensile stress of 10 GPa in [0 1 (1) over tilde 0] nanowires. The BCT-4 phase can be stabilized at tensile stresses above 7 GPa in [0 0 0 1] nanowires. The RS phase is stable at hydrostatic pressures above 8.2 GPa. The identification and characterization of these phase transformations reveal a more extensive polymorphism of ZnO than previously known. A crystalline structure-load triaxiality map is developed to summarize the new understanding. (c) 2007 Elsevier Ltd. All rights reserved.

Comparison Of The Quasi-Static Method And The Dynamic Method For Simulating Fracture Processes In Concrete

Concrete is heterogeneous and usually described as a three-phase material, where matrix, aggregate and interface are distinguished. To take this heterogeneity into consideration, the Generalized Beam (GB) lattice model is adopted. The GB lattice model is much more computationally efficient than the beam lattice model. Numerical procedures of both quasi-static method and dynamic method are developed to simulate fracture processes in uniaxial tensile tests conducted on a concrete panel. Cases of different loading rates are compared with the quasi-static case. It is found that the inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, an unrealistic result will be obtained if a fracture process including unstable cracking is simulated by the quasi-static procedure.

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.

Effect of dentin tubules on the mechanical properties of dentin - part III: Numerical, analysis

Imitating a real tooth and the periodontal supporting tissues, we have established a 2D finite element model and carried out a numerical analysis based on the inhomogeneous and anisotropic (IA) stress-strain relation and strength model of dentin proposed in the preceding Parts I and II, and the conventional homogeneous and isotropic (III) model, respectively. Quite a few cases of loadings for a non-defected and a defected tooth are considered. The numerical results show that the stress level predicted by the IA model is remarkably higher than that by the III model, revealing that the effect of the dentin tubules should be taken into a serious consideration from the viewpoint of biomechanics.

Regularity of strain distribution in short-fiber/whisker reinforced composites

Based on studies on the strain distribution in short-fiber/whisker reinforced metal matrix composites, a deformation characteristic parameter, lambda is defined as a ratio of root-mean-square strain of the reinforcers identically oriented to the macro-linear strain along the same direction. Quantitative relation between lambda and microstructure parameters of composites is obtained. By using lambda, the stiffness moduli of composites with arbitrary reinforcer orientation density function and under arbitrary loading condition are derived. The upper-bound and lower-bound of the present prediction are the same as those from the equal-strain theory and equal-stress theory, respectively. The present theory provides a physical explanation and theoretical base for the present commonly-used empirical formulae. Compared with the microscopic mechanical theories, the present theory is competent for stiffness modulus prediction of practical engineering composites in accuracy and simplicity.

Damage evaluation and damage localization of rock

Knowledge of damage accumulation and corresponding failure evolution are prerequisite for effective maintenance of civil engineering so as to avoid disaster. Based on statistical mesoscopic damage mechanics, it was revealed that there are three stages in the process of deformation, damage and failure of multiscale heterogeneous elastic-brittle medium. These are uniformly distributed damage, localized damage and catastrophic failure. In order to identify the transitions from scattering damage to macroscopically localized one, a condition for damage localization was given. The experiments of rock under uniaxial compression with the aid of observations of acoustic emission and speckle correlation do support the concept of localization. This provides a potential approach to properly evaluate damage accumulation in practice. In addition, it is found in the experiments that catastrophic failure displays critical sensitivity. This gives a helpful clue to the prediction of catastrophic failure. (C) 2004 Elsevier Ltd. All rights reserved.

Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.