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.
- NSF-PAR ID:
- 10418018
- Date Published:
- Journal Name:
- Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Volume:
- 479
- Issue:
- 2271
- ISSN:
- 1364-5021
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
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
-
null (Ed.)more » « less
Abstract 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.
-
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
-
null (Ed.)During this study, full-size wood composite sandwich panels, 1.2 m by 2.4 m (4 ft by 8 ft), with a biaxial corrugated core were evaluated as a building construction material. Considering the applications of this new building material, including roof, floor, and wall paneling, sandwich panels with one and two corrugated core(s) were fabricated and experimentally evaluated. Since primary loads applied on these sandwich panels during their service life are live load, snow load, wind, and gravity loads, their bending and compression behavior were investigated. To improve the thermal characteristics, the cavities within the sandwich panels created by the corrugated geometry of the core were filled with a closed-cell foam. The R-values of the sandwich panels were measured to evaluate their energy performance. Comparison of the weight indicated that fabrication of a corrugated panel needs 74% less strands and, as a result, less resin compared to a strand-based composite panel, such as oriented strand board (OSB), of the same size and same density. Bending results revealed that one-layer core sandwich panels with floor applications under a 4.79 kPa (100 psf) bending load are able to meet the smallest deflection limit of L/360 when the span length (L) is 137.16 cm (54 in) or less. The ultimate capacity of two-layered core sandwich panels as a wall member was 94% and 158% higher than the traditional walls with studs under bending and axial compressive loads, respectively. Two-layered core sandwich panels also showed a higher ultimate capacity compared to structural insulated panels (SIP), at 470% and 235% more in bending and axial compression, respectively. Furthermore, normalized R-values, the thermal resistance, of these sandwich panels, even with the presence of thermal bridging due to the core geometry, was about 114% and 109% higher than plywood and oriented strand board, respectively.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1098/rspa.2022.0754
