Elastic Solution of a Polygon-Shaped Inclusion With a Polynomial Eigenstrain

2021 ◽  
Vol 88 (6) ◽  
Author(s):  
Chunlin Wu ◽  
Huiming Yin

Abstract This paper presents the Eshelby’s tensor of a polygonal inclusion with a polynomial eigenstrain, which can provide an elastic solution to an arbitrary, convex inclusion with a continuously distributed eigenstrain by the Taylor series approximation. The Eshelby’s tensor for plane strain problem is derived from the fundamental solution of isotropic Green’s function with the Hadmard regularization, which is composed of the integrals of the derivatives of the harmonic and biharmonic potentials over the source domain. Using the Green’s theorem, they are converted to two line (contour) integrals over the polygonal cross section. This paper evaluates them by direct analytical integrals. Following Mura’s work, this paper formulates the method to derive linear, quadratic, and higher order of the Eshelby’s tensor in the polynomial form for arbitrary, convex polygonal shapes of inclusions. Numerical case studies were performed to verify the analytic results with the original Eshelby’s solution for a uniform eigenstrain in an ellipsoidal domain. It is of significance to consider higher order terms of eigenstrain for the polygon-shape inclusion problem because the eigenstrain distribution is generally non-uniform when Eshelby’s equivalent inclusion method is used. The stress disturbance due to a triangle particle in an infinite domain is demonstrated by comparison with the results of the finite element method (FEM). The present solution paves the way to accurately simulate the particle-particle, partial-boundary interactions of polygon-shape particles.

2021 ◽  
pp. 1-35
Author(s):  
Chunlin Wu ◽  
Liangliang Zhang ◽  
Huiming Yin

Abstract The paper extends the recent work (JAM, 88, 061002, 2021) of the Eshelby's tensors for polynomial eigenstrains from a two dimensional (2D) to three dimensional (3D) domain, which provides the solution to the elastic field with continuously distributed eigenstrain on a polyhedral inclusion approximated by the Taylor series of polynomials. Similarly, the polynomial eigenstrain is expanded at the centroid of the polyhedral inclusion with uniform, linear and quadratic order terms, which provides tailorable accuracy of the elastic solutions of polyhedral inhomogeneity by using Eshelby's equivalent inclusion method. However, for both 2D and 3D cases, the stress distribution in the inhomogeneity exhibits a certain discrepancy from the finite element results at the neighborhood of the vertices due to the singularity of Eshelby's tensors, which makes it inaccurate to use the Taylor series of polynomials at the centroid to catch the eigenstrain at the vertices. This paper formulates the domain discretization with tetrahedral elements to accurately solve for eigenstrain distribution and predict the stress field. With the eigenstrain determined at each node, the elastic field can be predicted with the closed-form domain integral of Green's function. The parametric analysis shows the performance difference between the polynomial eigenstrain by the Taylor expansion at the centroid and the 𝐶0 continuous eigenstrain by particle discretization. Because the stress singularity is evaluated by the analytical form of the Eshelby's tensor, the elastic analysis is robust, stable and efficient.


2001 ◽  
Vol 68 (5) ◽  
pp. 740-750 ◽  
Author(s):  
Q. Yang ◽  
W. Y. Zhou ◽  
G. Swoboda

In this paper, a three-dimensional penny-shaped isotropic inhomogeneity surrounded by unbounded isotropic matrix in a uniform stress field is studied based on Eshelby’s equivalent inclusion method. The solution including the deduced equivalent eigenstrain and its asymptotic expressions is presented in tensorial form. The so-called energy-based equivalent inclusion method is introduced to remove the singularities of the size and eigenstrain of the Eshelby’s equivalent inclusion of the penny-shaped inhomogeneity, and yield the same energy disturbance. The size of the energy-based equivalent inclusion can be used as a generic damage measurement.


Author(s):  
Y.-G. Lee ◽  
W.-N. Zou ◽  
E. Pan

This paper presents a closed-form solution for the arbitrary polygonal inclusion problem with polynomial eigenstrains of arbitrary order in an anisotropic magneto-electro-elastic full plane. The additional displacements or eigendisplacements, instead of the eigenstrains, are assumed to be a polynomial with general terms of order M + N . By virtue of the extended Stroh formulism, the induced fields are expressed in terms of a group of basic functions which involve boundary integrals of the inclusion domain. For the special case of polygonal inclusions, the boundary integrals are carried out explicitly, and their averages over the inclusion are also obtained. The induced fields under quadratic eigenstrains are mostly analysed in terms of figures and tables, as well as those under the linear and cubic eigenstrains. The connection between the present solution and the solution via the Green's function method is established and numerically verified. The singularity at the vertices of the arbitrary polygon is further analysed via the basic functions. The general solution and the numerical results for the constant, linear, quadratic and cubic eigenstrains presented in this paper enable us to investigate the features of the inclusion and inhomogeneity problem concerning polynomial eigenstrains in semiconductors and advanced composites, while the results can further serve as benchmarks for future analyses of Eshelby's inclusion problem.


