Cauchy singular integral equation method for transient antiplane dynamic problems

In this paper, by use of the boundary integral equation method and the techniques of Green basic solution and singularity analysis, the dynamic problem of antiplane is investigated. The problem is reduced to solving a Cauchy singular integral equation in Laplace transform space. This equation is strictly proved to be equivalent to the dual integral equations obtained by Sih [Mechanics of Fracture, Vol. 4. Noordhoff, Leyden (1977)]. On this basis, the dynamic influence between two parallel cracks is also investigated. By use of the high precision numerical method for the singular integral equation and Laplace numerical inversion, the dynamic stress intensity factors of several typical problems are calculated in this paper. The related numerical results are compared to be consistent with those of Sih. It shows that the method of this paper is successful and can be used to solve more complicated problems. Copyright (C) 1996 Elsevier Science Ltd

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突加载荷下，裂纹尖端动应力强度因子的时间响应。首先，运用积分变换方法将黏弹性混合边值问题化成变换域上的对偶积分方程。通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程。运用分片连续函数法数值求解奇异积分方程，得到变换域内的动应力强度因子。再用Laplace积分变换数值反演方法，将变换域的解反演到时间域内，最终求得动应力强度因子的时间响应，并对黏弹性参数的影响进行了分析。

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.

A Study of Size-Dependent Microindentation

We recently proposed a strain gradient theory to account for the size dependence of plastic deformation at micron and submicron length scales. The strain gradient theory includes the effects of both rotation gradient and stretch gradient such that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the stretch gradient measures explicitly enter the constitutive relations through the instantaneous tangent modulus. Indentation tests at scales on the order of one micron have shown that measured hardness increases significantly with decreasing indent size. In the present paper, the strain gradient theory is used to model materials undergoing small-scale indentations. A strong effect of including strain gradients in the constitutive description is found with hardness increasing by a factor of two or more over the relevant range behavior. Comparisons with the experimental data for polycrystalline copper and single crystal copper indeed show an approximately linear dependence of the square of the hardness, H 2, on the inverse of the indentation depth, 1/h, I.e., H-2 proportional to 1/h, which provides an important self-consistent check of the strain gradient theory proposed by the authors earlier.

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.

A New Deformation Theory with Strain Gradient Effects

A new phenomenological deformation theory with strain gradient effects is proposed. This theory, which belongs to nonlinear elasticity, fits within the framework of general couple stress theory and involves a single material length scale l. In the present theory three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom u(i). omega(i) has no direct dependence upon ui and is called the micro-rotation, i.e. the material rotation theta(i) plus the particle relative rotation. The strain energy density is assumed to only be a function of the strain tensor and the overall curvature tensor, which results in symmetric Cauchy stresses. Minimum potential principle is developed for the strain gradient deformation theory version. In the limit of vanishing 1, it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in details. Comparisons between the present theory and the theory proposed by Shizawa and Zbib (Shizawa, K., Zbib, H.M., 1999. A thermodynamical theory gradient elastoplasticity with dislocation density Censor: fundamentals. Int. J. Plast. 15, 899) are given. With the same hardening law as Fleck et al. (Fleck, N.A., Muller, G.H., Ashby, M.F., Hutchinson, JW., 1994 Strain gradient plasticity: theory and experiment. Acta Metall. Mater 42, 475), the new strain gradient deformation theory is used to investigate two typical examples, i.e. thin metallic wire torsion and ultra-thin metallic beam bend. The results are compared with those given by Fleck et al, 1994 and Stolken and Evans (Stolken, J.S., Evans, A.G., 1998. A microbend test method for measuring the plasticity length scale. Acta Mater. 46, 5109). In addition, it is explained for a unit cell that the overall curvature tensor produced by the overall rotation vector is the work conjugate of the overall couple stress tensor. (C) 2002 Elsevier Science Ltd. All rights reserved.

On a Degenerate Parabolic Problem in Holder Spaces

In this paper, we study some degenerate parabolic equation with Cauchy-Dirichlet boundary conditions. This problem is considered in little Holder spaces. The optimal regularity of the solution v is obtained and is specified in terms of those of the second member when some conditions upon the Holder exponent with respect to the degeneracy are satisfied. The proofs mainly use the sum theory of linear operators with or without density of domains and the results of smoothness obtained in the study of some abstract linear differential equations of elliptic type.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.