Exact solution for mixed boundary value problems at anisotropic piezoelectric bimaterial interface and unification of various interface defects

The special mixed boundary value problem in which a debonded conducting rigid line inclusion is embedded at the interface of two piezoelectric half planes is solved analytically by employing the 8-D Stroh formalism. Different from existing interface insulating crack model and interface conducting rigid line inclusion model, the presently analyzed model is based on the assumption that all of the physical quantities, i.e., tractions, displacements, normal component of electric displacements and electric potential, are discontinuous across the interface defect. Explicit solutions for stress singularities at the tips of debonded conducting rigid line inclusion are obtained. Closed form solutions for the distribution of tractions on the interface, surface opening displacements and jump in electric potential on the debonded inclusion are also obtained, in addition real form solutions for these physical quantities are derived. Various forms of interface defect problems encountered in practice are solved within a unified framework and the stress singularities induced by those interface defects are discussed in detail. Particularly, we find that the analysis of interface cracks between the embedded electrode layer and piezoelectric ceramics can also be carried out within the unified framework.     

Analysis of a Partially Debonded Conducting Rigid Elliptical Inclusion in a Piezoelectric Matrix

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｅｔｏｔｈｅｉｒｉｎｔｒｉｎｓｉｃｅｌｅｃｔｒｏｍｅｃｈａｎｉｃａｌｃｏｕｐｌｉｎｇｐｒｏｐｅｒｔｉｅｓ,ｐｉｅｚｏｅｌｅｃｔｒｉｃｃｅｒａｍｉｃｓｈａｖｅｂｅｅｎｅｘｔｅｎｓｉｖｅｌｙｕｓｅｄｉｎｄｅｓｉｇｎｏｆｖａｒｉｏｕｓｅｌｅｃｔｒｏｎｉｃａｎｄｅｌｅｃｔｒｏｍｅｃｈａｎｉｃａｌｄｅｖｉｃｅｓｓｕｃｈａｓｓｅｎｓｏｒｓａｎｄａｃｔｕａｔｏｒｓ.Ｉｎｒｅｃｅｎｔｙｅａｒｓ,ｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｏｆｄｉｓｌｏｃａｔｉｏｎｓ ,ｃｒａｃｋｓ,ｃａｖｉｔｉｅ…     

Loosening of elastic inclusions

A numerical method is presented for simulating the occurrence of localized slip and separation along the interfaces of multiple, randomly distributed, circular elastic inclusions in an infinite elastic plane. The method is an extension of a direct boundary integral approach previously described elsewhere for solving problems involving perfectly bonded circular inclusions. Here, we allow displacement discontinuities to develop along the inclusion/matrix interfaces in accordance with a linear Mohr–Coulomb yield condition combined with a tensile strength cut-off. The displacements, tractions, and displacement discontinuities on the inclusion boundaries are all represented by truncated Fourier series, and an explicit iterative algorithm is adopted to determine zones of slip and separation under the prevailing loading conditions. Several examples are given to demonstrate the accuracy and generality of the approach.     

Three-dimensional solution of smart laminated anisotropic circular cylindrical shells with imperfect bonding

This paper presents an exact solution for a simply-supported and laminated anisotropic cylindrical shell strip with imperfect bonding at the off-axis elastic layer interfaces and with attached anisotropic piezoelectric actuator and sensor subjected to transverse loading. In this research, the imperfect interface conditions are described in terms of linear relations between the interface tractions in the normal and tangential directions, and the respective discontinuities in displacements. The solution for an elastic (or piezoelectric) layer of the smart laminated cylindrical shell strip is obtained in terms of the six-dimensional (or eight-dimensional) pseudo-Stroh formalism, solution for multilayered system is then derived based on the transfer matrix method. Finally, a numerical example is presented to demonstrate the effect of imperfect interface on the static response of the smart laminated cylindrical shell. The derived solutions can serve as benchmark results to assess various approximate shell theories and numerical methods.     

Uniform Stress Inside an Anisotropic Elliptic Inclusion with Imperfect Interface Bonding

