A Numerical Study of Indentation Using Indenters of Different Geometry

Finite element simulation of the Berkovich, Vickers, Knoop, and cone indenters was carried out for the indentation of elastic-plastic material. To fix the semiapex angle of the cone, several rules of equivalence were used and examined. Despite the asymmetry and differences in the stress and strain fields, it was established that for the Berkovich and Vickers indenters, the load-displacement relation can closely be simulated by a single cone indenter having a semiapex angle equal to 70.3degrees in accordance with the rule of the volume equivalence. On the other hand, none of the rules is applicable to the Knoop indenter owing to its great asymmetry. The finite element method developed here is also applicable to layered or gradient materials with slight modifications.

A Dynamic Comprehensive Method for Landslide Control

A slope failure is developed due to progressive external loads and deteriorations of slope geomaterials, thus forming a progressive and dynamic development and occurrence of landslides. Site geological properties and other active factors such as hydrodynamic load and human activities are complex and usually unknown, thus this dynamic development and occurrence of landslides can only be understood through the progressive accumulation of knowledge on the landslides. For such a progressive process, this paper proposes a dynamic comprehensive control method for landslide control. This control method takes full advantage of updated monitoring data and site investigations of landslides, and emphasizes the implementation of possible measures for landslide control at reasonable stages and in different groups. These measures are to prevent the occurrence of a landslide disaster. As a case study, a landslide project at the Panluo open-pit iron mine is analyzed to illustrate this dynamic comprehensive control method.

Thermoelastic Stresses in SiC Single Crystals Grown by the Physical Vapor Transport Method

A finite element-based thermoelastic anisotropic stress model for hexagonal silicon carbide polytype is developed for the calculation of thermal stresses in SiC crystals grown by the physical vapor transport method. The composite structure of the growing SiC crystal and graphite lid is considered in the model. The thermal expansion match between the crucible lid and SiC crystal is studied for the first time. The influence of thermal stress on the dislocation density and crystal quality is discussed.

An Effective Method of Determining the Residual Stress Gradients in a Micro-Cantilever

A method of determining the micro-cantilever residual stress gradients by studying its deflection and curvature is presented. The stress gradients contribute to both axial load and bending moment, which, in prebuckling regime, cause the structural stiffness change and curving up/down, respectively. As the axial load corresponds to the even polynomial terms of stress gradients and bending moment corresponds to the odd polynomial terms, the deflection itself is not enough to determine the axial load and bending moment. Curvature together with the deflection can uniquely determine these two parameters. Both linear analysis and nonlinear analysis of micro-cantilever deflection under axial load and bending moment are presented. Because of the stiffening effect due to the nonlinearity of (large) deformation, the difference between linear and nonlinear analyses enlarges as the micro-cantilever deflection increases. The model developed in this paper determines the resultant axial load and bending moment due to the stress gradients. Under proper assumptions, the stress gradients profile is obtained through the resultant axial load and bending moment.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

Finite Element Analysis of Deformed Legs of Offshore Platform Structures

The element stiffness matrix of the equivalent beam or pipe element of the deformed leg of the platform is derived by the finite element method. The stresses and displacements of some damaged components are calculated, and the numeri-cal solutions agree well with those obtained by the fine mesh finite element method. Finally, as an application of this method, the stresses of some platform structures are calculated and analyzed.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

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.

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.

Experimental investigation of a control method for SI-HCCI-SI transition in a multi-cylinder gasoline engine

In HCCI engines, the Air/Fuel Ratio (AFR) and Residual Gas Fraction (RGF) are difficult to control during the SI-HCCI-SI transition, and this may result in incomplete combustion and/or high pressure raise rates. As a result, there may be undesirably high engine load fluctuations. The objectives of this work are to further understand this process and develop control methods to minimize these load fluctuations. This paper presents data on instantaneous AFR and RGF measurements, both taken by novel experimental techniques. The data provides an insight into the cyclic AFR and RGF fluctuations during the switch. These results suggest that the relatively slow change in the intake Manifold Air Pressure (MAP) and actuation time of the Variable Valve Timing (VVT) are the main causes of undesired AFR and RGF fluctuations, and hence an unacceptable Net IMEP (NIMEP) fluctuation. We also found large cylinder-to-cylinder AFR variations during the transition. Therefore, besides throttle opening control and VVT shifting, cyclic and individual cylinder fuel injection control is necessary to achieve a smooth transition. The control method was developed and implemented in a test engine, and the result was a considerably reduced NIMEP fluctuation during the mode switch. The instantaneous AFR and RGF measurements could furthermore be adopted to develop more sophisticated control methods for SI-HCCI-SI transitions. © 2010 SAE International.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Asymptotic analysis based on K-S constitutive law for a rubber wedge compressed by a line load at its tip

In the present paper, a rubber wedge compressed by a line load at its tip is asymptotically analyzed using a special constitutive law proposed by Knowles and Sternberg (K-S elastic law) [J. Elasticity 3 (1973) 67]. The method of dividing sectors proposed by Gao [Theoret. Appl. Fract, Mech. 14 (1990) 219] is used. Domain near the wedge tip can be divided into one expanding sector and two narrowing sectors. Asymptotic equations of the strain-stress field near the wedge tip are derived and solved numerically. The deformation pattern near a wedge tip is completely revealed. A special case. i.e. a half space compressed by a line load is solved while the wedge angle is pi.

A general design method for electric machines using magnetic circuit model considering the flux saturation problem

This paper proposes an analytical approach that is generalized for the design of various types of electric machines based on a physical magnetic circuit model. Conventional approaches have been used to predict the behavior of electric machines but have limitations in accurate flux saturation analysis and hence machine dimensioning at the initial design stage. In particular, magnetic saturation is generally ignored or compensated by correction factors in simplified models since it is difficult to determine the flux in each stator tooth for machines with any slot-pole combinations. In this paper, the flux produced by stator winding currents can be calculated accurately and rapidly for each stator tooth using the developed model, taking saturation into account. This aids machine dimensioning without the need for a computationally expensive finite element analysis (FEA). A 48-slot machine operated in induction and doubly-fed modes is used to demonstrate the proposed model. FEA is employed for verification.

The Newest Developments Of Load-Unload Response Ratio (Lurr)

The Load Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It is inspiring that its predictions using LURR have been improving. Since 2004 we have made a major breakthrough in intermediate-term earthquake forecasting of the strong earthquakes on the Chinese mainland using LURR and successfully predicted the Pakistan earthquake with magnitude M 7.6 on October 8, 2005. The causes for improving the prediction in terms of LURR have been discussed in the present paper.