A point heat source on the surface of a semi-infinite transversely isotropic electro-magneto-thermo-elastic material

We use the compact mono-harmonic general solutions of transversely isotropic electro-magneto-thermo-elastic material to construct the three-dimensional Green’s function for a steady point heat source on the surface of a semi-infinite transversely isotropic electro-magneto-thermo-elastic material by five newly introduced mono-harmonic functions. All components of coupled field are expressed in terms of elementary functions and are convenient to use. Numerical results are given graphically by contours.     

Three‐dimensional fundamental solution for transversely isotropic piezothermoelastic material

Fundamental solutions play an important role in electroelastic analyses and numerical methods of piezoelectric material. However, most works available on this topic are on the case of identical temperature. We use the compact mono‐harmonic general solutions of transversely isotropic piezothermoelastic material to construct the three‐dimensional fundamental solution of a steady point heat source in an infinite piezothermoelastic material by four newly introduced mono‐harmonic functions. All components of coupled field are expressed in terms of elementary functions and are convenient to use. Numerical results for cadmium selenide are given graphically by contours. Copyright © 2008 John Wiley & Sons, Ltd.     

Green’s functions for transversely isotropic magnetoelectroelastic media

Based on the governing equations of transversely isotropic magnetoelectroelastic media, four general solutions on the cases of distinct eigenvalues and multiple eigenvalues are given and expressed in five mono-harmonic displacement functions. Then, based on these general solutions, employing the trial-and-error method, the three-dimensional Green’s functions of infinite, two-phase and semi-infinite magnetoelectroelastic media under point forces, point charge and magnetic monopole are all presented in terms of elementary functions for all cases of distinct eigenvalues and multiple eigenvalues. Numerical results are also presented.     

Green's functions for the isotropic or transversely isotropic space containing a circular crack

Summary Green's functions for an infinite three-dimensional elastic solid containing a circular crack are derived in terms of integrals of elementary functions. The solid is assumed to be either isotropic or transversely isotropic (with the crack being parallel to the plane isotropy).     

Three-dimensional analysis of piezoelectric/piezomagnetic elastic media

This paper presents an exact three-dimensional method of solution for a transversely isotropic piezoelectric/piezomagnetic elastic media. The control partial differential equation set that is expressed by the elastic displacement, electric potential and magnetism potential function are established. By means of drawing two-displacement potential functions the general solution for three-dimensional transversely isotropic magneto-electro-elastic is derived. As an illustrative example, the analysis solutions of a half space body of magneto-electro-elastic acted by a point force at the origin directed along the z-axis and a point charge and a point electric current at the origin are obtained.     

2D general solution and fundamental solution for orthotropic thermoelastic materials

The 2D general solution for the plane problem of thermoelastic materials is derived in terms of three harmonic functions using strict differential operator theory for the case of distinct eigenvalues. Based on the obtained general solution, the 2D fundamental solution for a steady line heat source in an infinite and a semi-infinite thermoelastic plane is obtained by three newly introduced harmonic functions. All components of coupled fields are expressed in terms of elementary functions and they are convenient to be used.     

Analysis of elliptical cracks perpendicular to the interface of two joined transversely isotropic solids

This paper presents a boundary element analysis of elliptical cracks in two joined transversely isotropic solids. The boundary element method is developed by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of stress intensity factors (SIFs) are obtained by using crack opening displacements. The results of the proposed method compare well with the existing exact solutions for an elliptical crack parallel to the isotropic plane of a transversely isotropic solid of infinite extent. Elliptical cracks perpendicular to the interface of transversely isotropic bi-material solids of either infinite extent or occupying a cubic region are further examined in detail. The crack surfaces are subject to the uniform normal tractions. The stress intensity factor values of the elliptical cracks of the two types are analyzed and compared. Numerical results have shown that the stress intensity factors are strongly affected by the anisotropy and the combination of the two joined solids.     

Boundary element analysis of three-dimensional crack problems in two joined transversely isotropic solids

