GREEN''''S FUNCTIONS FOR 2D PROBLEMS OF ANISOTROPIC PIEZOELECTRIC MATERIAL WITH A STRIP REGION

By using Stroh's formalism and the conformal mapping technique, we derive the simple explicit elastic fields of a generalized line dislocation and a generalized line force in a general anisotropic piezoelectric strip with fixed surfaces, which are two fixed conductor electrodes. The solutions obtained are usually considered as Green's functions which play important roles in the boundary element methods. The Coulomb forces of the distributed charges along the region boundaries on the line chargeq atz ⁰ are analysed in detail. The results are valid not only for plane and antiplane problems but also for the coupled problems between inplane and outplane deformations.     

Interaction between elliptical and ellipsoidal inclusions under bending stress fields

Summary  This paper deals with interaction problems of elliptical and ellipsoidal inclusions under bending, using singular integral equations of the body force method. The problems are formulated as a system of singular integral equations with Cauchy-type or logarithmic-type singularities, where unknown functions are densities of body forces distributed in the x,y and r,θ,z directions in infinite bodies having the same elastic constants as those of the matrix and inclusions. In order to satisfy the boundary conditions along the elliptical and the ellipsoidal boundaries, the unknown functions are approximated by a linear combination of fundamental density functions and polynomials. The present method is found to yield the exact solutions for a single elliptical or spherical inclusion under a bending stress field. It yields rapidly converging numerical results for interface stresses in the interaction of inclusions. Received 9 September 1999; accepted for publication 15 January 2000     

Three-dimensional solutions for general anisotropy

The Stroh formalism is extended to provide a new class of three-dimensional solutions for the generally anisotropic elastic material that have polynomial dependence on x₃, but which have quite general form in x₁,x₂. The solutions are obtained by a sequence of partial integrations with respect to x₃, starting from Stroh's two-dimensional solution. At each stage, certain special functions have to be introduced in order to satisfy the equilibrium equation. The method provides a general analytical technique for the solution of the problem of the prismatic bar with tractions or displacements prescribed on its lateral surfaces. It also provides a particularly efficient solution for three-dimensional boundary-value problems for the half-space. The method is illustrated by the example of a half-space loaded by a linearly varying line force.     

Explicit expressions of the Barnett-Lothe tensors for anisotropic materials

The three Barnett-Lothe tensorsS, H, andL appear frequently in the real form solutions of two-dimensional anisotropic elasticity problems. Explicit expressions for the components of these tensors are presented for general anisotropic materials. The special cases of monoclinic materials with the plane of material symmetry at x₃=0, x₂=0, and x₁=0 are then deduced. For monoclinic materials with the symmetry plane at x₂=0 or x₁=0, the locations of image singularities for the Green's functions for a half-space have a special geometry.     

Elastic fields of dislocation loops in three-dimensional anisotropic bimaterials

By applying semi-analytical point-force Green's functions obtained via the Stroh formulism, we derive simple line integrals to calculate the elastic displacement and stress fields for a three-dimensional dislocation loop in an anisotropic bimaterial system. The solutions for the case of anisotropy are more convenient for treating an arbitrary dislocation loop compared with traditional area integration. With this new formulation, we numerically examine the displacement, stress, and energy due to the interaction between a dislocation loop and the bimaterial interface in an Al–Cu system. The interactive image energy due to the elastic moduli mismatch across the interface is then numerically evaluated. The result shows that a dislocation loop is subjected to an attractive force by the interface when it lies in the stiff material, and a repulsive force when it lies in the soft material. Moreover, the dependence of the interactive image energy of a dislocation loop on the position and size of the dislocation loop are also demonstrated and discussed. Significantly, it is found that the interactive image energy for a dislocation loop depends only on the ratio d/a, where a is the loop diameter and d is its distance to the interface. The examples studied provide benchmark solutions for anisotropic bimaterial dislocation problems.     

Thermo-elastic crack branching in general anisotropic media

Thermal fields may exist in addition to mechanical loading, for example, due to short term exposure to fire. In this paper, the branching of cracks in the presence of combined thermal and mechanical loads is investigated for general anisotropic media by employing the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. A general solution to the thermo-elastic crack problem for an anisotropic material under arbitrary loading is obtained in a compact form. Green’s functions are also presented for a thermal dislocation (heat vortex) and a conventional dislocation (or, referred as mechanical dislocation), which are formulated considering the cuts located at an arbitrary angle with respect to the x₁ axis of the coordinate system (x₁, x₂, x₃). Using the derived compact expressions, the interaction between the crack and the dislocation is studied and a closed form solution for this interaction is obtained. The branching portion of the thermo-elastic crack is modelled as a continuous distribution of dislocations. This problem is then converted into a set of singular integral equations. Numerical results are presented to illustrate the possible effects of thermal loading on the propagation of the branched crack.     

条带域内二维各向异性压电介质的Green函数

以压电各向异性弹性介质广义平面变形的Stroh一般解为基础，采用复变函数方法（即保角变换技术），研究了条带域介质内物理场的封闭形式解，求得了介质内某一点同时存在广义线位错和广义线力作用时的简单明确解，它就是边界元法中的Green函数，还分析了极化介质表面的电荷分布情况，并进而讨论了线电荷q与边界分布电荷间的库仑力问题，文中结果不仅适用于平面或反平面变形问题，而且也适用于两者耦合的二维变形问题。     

Rectilinear dislocation in an anisotropic plate

A solution has been found to the problem of calculating the stress and displacement fields caused by a rectilinear dislocation in an anisotropic elastic plate. Special cases of anisotropy have been found with solutions represented by elementary functions.Certain problems in describing crystal plastic deformation phenomena make it vital to know the fields of the elastic stresses and displacements caused by an individual dislocation in a bounded crystal. It is interesting to study the effect of crystal boundaries on these fields with a simple model which approximates fairly closely to experimental conditions.The model selected is shown in Fig, 1. A dislocation with a Burgers vector (b₁, b₂, b₃) is situated in an infinite elastic anisotropic plate of thickness 2h. The dislocation line is parallel to the plate boundaries. The following restriction is introduced in relation to the plate's elastic properties: the medium has a plane of elastic symmetry perpendicular to the dislocation line. The selection of the coordinate system and position of the dislocation are shown in Fig. 1. The requirement is to find the stresses and displacements at an arbitrary point in the plate.One limited special form of this problem has been solved by Kroupa [1]. The limitations which he introduced are as follows: the medium is isotropic, the dislocation is at the precise center of the band and the Burgers vector has only one component b₂ differing from zero (the same coordinates were chosen in [1] as in Fig. 1).Thus Kroupa's results can be obtained from the results of the present work as a special case. Other special cases arising from this problem are those concerning the elastic stress and displacement fields caused by a dislocation in anisotropic semi-bounded [2] and bounded [3] media.It is immediately apparent that the problem is a plane one, in the sense that the fields to be found do not depend on coordinate z. Since the medium has a plane of elastic symmetry perpendicular to the dislocation line, it is clear from [4] that the system of stresses and strains in such a medium can be divided into two independent subsystems. The first of these is plane deformation with stress components _xx, _yy and _xy differing from zero and displacement vector components u_x and u_y, the second is antiplane deformation with stress components _xz and _yz differing from zero and the displacement vector component u_z.In the case under examination, the plane deformation is caused by the Burgers vector edge components b_x and b_y and the antiplane deformation by the screw component b_z. The solution is therefore divided into two stages, corresponding to edge and screw dislocations.In conclusion, I wish to thank A. M. Kosevich for his valuable advice and L. A. Pastur for his constant vigilance and assistance in the work.     

Anisotropic Elastic Materials that Uncouple All Three Displacement Components, and Existence of One-Displacement Green's Function

It was shown in an earlier paper that, under a two-dimensional deformation, there are anisotropic elastic materials for which the antiplane displacement u ₃ and the inplane displacements u ₁, u ₂ are uncoupled but the antiplane stresses σ₃₁, σ₃₂ and the inplane stresses σ₁₁, σ₁₂, σ₂₂ remain coupled. The conditions for this to be possible were derived, but they have a complicated expression. In this paper new and simpler conditions are obtained, and a general anisotropic elastic material that satisfies the conditions is presented. For this material, and for certain monoclinic materials with the symmetry plane at x ₃ = 0, we show that the unnormalized Stroh eigenvectors a _k for k = 1, 2, 3 are all real. The matrix A =[a ₁, a ₂, a ₃] is a unit matrix when the material has a symmetry plane at x ₂ = 0. Thus any one of the u ₁, u ₂, u ₃ can be the only nonzero displacement, and the solution is a one-displacement field. Application to the Green's function due to a line of concentrated force f and a line dislocation with Burgers vector v in the infinite space, the half-space with a rigid boundary, and the infinite space with an elliptic rigid inclusion shows that one can indeed have a one-displacement field u ₁, u ₂ or u ₃. One can also have a two-displacement field polarized on a plane other than the (x ₁, x ₂)-plane. The material that uncouples u ₁, u ₂, u ₃ is not as restrictive as one might have thought. It can be triclinic, monoclinic, orthotropic, tetragonal, transversely isotropic, or cubic. However, it cannot be isotropic. This revised version was published online in August 2006 with corrections to the Cover Date.     

Pointwise and Other Decay Estimates for an Isotropic Elastic Strip

This paper is concerned with an homogeneous isotropic linear elastic strip, in plane strain. It is supposed that its lateral boundaries are displacement-free and that the deformation is generated by actions on the ends. A cross-sectional measure of deformation, complementing that of a previous paper, is defined and shown to satisfy a generalised convexity condition in the axial variable x ₁, for materials with negative Poisson's ratio σ. An enhanced measure is subsequently defined, and, in the case of a semi-infinite strip, is shown to yield pointwise exponential decay estimates for both the axial and the transverse displacement components for materials with positive Poisson's ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

A moving dislocation in a magneto–electro-elastic solid

This work considers the generalized plane problem of a moving dislocation in an anisotropic elastic medium with piezoelectric, piezomagnetic and magnetoelectric effects. The closed-form expressions for the elastic, electric and magnetic fields are obtained using the extended Stroh formalism for steady-state motion. The radial components, E_rand H_r, of the electric and magnetic fields as well as the hoop components, D_θ and B_θ, of electric displacement and magnetic flux density are found to be independent of θ in a polar coordinate system. This interesting phenomenon is proven to be is a consequence of the electric and magnetic fields, electric displacement and magnetic flux density that exhibit the singularity r⁻¹ near the dislocation core. As an illustrative example, the more explicit results for a moving dislocation in a transversely isotropic magneto–electro-elastic medium are provided and the behavior of the coupled fields is analyzed in detail.     

A piezoelectric screw dislocation interacting with a half-plane trimaterial composite

The electro-elastic stress investigation on the interaction between a screw dislocation and a half-plane trimaterial composite composed of three bonded dissimilar transversely isotropic piezoelectric materials is analyzed in the framework of linear piezoelectricity. Each layer is assumed to have the same material orientation with x ₃ in the poling direction. The dislocations are characterized by a discontinuous displacement and electric potential across the slip plane and are subjected to a line force and a line charge at the core. Based on the complex variable and the method of alternating technique, the solution of electric field and displacement field is expressed in terms of explicit series form. The solutions derived here can be applied to a variety of problems, for example, a half-plane bimaterial, a quarter-plane bimaterial, a quarter-plane material and a rectangular strip etc. Numerical results are provided to show the influences of the material combinations and geometric configurations on the electro-elastic fields and image force calculated through the generalized Peach-Koehler formula. The solutions proposed here can be served as Green??s functions for the analyses corresponding piezoelectric cracking problems.     

A Crack with an Electric Displacement Saturation Zone in an Electrostrictive Material

A crack with an electric displacement saturation zone in an electrostrictive material under purely electric loading is analyzed. A strip saturation model is here employed to investigate the effect of the electrical polarization saturation on electric fields and elastic fields. A closed form solution of electric fields and elastic fields for the crack with the strip saturation zone is obtained by using the complex function theory. It is found that the K _I -dominant region is very small compared to the strip saturation zone. The generalized Dugdale zone model is also employed in order to investigate the effect of the saturation zone shape on the stress intensity factor. Using the body force analogy, the stress intensity factor for the asymptotic problem of a crack with an elliptical saturation zone is evaluated numerically.     

Surface waves in an exponentially graded, general anisotropic elastic material under the influence of gravity

In a recent paper Destrade [1] studied surface waves in an exponentially graded orthotropic elastic material. He showed that the quartic equation for the Stroh eigenvalue p is, after properly modified, a quadratic equation in p² with real coefficients. He also showed that the displacement and the stress decay at different rates with the depth x₂ of the half-space. Vinh and Seriani [2] considered the same problem and added the influence of gravity on surface waves. In this paper we generalize the problem to exponentially graded general anisotropic elastic materials. We prove that the coefficients of the sextic equation for p remain real and that the different decay rates for the displacement and the stress hold also for general anisotropic materials. A surface wave exists in the graded material under the influence of gravity if a surface wave can propagate in the homogeneous material without the influence of gravity in which the material parameters are taken at the surface of the graded half-space. As the wave number k → ∞, the surface wave speed approaches the surface wave speed for the homogeneous material. A new matrix differential equation for surface waves in an arbitrarily graded anisotropic elastic material under the influence of gravity is presented. Finally we discuss the existence of one-component surface waves in the exponentially graded anisotropic elastic material with or without the influence of gravity.     

A layered composite containing a crack in its lower layer loaded by a rigid stamp

The elastostatic plane problem of a layered composite containing an internal or edge crack perpendicular to its boundaries in its lower layer is considered. The layered composite consists of two elastic layers having different elastic constants and heights and rests on two simple supports. Solution of the problem is obtained by superposition of solutions for the following two problems: The layered composite subjected to a concentrated load through a rigid rectangular stamp without a crack and the layered composite having a crack whose surface is subjected to the opposite of the stress distribution obtained from the solution of the first problem. Using theory of Elasticity and Fourier transform technique, the problem is formulated in terms of two singular integral equations. Solving these integral equations numerically by making use of Gauss–Chebyshev integration, numerical results related to the normal stress σ_x(0,y), the stress-intensity factors, and the crack opening displacements are presented and shown graphically for various dimensionless quantities.     

Three-dimensional elastic displacements induced by a dislocation of polygonal shape in anisotropic elastic crystals

Dislocations and the elastic fields they induce in anisotropic elastic crystals are basic for understanding and modeling the mechanical properties of crystalline solids. Unlike previous solutions that provide the strain and/or stress fields induced by dislocation loops, in this paper, we develop, for the first time, an approach to solve the more fundamental problem—the anisotropic elastic dislocation displacement field. By applying the point-force Green’s function for a three-dimensional anisotropic elastic material, the elastic displacement induced by a dislocation of polygonal shape is derived in terms of a simple line integral. It is shown that the singularities in the integrand of this integral are all removable. The proposed expression is applied to calculate the elastic displacements of dislocations of two different fundamental shapes, i.e. triangular and hexagonal. The results show that the displacement jump across the dislocation loop surface exactly equals the assigned Burgers vector, demonstrating that the proposed approach is accurate. The dislocation-induced displacement contours are also presented, which could be used as benchmarks for future numerical studies.     

Energetics of misfit dislocation dipoles in anisotropic epitaxial films with nanoscale compositional modulation

The elastic strain and stress fields associated with nanoscale compositional modulation in an anisotropic epitaxial film on an anisotropic substrate are obtained by using Stroh formalism and the Eshelby-type inclusion method. The composition of the epitaxial film is considered to periodically fluctuate in a surface soft mode, with the amplitude of the composition modulation maximal near the growing surface and decreasing exponentially into the film. It has been experimentally observed that the composition modulation affects the formation of a new type of crystal defects, i.e., misfit dislocation dipoles, in III–V compound semiconductor materials. The formation energy of a misfit dislocation dipole under the elastic fields due to the composition modulation is calculated in this study. It is composed of the core and self energies of two dislocations, the interaction energy between two dislocations, and the interaction energies between the composition modulation and two dislocations. Numerical calculations are performed for a dislocation dipole in a lattice-matched Ga_0.5In_0.5P film on a GaAs substrate.     

End effects for plane deformations of an elastic anisotropic semi-infinite strip

In the linear theory of elasticity, Saint-Venant's principle is used to justify the neglect of edge effects when determining stresses in a body. For isotropic materials, the validity of this is well established. However for anisotropic and composite materials, experimental results have shown that edge effects may persist much farther into the material than for isotropic materials and as a result cannot be neglected. This paper further examines the effects of material anisotropy on the exponential decay rate for stresses in a semi-infinite elastic strip. A linearly elastic semi-infinite strip in a state of plane stress/strain subject to a self-equilibrated end load is considered first for a specially orthotropic material and then for the general anisotropic material. The problem is governed by a fourth-order elliptic partial differential equation with constant coefficients. In the former case, just a single dimensionless material parameter appears, while in the latter, only three dimensionless parameters are required. Energy methods are used to establish lower bounds on the actual stress decay rate. Both analytic and numerical estimates are obtained in terms of the elastic constants of the material and results are shown for several contemporary engineering materials. When compared with the exact stress decay rate computed numerically from the eigenvalues of a fourth-order ordinary differential equation, the results in some cases show a high degree of accuracy. In particular, for strongly orthotropic materials, an asymptotic estimate provides extremely accurate estimates for the decay rate. Results of the type obtained here have several important practical applications. For example, they provide physical insight into the mechanical testing of anisotropic and laminated composite structures (including the off-axis tension test), are useful in assessing the influence of fasteners, joints, etc. on the behavior of composite structures and allow for tailoring a material with specific properties to ensure that local stresses attenuate at a desired rate.     

Explicit expressions of S(v), H(v) and L(v) for anisotropic elastic materials

The three matrices L(v), S(v) and H(v), appearing frequently in the investigations of the two-dimensional steady state motions of elastic solids, are expressed explicitly in terms of the elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with a plane of symmetry at x₃ = 0, x₁ = 0, and x₂ = 0 are all deduced. Results for orthotropic materials appearing in the literature may be recovered from the present explicit expressions.     

Magnetic field and magneto elastic stress in an infinite plate containing an elliptical hole with an edge crack under uniform electric current

Two-dimensional solutions of the electric current, magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack under uniform electric current. Using a rational mapping function, the each solution is obtained as a closed form. The linear constitutive equation is used for the magnetic field and the stress analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate which raises a plane stress state for a thin plate and the deformation of the plate thickness. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, electric current, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that the stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solving the present magneto elastic stress problem, dislocation and rotation terms appear, which makes the present problem complicate. Solutions of the magneto elastic stress are nonlinear for the direction of electric current. Stresses in the direction of the plate thickness are caused and the solution is also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length and the electric current direction.