In this paper we study the two-dimensional deformation of an anisotropic elliptic inclusion embedded in an infinite dissimilar anisotropic matrix subject to a uniform loading at infinity. The interface is assumed to be imperfectly bonded. The surface traction is continuous across the interface while the displacement is discontinuous. The interface function that relates the surface traction and the displacement discontinuity across the interface is a tensor function, not a scalar function as employed by most work in the literature. We choose the interface function such that the stress inside the elliptic inclusion is uniform. Explicit solution for the inclusion and the matrix is presented. The materials in the inclusion and in the matrix are general anisotropic elastic materials so that the antiplane and inplane displacements are coupled regardless of the applied loading at infinity. T.C.T. Ting is Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

On the Elastic Field of a Shpherical Inhomogeneity with an Imperfectly Bonded Interface

This study is devoted to the development of a unified and explicit elastic solution to the problem of a spherical inhomogeneity with an imperfectly bonded interface. Both tangential and normal displacement discontinuities at the interface are considered and a linear interfacial condition, which assumes that the tangential and the normal displacement jumps are proportional to the associated tractions, is adopted. The elastic disturbance due to the presence of an imperfectly bonded inhomogeneity is decomposed into two parts: the first is formulated in terms of an equivalent nonuniform eigenstrain distributed over a perfectly bonded spherical inclusion, while the second is formulated in terms of an imaginary Somigliana dislocation field which models the interfacial sliding and normal separation. The exact form of the equivalent nonuniform eigenstrain and the imaginary Somigliana dislocation are fully determined in this paper.     

直角平面区域内固定圆形刚性夹杂问题的Green函数解

利用复变函数法、多极坐标移动技术研究了直角平面区域内含有固定圆形夹杂时的反平面问题Green函数解.首先构造出不含夹杂的完整直角平面区域内满足边界应力条件的入射位移场;其次,建立直角平面区域内固定圆形夹杂对该入射场产生的满足直角边界应力自由条件的散射波解,并由叠加原理得到介质内的总波场.最后利用夹杂边界处的位移条件确定出散射波解中的未知系数,最终得到问题的Green函数解,还通过算例讨论了夹杂边界处的径向应力和环向应力随不同波数、角度和不同夹杂位置及不同点源位置的变化情况.算例结果表明了该文方法的有效实用性.     

Scattering far field solution of SH-wave by movable rigid cylindrical interface inclusion

The Green's function is used to solve the scattering far field solution of SH-wave by a movable rigid cylindrical interface inclusion in a linear elastic body. First, a suitable Green's function is developed, which is the fundamental displacement solution of an elastic half space with a movable rigid half-cylindrical inclusion impacted by out-of-plane harmonic line source loaded at any point of its horizontal surface. By using the Green's function, a series of Fredholm integral equations of the first kind which determine the scattering far field can be set up. Then the paper gives the expressions on the far field including the displacement mode of scattering wave and the solution of scattering cross-section. Finally, some examples and numerical results are discussed to analyze the influence of the combination of different media parameters on the answer of far field.     

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.     

Indentation of a spherical cavity in an elastic body by a rigid spherical inclusion: influence of non-classical interface conditions

This paper examines the indentation of an elastic body by a rigid spherical inclusion. In contrast to conventional treatments where the contact between a rigid inclusion and the elastic medium is regarded as being perfectly bonded, we examine the influence of non-classical interface conditions including frictionless bilateral contact, separation and Coulomb friction on the load–displacement behaviour of the spherical rigid inclusion. Both analytical methods and boundary element techniques are used to examine the inclusion/elastic medium interaction problems. This paper also provides a comprehensive review of non-classical interface conditions between inclusions and the surrounding elastic media.     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法, 在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路, 从Somigliana等式出发, 利用格林函数性质,得到了一种边界积分法, 使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到. 应用此新方法, 求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接, 将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设, 减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法, 求得界面处的位移与应力的值. 然后再求解域内位移与应力.得到了问题的精确解析解, 当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解. 求解过程表明,若问题的求解区域包含无穷远处时, 所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点, 试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.      