The present paper presents a boundary element analysis of penny-shaped crack problems in two joined transversely isotropic solids. The boundary element analysis is carried out by incorporating the fundamental singular solution for a concentrated point load in a transversely isotropic bi-material solid of infinite space into the conventional displacement boundary integral equations. The conventional multi-region method is used to analyze the crack problems. The traction-singular elements are employed to capture the singularity around the crack front. The values of the stress intensity factors are obtained by using crack opening displacements. The numerical scheme results are verified with the closed-form solutions available in the literature for a penny-shaped crack parallel to the plane of the isotropy of a homogeneous and transversely isotropic solid of infinite extent. The new problem of a penny-shaped crack perpendicular to the interface of a transversely isotropic bi-material solid is then examined in detail. The crack surfaces are subject to the three normal tractions and the uniform shear traction. The associated stress intensity factor values are obtained and analyzed. The present results can be used for the prediction of the stability of composite structures and the hydraulic fracturing in deep rock strata and reservoir engineering.     

Stress intensity factor for transient thermal stresses in a transversely isotropic infinite body with an external circular crack

Summary The present work deals with the transient thermal stress in a transversely isotropic infinite body with an external circular crack. The surface cooling of the crack depends on position and time. Since it it usually very difficult to obtain an analytical solution for the temperature field, a finite difference formulation with respect to a tive variable is introduced. In the first step, applying this method to the general heat conduction equation in an orthotropic body, a very compact difference equation with respect to the spatial variables is obtained. In the second step, this method is applied to the transient thermoelastic problem in a transversely isotropic infinite body with an external circular crack subjected to heat exchange on the crack surface. Thermal stresses are analyzed by means of the transversely isotropic potential functions method.With 7 Figures     

Analytical solutions for an infinite transversely isotropic functionally graded sectorial plate subjected to a concentrated force or couple at the tip

This paper considers the elastic responses of an infinite sectorial plate made of transversely isotropic functionally graded material (FGM), which is subjected to a concentrated force or couple at the tip. There is no load acting on the upper and lower surfaces, and the elastic coefficients can vary arbitrarily through the plate thickness. No constraint is required on the symmetry of the plate in the thickness direction. Based on the displacement assumption for the bending of an FGM plate and by using the complex variable method, this paper presents the general solutions to the basic equations governing transversely isotropic FGM plates, which are expressed in terms of four analytical functions (or complex potentials). The boundary conditions are a combination of those from the plane elasticity and those from the classical plate theory. For a particular boundary value problem, such as the ones considered here for a sectorial plate, with the specific conditions for determining solutions, the four analytic functions can be assumed in appropriate forms, which contain only some unknown constants. Once these constants are determined from the specific conditions, the complete solutions are readily derived too. Among the solutions presented here, the solutions for the infinite FGM sectorial plate under a concentrated couple are absolutely new to the literature, and they are also applicable to isotropic FGM sectorial plates. The solutions degenerate into the ones for a homogeneous sectorial plate, which coincide with the available solutions from the plane elasticity theory. There are three-dimensional correction terms in the mid-plane displacements.     

Boundary integral equation method for conductive cracks in two and three-dimensional transversely isotropic piezoelectric media

Using the fundamental solutions and the Somigliana identity of piezoelectric medium, the boundary integral equations are obtained for a conductive planar crack of arbitrary shape in three-dimensional transversely isotropic piezoelectric medium. The singular behaviors near the crack edge are studied by boundary integral equation approach, and the intensity factors are derived in terms of the displacement discontinuity and the electric displacement boundary value sum near the crack edge on crack faces. The boundary integral equations for two dimensional crack problems are deduced as a special case of infinite strip planar crack. Based on the analogy of the obtained boundary integral equations and those for cracks in conventional isotropic elastic material and for contact problem of half-space under the action of a rigid punch, an analysis method is proposed. As an example, the solution to conductive Griffith crack is derived.     

3D point force solution for a permeable penny-shaped crack embedded in an infinite transversely isotropic piezoelectric medium

