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
This content will become publicly available on November 25, 2024
The effects of boundary and inhomogeneities on the delamination of a bi-layered material system
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- NSF-PAR ID:
- 10475615
- Publisher / Repository:
- SAGE Publications
- Date Published:
- Journal Name:
- International Journal of Damage Mechanics
- ISSN:
- 1056-7895
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Hydrogels are a class of soft, highly deformable materials formed by swelling a network of polymer chains in water. With mechanical properties that mimic biological materials, hydrogels are often proposed for load bearing biomedical or other applications in which their deformation and failure properties will be important. To study the failure of such materials a means for the measurement of deformation fields beyond simple uniaxial tension tests is required. As a non-contact, full-field deformation measurement method, Digital Image Correlation (DIC) is a good candidate for such studies. The application of DIC to hydrogels is studied here with the goal of establishing the accuracy of DIC when applied to hydrogels in the presence of large strains and large strain gradients. Experimental details such as how to form a durable speckle pattern on a material that is 90% water are discussed. DIC is used to measure the strain field in tension loaded samples containing a central hole, a circular edge notch and a sharp crack. Using a nonlinear, large deformation constitutive model, these experiments are modeled using the finite element method (FEM). Excellent agreement between FEM and DIC results for all three geometries shows that the DIC measurements are accurate up to strains of over 10, even in the presence of very high strain gradients near a crack tip. The method is then applied to verify a theoretical prediction that the deformation field in a cracked sample under relaxation loading, i.e. constant applied boundary displacement, is stationary in time even as the stress relaxes by a factor of three.more » « less
-
more » « less
Summary The discrete crack mechanics (DCM) method is a dislocation‐based crack modeling technique where cracks are constructed using Volterra dislocation loops. The method allows for the natural introduction of displacement discontinuities, avoiding numerically expensive techniques. Mesh dependence in existing computational modeling of crack growth is eliminated by utilizing a superposition procedure. The elastic field of cracks in finite bodies is separated into two parts: the infinite‐medium solution of discrete dislocations and an finite element method solution of a correction problem that satisfies external boundary conditions. In the DCM, a crack is represented by a dislocation array with a fixed outer loop determining the crack tip position encompassing additional concentric loops free to expand or contract. Solving for the equilibrium positions of the inner loops gives the crack shape and stress field. The equation of motion governing the crack tip is developed for quasi‐static growth problems. Convergence and accuracy of the DCM method are verified with two‐ and three‐dimensional problems with well‐known solutions. Crack growth is simulated under load and displacement (rotation) control. In the latter case, a semicircular surface crack in a bent prismatic beam is shown to change shape as it propagates inward, stopping as the imposed rotation is accommodated.
-
This article presents a numerical strategy for actively manipulating electromagnetic (EM) fields in layered media. In particular, we develop a scheme to characterize an EM source that will generate some predetermined field patterns in prescribed disjoint exterior regions in layered media. The proposed question of specifying such an EM source is not an inverse source problem (ISP) since the existence of a solution is not guaranteed. Moreover, our problem allows for the possibility of prescribing different EM fields in mutually disjoint exterior regions. This question involves a linear inverse problem that requires solving a severely ill-posed optimization problem (i.e. suffering from possible non-existence or non-uniqueness of a solution). The forward operator is defined by expressing the EM fields as a function of the current at the source using the layered media Green’s function (LMGF), accounting for the physical parameters of the layered media. This results to integral equations that are then discretized using the method of moments (MoM), yielding an illposed system of linear equations. Unlike in ISPs, stability with respect to data is not an issue here since no data is measured. Rather, stability with respect to input current approximation is important. To get such stable solutions, we applied two regularization methods, namely, the truncated singular value decomposition (TSVD) method and the Tikhonov regularization method with the Morozov Discrepancy Principle. We performed several numerical simulations to support the theoretical framework and analyzes, and to demonstrate the accuracy and feasibility of the proposed numerical algorithms.more » « less
