Elastic study on singularities interacting with interfaces using alternating technique: Part I. Anisotropic trimaterial

Schwarz–Neumann's alternating technique is applied to singularity problems in an anisotropic `trimaterial', which denotes an infinite body composed of three dissimilar materials bonded along two parallel interfaces. Linear elastic materials under general plane deformations are assumed, in which the plane of deformation is perpendicular to the two parallel interface planes.It is well known that if the solution is known for singularities in a homogeneous anisotropic medium, the solution for the same singularities in an anisotropic bimaterial can be constructed by the method of analytic continuation. It is shown here that the solution for singularities in a homogeneous medium may also be used as a base of the solution for the same singularities in a trimaterial. The alternating technique is applied to derive the trimaterial solution in a series form, whose convergence is guaranteed. The solution procedure is universal in the sense that no specific information about the singularity is needed. The energetic forces exerted on a dislocation due to interfaces are also evaluated from the trimaterial solution. The trimaterial solution studied here can be applied to a variety of problems, e.g. a bimaterial (including a half-plane problem), a finite thin film on semi-infinite substrate, and a finite strip of thin film, etc. Some examples are presented to verify the usefulness of the obtained solutions.     

Fracture Analysis for a Frictional Fastener Hole with Edge Cracks in an Anisotropic Plate∗

Abstract

