skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]

Abstract

Here, we present a numerical method for modeling the mechanical effects of nonlinearly-compliant joints in elasto-plastic media. The method uses a series of strain-rate and stress update algorithms to determine joint closure, slip, and solid stress within computational cells containing multiple “embedded” joints. This work facilitates efficient modeling of nonlinear wave propagation in large spatial domains containing a large number of joints that affect bulk mechanical properties. We implement the method within the massively parallel Lagrangian code GEODYN-L and provide verification and examples. We highlight the ability of our algorithms to capture joint interactions and multiple weakness planes within individual computational cells, as well as its computational efficiency. We also discuss the motivation for developing the proposed technique: to simulate large-scale wave propagation during the Source Physics Experiments (SPE), a series of underground explosions conducted at the Nevada National Security Site (NNSS).

Authors:
 [1];  [1];  [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1406448
Report Number(s):
LLNL-JRNL-692181
Journal ID: ISSN 0178-7675; TRN: US1702910
Grant/Contract Number:
AC52-07NA27344
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Computational Mechanics
Additional Journal Information:
Journal Volume: 60; Journal Issue: 2; Journal ID: ISSN 0178-7675
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; 58 GEOSCIENCES; Discontinuities; Elastic-plastic solids; Rock mechanics; Wave

Citation Formats

Hurley, R. C., Vorobiev, O. Y., and Ezzedine, S. M. An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]. United States: N. p., 2017. Web. doi:10.1007/s00466-017-1403-6.
Hurley, R. C., Vorobiev, O. Y., & Ezzedine, S. M. An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]. United States. doi:10.1007/s00466-017-1403-6.
Hurley, R. C., Vorobiev, O. Y., and Ezzedine, S. M. Thu . "An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]". United States. doi:10.1007/s00466-017-1403-6. https://www.osti.gov/servlets/purl/1406448.
@article{osti_1406448,
title = {An algorithm for continuum modeling of rocks with multiple embedded nonlinearly-compliant joints [Continuum modeling of elasto-plastic media with multiple embedded nonlinearly-compliant joints]},
author = {Hurley, R. C. and Vorobiev, O. Y. and Ezzedine, S. M.},
abstractNote = {Here, we present a numerical method for modeling the mechanical effects of nonlinearly-compliant joints in elasto-plastic media. The method uses a series of strain-rate and stress update algorithms to determine joint closure, slip, and solid stress within computational cells containing multiple “embedded” joints. This work facilitates efficient modeling of nonlinear wave propagation in large spatial domains containing a large number of joints that affect bulk mechanical properties. We implement the method within the massively parallel Lagrangian code GEODYN-L and provide verification and examples. We highlight the ability of our algorithms to capture joint interactions and multiple weakness planes within individual computational cells, as well as its computational efficiency. We also discuss the motivation for developing the proposed technique: to simulate large-scale wave propagation during the Source Physics Experiments (SPE), a series of underground explosions conducted at the Nevada National Security Site (NNSS).},
doi = {10.1007/s00466-017-1403-6},
journal = {Computational Mechanics},
number = 2,
volume = 60,
place = {United States},
year = {Thu Apr 06 00:00:00 EDT 2017},
month = {Thu Apr 06 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 2works
Citation information provided by
Web of Science

Save / Share:
  • A mathematical model, based on the finite-element technique and incorporating thermo-elasto-plastic behavior during the water spray quenching of steel, has been developed. In the model, the kinetics of diffusion-dependent phase transformation and martensitic transformation have been coupled with the transient heat flow to predict the microstructural evolution of the steel. Furthermore, an elasto-plastic constitutive relation has been applied to calculate internal stresses resulting from phase changes as well as temperature variation. The computer code has been verified for internal consistency with previously published results for pure iron bars. The model has been applied to the water spray quenching of twomore » grades of steel bars, 1035 carbon and nickel-chromium alloyed steel; the calculated temperature, hardness, distortion, and residual stresses in the bars agreed well with experimental measurements. The results show that the phase changes occurring during this process affect the internal stresses significantly and must be included in the thermomechanical model.« less
  • The molar volume difference between the matrix and the precipitate phases in the case of solid state phase transformations results in the creation of stain energy in the system due to the misfit strains. A finite element model based on the initial strain approach is proposed to evaluate elasto-plastic accommodation energies during solid state transformation. The three-dimensional axisymmetric model has been used to evaluate energies as a function of transformation for {alpha}-{beta} hydrogen transformations in the Nb-H system. The transformation has been analyzed for the cases of transformation progressing both from the center to surface and from the surface tomore » center of the system. The effect of plastic deformation has been introduced to make the model realistic, specifically to the Nb-NbH phase transformation which involves a 4% linear misfit strain. It has been observed that plastic deformation reduces the strain energies compared to the linear elastic analysis.« less
  • The elasto-plastic deformation due to thermal and mechanical loading of layered metal-ceramic composites with compositionally graded interfaces is analyzed using detailed finite element models. The model material considered is a Ni-Al{sub 2}O{sub 3} layered system with a compositionally graded interface. The analyzes consider planar geometries with perfectly periodic arrangements of the constituent phases, by recourse of new classes of square-packing and hexagonal-packing unit cell formulations for the graded material. Also considered are graded phase arrangements within which large numbers of microstructural units of the two phases are randomly placed within the unit cell. It is found that square-packing arrangements providemore » the best possible bounds for the thermal strains and coefficient of thermal expansion (CTE) of the graded multi-layer, among the different unit cell models examined; however, no unique bounds could be identified for mechanical loading. The numerical predictions of thermal and mechanical response are compared with those provided by the mean-field approach involving an incremental Mori-Tanaka analysis and by the simple rule-of-mixture approximations. The former method provides a stiffer mechanical response than the finite-element unit cell models. The finite element predictions of bending CTE due to thermal excursions match the overall trends observed experimentally for the Ni-Al{sub 2}O{sub 3} graded system, and further provide a quantitative prediction of the temperature for the onset of plastic flow in the graded material.« less
  • Elastic and elasto-plastic modelling of indentation in CFRP cross-ply laminates has been performed. Detailed knowledge of the field solutions in the volume below the indentor forms the basis for the reported micromechanical interpretation of the observed damage in test specimens. The analysis shows that matrix cracks originate at sites of maximum tensile stress perpendicular to fibers. The predicted stress fields due to indentation show that stress concentrations occur in the interface between alternating plies. It is found that microcracking in this zone is a precursor to the observed failure. This analysis is supported by in-situ scanning electron microscopy during loadingmore » by a cylindrical indentor onto the laminate supported on a rigid substrate. The microscopy reveals microdamage in the region of interfacial tensile stress concentrations. The onset of indentation failure in these layered composites suggests that plastic interleaves would delay failure. It is shown numerically that plastic deformation of the interleaves redistributes stresses and thereby weakens the tensile stress concentrations which arise during indentation. Experimentally it is shown that aluminium interleaves affect the formation of indentation failure. In a cross-ply laminate, where alternating ply groups are separated by aluminium sheets, matrix cracking and delamination failures are suppressed by the occurrence of plastic deformation. Since the aluminium is likely to be weakly bonded to the plies, it is seen that weak interlaminar fracture toughness does not necessarily cause delaminations, nor lead to a lower indentation strength. High indentation strength and delamination resistance are complex qualities which, among others, seem to be achieved in laminate geometries which have a minimum of stress concentration at interfaces between ply groups of different orientation.« less