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
Oscillatory and tip-splitting instabilities in 2D dynamic fracture: The roles of intrinsic material length and time scales
Recent theoretical and computational progress has led to unprecedented understanding of symmetry-breaking instabilities in 2D dynamic fracture. At the heart of this progress resides the identification of two intrinsic, near crack tip length scales — a nonlinear elastic length scale ℓ and a dissipation length scale ξ — that do not exist in Linear Elastic Fracture Mechanics (LEFM), the classical theory of cracks. In particular, it has been shown that at a propagation velocity v of about 90% of the shear wave-speed, cracks in 2D brittle materials undergo an oscillatory instability whose wavelength varies linearly with ℓ, and at larger loading levels (corresponding to yet higher propagation velocities), a tip-splitting instability emerges, both in agreements with experiments. In this paper, using phase-field models of brittle fracture, we demonstrate the following properties of the oscillatory instability: (i) It exists also in the absence of near-tip elastic nonlinearity, i.e. in the limit ℓ→0, with a wavelength determined by the dissipation length scale ξ. This result shows that the instability crucially depends on the existence of an intrinsic length scale associated with the breakdown of linear elasticity near crack tips, independently of whether the latter is related to nonlinear elasticity or to dissipation. (ii) It is a supercritical Hopf bifurcation, featuring a vanishing oscillations amplitude at onset. (iii) It is largely independent of the phenomenological forms of the degradation functions assumed in the phase-field framework to describe the cohesive zone, and of the velocity-dependence of the fracture energy Γ(v) that is controlled by the dissipation time scale in the Ginzburg-Landau-type evolution equation for the phase-field. These results substantiate the universal nature of the oscillatory instability in 2D. In addition, we provide evidence indicating that the tip-splitting instability is controlled by the limiting rate of elastic energy transport inside the crack tip region. The latter is sensitive to the wave-speed inside the dissipation zone, which can be systematically varied within the phase-field approach. Finally, we describe in detail the numerical implementation scheme of the employed phase-field fracture approach, allowing its application in a broad range of materials failure problems.
more »
« less
- Award ID(s):
- 1827343
- NSF-PAR ID:
- 10216227
- Date Published:
- Journal Name:
- Journal of the mechanics and physics of solids
- ISSN:
- 0022-5096
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Double-network gels are a class of tough soft materials comprising two elastic networks with contrasting structures. The formation of a large internal damage zone ahead of the crack tip by the rupturing of the brittle network accounts for the large crack resistance of the materials. Understanding what determines the damage zone is the central question of the fracture mechanics of double-network gels. In this work, we found that at the onset of crack propagation, the size of necking zone, in which the brittle network breaks into fragments and the stretchable network is highly stretched, distinctly decreases with the increase of the solvent viscosity, resulting in a reduction in the fracture toughness of the material. This is in sharp contrast to the tensile behavior of the material that does not change with the solvent viscosity. This result suggests that the dynamics of stretchable network strands, triggered by the rupture of the brittle network, plays a role. To account for this solvent viscosity effect on the crack initiation, a delayed blunting mechanism regarding the polymer dynamics effect is proposed. The discovery on the role of the polymer dynamic adds an important missing piece to the fracture mechanism of this unique material.more » « less
-
- (Ed.)Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.more » « less
-
Prior to fracture, a polyacrylamide hydrogel has a stress-stretch curve of nearly perfect elasticity, but it has been suggested that an inelastic zone exists around a crack tip. This inelastic zone, however, has never been observed directly in a polyacrylamide hydrogel. Here we identify the inelastic zone using digital image correlation (DIC). We prepare a polyacrylamide hydrogel with a precut crack. While a sample of the hydrogel is stretched, the speckle patterns are recorded using a microscope or a camera, with pixel size 2.3 μm and 22.7 μm, respectively. The speckle patterns recorded by the microscope and camera are processed using the DIC software, and merged to provide the deformation field over the entire sample. The measured field of deformation is used to calculate the field of energy density according to the neo-Hookean model. When the body is perfectly elastic, the field of energy density around the crack tip is inversely proportional to the distance from the crack tip. The difference between the measured field and the predicted elastic field identifies the inelastic zone. The measured size of the inelastic zone is ∼ 0.6 mm. We further confirm that, when a sample is much larger than the inelastic zone, an annulus exists, in which the elastic crack tip field prevails.more » « less
-
Energy dissipation around a propagating crack is the primary mechanism for the enhanced fracture toughness in viscoelastic solids. Such dissipation is spatially non-uniform and is highly coupled to the crack propagation process due to the history-dependent nature of viscoelasticity. We present an experimental approach to map the dissipation field during crack propagation in soft viscoelastic solid. Specifically, we track randomly distributed tracer particles to measure the evolving deformation field. The measured deformation field is then put into a nonlinear constitutive model to determine the dissipation field. Our methodology was used to investigate the deformation and dissipation fields around a propagating crack in a Polyampholyte (PA) hydrogel. The deformation field measurements allowed us to assess whether the commonly assumed translational invariance in viscoelastic fracture theories holds true in practical experiments. Furthermore, by combining the obtained deformation fields with a nonlinear viscoelastic model, we captured the complete history of the dissipation field during crack propagation. We found that dissipation occurred even at material points that are a few millimeters away from the crack tip. The mapped dissipation field also enabled the separate determination of the intrinsic and dissipative components of fracture toughness for the viscoelastic hydrogel.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/104372
