Ab initio investigation of the elasticity and stability of aluminium

On the basis of the pseudopotential plane-wave (PP-PW) method in combination with the local density functional theory (LDFT), complete stress-strain curves for the uniaxial loading and uniaxial deformation along the [001] and [111] directions, and the biaxial proportional extension along [010] and [001] for aluminium are obtained. During the uniaxial loading, certain general behaviours of the energy versus the stretch and the load versus the stretch are confirmed; in each case, there exist three special unstressed structures: f.c.c., b.c.c., and f.c.t. for [001]; f.c.c., s.c., and b.c.c. for [111]. Using stability criteria, we find that all of these states are unstable, and always occur together with shear instability, except the natural f.c.c. structure. A Pain transformation from the stable f.c.c. structure to the stable b.c.c. configuration cannot be obtained by uniaxial compression along any equivalent [001] and [111] direction. The tensile strengths are similar for the two directions. For the higher energy barrier of the [111] direction, the compressive strength is greater than that for the [001] direction. With increase in the ratio of the biaxial proportional extension, the stress and tensile strength increase; however, the critical strain does not change significantly. Our results add to the existing ab initio database for use in fitting and testing interatomic potentials.

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.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

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.

Time-Frequency Analysis of Nonlinear Systems: The Skeleton Linear Model and the Skeleton Curves

The joint time-frequency analysis method is adopted to study the nonlinear behavior varying with the instantaneous response for a class of S.D.O.F nonlinear system. A time-frequency masking operator, together with the conception of effective time-frequency region of the asymptotic signal are defined here. Based on these mathematical foundations, a so-called skeleton linear model (SLM) is constructed which has similar nonlinear characteristics with the nonlinear system. Two skeleton curves are deduced which can indicate the stiffness and damping in the nonlinear system. The relationship between the SLM and the nonlinear system, both parameters and solutions, is clarified. Based on this work a new identification technique of nonlinear systems using the nonstationary vibration data will be proposed through time-frequency filtering technique and wavelet transform in the following paper.

Closed form solutions of T-stress in plane elasticity crack problems

The T-stress is considered as an important parameter in linear elastic fracture mechanics. In this paper, several closed form solutions of T-stress in plane elasticity crack problems in an infinite plate are investigated using the complex potential theory. In the line crack case, if the applied loading is the remote stress or the concentrated forces, the T-stress can be derived from the basic field. Here, the basic field is defined as the field caused by the applied loading in the infinite plate without the crack. For the circular are crack, the T-stress can be abstracted from a known solution. For the cusp crack problems, the T-stress can be separated from the obtained stress solution for which the conformal mapping technique is used.

Flow-Induced Vibrations Of A Cylinder With Two Degrees Of Freedom Near Rigid Plane Boundary

Based on similarity analyses, the flow-induced vibrations of a near-wall cylinder with 2 degrees of freedom are investigated experimentally by employing a hydroelastic apparatus in conjunction with a flume. The cylinder's vibration amplitude, vibration frequency and vortex shedding frequency were measured and analyzed. The effects of gap-to-diameter ratio (e,ID) upon the vibration responses are further investigated. The experimental results indicate that, when the reduced velocity (Vr) is small (e.g. Vr = 1.2 similar to 2.6), only streamwise vibration occurs, and its frequency is quite close to its natural frequency in still water. When increasing Vr (e.g. Vr > 3.4), both streamwise and transverse vibrations of the near-wall cylinder may occur. In the examined range of gap-to-diameter ratio (0.42 < e(0)/D < 2.68), 2 vibration stages (in terms of Vr) of streamwise vibrations usually exist: First Streamwise Vibration (FSV) and Second Streamwise Vibration (SSV). In the SSV stage, the vortex shedding frequency may either undergo a jump to that of the streamwise vibration, or stay consistent with that of the transverse vibration. The amplitudes of transverse vibration are usually much larger than those of streamwise vibration for the same value of e(0)/D. The maximum amplitudes of both streamwise and transverse vibration get larger with the increase of e(0)/D (0.42 < e(0)/D < 2.68).

Instability Of Plane Poiseuille Flow In A Fluid-Porous System

The instability of Poiseuille flow in a fluid-porous system is investigated. The system consists of a fluid layer overlying porous media and is subjected to a horizontal plane Poiseuille flow. We use Brinkman's model instead of Darcy's law to describe the porous layer. The eigenvalue problem is solved by means of a Chebyshev collocation method. We study the influence of the depth ratio (d) over cap and the Darcy number delta on the instability of the system. We compare systematically the instability of Brinkman's model with the results of Darcy's model. Our results show that no satisfactory agreement between Brinkman's model and Darcy's model is obtained for the instability of a fluid-porous system. We also examine the instability of Darcy's model. A particular comparison with early work is made. We find that a multivalued region may present in the (k, Re) plane, which was neglected in previous work. Here k is the dimensionless wavenumber and Re is the Reynolds number. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.3000643]

Plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant materials

The plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant material are investigated in this paper. Applying the stress-strain relation for the pressure-sensitive dilatant material, we have obtained an exact asymptotic solution for the plane strain tip fields for two types of cracks, one of which lies in the pressure-sensitive dilatant material and the other in the elastic material and their tips touch both the bimaterial interface. In cases, numerical results show that the singularity and the angular variations of the fields obtained depend on the material hardening exponent n, the pressure sensitivity parameter mu and geometrical parameter lambda.

Asymptotic analysis of plane-strain mode I steady-state crack growth in transformation toughening ceramics (II)

Based on a constitutive law which includes the shear components of transformation plasticity, the asymptotic solutions to near-tip fields of plane-strain mode I steadity propagating cracks in transformed ceramics are obtained for the case of linear isotropic hardening. The stress singularity, the distributions of stresses and velocities at the crack tip are determined for various material parameters. The factors influencing the near-tip fields are discussed in detail.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.