An official website of the United States government
Toughness of soft materials such as elastomers and gels depends on their ability to dissipate energy and to reduce stress concentration at the crack tip. The primary energy dissipation mechanism is viscoelasticity. Most analyses and models of fracture are based on linear viscoelastic theory (LVT) where strains are assumed to be small and the relaxation mechanisms are independent of stress or strain history. A well-known paradox is that the size of the dissipative zone predicted by LVT is unrealistically small. Here we use a physically based nonlinear viscoelastic model to illustrate why the linear theory breaks down. Using this nonlinear model and analogs of crack problems, we give a plausible resolution to this paradox. In our model, viscoelasticity arises from the breaking and healing of physical cross-links in the polymer network. When the deformation is small, the kinetics of bond breaking and healing are independent of the strain/stress history and the model reduces to the standard linear theory. For large deformations, localized bond breaking damages the material near the crack tip, reducing stress concentration and dissipating energy at the same time. The damage zone size is a new length scale which depends on the strain required to accelerate bond breaking kinetics. These effects are illustrated by considering two cases with stress concentrations: the evolution of spherical damage in a viscoelastic body subjected to internal pressure, and a zero degree peel test.
more »
« less
Nanoscale crack propagation in clay with water adsorption through reactive MD modeling
Abstract The atomic‐scale cracking mechanism in clay is vital in discovering the cracking mechanism of clay at the continuum scale in that clay is a nanomaterial. In this article, we investigate mechanisms of modes I and II crack propagations in pyrophyllite and Ca‐montmorillonite with water adsorption through reactive molecular dynamics (MD) with a bond‐order force field. Clay water adsorption is considered by adding water molecules to the clay surface. During the equilibration stage, water adsorption could cause bending deformation of the predefined edge crack region. The relatively small orientating angle of water molecules indicates the formation of hydrogen bonds in the crack propagation process. The peak number density of adsorbed water decreases with the increasing strains. The atomistic structure evolution of the crack tip under loading is analyzed to interpret the nanoscale crack propagation mechanism. The numerical results show that the crack tip first gets blunted with a significant increase in the radius of the curvature of the crack tip and a slight change in crack length. The crack tip blunting is studied by tracking the crack tip opening distance and O–Si–O angle in the tetrahedral Si–O cell in modes I and II cracks. We compare bond‐breaking behaviors between Al–O and Si–O. It is found that Si–O bond breaking is primarily responsible for crack propagation. The critical stress intensity factor and critical energy release rate are determined from MD simulation results.
more »
« less
- PAR ID:
- 10411042
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical and Analytical Methods in Geomechanics
- Volume:
- 47
- Issue:
- 7
- ISSN:
- 0363-9061
- Page Range / eLocation ID:
- p. 1103-1133
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Objective. The insertion of penetrating neural probes into the brain is crucial for advancing neuroscience, yet it involves various inherent risks. Prototype probes are typically inserted into hydrogel-based brain phantoms and the mechanical responses are analyzed in order to inform the insertion mechanics duringin vivoimplantation. However, the underlying mechanism of the insertion dynamics of neural probes in hydrogel brain phantoms, particularly the phenomenon of cracking, remains insufficiently understood. This knowledge gap leads to misinterpretations and discrepancies when comparing results obtained from phantom studies to those observed under thein vivoconditions. This study aims to elucidate the impact of probe sharpness and dimensions on the cracking mechanisms and insertion dynamics characterized during the insertion of probes in hydrogel phantoms.Approach. The insertion of dummy probes with different shank shapes defined by the tip angle, width, and thickness is systematically studied. The insertion-induced cracks in the transparent hydrogel were accentuated by an immiscible dye, tracked byin situimaging, and the corresponding insertion force was recorded. Three-dimensional finite element analysis models were developed to obtain the contact stress between the probe tip and the phantom.Main results. The findings reveal a dual pattern: for sharp, slender probes, the insertion forces remain consistently low during the insertion process, owing to continuously propagating straight cracks that align with the insertion direction. In contrast, blunt, thick probes induce large forces that increase rapidly with escalating insertion depth, mainly due to the formation of branched crack with a conical cracking surface, and the subsequent internal compression. This interpretation challenges the traditional understanding that neglects the difference in the cracking modes and regards increased frictional force as the sole factor contributing to higher insertion forces. The critical probe sharpness factors separating straight and branched cracking is identified experimentally, and a preliminary explanation of the transition between the two cracking modes is derived from three-dimensional finite element analysis.Significance. This study presents, for the first time, the mechanism underlying two distinct cracking modes during the insertion of neural probes into hydrogel brain phantoms. The correlations between the cracking modes and the insertion force dynamics, as well as the effects of the probe sharpness were established, offering insights into the design of neural probes via phantom studies and informing future investigations into cracking phenomena in brain tissue during probe implantations.more » « less
-
This paper presents the development, implementation and experimental evaluation of a new crack detection mechanism for centrifuge modelling. The proposed mechanism is a brittle conductor bonded to cement providing a binary indication of if, and when, a sensor is cracked. The results of a pair of large centrifuge tests were used to evaluate the effectiveness of the proposed crack detection mechanism. Each test model included a soil profile consisting of a 23 m thick layer of lightly over-consolidated clay, underlain and overlain by thin layers of dense sand. The centrifuge models had two separate zones, a zone without reinforcement and a zone with an ‘embedded’ soil–cement grid, which had a unit cell area replacement ratio A r = 24%. Models were subjected to 13 different shaking events with peak base accelerations ranging from 0·01 to 0·55g. The performance of the proposed crack detection mechanism was examined using (i) post-test crack mapping in the soil–cement grids, (ii) results of the crack detection system and (iii) time series of accelerations, displacements and footing rotation. The results from the centrifuge test showed that the proposed crack detection method accurately captured if, and when, cracking occurred in the soil–cement grid at the locations of the sensors.more » « less
-
null (Ed.)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
-
Summary The mechanical and hydraulic properties of unsaturated clay under nonisothermal conditions have practical implications in geotechnical engineering applications such as geothermal energy harvest, landfill cover design, and nuclear waste disposal facilities. The water menisci among clay particles impact the mechanical and hydraulic properties of unsaturated clay. Molecular dynamics (MD) modeling has been proven to be an effective method in investigating clay structures and their hydromechanical behavior at the atomic scale. In this study, we examine the impact of temperature increase on the capillary force and capillary pressure of the partially saturated clay‐water system through high‐performance computing. The water meniscus formed between two parallel clay particles is studied via a full‐scale MD modeling at different elevated temperatures. The numerical results have shown that the temperature increase impacts the capillary force, capillary pressure, and contact angle at the atomic scale. The capillary force on the clay particle obtained from MD simulations is also compared with the results from the macroscopic theory. The full‐scale MD simulation of the partially saturated clay‐water system can not only provide a fundamental understanding of the impact of temperature on the interface physics of such system at the atomic scale, but also has practical implication in formulating physics‐based multiscale models for unsaturated soils by providing interface physical properties of such materials directly through high‐performance computing.more » « less
