skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 10 until 2:00 AM ET on Friday, October 11 due to maintenance. We apologize for the inconvenience.


Title: Controlling the fracture response of structures via topology optimization: From delaying fracture nucleation to maximizing toughness
It is now a well-established fact that even simple topology variations can drastically change the fracture response of structures. With the objective of gaining quantitative insight into this phenomenon, this paper puts forth a density-based topology optimization framework for the fracture response of structures subjected to quasistatic mechanical loads. One of the two key features of the proposed framework is that it makes use of a complete phase-field fracture theory that has been recently shown capable of accurately describing the nucleation and propagation of brittle fracture in a wide range of nominally elastic materials under a wide range of loading conditions. The other key feature is that the framework is based on a multi-objective function that allows optimizing in a weighted manner: ( ) the initial stiffness of the structure, ( ) the first instance at which fracture nucleates, and ( ) the energy dissipated by fracture propagation once fracture nucleation has occurred. The focus is on the basic case of structures made of a single homogeneous material featuring an isotropic linear elastic behavior alongside an isotropic strength surface and toughness. Novel interpolation rules are proposed for each of these three types of material properties. As a first effort to gain quantitative insight, the framework is deployed to optimize the fracture response of 2D structures wherein the fracture is bound to nucleate in three different types of regions: within the bulk, from geometric singularities (pre-existing cracks and sharp corners), and from smooth parts of the boundary. The obtained optimized structures are shown to exhibit significantly enhanced fracture behaviors compared to those of structures that are optimized according to conventional stiffness maximization. Furthermore, the results serve to reveal a variety of strengthening and toughening mechanisms. These include the promotion of highly porous structures, the formation of tension-compression asymmetric regions, and the removal of cracks and sharp corners. The particular mechanism that is preferred by a given structure, not surprisingly, correlates directly to the elastic, strength, and toughness properties of the material that is made of.  more » « less
Award ID(s):
2127134
NSF-PAR ID:
10492195
Author(s) / Creator(s):
; ;
Publisher / Repository:
ScienceDirect
Date Published:
Journal Name:
Journal of the Mechanics and Physics of Solids
Volume:
173
Issue:
C
ISSN:
0022-5096
Page Range / eLocation ID:
105227
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. An experimental and theoretical study of delayed fracture of polydimethlsiloxane (PDMS) is presented. Previous works have demonstrated that delayed fracture in single edge notch specimens is caused by time dependent damage due to chain scission. Here we study the interactions between damage and the elastic field using different specimens and crack geometries with blunt and sharp cracks. Our experiments show that initial toughness is not well defined, as stable slow crack growth can occur over a range of applied loads. Our experiments demonstrate that there is a universal relation between crack growth rate and applied energy release rate. A model coupling the nonlinear elastic deformation and rate dependent bond scission is proposed and is in good agreement with experimental data. 
    more » « less
  2. Anisotropy in additive manufacturing (AM), particularly in the material extrusion process, plays a crucial role in determining the actual structural performance, including the stiffness and strength of the printed parts. Unless accounted for, anisotropy can compromise the objective performance of topology-optimized structures and allow premature failures for stress-sensitive design domains. This study harnesses process-induced anisotropy in material extrusion-based 3D printing to design and fabricate stiff, strong, and lightweight structures using a two-step framework. First, an AM-oriented anisotropic strength-based topology optimization formulation optimizes the structural geometry and infill orientations, while assuming both anisotropic (i.e., transversely isotropic) and isotropic infill types as candidate material phases. The dissimilar stiffness and strength interpolation schemes in the formulation allow for the optimized allocation of anisotropic and isotropic material phases in the design domain while satisfying their respective Tsai–Wu and von Mises stress constraints. Second, a suitable fabrication methodology realizes anisotropic and isotropic material phases with appropriate infill density, controlled print path (i.e., infill directions), and strong interfaces of dissimilar material phases. Experimental investigations show up to 37% improved stiffness and 100% improved strength per mass for the optimized and fabricated structures. The anisotropic strength-based optimization improves load-carrying capacity by simultaneous infill alignment along the stress paths and topological adaptation in response to high stress concentration. The adopted interface fabrication methodology strengthens comparatively weaker anisotropic joints with minimal additional material usage and multi-axial infill patterns. Furthermore, numerically predicted failure locations agree with experimental observations. The demonstrated framework is general and can potentially be adopted for other additive manufacturing processes that exhibit anisotropy, such as fiber composites. 
    more » « less
  3. Abstract

    Brittle fracture propagation in rocks is a complex process due to significant grain‐scale heterogeneity and evolving stress states under dynamic loading conditions. In this work, we use digital image correlation and linear elastic fracture mechanics to make instantaneous measurements of the opening (mode I) and in plane shear (mode II) components of the stress intensity field during dynamic mixed mode crack initiation and propagation in crystalline and granular rocks. Both rock types display some similar fracture behaviors as observed in engineered materials, including rate dependent fracture initiation toughness and a direct relationship between propagation toughness and crack velocity; however, measured propagation toughness is higher than quasi‐static values at crack velocities well below the branching velocity in both rocks. Additionally, due to grain scale controls on the fracture process, mixed mode crack propagation is fundamentally different between these two rock types. Mixed mode propagation is energetically more favorable than pure opening mode propagation in sandstone, while the opposite is true in granite. Furthermore, following initiation, propagation in granite occurs so as to minimize the mode II contribution, irrespective of the initiation conditions, while fractures in sandstone maintain a non‐negligible mode II contribution during propagation across the sample.

     
    more » « less
  4. 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
  5. Abstract Fracture dictates the service limits of metallic structures. Damage tolerance of materials may be characterized by fracture toughness rigorously developed from fracture mechanics, or less rigorous yet more easily obtained impact toughness (or impact energy as a variant). Given the promise of high-entropy alloys (HEAs) in structural and damage-tolerance applications, we compiled a dataset of fracture toughness and impact toughness/energy from the literature till the end of the 2022 calendar year. The dataset is subdivided into three categories, i.e., fracture toughness, impact toughness, and impact energy, which contain 153, 14, and 78 distinct data records, respectively. On top of the alloy chemistry and measured fracture quantities, each data record also documents the factors influential to fracture. Examples are material-processing history, phase structures, grain sizes, uniaxial tensile properties, such as yield strength and elongation, and testing conditions. Data records with comparable conditions are graphically visualized by plots. The dataset is hosted in Materials Cloud, an open data repository. 
    more » « less