2007 ◽  
Vol 336-338 ◽  
pp. 2513-2516
Author(s):  
Hua Jian Chang ◽  
Shu Wen Zhan

A micromechanical approach is developed to investigate the behavior of composite materials, which undergo interfacial delamination. The main objective of this approach is to build a bridge between the intricate theories and the engineering applications. On the basis of the spring-layer model, which is useful to treat the interfacial debonding and sliding, the present paper proposes a convenient method to assess the effects of delamination on the overall properties of composites. By applying the Equivalent Inclusion Method (EIM), two fundamental tensors are derived in the present model, the modified Eshelby tensor, and the compliance tensor (or stiffness tensor) of the weakened inclusions. Both of them are the fundamental tensors for constructing the overall constitutive law of composite materials. By simply substituting these tensors into an existing constitutive model, for instance, the Mori-Tanaka model, one can easily evaluate the effects of interfacial delamination on the overall properties of composite materials. Therefore, the present method offers a pretty convenient tool. Some numerical results are carried out in order to demonstrate the performance of this model.


2007 ◽  
Vol 334-335 ◽  
pp. 33-36 ◽  
Author(s):  
Akihiro Wada ◽  
Yusuke Nagata ◽  
Shi Nya Motogi

In this study, partially debonded spherical particles in a particulate composite are analyzed by three-dimensional finite element method to investigate their load carrying capacities, and the way to replace a debonded particle with an equivalent inclusion is examined. The variation in Young’s modulus and Poisson’s ratio of a composite with the debonded angle was evaluated for different particle arrangements and particle volume fractions, which in turn compared with the results derived from the equivalent inclusion method. Consequently, it was found that by replacing a debonded particle with an equivalent orthotropic one, the macroscopic behavior of the damaged composite could be reproduced so long as the interaction between neighboring particles is negligible.


1995 ◽  
Vol 62 (3) ◽  
pp. 590-594
Author(s):  
N. Chandra ◽  
Zhiyum Xie

A pair of two new tensors called GPS tensors S and D is proposed for the concentric cylindrical inclusion problem. GPS tensor S relates the strain in the inclusion constrained by the matrix of finite radius to the uniform transformation strain (eigenstrain), whereas tensor D relates the strain in the matrix to the same eigenstrain. When the cylindrical matrix is of infinite radius, tensor S reduces to the appropriate Eshelby’s tensor. Explicit expressions to evaluate thermal residual stresses σr, σθ and σz in the matrix and the fiber using tensor D and tensor S, respectively, are developed. Since the geometry of the present problem is of finite radius, the effect of fiber volume fraction on the stress distribution can be easily studied. Results for the thermal residual stress distributions are compared with Eshelby’s infinite domain solution and finite element results for a specified fiber volume fraction.


2000 ◽  
Vol 68 (1) ◽  
pp. 3-10 ◽  
Author(s):  
H. M. Shodja ◽  
A. S. Sarvestani

Consider a double-inhomogeneity system whose microstructural configuration is composed of an ellipsoidal inhomogeneity of arbitrary elastic constants, size, and orientation encapsulated in another ellipsoidal inhomogeneity, which in turn is surrounded by an infinite medium. Each of these three constituents in general possesses elastic constants different from one another. The double-inhomogeneity system under consideration is subjected to far-field strain (stress). Using the equivalent inclusion method (EIM), the double inhomogeneity is replaced by an equivalent double-inclusion (EDI) problem with proper polynomial eigenstrains. The double inclusion is subsequently broken down to single-inclusion problems by means of superposition. The present theory is the first to obtain the actual distribution rather than the averages of the field quantities over the double inhomogeneity using Eshelby’s EIM. The present method is precise and is valid for thin as well as thick layers of coatings, and accommodates eccentric heterogeneity of arbitrary size and orientation. To establish the accuracy and robustness of the present method and for the sake of comparison, results on some of the previously reported problems, which are special cases encompassed by the present theory, will be re-examined. The formulations are easily extended to treat multi-inhomogeneity cases, where an inhomogeneity is surrounded by many layers of coatings. Employing an averaging scheme to the present theory, the average consistency conditions reported by Hori and Nemat-Nasser for the evaluation of average strains and stresses are recovered.


Sign in / Sign up

Export Citation Format

Share Document