An official website of the United States government
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
A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials
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
- Award ID(s):
- 1911836
- PAR ID:
- 10494470
- Editor(s):
- -
- Publisher / Repository:
- Comput. Methods Appl. Mech. Engrg.
- Date Published:
- Journal Name:
- Computer Methods in Applied Mechanics and Engineering
- Volume:
- 410
- Issue:
- C
- ISSN:
- 0045-7825
- Page Range / eLocation ID:
- 115982
- Subject(s) / Keyword(s):
- Phase field fracture approach Non-affine microsphere model Multiscale model Fracture propagation Rubber-like materials
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Molecular dynamics simulations are used to investigate the conformations of a single polymer chain, represented by the Kremer-Grest bead-spring model, in a solution with a Lennard-Jones liquid as the solvent when the interaction strength between the polymer and solvent is varied. Results show that when the polymer-solvent interaction is unfavorable, the chain collapses as one would expect in a poor solvent. For more attractive polymer-solvent interactions, the solvent quality improves and the chain is increasingly solvated and exhibits ideal and then swollen conformations. However, as the polymer-solvent interaction strength is increased further to be more than about twice the strength of the polymer-polymer and solvent-solvent interactions, the chain exhibits an unexpected collapsing behavior. Correspondingly, for strong polymer-solvent attractions, phase separation is observed in the solutions of multiple chains. These results indicate that the solvent becomes effectively poor again at very attractive polymer-solvent interactions. Nonetheless, the mechanism of chain collapsing and phase separation in this limit differs from the case with a poor solvent rendered by unfavorable polymer-solvent interactions. In the latter, the solvent is excluded from the domain of the collapsed chains while in the former, the solvent is still present in the pervaded volume of a collapsed chain or in the polymer-rich domain that phase separates from the pure solvent. In the limit of strong polymer-solvent attractions, the solvent behaves as a glue to stick monomers together, causing a single chain to collapse and multiple chains to aggregate and phase separate.more » « less
-
ABSTRACT Strongly anisotropic geomaterials, such as layered shales, have been observed to undergo fracture under compressive loading. This paper applies a phase‐field fracture model to study this fracture process. While phase‐field fracture models have several advantages—primarily that the fracture path is not predetermined but arises naturally from the evolution of a smooth non‐singular damage field—they provide unphysical predictions when the stress state is complex and includes compression that can cause crack faces to contact. Building on a recently developed phase‐field model that accounts for compressive traction across the crack face, this paper extends the model to the setting of anisotropic fracture. The key features of the model include the following: (1) a homogenized anisotropic elastic response and strongly anisotropic model for the work to fracture; (2) an effective damage response that accounts consistently for compressive traction across the crack face, that is derived from the anisotropic elastic response; (3) a regularized crack normal field that overcomes the shortcomings of the isotropic setting, and enables the correct crack response, both across and transverse to the crack face. To test the model, we first compare the predictions to phase‐field fracture evolution calculations in a fully resolved layered specimen with spatial inhomogeneity, and show that it captures the overall patterns of crack growth. We then apply the model to previously reported experimental observations of fracture evolution in laboratory specimens of shales under compression with confinement, and find that it predicts well the observed crack patterns in a broad range of loading conditions. We further apply the model to predict the growth of wing cracks under compression and confinement. Prior approaches to simulate wing cracks have treated the initial cracks as an external boundary, which makes them difficult to apply to general settings. Here, the effective crack response model enables us to treat the initial crack simply as a nonsingular damaged zone within the computational domain, thereby allowing for easy and general computations.more » « less
-
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
-
Histogram-reweighting grand canonical Monte Carlo simulations are used to obtain the critical properties of lattice chains composed of solvophilic and solvophobic monomers. The model is a modification of one proposed by Larson et al. [J. Chem. Phys. 83, 2411 (1985)], lowering the “contrast” between beads of different types to prevent aggregation into finite-size micelles that would mask true phase separation between bulk high- and low-density phases. Oligomeric chains of lengths between 5 and 24 beads are studied. Mixed-field finite-size scaling methods are used to obtain the critical properties with typical relative accuracies of better than 10−4 for the critical temperature and 10−3 for the critical volume fraction. Diblock chains are found to have lower critical temperatures and volume fractions relative to the corresponding homopolymers. The addition of solvophilic blocks of increasing length to a fixed-length solvophobic segment results in a decrease of both the critical temperature and the critical volume fraction, with an eventual slow asymptotic approach to the long-chain limiting behavior. Moving a single solvophobic or solvophilic bead along a chain leads to a minimum or maximum in the critical temperature, with no change in the critical volume fraction. Chains of identical length and composition have a significant spread in their critical properties, depending on their precise sequence. The present study has implications for understanding biomolecular phase separation and for developing design rules for synthetic polymers with specific phase separation properties. It also provides data potentially useful for the further development of theoretical models for polymer and surfactant phase behavior.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.cma.2023.115982
