skip to main content

Explore Research Products in the NSF-PAR

Title: Elastic Solution of a Polygon-Shaped Inclusion With a Polynomial Eigenstrain
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.  more » « less
Award ID(s):
1762891 1738802
NSF-PAR ID:
10276154
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Applied Mechanics
Volume:
88
Issue:
6
ISSN:
0021-8936
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Journal of Applied Mechanics)
    null (Ed.)
    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. 
    more » « less
  2. ; ; ( , Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)
    This paper investigates the mechanical behaviour of a bi-layered panel containing many particles in one layer and demonstrates the size effect of particles on the deflection. The inclusion-based boundary element method (iBEM) considers a fully bounded bi-material system. The fundamental solution for two-jointed half spaces has been used to acquire elastic fields resulting from source fields over inclusions and boundary-avoiding multi-domain integral along the interface. Eshelby’s equivalent inclusion method is used to simulate the material mismatch with a continuously distributed eigenstrain field over the equivalent inclusion. The eigenstrain is expanded at the centre of the inclusion, which provides tailorable accuracy based on the order of the polynomial of the eigenstrain. As a single-domain approach, the iBEM algorithm is particularly suitable for conducting virtual experiments of bi-layered composites with many defects or reinforcements for both local analysis and homogenization purposes. The maximum deflection of solar panel coupons is studied under uniform vertical loading merged with inhomogeneities of different material properties, dimensions and volume fractions. The size of defects or reinforcements plays a significant role in the deflection of the panel, even with the same volume fraction, as the substrate is relatively thin. 
    more » « less
  3. ; ; ( , IMA Journal of Numerical Analysis)
    Abstract In this paper we study the biharmonic equation with Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem as a system of Poisson equations. Our method differs from the naive mixed method that leads to two Poisson problems but only applies to convex domains; our decomposition involves a third Poisson equation to confine the solution in the correct function space, and therefore can be used in both convex and nonconvex domains. A $C^0$ finite element algorithm is in turn proposed to solve the resulting system. In addition, we derive optimal error estimates for the numerical solution on both quasi-uniform meshes and graded meshes. Numerical test results are presented to justify the theoretical findings. 
    more » « less
  4. ; ; ; ; ( , ESAIM: Mathematical Modelling and Numerical Analysis)
    We consider an unconstrained tangential Dirichlet boundary control problem for the Stokes equations with an $ L^2 $ penalty on the boundary control.  The contribution of this paper is twofold.  First, we obtain well-posedness and regularity results for the tangential Dirichlet control problem on a convex polygonal domain.  The analysis contains new features not found in similar Dirichlet control problems for the Poisson equation; an interesting result is that the optimal control has higher local regularity on the individual edges of the domain compared to the global regularity on the entire boundary.  Second, we propose and analyze a hybridizable discontinuous Galerkin (HDG) method to approximate the solution.  For convex polygonal domains, our theoretical convergence rate for the control is optimal with respect to the global regularity on the entire boundary.  We present numerical experiments to demonstrate the performance of the HDG method. 
    more » « less
  5. ; ; ; ( , Journal of Applied Mechanics)
    Abstract The Green’s function of a bimaterial infinite domain with a plane interface is applied to thermal analysis of a spherical underground heat storage tank. The heat transfer from a spherical source is derived from the integral of the Green’s function over the spherical domain. Because the thermal conductivity of the tank is generally different from soil, the Eshelby’s equivalent inclusion method (EIM) is used to simulate the thermal conductivity mismatch of the tank from the soil. For simplicity, the ground with an approximately uniform temperature on the surface is simulated by a bimaterial infinite domain, which is perfectly conductive above the ground. The heat conduction in the ground is investigated for two scenarios: First, a steady-state uniform heat flux from surface into the ground is considered, and the heat flux is disturbed by the existence of the tank due to the conductivity mismatch. A prescribed temperature gradient, or an eigen-temperature gradient, is introduced to investigate the local temperature field in the neighborhood of the tank. Second, when a temperature difference exists between the water in the tank and soil, the heat transfer between the tank and soil depends on the tank size, conductivity, and temperature difference, which provide a guideline for heat exchange design for the tank size. The modeling framework can be extended to two-dimensional cases, periodic, or transient heat transfer problems for geothermal well operations. The corresponding Green’s functions are provided for those applications. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email