Stress intensity factors are evaluated for a singly or doubly cracked fastener hole with frictional traction in an anisotropic plate, using a special kernel boundary integral equation (BIE) approach. The integration kernel (Green's function) used in this BIE approach has already taken the presence of the crack (or cracks) into account, thus.avoiding the need for element discretization to model the stress singularity at the crack tip. The Green's function employed is that of an infinite anisotropic plate containing an elliptical hole or crack, and subjected to an arbitrarily positioned point force. Several types of normal and shear traction conditions at the pinhole interface are considered. Numerical results are obtained for various geometrical and loading conditions and are compared with known solutions, where available, for their isotropic counterparts.     

Elastic study on singularities interacting with interfaces using alternating technique: Part II. Isotropic trimaterial

Singularity problems in an isotropic trimaterial are analyzed by the same procedure as in an anisotropic trimaterial of Part I [Int. J. Solids Struct. 39, 943–957]. `Trimaterial' denotes an infinite body composed of three dissimilar materials bonded along two parallel interfaces. Linear elastic isotropic materials under plane deformations are assumed, in which the plane of deformation is perpendicular to the two parallel interface planes, and thus Muskhelishvili's complex potentials are used. The method of analytic continuation is alternatively applied across the two parallel interfaces in order to derive the trimaterial solution in a series form from the corresponding homogeneous solution. A variety of problems, e.g. a bimaterial (including a half-plane problem), a finite thin film on semi-infinite substrate, and a finite strip of thin film, etc, can be analyzed as special cases of the present study. A film/substrate structure with a dislocation is exemplified to verify the usefulness of the solutions obtained.     

Anti-plane time-harmonic Green's functions for a circular inhomogeneity in piezoelectric medium with a spring- or membrane-type interface

The present work focuses on the two-dimensional anti-plane time-harmonic dynamic Green's functions for a circular inhomogeneity immersed in an infinite matrix with an imperfect interface, where both the inhomogeneity and matrix are assumed to be piezoelectric and transversely isotropic. Two types of imperfect interface, the spring-type interface with electromechanical coupling and the membrane-type interface, are considered. The former type is often used to model the electromechanical damage along the interface while the latter is largely employed to simulate surface/interface effect of nano-sized inhomogeneity. By using the Bessel function expansions, explicit solutions for the electromechanical fields induced by a time-harmonic anti-plane line force and line charge located in an unbounded matrix as well as the circular inhomogeneity are respectively derived. The present solutions can recover the anti-plane Green's functions for some special cases, such as the dynamic or quasi-static Green's functions of piezoelectricity with perfect interface as well as the dynamic or quasi-static Green's functions of pure elasticity with imperfect interface. For detailed discussions, selected analytical results are presented. Dependence of the electromechanical fields on circular frequency as well as interface properties is examined. The size effect related to interfacial imperfection is also discussed.     

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.     

Time-dependent Green's functions for an anisotropic bimaterial with viscous interface

By virtue of the Stroh formalism, we derive the exact closed-form solutions for the time-dependent two-dimensional Green's functions due to a line force and line dislocation in an anisotropic bimaterial with a viscous interface. We first reduce the boundary value problem to two coupled homogeneous first-order partial differential equations, which can be solved using a decoupling technique. The full-field expressions of the time-dependent displacements and stresses due to the line force and line dislocation interacting with the viscous interface are obtained.     

Green's function solution and displacement potentials for Transformed Transversely Isotropic materials

A totally non-degenerate expression for the Green's function of infinite Transversely Isotropic (TI) materials is first deduced from the solutions given by Pan and Chou [Pan, Y.-C., Chou, T.-W., 1976. Point force solution for an infinite transversely isotropic solid. Trans. ASME, J. Appl. Mech. 43 (4), 608–612]. Then this solution and also the displacement potentials for TI materials are extended by a linear transformation to a larger family of anisotropic materials (Transformed Transversely Isotropic or TraTI materials). This family depends on 12 independent parameters and contains non-orthotropic materials and in this way a first explicit analytical solution for the Green's function for a non-orthotropic material is obtained. The TraTI materials which have orthotropic Symmetry (StraTI materials) constitute a sub-family depending on 6 independent parameters in the symmetry basis of the material. These materials present a 3D anisotropy (different stiffnesses in three orthogonal directions). General displacement potentials and the Green's function solution for STraTI materials can be deduced by a simple change and introducing one additional parameter in the well-known TI solutions.     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

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.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

A moving screw dislocation near interfacial rigid lines in two dissimilar anisotropic media

This paper attempts to investigate the problem for the interaction between a uniformly moving screw dislocation and interface rigid lines in two dissimilar.anisotropic. materials. Integrating Riemann-Schwarz＇s symmetry principle with the analysis singularity of complex functions, we present the general elastic solutions of this problem and the closed form solutions for interfaces containing one and two rigid lines. The expressions of stress intensity factors, at the rigid line tips and image force acting on moving dislocation are derived explicitly. The results show that dislocation velocity has an antishielding effect on the rigid line tip and a larger dislocation velocity leads to the equilibrium position of dislocation closing with the rigid line. The presented solutions contain previously known results as the special cases.     

On the solutions in the linear theory of micropolar viscoelasticity

In the present paper the linear theory of micropolar viscoelasticity is considered. The explicit expression of fundamental solution of the system of equations of steady vibrations is constructed by means of elementary functions and its basic properties are established. The Green's formulas in the considered theory are obtained. The formulas of integral representations of Somigliana-type of regular vector and regular (classical) solution are presented. The representation formulas of Galerkin-type solution of the system of nonhomogeneous equations and of the general solution of the system of homogeneous equations by means of eight metaharmonic functions are presented. The completeness of these solutions is proved.     

Magnetoelectric Green's functions and their application to the inclusion and inhomogeneity problems

Explicit expressions of magnetoelectric Green's functions are obtained for a transversely isotropic medium exhibiting coupling between the static electric and magnetic fields utilizing the contour integral representation. Four Green's functions exist which represent the coupled static electric and magnetic response to a unit point electric or magnetic charge. The Green's functions are applied to analyze the inclusion and inhomogeneity problems in an infinite magnetoelectric medium, and explicit, closed form expressions are obtained for the Eshelby type tensors. The magnetoelectric Eshelby's tensors can be readily used in the solution of numerous problems in the mechanics and physics of magnetoelectric solids.     

Non-singular direct formulation of boundary integral equations for potential flows

This paper presents a general direct integral formulation for potential flows. The singularities of Green's functions are desingularized theoretically, using a subtracting and adding back technique, so that Gaussian quadrature or any other numerical integration methods can be applied directly to evaluate all the integrals without any difficulty. When high-order quadrature formulas are applied globally, the number of unknowns can be reduced. Interpolation functions are not necessary for unknown variables in the present paper. Therefore, the present method is much simpler and more efficient than the conventional one. Several numerical examples are calculated and compared satisfactorily with analytical solutions or published results. © 1998 John Wiley & Sons, Ltd.     

A point heat source on the apex of a transversely isotropic piezothermoelastic cone

We use the compact harmonic general solutions of transversely isotropic piezothermoelastic materials to construct the three-dimensional solutions of a steady point heat source on the apex of a transversely isotropic piezothermoelastic cone by four newly introduced harmonic functions. All components of coupled field are expressed in terms of elementary functions and are convenient to use. When the apex angle 2α equals to π, the solutions reduce to those of the semi-infinite body problem, which is called Green's function. Numerical results are given graphically by contours.     

Uniformly moving screw dislocation interacting with interface cracks in anisotropic bimaterials

This paper attempts to investigate the problem for the interaction between a uniformly subsonic moving screw dislocation and interface cracks in two dissimilar anisotropic materials. Using Riemann–Schwarz’s symmetry principle integrated with the analysis singularity of complex functions, we present the general elastic solutions of this problem and the closed form solutions for interface containing one and two cracks. The expressions of stress intensity factors at the crack tips and image force acting on moving dislocation are derived explicitly. The results show that the stress intensity factors at the crack tips decrease with increasing velocity of dislocation, and larger dislocation velocity leads to the equilibrium position of dislocation leaving from crack tips. The presented solutions contain previously known results as the special cases.     

Some basic problems in anisotropic bimaterials with a viscoelastic interface

This research is devoted to the study of anisotropic bimaterials with Kelvin-type viscoelastic interface under antiplane deformations. First we derive the Green’s function for a bimaterial with a Kelvin-type viscoelastic interface subjected to an antiplane force and a screw dislocation by means of the complex variable method. Explicit expressions are derived for the time-dependent stress field induced by the antiplane force and screw dislocation. Also presented is the time-dependent image force acting on the screw dislocation due to its interaction with the Kelvin-type viscoelastic interface. Second we investigate a rectangular inclusion with uniform antiplane eigenstrains embedded in one of the two bonded anisotropic half-planes by virtue of the derived Green’s function for a line force. The explicit expressions for the time-dependent stress field induced by the rectangular inclusion are obtained in terms of the simple logarithmic and exponential integral functions. It is observed that in general the stresses exhibit the logarithmic singularity at the four corners of the rectangular inclusion. Our results also show that when one side of the rectangular inclusion lies on the viscoelastic interface, the interfacial tractions are still regular at the two corners of the inclusion which are located on the interface. Last we address a finite Griffith crack normal to the viscoelastic interface by means of the obtained Green’s function for a screw dislocation. The crack problem is formulated in terms of a resulting singular integral equation which is solved numerically. The time-dependent stress intensity factors at the two crack tips are obtained and some interesting features are discussed.     

A magnetoelectric screw dislocation interacting with a circular layered inclusion

Within the framework of the linear theory of magnetoelectroelasticity, the problem of a circular layered inclusion interacting with a generalized screw dislocation under remote anti-plane shear stress and in-plane magnetoelectric loads is investigated in this paper. The generalized dislocation can be located either in the matrix or in the circular layered inclusion. The layers are coaxial cylinders of annular cross-sections with arbitrary radii and different material properties. Using complex variable theory and the alternating technique, the solution of the present problem is expressed in terms of the solution of the corresponding homogeneous medium problem subjected to the same loading. Some numerical results are provided to investigate the influence of material combinations on the shear stress, electric field, magnetic and image force. These solutions can be used as Green's functions for the analysis of the corresponding magnetoelectric crack problem.     

High-order derivatives of Green’s functions in magneto-electro-elastic materials

Three-dimensional Green’s functions and their arbitrary order derivatives in general anisotropic magneto-electro-elastic materials are derived by using Fourier transform. They are analytical solutions expressed in line integral forms, and can be evaluated by a standard numerical integration method. With this method, we can obtain results with high accuracy. Besides, a numerical finite difference method is also given to evaluate the second-order derivatives quickly. When setting the appropriate material coefficients to zero, the piezoelectric, piezomagnetic, and purely anisotropic elastic Green’s functions and their derivatives can all be obtained from the current solutions.     

Two-dimensional Green's functions in anisotropic multiferroic bimaterials with a viscous interface

We derive, by virtue of the unified Stroh formalism, the extremely concise and elegant solutions for two-dimensional and (quasi-static) time-dependent Green's functions in anisotropic magnetoelectroelastic multiferroic bimaterials with a viscous interface subjected to an extended line force and an extended line dislocation located in the upper half-plane. It is found for the first time that, in the multiferroic bimaterial Green's functions, there are 25 static image singularities and 50 moving image singularities in the form of the extended line force and extended line dislocation in the upper or lower half-plane. It is further observed that, as time evolves, the moving image singularities, which originate from the locations of the static image singularities, will move further away from the viscous interface with explicit time-dependent locations. Moreover, explicit expression of the time-dependent image force on the extended line dislocation due to its interaction with the viscous interface is derived, which is also valid for mathematically degenerate materials. Several special cases are discussed in detail for the image force expression to illustrate the influence of the viscous interface on the mobility of the extended line dislocation, and various interesting features are observed. These Green's functions can not only be directly applied to the study of dislocation mobility in the novel multiferroic bimaterials, they can also be utilized as kernel functions in a boundary integral formulation to investigate more complicated boundary value problems where multiferroic materials/composites are involved.