ABSTRACT A recently developed Timoshenko‐based peridynamic model with a variable micropolar shear influence factor is extended to study the behavior of dynamic crack propagation in functionally graded materials (FGMs). To this end, first, the proposed model is validated against two experimental three‐point bending benchmark problems with different material functions as well as varying loading rates and durations. Then, numerous additional cases with different boundary conditions and material distribution are studied to predict crack initiation and propagation in such mediums. The examples consist of three‐point bending and Kalthoff–Winkler specimens with various material functions under dynamic loads. Finally, the effects of material anisotropy induced by functionally varying material properties on crack propagation path are addressed. It is shown that this new model is advantageous because of its capability to account for shear deformation effects in the bonds previously ignored by the original bond‐based peridynamic models. Moreover, comparing the proposed modified bond‐based model to more complex methods, such as state‐based peridynamics, reveals that the simplicity of the current approach results in lower computational costs while still achieving comparable results.
more »
« less
A modified mixed-mode Timoshenko-based peridynamics model considering shear deformation
The original two-dimensional bond-based peridynamic (BBPD) framework, which only considers the pairwise forces (compression and tension) between two material points, is extended by incorporating the effect of shear deformation in the calculations and its influence on the failure of the bonds. To this end, each bond is considered as a short Timoshenko beam, and by doing so, the traditional BBPD is enhanced into a more comprehensive model known as multi-polar peridynamic (MPPD). The proposed novel approach explicitly considers the shear influence factor used in Timoshenko beams and introduces a strain-based shear deformation failure criterion. The model is then validated against two benchmark experimental tests (i.e., a standard pure mode I edge crack, and a Kalthoff-Winkler configuration) reported in the literature under in-plane dynamic loading and plane stress conditions. In most cases, the developed model is shown to be more accurate in predicting the crack paths obtained from the experimental results when compared to other theoretical methods delineated in the literature. Furthermore, a noticeable change in crack branching and crack path is observed in a study on the effects of Poisson’s ratio and the loading rate. This investigation also demonstrated that the proposed MPPD model can accommodate materials with Poisson’s ratios up to 1/3, expanding the range beyond the traditional BBPD limitations.
more »
« less
- Award ID(s):
- 2317406
- PAR ID:
- 10553178
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- International Journal of Mechanical Sciences
- Volume:
- 285
- Issue:
- C
- ISSN:
- 0020-7403
- Page Range / eLocation ID:
- 109802
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
For the first time, the mixed-mode dynamic fracture in anisotropic functionally varying microcellular structures is investigated herein. To this end, a recently developed homogenization MATLAB implementation capable of considering material and geometry-induced anisotropy is used, and a continuous medium with equivalent functionality distributed mechanical properties to the original microcellular domain is obtained. Then, the resulting material domain is subjected to dynamic loads, and the crack propagation is predicted by using a novel Timoshenko-based peridynamic model. This innovative method unprecedentedly accounts for a bond-length dependent shear influence factor and a shear strain-based failure criterion. Finally, numerous cases consisting of compact-tension (CT) and Kalthoff-Winkler specimens with several void sizes, shapes, and distribution patterns are numerically solved. The results demonstrate that the crack path is significantly influenced by the void distribution pattern near the crack tip, providing a foundation for engineering crack propagation to prevent it from reaching critical areas of a structure.more » « less
-
Abstract The recently conceived gap test and its simulation revealed that the fracture energy Gf (or Kc, Jcr) of concrete, plastic-hardening metals, composites, and probably most materials can change by ±100%, depending on the crack-parallel stresses σxx, σzz, and their history. Therefore, one must consider not only a finite length but also a finite width of the fracture process zone, along with its tensorial damage behavior. The data from this test, along with ten other classical tests important for fracture problems (nine on concrete, one on sandstone), are optimally fitted to evaluate the performance of the state-of-art phase-field, peridynamic, and crack band models. Thanks to its realistic boundary and crack-face conditions as well as its tensorial nature, the crack band model, combined with the microplane damage constitutive law in its latest version M7, is found to fit all data well. On the contrary, the phase-field models perform poorly. Peridynamic models (both bond based and state based) perform even worse. The recent correction in the bond-associated deformation gradient helps to improve the predictions in some experiments, but not all. This confirms the previous strictly theoretical critique (JAM 2016), which showed that peridynamics of all kinds suffers from several conceptual faults: (1) It implies a lattice microstructure; (2) its particle–skipping interactions are a fiction; (4) it ignores shear-resisted particle rotations (which are what lends the lattice discrete particle model (LDPM) its superior performance); (3) its representation of the boundaries, especially the crack and fracture process zone faces, is physically unrealistic; and (5) it cannot reproduce the transitional size effect—a quintessential characteristic of quasibrittleness. The misleading practice of “verifying” a model with only one or two simple tests matchable by many different models, or showcasing an ad hoc improvement for one type of test while ignoring misfits of others, is pointed out. In closing, the ubiquity of crack-parallel stresses in practical problems of concrete, shale, fiber composites, plastic-hardening metals, and materials on submicrometer scale is emphasized.more » « less
-
ABSTRACT:Due to rock mass being commonly subjected to compressive or shear loading, the mode II fracture toughness is an important material parameter for rocks. Fracturing in rocks is governed by the behavior of a nonlinear region surrounding the crack tip called the fracture process zone (FPZ). However, the characteristics of mode II fracture are still determined based on the linear elastic fracture mechanics (LEFM), which assumes that a pure mode II loading results in a pure mode II fracture. In this study, the FPZ development in Barre granite specimens under mode II loading was investigated using the short beam compression (SBC) test. Additionally, the influence of lateral confinement on various characteristics of mode II fracture was studied. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique (2D-DIC) to identify fracture mode and characterize the FPZ evolution in Barre granite specimens. The 2D-DIC analysis showed a dominant mixed-mode I/II fracture in the ligament between two notches, irrespective of confinement level on the SBC specimens. The influence of confinement on the SBC specimens was assessed by analyzing the evolution of crack displacement and changes in value of mode II fracture toughness. Larger levels of damage in confined specimens were observed prior to the failure than the unconfined specimens, indicating an increase in the fracture resistance and therefore mode II fracture toughness with the confining stress. 1. INTRODUCTIONThe fracturing in laboratory-scale rock specimens is often characterized by the deformation of the inelastic region surrounding the crack tips, also known as the fracture process zone (FPZ) (Backers et al., 2005; Ghamgosar and Erarslan, 2016). While the influence of the FPZ on mode I fracture in rocks has been extensively investigated, there are limited studies on FPZ development in rocks under pure mode II loading (Ji et al., 2016; Lin et al., 2020; Garg et al., 2021; Li et al., 2021).more » « less
-
A phase-field model for thermomechanically-induced fracture in NiTi at the single crystal level, i.e., fracture under loading paths that may take advantage of either of the functional properties of NiTi–superelasticity or shape memory effect–, is presented, formulated within the kinematically linear regime. The model accounts for reversible phase transformation from austenite to martensite habit plane variants and plastic deformation in the austenite phase. Transformation-induced plastic deformation is viewed as a mechanism for accommodation of the local deformation incompatibility at the austenite–martensite interfaces and is accounted for by introducing an interaction term in the free energy derived based on the Mori–Tanaka and Kröner micromechanical assumptions and the hypothesis of martensite instantaneous growth within austenite. Based on experimental observations suggesting that NiTi fractures in a stress-controlled manner, damage is assumed to be driven by the elastic energy, i.e., phase transformation and plastic deformation are assumed to contribute in crack formation and growth indirectly through stress redistribution. The model is restricted to quasistatic mechanical loading (no latent heat effects), thermal loading sufficiently slow with respect to the time rate of heat transfer by conduction (no thermal gradients), and a temperature range below 𝑀𝑑, which is the temperature above which the austenite phase is stable, i.e., stress-induced martensitic transformation is suppressed. The numerical implementation of the model is based on an efficient scheme of viscous regularization in both phase transformation and plastic deformation, an explicit numerical integration via a tangent modulus method, and a staggered scheme for the coupling of the unknown fields. The model is shown able to capture transformation-induced toughening, i.e., stable crack advance attributed to the shielding effect of inelastic deformation left in the wake of the growing crack under nominal isothermal loading, actuation-induced fracture under a constant bias load, and crystallographic dependence on crack pattern.more » « less
An official website of the United States government