This paper considers the non-axisymmetric three-dimensional problem of a penny-shaped crack with permeable electric conditions imposed on the crack surfaces, subjected to a pair of point normal forces applied symmetrically with respect to the crack plane. The crack is embedded in an infinite transversely isotropic piezoelectric body with the crack face perpendicular to the axis of material symmetry. Applying the symmetry of the problem under consideration then leads to a mixed–mixed boundary value problem of a half-space, for which potential theory method is employed for the purpose of analysis. The cases of equal eigenvalues are also discussed. Although the treatment differs from that for an impermeable crack reported in literature, the resulting governing equation still has a familiar structure. For the case of a point force, exact expressions for the full-space electro-elastic field are derived in terms of elementary functions with explicit stress and electric displacement intensity factors presented. The exact solution for a uniform loading is also given.     

Three-dimensional exact solutions for free vibrations of simply supported magneto-electro-elastic cylindrical panels

This work investigates the free vibrations of magneto-electro-elastic cylindrical panels based on three-dimensional theory. Firstly, the general solutions for transversely isotropic magneto-electro-elastic materials are introduced and the displacement functions in the general solutions are expanded in trigonometric functions along the circumferential and axial directions. Then an ordinary differential equation of the displacement functions in radial direction is derived and solved. As a result, the frequency equations are obtained through the traction-free conditions on the cylindrical surfaces of the panel as well as the electric and magnetic conditions. For the torsion and thickness-shear modes, the frequency equations in simpler forms are presented. It is found that the magneto-electro-elastic coupling effects disappeared in torsion vibration. Meanwhile, the frequencies of pure elastic materials and magneto-electro-elastic materials have an explicit relation for the thickness-shear modes. The aforementioned solutions satisfy all the governing equations and boundary conditions point by point and they are three-dimensionally exact. Finally the numerical example demonstrates the present method and is compared with those from finite element method. Parametric investigation is also conducted to show the behavior of free vibrations of cylindrical panels.     

2D fundamental solution for orthotropic pyroelectric media

Fundamental solutions play an important role in electroelastic analyses of piezoelectric media. However, most works available are on the topic of identical temperature. On the basis of the compact general solution of orthotropic pyroelectric media, which is expressed in harmonic functions, and employing the trial-and-error method, the 2D fundamental solution for a steady point heat source in the interior of an infinite pyroelectric plane or on the surface of a semi-infinite pyroelectric plane are presented by four newly induced harmonic functions. All components of the coupled field are expressed in terms of elementary functions and are convenient to use.     

Hygrothermal stresses for a plane crack in a generally anisotropic plate

The solutions are presented for the hygrothermal stress field of a generally anisotropic plate under uniform heat flux and moisture concentration transfer obstructed by a hygrothermally insulated crack. For uncoupled diffusion of temperature and moisture, the solutions of both temperature and moisture are obtained directly from the Hilbert problem approach, and are treated as the particular solutions to a pair of nonhomogeneous partial differential equations for an uncoupled hygrothermoelastic system. The associated homogeneous solutions are expressed in terms of three stress functions based on the complex variable approach of Lekhnitskii. With some identities concerning the eigenvalues and eigenvectors, the general expressions of the stress and displacement fields can then be found in an explicit form. The stress intensity factors, crack opening displacements and energy release rate are expressed in terms of the heat flow, moisture concentration, material geometry, elastic and hygrothermal anisotropy. The simultaneous existence of mode I, II and III fracture is found to be due to material inherent anisotropy. Special cases for isotropic and orthotropic materials are also discussed.     

Application of finite-part integrals to three-dimensional fracture problems for piezoelectric media Part I: Hypersingular integral equation and theoretical analysis

