skip to main content

Title: Effect of defects and reinforcements on the mechanical behaviour of a bi-layered panel
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
Award ID(s):
1941244 1738802 1762891
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The inclusion-based boundary element method (iBEM) is developed to calculate the elastic fields of a bi-layered composite with inhomogeneities in one layer. The bi-material Green’s function has been applied to obtain the elastic field caused by the domain integral of the source fields on inclusions and the boundary integral of the applied loads on the surface. Using Eshelby’s equivalent inclusion method (EIM), the material mismatch between the particle and matrix phases is simulated with a continuously distributed source field, namely eigenstrain, on inhomogeneities so that the iBEM can calculate the local field. The stress singularity along the interface leads to the delamination of the bimaterials under a certain load. The crack’s energy release rate (  J) is obtained through the J-integral, which predicts the stability of the delamination. When the stiffness of one layer increases, the J-integral increases with a higher gradient, leading to lower stability. Particularly, the effect of the boundary and inhomogeneity on the J-integral is illustrated by changing the crack length and inhomogeneity configuration, which shows the crack is stable at the beginning stage and becomes unstable when the crack tip approaches the boundary; a stiffer inhomogeneity in the neighborhood of a crack tip decreases J and improves the fracture resistance. For the stable cracking phase, the J-integral increases with the volume fraction of inhomogeneity are evaluated. The model is applied to a dual-glass solar module with air bubbles in the encapsulant layer. The stress distribution is evaluated with the iBEM, and the J-integral is evaluated to predict the delamination process with the energy release rate, which shows that the bubbles significantly increase the J-integral. The effect of the bubble size, location, and number on the J-integral is also investigated. The present method provides a powerful tool for the design and analysis of layered materials and structures.

    more » « less
  2. null (Ed.)
    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
  3. 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
  4. We discuss a novel approach for imaging local faults inside an infinite bi-periodic layered medium in ℝ3 using acoustic measurements of scattered fields at the bottom or the top of the layer. The faulted area is represented by compactly supported perturbations with erroneous material properties. Our method reconstructs the support of perturbations without knowing or reconstructing the constitutive material parameters of healthy or faulty bi-period layer; only the size of the period is needed. This approach falls under the class of non-iterative imaging methods, known as the generalized linear sampling method with differential measurements, first introduced in [2] and adapted to periodic layers in [25]. The advantage of applying differential measurements to our inverse problem is that instead of comparing the measured data against measurements due to healthy structures, one makes use of periodicity of the layer where the data operator restricted to single Floquet-Bloch modes plays the role of the one corresponding to healthy material. This leads to a computationally efficient and mathematically rigorous reconstruction algorithm. We present numerical experiments that confirm the viability of the approach for various configurations of defects 
    more » « less
  5. null (Ed.)

    Composites can be tailored to specific applications by adjusting process variables. These variables include those related to composition, such as volume fraction of the constituents and those associated with processing methods, methods that can affect composite topology. In the case of particle matrix composites, orientation of the inclusions affects the resulting composite properties, particularly so in instances where the particles can be oriented and arranged into structures. In this work, we study the effects of coupled electric and magnetic field processing with externally applied fields on those structures, and consequently on the resulting material properties that arise. The ability to vary these processing conditions with the goal of generating microstructures that yield target material properties adds an additional level of control to the design of composite material properties. Moreover, while analytical models allow for the prediction of resulting composite properties from constituents and composite topology, these models do not build upward from process variables to make these predictions.

    This work couples simulation of the formation of microscale architectures, which result from coupled electric and magnetic field processing of particulate filled polymer matrix composites, with finite element analysis of those structures to provide a direct and explicit linkages between process, structure, and properties. This work demonstrates the utility of these method as a tool for determining composite properties from constituent and processing parameters. Initial particle dynamics simulation incorporating electromagnetic responses between particles and between the particles and the applied fields, including dielectrophoresis, are used to stochastically generate representative volume elements for a given set of process variables. Next, these RVEs are analyzed as periodic structures using FEA yielding bulk material properties. The results are shown to converge for simulation size and discretization, validating the RVE as an appropriate representation of the composite volume. Calculated material properties are compared to traditional effective medium theory models. Simulations allow for mapping of composite properties with respect to not only composition, but also fundamentally from processing simulations that yield varying particle configurations, a step not present in traditional or more modern effective medium theories such as the Halpin Tsai or double-inclusion theories.

    more » « less