An official website of the United States government
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 decreasesJand 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
This content will become publicly available on March 1, 2026
Application of J-Integral to a Random Elastic Medium
Abstract This study investigates the use of the J-integral to compute the statistics of the energy release rate of a random elastic medium. The spatial variability of the elastic modulus is modeled as a homogeneous lognormal random field. Within the framework of Monte Carlo simulation, a modified contour integral is applied to evaluate the first and second statistical moments of the energy release rate. These results are compared with the energy release rate calculated from the potential energy function. The comparison shows that, if the random field of elastic modulus is homogenous in space, the path independence of the classical J-integral remains valid for calculating the mean energy release rate. However, this path independence does not extend to the higher order statistical moments. The simulation further reveals the effect of the correlation length of the spatially varying elastic modulus on the energy release rate of the specimen.
more »
« less
- Award ID(s):
- 2151209
- PAR ID:
- 10621745
- Publisher / Repository:
- ASME
- Date Published:
- Journal Name:
- Journal of Applied Mechanics
- Volume:
- 92
- Issue:
- 3
- ISSN:
- 0021-8936
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In this work, the essential work of fracture (EWF) method is introduced for a peridynamic (PD) material model to characterize fracture toughness of ductile materials. First, an analytical derivation for the path-independence of the PD J -integral is provided. Thereafter, the classical J -integral and PD J -integral are computed on a number of analytical crack problems, for subsequent investigation on how it performs under large scale yielding of thin sheets. To represent a highly nonlinear elastic behavior, a new adaptive bond stiffness calibration and a modified bond-damage model with gradual softening are proposed. The model is employed for two different materials: a lower-ductility bainitic-martensitic steel and a higher-ductility bainitic steel. Up to the start of the softening phase, the PD model recovers the experimentally obtained stress–strain response of both materials. Due to the high failure sensitivity on the presence of defects for the lower-ductility material, the PD model could not recover the experimentally obtained EWF values. For the higher-ductility bainitic material, the PD model was able to match very well the experimentally obtained EWF values. Moreover, the J -integral value obtained from the PD model, at the absolute maximum specimen load, matched the corresponding EWF value.more » « less
-
null (Ed.)The weak localization or enhanced backscattering phenomenon has received a lot of attention in the literature. The enhanced backscattering cone refers to the situation that the wave backscattered by a random medium exhibits an enhanced intensity in a narrow cone around the incoming wave direction. This phenomenon can be analyzed by a formal path integral approach. Here a mathematical derivation of this result is given based on a system of equations that describes the second-order moments of the reflected wave. This system derives from a multiscale stochastic analysis of the wave field in the situation with high-frequency waves and propagation through a lossy medium with fine scale random microstructure. The theory identifies a duality relation between the spreading of the wave and the enhanced backscattering cone. It shows how the cone, its regularity and width relate to the statistical structure of the random medium. We discuss how this information in particular can be used to estimate the internal structure of the random medium based on observations of the reflected wave.more » « less
-
The Pure Shear (PS) crack specimen is widely employed to assess the fracture toughness of soft elastic materials. It serves as a valuable tool for investigating the behavior of crack growth in a steady-state manner following crack initiation. One of its advantages lies in the fact that the energy release rate (J) remains approximately constant for sufficiently long cracks, independent of crack length. Additionally, the PS specimen facilitates the easy evaluation of J for long cracks by means of a tension test conducted on an uncracked sample. However, the lack of a published expression for short cracks currently restricts the usefulness of this specimen. To overcome this limitation, we conducted a series of finite element (FE) simulations utilizing three different constitutive models, namely the neo-Hookean (NH), Arruda-Boyce (AB), and Mooney-Rivlin (MR) models. Our finite element analysis (FEA) encompassed practical crack lengths and strain levels. The results revealed that under a fixed applied displacement, the energy release rate (J) monotonically increases with the crack length for short cracks, reaches a steady-state value when the crack length exceeds the height of the specimen, and subsequently decreases as the crack approaches the end of the specimen. Drawing from these findings, we propose a simple closed-form expression for J that can be applied to most hyper-elastic models and is suitable for all practical crack lengths, particularly short cracks.more » « less
-
A bstract It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p -spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.more » « less
