Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

Synchronous Multi-Scale Observations on Rock Damage and Rupture

This paper reports a multi-scale study on damage evolution process and rupture of gabbro under uniaxial compression with several experimental techniques, including MTS810 testing machine, white digital speckle correlation method, and acoustic emission technique. In particular, the synchronization of the three experimental systems is realized for the study of relationship of deformation and damage at multiple scales. It is found that there are significant correlation between damage evolution at small and large length scales, and rupture at sample scale, especially it displays critical sensitivity at multiple scales and trans-scale fluctuations.

Modeling Study on the Characteristics of Laminar and Turbulent Argon Plasma Jets Impinging Normally upon a Flat Plate in Ambient Air

Modeling study is performed to compare the flow and heat transfer characteristics of laminar and turbulent argon thermal-plasma jets impinging normally upon a flat plate in ambient air. The combined-diffusion-coefficient method and the turbulence-enhanced combined-diffusion-coefficient method are employed to treat the diffusion of argon in the argon-air mixture for the laminar and the turbulent cases, respectively. Modeling results presented include the flow, temperature and argon concentration fields, the air mass flow-rates entrained into the impinging plasma jets, and the distributions of the heat flux density on the plate surface. It is found that the formation of a radial wall jet on the plate surface appreciably enhances the mass flow rate of the ambient air entrained into the laminar or turbulent plasma impinging-jet. When the plate standoff distance is comparatively small, there exists a significant difference between the laminar and turbulent plasma impinging-jets in their flow fields due to the occurrence of a large closed recirculation vortex in the turbulent plasma impinging-jet, and no appreciable difference is found between the two types of jets in their maximum values and distributions of the heat flux density at the plate surface. At larger plate standoff distances, the effect of the plate on the jet flow fields only appears in the region near the plate, and the axial decaying-rates of the plasma temperature, axial velocity and argon mass fraction along the axis of the laminar plasma impinging-jet become appreciably less than their turbulent counterparts.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Theory and Application of the Transient Injection Well Test

Based on the theory of the pumping well test, the transient injection well test was suggested in this paper. The design method and the scope of application are discussed in detail. The mathematical models are developed for the short-time and long-time transient injection test respectively. A double logarithm type curve matching method was introduced for analyzing the field transient injection test data. A set of methods for the transient injection test design, experiment performance and data analysis were established. Some field tests were analyzed, and the results show that the test model and method are suitable for the transient injection test and can be used to deal with the real engineering problems.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

Study of Growth Mechanism of Lysozyme Crystal by Batch Crystallization Method

The lysozyme crystals were made by batch crystallization method and the distribution of aggregate in solution were measured by dynamic light scattering. The results showed that the dimension of aggregate increased with the increase of the concentration of lysozyme and NaCl, lysozyme molecules aggregated gradually in solution and finally arrived at balance each other. The higher the concentrations of lysozyme and NaCl were, the faster the growth rate of (I 10) face was. The growth rates of lysozyme crystal were obtained by a Zeiss microscope, and the effective surface energy (a) of growing steps were calculated about 4.01 X 10(-8) J.cm(-2) according to the model of multiple two-dimensional nucleation mechanism.

A new method for the acquisition of ultrasonic strain image volumes.

This article presents a new method for acquiring three-dimensional (3-D) volumes of ultrasonic axial strain data. The method uses a mechanically-swept probe to sweep out a single volume while applying a continuously varying axial compression. Acquisition of a volume takes 15-20 s. A strain volume is then calculated by comparing frame pairs throughout the sequence. The method uses strain quality estimates to automatically pick out high quality frame pairs, and so does not require careful control of the axial compression. In a series of in vitro and in vivo experiments, we quantify the image quality of the new method and also assess its ease of use. Results are compared with those for the current best alternative, which calculates strain between two complete volumes. The volume pair approach can produce high quality data, but skillful scanning is required to acquire two volumes with appropriate relative strain. In the new method, the automatic quality-weighted selection of image pairs overcomes this difficulty and the method produces superior quality images with a relatively relaxed scanning technique.

Experiments on interaction between current-induced vibration and scour of submarine pipelines on sandy bottom

In order to understand the dynamic behavior of submarine pipelines exposed to current and the mechanism of the interaction between current-induced vibration and scour of pipelines on a sandy bottom, an experimental investigation is conducted with a small scale model A test model which can be tested in the flume is set up by taking into account the typical working conditions of the pipelines and by applying the similarity theory. The interactions between the shape of the scour hole and the behavior of the pipeline as well as the flow patterns of the current are detailed, and the interaction mechanism outlined. The effect of vibration of the pipeline on the development of dynamic scour at different stages is found out. The proposed experimental method and test results provide an effective means for design of marine pipelines against scouring.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.

Hill Instability Analysis Of Tlp Tether Subjected To Combined Platform Surge And Heave Motions

This paper presents the Hill instability analysis of Tension Leg Platform (TLP) tether it, deep sea. The 2-D nonlinear beam model which is Undergoing Coupled axial and transverse vibrations, is applied. The governing equations are reduced to nonlinear Hill equation by use of the Galerkin's method and the modes superposition principle. The Hill instability charted Lip to large parameters is obtained. An important parameter M is defined and can he expressed as the functions of tether length, the platform surge and heave motion amplitudes. Some example studies are performed for various environmental conditions. The results demonstrate that the nonlinear coupling between the axial and transverse vibrations has a significant effect on the response of structure.. It needs to be considered for the accurate dynamic analysis of long TLP tether subjected to the combined platform surge and heave motions.

Comparison Between Lurr And State Vector Analysis Before Strong Earthquakes In Southern California Since 1980

There are seven strong earthquakes with M >= 6.5 that occurred in southern California during the period from 1980 to 2005. In this paper, these earthquakes were studied by the LURR (Load/Unload Response Ratio) method and the State Vector method to detect if there are anomalies before them. The results show that LURR anomalies appeared before 6 earthquakes out of 7 and State Vector anomalies appeared before all 7 earthquakes. For the LURR method, the interval between maximum LURR value and the forthcoming earthquake is 1 to 19 months, and the dominant mean interval is about 10.7 months. For the State Vector method, the interval between the maximum modulus of increment State Vector and the forthcoming earthquake is from 3 to 27 months, but the dominant mean interval between the occurrence time of the maximum State Vector anomaly and the forthcoming earthquake is about 4.7 months. The results also show that the minimum valid space window scale for the LURR and the State Vector is a circle with a radius of 100 km and a square of 3 degrees 3 degrees, respectively. These results imply that the State Vector method is more effective for short-term earthquake prediction than the LURR method, however the LURR method is more effective for location prediction than the State Vector method.