This paper deals with some basic linear elastic fracture problems for an arbitrary-shaped planar crack in a three-dimensional infinite transversely isotropic piezoelectric media. The finite-part integral concept is used to derive hypersingular integral equations for the crack from the point force and charge solutions with distinct eigenvalues s _i(i=1,2,3) of an infinite transversely isotropic piezoelectric media. Investigations on the singularities and the singular stress fields and electric displacement fields in the vicinity of the crack are made by the dominant-part analysis of the two-dimensional integrals. Thereafter the stress and electric displacement intensity factor K-fields and the energy release rate G are exactly obtained by using the definitions of stress and electric displacement intensity factors and the principle of virtual work, respectively. The hypersingular integral equations under axially symmetric mechanical and electric loadings are solved analytically for the case of a penny-shaped crack.     

Sensitivities of Two-Dimensional Fracture Problems to the Near-tip Mesh Parameters

Analytical results for a penny-shaped crack with a plastic zone at the crack front are given. The crack is embedded in an infinite transversely isotropic elastic medium and is assumed to be subjected to two identical axisymmetric loads on the upper and lower crack faces. The size of the plastic zone at the crack front is determined by applying Dugdale hypothesis to the elasticity results for a penny-shaped crack. The size of the plastic zone is derived in terms of hyper-geometric functions. Expression of the normal stress outside the plastic zone is also given.     

Three-dimensional fracture analysis in transversely isotropic solids

In this paper a boundary element formulation for three-dimensional crack problems in transversely isotropic bodies is presented. Quarter-point and singular quarter-point elements are implemented in a quadratic isoparametric element context. The point load fundamental solution for transversely isotropic media is implemented. Numerical solutions to several three-dimensional crack problems are obtained. The accuracy and robustness of the present approach for the analysis of fracture mechanics problems in transversely isotropic bodies are shown by comparison of some of the results obtained with existing analytical solutions. The approach is shown to be a simple and useful tool for the evaluation of stress intensity factors in transversely isotropic media.     

Dynamic crack problems in three-dimensional transversely isotropic solids

In this paper, a general boundary element approach for three-dimensional dynamic crack problems in transversely isotropic bodies is presented for the first time. Quarter-point and singular quarter-point elements are implemented in a quadratic isoparametric element context. The procedure is based on the subdomain technique, the displacement integral representation for elastodynamic problems and the expressions of the time-harmonic point load fundamental solution for transversely isotropic media. Numerical results corresponding to cracks under the effects of impinging waves are presented. The accuracy of the present approach for the analysis of dynamic fracture mechanics problems in transversely isotropic solids is shown by comparison of the obtained results with existing solutions.     

Dynamic steady-state crack propagation in a transversely isotropic viscoelastic body

The problem considered herein is the dynamic, subsonic, steady-state propagation of a semi-infinite, generalized plane strain crack in an infinite, transversely isotropic, linear viscoelastic body. The corresponding boundary value problem is considered initially for a general anisotropic, linear viscoelastic body and reduced via transform methods to a matrix Riemann–Hilbert problem. The general problem does not readily yield explicit closed form solutions, so attention is addressed to the special case of a transversely isotropic viscoelastic body whose principal axis of material symmetry is parallel to the crack edge. For this special case, the out-of-plane shear (Mode III), in-plane shear (Mode II) and in-plane opening (Mode I) modes uncouple. Explicit expressions are then constructed for all three Stress Intensity Factors (SIF). The analysis is valid for quite general forms for the relevant viscoelastic relaxation functions subject only to the thermodynamic restriction that work done in closed cycles be non-negative. As a special case, an analytical solution of the Mode I problem for a general isotropic linear viscoelastic material is obtained without the usual assumption of a constant Poissons ratio or exponential decay of the bulk and shear relaxation functions. The Mode I SIF is then calculated for a generalized standard linear solid with unequal mean relaxation times in bulk and shear leading to a non-constant Poissons ratio. Numerical simulations are performed for both point loading on the crack faces and for a uniform traction applied to a compact portion of the crack faces. In both cases, it is observed that the SIF can vanish for crack speeds well below the glassy Rayleigh wave speed. This phenomenon is not seen for Mode I cracks in elastic material or for Mode III cracks in viscoelastic material.