A Finite Strain Analysis of Cavity Formation at a Rigid Inhomogeneity

The formation of a cavity by inclusion-matrix interfacial separation is examined by analyzing the response of a plane rigid inclusion embedded in an unbounded incompressible matrix subject to remote equibiaxial dead load traction. A vanishingly thin interfacial cohesive zone, characterized by normal and tangential interface force-separation constitutive relations, is assumed to govern separation behavior. Rotationally symmetric cavity shapes (circles) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation and the remote loading. Nonsymmetrical cavity formation, under rotationally symmetric conditions of geometry and loading, is treated within the theory of infinitesimal strain superimposed on a given finite strain state. Rotationally symmetric and nonsymmetric bifurcations are analyzed and detailed results, for the Mooney–Rivlin strain energy density and for an exponential interface force-separation law, are presented. For the nonsymmetric rigid body displacement mode, a simple formula for the critical load is presented. The effect on bifurcation behavior of interfacial shear stiffness and other interface parameters is treated as well. In particular we demonstrate that (i) for the smooth interface nonsymmetric bifurcation always precedes rotationally symmetric bifurcation, (ii) unlike rotationally symmetric bifurcation, there is no threshold value of interface parameter for which nonsymmetric bifurcation will not occur and (iii) interfacial shear may significantly delay the onset of nonsymmetric bifurcation. Also discussed is the range of validity of a nonlinear infinitesimal strain theory previously presented by the author (Levy [1]). This revised version was published online in July 2006 with corrections to the Cover Date.     

Shear deformation effect in plates stiffened by parallel beams

In this paper a general solution for the analysis of shear deformable stiffened plates subjected to arbitrary loading is presented. According to the proposed model, the arbitrarily placed parallel stiffening beams of arbitrary doubly symmetric cross section are isolated from the plate by sections in the lower outer surface of the plate, taking into account the arising tractions in all directions at the fictitious interfaces. These tractions are integrated with respect to each half of the interface width resulting two interface lines, along which the loading of the beams as well as the additional loading of the plate is defined. Their unknown distribution is established by applying continuity conditions in all directions at the interfaces. The utilization of two interface lines for each beam enables the nonuniform distribution of the interface transverse shear forces and the nonuniform torsional response of the beams to be taken into account. The analysis of both the plate and the beams is accomplished on their deformed shape taking into account second-order effects. The analysis of the plate is based on Reissner’s theory, which may be considered as the standard thick plate theory with which all others are compared, while the analysis of the beams is performed employing the linearized second order theory taking into account shear deformation effect. Six boundary value problems are formulated and solved using the analog equation method (AEM), a BEM based method. The solution of the aforementioned plate and beam problems, which are nonlinearly coupled, is achieved using iterative numerical methods. The adopted model permits the evaluation of the shear forces at the interfaces in both directions, the knowledge of which is very important in the design of prefabricated ribbed plates. The effectiveness, the range of applications of the proposed method and the influence of shear deformation effect are illustrated by working out numerical examples with great practical interest.     

滑动界面球形夹杂对平面压缩波的散射

非理想粘结界面对多相材料力学性能的影响日益受到重视。本文研究了无限各向同性基体中的滑介面球形单夹杂对平面压缩的散射问题。夹杂与基体间的界面为非理想粘结界面,在剪应力的作用下将出现界面两侧相对滑移。假定界面相对滑动位移与界面剪应力成正比,在这种线弹簧型滑动界面条件下,通过波函数的级数展开法,获得了夹杂在基体中反射波和折射波以及应力场的解析表达式,并讨论了界面自由滑动和刚性夹杂等特例。     

Bifurcation phenomena in the rigid inclusion power law matrix composite sphere

