Scaling approach to conical indentation in elastic-plastic solids with work hardening

We derive, using dimensional analysis and finite element calculations, several scaling relationships for conical indentation in elastic-plastic solids with work hardening. Using these scaling relationships, we examine the relationships between hardness, contact area, initial unloading slope, and mechanical properties of solids. The scaling relationships also provide new insights into the shape of indentation curves and form the basis for understanding indentation measurements, including nano- and micro-indentation techniques. They may also be helpful as a guide to numerical and finite element calculations of indentation problems.

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.

Analysis of indentation loading curves obtained using conical indenters

Using dimensional analysis and finite-element calculations we determine the functional form of indentation loading curves for a rigid conical indenter indenting into elastic-perfectly plastic solids. The new results are compared with the existing theories of indentation using conical indenters, including the slip-line theory for rigid-plastic solids, Sneddon's result for elastic solids, and Johnson's model for elastic-perfectly plastic solids. In the limit of small ratio of yield strength (Y) to Young's modulus (E), both the new results and Johnson's model approach that predicted by slip-line theory for rigid-plastic solids. In the limit of large Y/E, the new results agree with that for elastic solids. For a wide range of Y/E, some difference is found between Johnson's model-and the present result. This study also demonstrates the possibilities and limitations of using indentation loading curves to extract fundamental mechanical properties of solids.

Relationships Between Initial Unloading Slope, Contact Depth, and Mechanical Properties for Conical Indentation in Linear Viscoelastic Solids

Using analytical and finite element modeling, we examine the relationships between initial unloading slope, contact depth, and mechanical properties for spherical indentation in viscoelastic solids with either displacement or load as the independent variable. We then investigate whether the Oliver-Pharr method for determining the contact depth and contact radius, originally proposed for indentation in elastic and elastic-plastic solids, is applicable to spherical indentation in viscoelastic solids. Finally, the analytical and numerical results are used to answer questions raised in recent literature about measuring viscoelastic properties from instrumented spherical indentation experiments.

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.

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total 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. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for 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. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

Application of the strain energy density criterion to ductile fracture

The strain energy density criterion due to Sih is used to predict fracture loads of two thin plates subjected to large elastic-plastic deformation. The prediction is achieved with a finite element analysis which is based on Hill's variational principle for incremental deformations capable of solving gross yielding problems involving arbitrary amounts of deformation. The computed results are in excellent agreement with those obtained in Sih's earlier analysis and with an experimental observation.

Ductile crack propagation with the strain-energy density criterion

The strain energy density criterion is used to characterize subcritical crack growth in a thin aluminum alloy sheet undergoing general yielding. A finite element analysis which incorporates both material and geometrical nonlinear behaviors of the cracked sheets is developed to predict fracture loads at varying crack growth increments. The predicted results are in excellent agreement with those measured experimentally, thus confirming the validity of the strain energy density criterion for characterizing ductile crack propagation.

Transient Temperature Distribution in a Sheet Metal by Low-Energy Laser Scanning

Low-energy laser-heating techniques are widely used in engineering applications such as, thinfilm deposition, surface treatment, metal forming and micro-structural pattern formation. In this paper,under the conditions of ignoring the thermo-mechanical coupling, a numerical simulation on the spatialand temporal temperature distribution in a sheet metal produced by the laser beam scanning in virtue of thefinite element method is presented. Both the three-dimensional transient temperature field and thetemperature evolution as a function of heat penetrating depth in the metal sheet are calculated. Thetemperature dependence of material properties was taken into account. It was shown that, after taking thetemperature dependence of the material absorbance effect into consideration, the temperature change ratealong the scanning direction and the temperature maximum were both increased.

Adiabatic shear localization: occurrence, theories, and applications

激光二极管端面抽运Nd∶YVO_4板条激光器及其热效应

激光器中激光介质采用板条状几何结构可以极大地降低它的热效应,但仍然需要进一步分析其影响,进而优化激光器效率。利用有限元分析方法分析了部分端面抽运的混合腔板条激光器中激光介质的热效应,计算的热透镜焦距与实测结果基本相符。分析了热效应对模式匹配的影响,分析结果对于优化激光器效率、改进谐振腔设计具有一定的参考价值。并在分析的基础上进行了混合腔实验,抽运功率为110 W时,获得连续输出激光功率41.5 W,光-光转换效率约38%,斜效率达58.8%,M2因子为非稳腔方向M2x=1.59,稳定腔方向M2y=1.55。

Novel method for design of surface relief guided-mode resonant gratings at normal incidence

We demonstrate that the surface relief guided-mode resonant gratings with specified central wavelength and FWHM in the visible wavelength range can be designed by analyzing the complex poles of Reflectance and transmission coefficient matrix algorithm (RTCM), a variant of S-matrix propagation algorithm proposed for calculation of multilayer gratings. In addition, FWHM is computed with couple-mode (CM) theory of resonant gratings which is firstly extended by Norton et al. in calculation of waveguide grating. Furthermore, the side band reflections of the filter can be reduced to less than 5% in the visible wavelength with the antireflection (AR) design technique widely used in the thin-film field. (C) 2008 Elsevier B.V. All rights reserved.