Numerical analyses of deformability tests in jointed rock: Joint perturbation and no tension finite element solutions
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
In this paper, the Barton–Bandis model for rock joints is extended to cyclic loading conditions, without any new material parameter. Also developed herein is an implicit solution algorithm for the extended Barton–Bandis model, which can also be used for the original Barton–Bandis model for which an implicit algorithm has been unavailable. To this end, we first cast the Barton–Bandis model into an incremental elasto-plastic framework, deriving an expression for the elastic shear stiffness being consistent with the original model formulation. We then extend the model formulation to cyclic loading conditions, incorporating the dependence of shear stress and dilation on the joint position and the shearing direction. The extension is achieved by introducing a few state-dependent variables which can be calculated with the existing material parameters. For robust and accurate utilization of the model, we also develop an implicit algorithm based on return mapping, which is unconditionally stable and guarantees the satisfaction of the strength criterion. We verify that the proposed model formulation and algorithm produce virtually the same results as the original Barton–Bandis model under monotonic shearing conditions. Here, we then validate the extended Barton–Bandis model against experimental data on natural rock joints under cycling loading conditions. The present work thus enables the Barton–Bandis model, which has been exceptionally popular in research and practice, to be applicable to a wider range of problems in rock mechanics and rock engineering.
Three major investigations using LANDSAT and Skylab imagery concerned with analyses of color anomalies and linear features of central Colorado are discussed. The studies conducted are concerned with the geological significance of spectral and spatial features on the images. Color anomalies in Skylab photographs were analyzed and evaluated for locating indicators of mineralization. The relationships were determined of all linear features in a LANDSAT image to the rock joint systems and the detectable larger geologic structures; techniques for extracting that geologic information are indicated. Some anomalous megalinear features in LANDSAT and Skylab images are analyzed which transect major structures and, their associated geologic features are described.
The Cheney-Palouse tract of the channeled scabland is the largest continuous tract of scabland in eastern Washington. The tract is composed of a varied assortment of bedrock erosional forms, loess islands, and gravel bars. Prominent bedrock longitudinal grooves and inner channels formed by macroturbulent plucking erosion of the jointed rock. Loess island forms vary as a function of their position within the flow. The three major types (submerged, partially submerged, and subaerially exposed) created sedimentologic conditions and resulting bar forms distinct from one another. Other bar forms, notably expansion bars, account for most of the sedimentation in the tract.
A new technique using external pressure was developed to simulate the lithostatic pressure due to self-gravity of small bodies. A 13-in. diameter cylindrical test chamber with L/D of 1 was fabricated to accommodate firing explosive charges with gas overpressures of up to 6000 psi. The chamber was hydrotested to 9000 psi. The method allows much larger scale factors that can be obtained with existing centrifuges and has the correct spherical geometry of self gravity. A simulant for jointed rock to be used in this fixture was developed using weakly cemented basalt. Various strength/pressure scaling theories can now be examined and tested.
Our joint research effort was aimed at developing techniques for X-ray diffractometry that was being investigated by NASA as possible flight instrumentation for the exploration of Mars. SJSU would provide the use of in-house X-ray facilities for calibration of the instrumentation , and would provide technical expertise regarding interpretation of data acquired during both laboratory testing, and during field testing of instruments on the Marsokhod rover at Ames. Accomplishments are: (1) quantification of X-ray signals from rock surfaces using San Jose State University (SJSU) diffractometer; (2) development of criteria for fingerprinting rock samples using pattern recognition of diffraction spectra, and augmentation of diffraction data with X-ray fluorescence information; (3) calibration of NASA instrumentation using SJSU-generator data; and (4) assistance in the development, lab testing, and field deployment of the NASA instrument on the Russian Marsokhod roving vehicle designed for martian exploration.
The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.
Granitic inselbergs are enduring landforms that punctuate planation surfaces in tectonically stable settings, yet the mechanisms and timescales of their persistence remain debated. Here, in this study, we present a multinuclide cosmogenic dataset ( 10 Be, 26 Al and 21 Ne) from bedrock, fluvial sediments and planation surfaces in the semi-arid Central Ceará Domain, north-eastern Brazil, to quantify denudation rates and reconstruct the long-term dynamics of residual relief formation. This triple-nuclide approach enables independent validation of denudation contrasts and long-term exposure histories, providing quantitative constraints on inselberg evolution that a single-nuclide approach could not resolve. Bedrock denudation rates (1–4 m Myr −1 ) are systematically lower than basin-wide rates (8–14 m Myr −1 ), revealing a regional pattern of vertical differential denudation. When combined with summit elevations and surface exposure ages of up to 5 Myr, these contrasts suggest that cumulative relief growth has been operating over timescales of tens of millions of years (1–57 Myr). Morphological and structural observations further indicate localized mass wasting and heterogeneous erosion styles, modulated by rock fabric and joint density. Cosmogenic nuclide inventories in detrital sediments reveal partial reworking from older sediment stocks, including surface pebbles with exposure ages up to 5 Myr, which record multistage of burial and re-exposure histories and reflect episodic remobilization in low-connectivity, semi-arid catchments. Rather than static relicts of ancient surfaces, inselbergs in this cratonic setting emerge as dynamic landforms actively sustained by slow but persistent differential denudation. These results provide new constraints on long-term cratonic landscape evolution and underscore the interplay between lithological resistance, structural inheritance and sediment transfer processes in sustaining topographic relief. They also call for a critical reassessment of classical inselberg models, suggesting that hybrid mechanisms, rather than single-process paradigms, more accurately capture the complexity of residual landform development in cratonic interiors.
Cryoplanation terraces are step- or table like residual landforms consisting of a nearly horizontal bedrock surface covered by a thin veneer of rock debris and bounded by ascending or descending scarps or both. Among examples studied, rubble-covered scarps range in height from 3 to 76 m and slope from 9 deg to 32 deg; nearly vertical scarps exist where bedrock is exposed or thinly buried. Simple transverse nivation hollows, which are occupied by large seasonal snow banks, commonly indent the lower surfaces of sharply angular ascending scarps. Terrace treads slope from 1 deg to 10 deg and commonly cut across bedrock structures such as bedding, rock contacts, foliation, joints, faults, and shear zones. Debris on terrace treads is generally 0.8-2.5 m thick. Permafrost table is generally present from 0.5 to 2 m below the tread surface. Permafrost is shallowest in the floors of nivation hollows and deepest in the well-drained margins of terrace treads. Side slopes of cryoplanation terraces are shallowly buried bedrock surfaces that are littered with a variety of mass-movement deposits.
In reexamining the accumulated magnetic data on lunar rocks, several common patterns of magnetic behavior are recognized. Their joint occurrence strongly suggests a new model of lunar rock magnetism, which is based on partial preferred textural alignment of the spontaneous moments of magnetic grains without requiring the existence of ancient lunar magnetic fields. This magnetic fabric, mimetic to locally oriented petrofabric, gives rise to an apparent 'textural remanent magnetization'. In order to account for the observed intensity of 'stable remanence' in lunar rocks, only a minute fraction (0.001 to 0.00001) of the single-domain iron grains present need be preferentially aligned. Several mechanisms operating on the lunar surface, including shock and diurnal thermal cycling, appear adequate for producing the required type and degree of magnetic alignment in all lunar rock classes. The model is supported by a wide variety of direct and indirect evidence, and its predictions (e.g., regarding anisotropic susceptibility and remanence acquisition) can be experimentally tested.
Abstract not provided.
Here we present our analysis of the near- and mid-infrared spectral properties of pyroclastic deposits within the floor fractured Oppenheimer Crater that are hypothesized to be Vulcanian in origin. These are the first results of our global study of lunar pyroclastic deposits aimed at constraining the range of eruption processes on the Moon. In the near-infrared, we have employed a new method of spectral analysis developed in Horgan et al. (2013) of the 1 m iron absorption band in Chandrayaan-1 Moon Mineralogy Mapper (M3) spectra. By analyzing both the position and shape of the 1 m band we can detect and map the distribution of minerals, glasses, and mixtures of these phases in pyroclastic deposits. We are also using mid-infrared spectra from the Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment to develop ~200 m/pixel Christiansen Feature (CF) maps, which correlate with silica abundance. One of the benefits of using CF maps for analysis of pyroclastic deposits is that they can be used to detect silicic country rock that may have been emplaced by Vulcanian-style eruptions, and are sensitive to iron abundance in glasses, neither of which is possible in the near-infrared. M3 analysis reveals that the primary spectral endmembers are low-calcium pyroxene and iron-bearing glass, with only minor high-calcium pyroxene, and no detectable olivine. The large deposit in the south shows higher and more extensive glass concentrations than the surrounding deposits. We interpret the M3 spectra of the pyroclastic deposits as indicating a mixture of low-calcium pyroxene country rock and juvenile glass, and no significant olivine. Analysis of Diviner CF maps of the Oppenheimer crater floor indicates an average CF value of 8.16, consistent with a mixture of primarily plagioclase and some pyroxene. The average CF values of the pyroclastic deposits range from 8.31 in the SW to 8.24 in the SE. Since CF values within the deposits are as high as 8.49, the lower average CF values of the deposits suggest that each deposit is a mixture of crater floor material and highly mafic juvenile material consistent with either olivine or Fe-bearing pyroclastic glass. Synthesizing our M3 and Diviner results indicates that the crater floor consists of plagioclase with some pyroxene, and the pyroclastic deposits are a mix of this substrate and a glass-rich juvenile material. While we cannot determine the iron content of the glass from M3 spectra alone, the high Diviner CF values suggest that the glass is relatively iron-rich. Indeed, FeO abundances inferred from CF values using the method of Allen et al. (2012) imply that the large southern deposit exhibits a significant enhancement in iron content. This supports our hypothesis that the glass in this deposit is relatively iron-rich.
Results are presented of three early applications of imagery from the NASA Earth Resources Technology Satellite to geologic studies in California. In the Coast Ranges near Monterey Bay, numerous linear drainage features possibly indicating unmapped fracture zones were mapped within one week after launch of the satellite. A similar study of the Sierra Nevada near Lake Tahoe revealed many drainage features probably formed along unmapped joint or faults in granitic rocks. The third study, in the Peninsular Ranges, confirmed existence of several major faults not shown on published maps. One of these, in the Sawtooth Range, crosses in Elsinore fault without lateral offset; associated Mid-Cretaceous structures have also been traced continuously across the fault without offset. It therefore appears that displacement along the Elsinore fault has been primarily of a dip-slip nature, at least in this area, despite evidence for lateral displacement elsewhere.
Igneous rocks are the primary building blocks of planetary crusts. Most igneous rocks originate via decompression melting and/or wet melting of protolith lithologies within planetary interiors and their classification and compositional, petrographic, and textural characteristics, arewell-studied. As our exploration of the Solar System continues, so too does the inventory of intrusive and extrusive igneous rocks, settings, and processes. The results of planetary exploration have also clearly demonstrated that impact cratering is a ubiquitous geological process that has affected, andwill continue to affect, all planetary objectswith a solid surface,whether that be rock or ice. It is nowrecognized that the production of igneous rocks is a fundamental outcome of hypervelocity impact. The goal of this review is to provide an up-to-date synthesis of our knowledge and understanding of igneous rocks formed by hypervelocity impact. Following a brief overview of the basics of the impact process, we describe howand whymelts are generated during impact events and howimpactmelting differs fromendogenic igneous processes.While the processmay differ,we showthat the products of hypervelocity impact can share close similarities with volcanic and shallow intrusive igneous rocks of endogenic origin. Such impact melt rocks, as they are termed, can display lobate margins and cooling cracks, columnar joints and at the hand specimen andmicroscopic scale, such rocks can displaymineral textures that are typical of volcanic rocks, such as quench crystallites, ophitic, porphyritic, aswell as features such as vesicles, flowtextures, and so on.Historically, these similarities led to themisidentification of some igneous rocks now known to be impact melt rocks as being of endogenic origin. This raises the question as to how to distinguish between an impact versus an endogenic origin for igneous-like rocks on other planetary bodies where fieldwork and sample analysis may not be possible and all that may be available is remote sensing data.While the interpretation of some impact melt rocks may be relatively straightforward (e.g., for clast-rich varieties and those with clear projectile contamination) we conclude that distinguishing between impact and endogenic igneous rocks is a non-trivial task that ultimately may require sample investigation and analysis to be conducted. Caution is, therefore, urged in the interpretation of igneous rocks on planetary surfaces.
Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.
The importance of the effects of structural joints on the damping of built-up structures is pointed out, and the energy dissipation mechanisms associated with squeezing, rocking, and shearing motions are discussed for simple joints that are dry, lubricated, or provided with viscoelastic inserts. The damping mechanisms and behaviors of built-up beams and of skin-stringer structures are discussed as far as they are currently understood, and available damping estimation methods are summarized. Difficulties in defining and measuring the damping of skin-stringer structures are indicated, and it is pointed out that particularly the high-frequency damping of built-up beams and the low-frequency damping of skin-stringer configurations require further investigation.
This report describes a three-month research program undertook jointly by the Robotics Institute at Carnegie Mellon University and Ames Research Center as part of the Ames' Joint Research Initiative (JRI.) The work was conducted at the Ames Research Center by Mr. Liam Pedersen, a graduate student in the CMU Ph.D. program in Robotics under the supervision Dr. Ted Roush at the Space Science Division of the Ames Research Center from May 15 1999 to August 15, 1999. Dr. Martial Hebert is Mr. Pedersen's research adviser at CMU and is Principal Investigator of this Grant. The goal of this project is to investigate and implement methods suitable for a robotic rover to autonomously identify rocks and minerals in its vicinity, and to statistically characterize the local geological environment. Although primary sensors for these tasks are a reflection spectrometer and color camera, the goal is to create a framework under which data from multiple sensors, and multiple readings on the same object, can be combined in a principled manner. Furthermore, it is envisioned that knowledge of the local area, either a priori or gathered by the robot, will be used to improve classification accuracy. The key results obtained during this project are: The continuation of the development of a rock classifier; development of theoretical statistical methods; development of methods for evaluating and selecting sensors; and experimentation with data mining techniques on the Ames spectral library. The results of this work are being applied at CMU, in particular in the context of the Winter 99 Antarctica expedition in which the classification techniques will be used on the Nomad robot. Conversely, the software developed based on those techniques will continue to be made available to NASA Ames and the data collected from the Nomad experiments will also be made available.
Recent analysis has shown that by 2050, geothermal energy can supply up to 90 GW of power to the U.S. electrical grid. Reducing the lifecycle cost of geothermal wells is critical to achieving this level of adoption. While time-related costs of drilling remain an important issue, there is also a need to focus on the costs associated with casing and cementing of geothermal wells. This report provides a roadmap of R&D activities needed to reduce the life cycle cost of casing and cementing while improving life-of-well performance. The purpose of the roadmap is to develop priorities for R&D in the areas of (1) high-performance and cost-effective materials for target geothermal well conditions, (2) well construction methods and techniques to reduce construction costs, and (3) methods and techniques to decrease to decrease long-term operating costs. The roadmap sets targets for the next 10-year period and purposely excludes R&D efforts into improved rock reduction. The roadmap was developed by a joint working group of experts from Oak Ridge National Laboratory, Brookhaven National Laboratory, and Sandia National Laboratories as well as 35 experts from industry and academia. Input from a series of information gathering sessions was used to inform and guide the development of the roadmap.