Bifurcation of interface separation related to cavity nucleation is analyzed for a radially loaded composite sphere consisting of a rigid inclusion separated from a power law matrix by a uniform, non-linear cohesive zone. Equations for the spherically symmetric and non-symmetric problems are obtained from a hyperelastic finite strain theory by a limiting process that preserves non-linear matrix and interface response at infinitesimal strain. A complete solution to the symmetric problem is presented including bifurcation load, stresses, and evolution of elasto-plastic boundary and interface separation. An analysis of non-symmetric bifurcation, under symmetric conditions of geometry and loading, yields the bifurcation load and first non-symmetric mode shape associated with rigid inclusion displacement. An energy analysis is carried out for both symmetric and non-symmetric problems in order to assess stability of spherically symmetric states to spherically symmetric and non-symmetric “rigid body mode” perturbations.Results are provided for an interface force law that captures interface failure in normal mode and linear response in shear mode. For the symmetric problem, (i) there are threshold parameter values above which bifurcation will generally not occur, (ii) threshold values below which there do not exist equilibria in the post bifurcation regime, (iii) bifurcation occurs after attainment of the maximum interface strength. For the non-symmetric problem, (i) bifurcation always occurs, although it can be delayed by interfacial shear, (ii) for the smooth interface, non-symmetric bifurcation occurs after attainment of the maximum interface strength and always precedes symmetric bifurcation.     

Perturbation analysis of Mode III interfacial cracks advancing in a dilute heterogeneous material

The paper addresses the problem of a Mode III interfacial crack advancing quasi-statically in a heterogeneous composite material, that is a two-phase material containing elastic inclusions, both soft and stiff, and defects, such as microcracks, rigid line inclusions and voids. It is assumed that the bonding between dissimilar elastic materials is weak so that the interface is a preferential path for the crack. The perturbation analysis is made possible by means of the fundamental solutions (symmetric and skew-symmetric weight functions) derived in Piccolroaz et al. (2009). We derive the dipole matrices of the defects in question and use the corresponding dipole fields to evaluate “effective” tractions along the crack faces and interface to describe the interaction between the main interfacial crack and the defects. For a stable propagation of the crack, the perturbation of the stress intensity factor induced by the defects is then balanced by the elongation of the crack along the interface, thus giving an explicit asymptotic formula for the calculation of the crack advance. The method is general and applicable to interfacial cracks with general distributed loading on the crack faces, taking into account possible asymmetry in the boundary conditions.The analytical results are used to analyse the shielding and amplification effects of various types of defects in different configurations. Numerical computations based on the explicit analytical formulae allows for the analysis of crack propagation and arrest.     

Analytical solutions of uniform extended dislocations and tractions over a circular area in anisotropic magnetoelectroelastic bimaterials

In this paper, we derive the analytical solutions in a three-dimensional anisotropic magnetoelectroelastic bimaterial space subject to uniform extended dislocations and tractions within a horizontal circular area. By virtue of the Stroh formalism and Fourier transformation, the final expression of solutions in the physical domain contains only line integrals over [0,2π] rather than infinite integrals. As the reduced cases, the half-space and homogeneous full-space solutions can be directly derived from the present solutions. Also, in terms of material domains, the present solutions can be reduced to the piezoelectric, piezomagnetic, purely elastic materials with different symmetries of material property. To carry out numerical calculations, Gauss quadrature is adopted. In the numerical examples, the effect of different loading locations on the response at the interface is analyzed. It is shown that, when the magnetic traction or electric dislocation is applied, the physical quantities on the interface may not decrease monotonically as the loading area moves away from the interface. The distributions of different in-plane physical quantities on the upper and lower interfaces under various extended horizontal loadings are compared and the differences are discussed. The work presented in this paper can serve as benchmarks for future numerical studies in related research fields.     

The in-plane loading of rigid disc inclusion embedded in a crack

The paper examines the in-plane loading of a disc shaped rigid disc inclusion which is embedded in bonded contact with the plane surfaces of a penny-shaped crack. The mixed boundary value problem governing the elastostatic problem is reduced to the solution of a system of coupled integral equations, which are solved numerically to determine results of engineering interest. These results include the in-plane stiffness of the disc inclusion and the crack opening mode stress intensity factor at the boundary of the penny-shaped crack.