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.

Continuous spectrum growth and modal instability in swirling duct flow

We present results on the stability of compressible inviscid swirling flows in an annular duct. Such flows are present in aeroengines, for example in the by-pass duct, and there are also similar flows in many aeroacoustic or aeronautical applications. The linearised Euler equations have a ('critical layer') singularity associated with pure convection of the unsteady disturbance by the mean flow, and we focus our attention on this region of the spectrum. By considering the critical layer singularity, we identify the continuous spectrum of the problem and describe how it contributes to the unsteady field. We find a very generic family of instability modes near to the continuous spectrum, whose eigenvalue wavenumbers form an infinite set and accumulate to a point in the complex plane. We study this accumulation process asymptotically, and find conditions on the flow to support such instabilities. It is also found that the continuous spectrum can cause a new type of instability, leading to algebraic growth with an exponent determined by the mean flow, given in the analysis. The exponent of algebraic growth can be arbitrarily large. Numerical demonstrations of the continuous spectrum instability, and also the modal instabilities are presented.

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.

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.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

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.