Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mechanical loading”

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.

At least 649 records · Page 36

The mechanism for acetate formation in electrochemical CO (2) reduction on Cu: selectivity with potential, pH, and nanostructuring

Nanostructured Cu catalysts have increased the selectivities and geometric activities for high value C–C coupled (C 2 ) products (ethylene, ethanol, and acetate) in the electrochemical CO (2) reduction reaction (CO (2) RR). The selectivity among the high-value C 2 products is also altered, where for instance the yield of acetate increases with alkalinity and is dependent on the catalyst morphology. The reaction mechanisms behind the selectivity towards acetate vs. other C 2 products remain controversial. In this work, we elucidate the reaction mechanism for acetate formation by using ab initio simulations, a coupled kinetic-transport model, and loading dependent experiments. We find that trends in acetate selectivity can be rationalized from variations in electrolyte pH and the local mass transport properties of the catalyst and not from changes in Cu's intrinsic activity. The selectivity mechanism originates from the transport of ketene, a stable (closed shell) intermediate, away from the catalyst surface into solution where it reacts to form acetate. While this type of mechanism has not yet been discussed in the CO (2) RR, variants of it may explain similar selectivity fluctuations observed for other stable intermediates like CO and acetaldehyde. Our proposed mechanism suggests that acetate selectivity increases with increasing pH, decreasing catalyst roughness and significantly varies with the applied potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An approach towards plastic scintillators from thermally activated delayed fluorescent dyes and cross-linkable bismuth compounds

Cross-linked polyvinyltoluene (PVT) based plastic scintillators loaded with dyes exhibiting thermally activated delayed fluorescence (TADF) and various amounts of (caproate)di-(methacrylate) bismuth (CMB) were developed. The co-polymerization of vinyl toluene monomer loaded with different weight ratios of CMB and 1 wt% of TADF dye yielded highly transparent cross-linked plastic scintillators with intense green photoluminescence. The thermal and mechanical properties of the scintillators were investigated using thermogravimetric analysis (TGA), differential scanning calorimetry (DSC) and micro indentation techniques. Here, their performance for radiation detection was assessed by performing gamma spectroscopy using a Cs-137 source and benchmarked to commercial plastic scintillators. The loading of CMB was guided by Monte Carlo N-Particle (MCNP) simulations. The results of this study show that the cross-linking approach and the use of TADF dyes yield scintillators with good mechanical properties with a uniaxial yield strength up to ~66 MPa for compositions loaded with 40 wt% of CMB and high optical quality. However, the light yield in gamma spectroscopy experiments is found to be reduced with the addition of CMB, suggesting a trade-off between mechanical yield strength and light yield in gamma spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metaproteomics reveals enzymatic strategies deployed by anaerobic microbiomes to maintain lignocellulose deconstruction at high solids

Economically viable production of cellulosic biofuels requires operation at high solids loadings—on the order of 15 wt%. To this end we characterize Nature’s ability to deconstruct and utilize mid-season switchgrass at increasing solid loadings using an anaerobic methanogenic microbiome. This community exhibits undiminished fractional carbohydrate solubilization at loadings ranging from 30 g/L to 150 g/L. Metaproteomic interrogation reveals marked increases in the abundance of specific carbohydrate-active enzyme classes. Significant enrichment of auxiliary activity family 6 enzymes at higher solids suggests a role for Fenton chemistry. Stress-response proteins accompanying these reactions are similarly upregulated at higher solids, as are β-glucosidases, xylosidases, carbohydrate-debranching, and pectin-acting enzymes—all of which indicate that removal of deconstruction inhibitors is important for observed undiminished solubilization. Our work provides insights into the mechanisms by which natural microbiomes effectively deconstruct and utilize lignocellulose at high solids loadings, informing the future development of defined cultures for efficient bioconversion.

09 BIOMASS FUELS↗

TEAMER Technical Support for Aquantis (Materials) (Final Report)

The Aquantis Tidal Power Tug is a unique synthesis of best-available technologies and materials configured as a novel spar vessel to create an optimal platform for tidal stream energy conversion. The Power Tug utilizes an upstream-facing horizontal, 2-bladed rotor. To drive down capital costs and extend life, Aquantis proposes to employ new materials for the blades that are both less expensive than current state-of-the-art materials and are potentially better suited for survival in a submerged, seawater environment. Coupons from a novel geopolymer material will be tested in a load frame to characterize material properties. The Requestor is requesting technical assistance for evaluating mechanical properties of a novel geopolymer material for application in tidal stream energy conversion systems. This work will establish characteristic properties of a geopolymer material under dry and saturated conditions through material coupon testing. Material coupons will be tested to characterize properties including: a) Tensile strength, b) Compression strength, c) Flexural strength, d) Modulus of elasticity, e) Fatigue strength at varying load levels. The types of material tests performed will be guided by characteristics required to validate material properties for use in full-scale structures. Characterization of material properties that consider the demands of a full-scale structure is a necessary step in research and development that can enable the use of geopolymer materials in tidal stream energy conversion systems.

16 TIDAL AND WAVE POWER↗

A Mechanical Test Frame for Property Evaluations at Cryogenic Temperature

On-board storage of hydrogen fuel has been designated as a limiting factor in the advancement of fuel cell technologies in the automotive industry. The use of cryo-compressed Type III pressure vessels presents one option to overcome this barrier. These multi-material vessels will be expected to perform at high pressure and extreme low temperatures, which creates a complex engineering design challenge. Many materials exhibit highly temperature-dependent properties, which must be considered for the efficient design of cryo-compressed pressure vessels. Complicating this issue, mechanical property data at cryogenic temperatures is sparse. This technical paper provides a description of a test apparatus commissioned specifically to help address this shortcoming. A mechanical test frame was retrofitted with a continuous flow cryostat capable to evaluate various mechanical properties throughout a broad temperature range of 25 °C to -269 °C with a load limit of 10 kN. An investigation of the thermomechanical properties of an epoxy resin was carried out as a demonstration.

Mechanical characterization, cryogenic pressure ve↗

High-speed synchrotron X-ray phase-contrast imaging for evaluating microscale damage mechanisms and tracking cracking behaviors inside cross-ply GFRCs

In this study, we integrated the high-speed synchrotron X-ray phase-contrast imaging (PCI) with a modified Kolsky compression bar loading platform to visualize the dynamic failure processes of cross-ply glass fiber reinforced composites (GFRCs). Four S-2 glass/SC-15 composite specimens were prepared, having similar thicknesses but different stacking sequences, namely as [0 12 /90 12 ], [90 12 /0 12 ], [0 8 /90 8 /0 8 ], and [90 8 /0 8 /90 8 ]. Three-dimensional synchrotron X-ray computed tomography and scanning electron microscopy (SEM) were employed to examine the microstructures and quantify the fiber volume fractions. Each specimen was notched and subjected to a dynamic three-point flexural loading. The onset of cracking close to the notch tip, crack propagation in 0° or 90° plies and their interface, crack opening, and ultimately failure of the specimen were captured by high-speed synchrotron X-ray PCI. Additional dynamic experiments were performed to determine the average time when the stress wave propagated through the specimen and correlate the X-ray images with the specimen's force-deflection response. Finally, the surface morphology of each specimen after the dynamic loading was imaged by SEM. Comparison between real-time X-ray images and post-fracture SEM images demonstrated the capability of the X-ray method to record damage evolution inside composites. Furthermore, the high spatial and temporal resolutions of the X-ray setup and edge enhancement by PCI enabled the identification of microscale damage features within 1 μs. Two damaging processes were identified, crack growth was quantified, and fracture toughness of the composites was evaluated. The method is deemed useful to reveal microscale damage mechanisms and track cracking behaviors inside cross-ply GFRCs under dynamic loading in real time.

42 ENGINEERING↗

Understanding Geometry and Microstructure Interactions In LPBF By Linking Lab Scale Studies To Real Part Case Examples

The mechanical properties derived from simple tensile tests are uniquely important in the design and qualification of load bearing metallic components. The tensile properties in standards for wrought metallic materials are generally accessible for design and simulation engineers, but this is not the case for additively manufactured metals. This work investigates the influence of tensile specimen geometry on the mechanical properties of laser powder bed fusion (LPBF) additively manufactured Ti-6Al-4V, particularly in relation to properties measured from specimens excised from an exemplar part. A comprehensive analysis was conducted across various specimen geometries, considering factors such as thermal history, cross-sectional area, shape, and surface roughness. Key findings reveal that sample size, particularly cross-sectional area, significantly affects reported tensile properties, with ASTM E8 flat and round specimens exhibiting differences in elastic modulus (13%), yield strength (17%), ultimate tensile strength (UTS) (14%), and elongation to failure (11%). Additionally, surface roughness was found to have a limited impact on ASTM round geometries, but it becomes critical in thin-walled and small high-throughput samples. The ASTM flat samples closely matched the properties of samples excised from the exemplar part, underscoring the necessity of testing witness samples that accurately represent end-use conditions. Furthermore, mechanical properties of samples built within the stitch zone of the exemplar part demonstrated reduced ductility, highlighting the importance of laser alignment and qualification in multi-laser additive manufacturing processes. Overall, this study emphasizes the need for careful consideration of specimen geometry and testing conditions to ensure accurate reporting of mechanical properties in LPBF Ti-6Al-4V.

36 MATERIALS SCIENCE↗

Energy absorption of architectured PrintCast interpenetrating composites in tension

Additively manufactured (AM) metal-metal composites consisting of PrintCasted 316 L austenitic stainless-steel lattice structures infiltrated with A356 casting alloy, have recently been developed for use in high energy absorption systems with potential applications ranging from static load bearing to dynamic blast containment structures. This system has a unique mechanical behavior as the volume fraction of lattice increases showing a transition from localized to de-localized failure and dramatic increase in energy absorption capability. In this work, PrintCast A356/316 L composite tensile specimens were produced with lattice volume fractions ranging from 20 % to 50 % to capture the range of this behavior. Finite element simulations support neutron diffraction measurements of stress state. Results illustrate that in tension, the reinforcement material is in tension while the matrix support material is in compression, information offering significant insight into the transition to de-localized failure. Moreover, the simulation results provide further insight into how interfacial bonding (or lack of bonding) affects the energy absorption capabilities of the PrintCast composites.

36 MATERIALS SCIENCE↗

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗

ATF Cladding Mechanical Properties Report: Capability Demonstration

This report documents mechanical testing capability demonstration activities performed in fiscal year (FY)2025 at Oak Ridge National Laboratory(ORNL)on chromium-coated (Cr-coated) and uncoated advanced zirconium alloy claddings irradiated in the High Flux Isotope Reactor (HFIR) to approximately 4 displacements per atom (dpa), corresponding to ~13GWd/t burnup. Specimens were prepared in axial tension (ATT) and ring tension (RTT) geometries, and passive silicon carbide(SiC)thermometry (TM) was employed to determine irradiation temperatures, which averaged 38–43 °C below the 330 °C design target. Mechanical testing at ambient temperature demonstrated the expected irradiation-induced hardening, with yield strength(YS)and ultimate tensile strength(UTS)values increasing substantially relative to unirradiated counterparts. However, this strengthening was accompanied by a reduction in ductility, as indicated by lower uniform and total elongations(UE and TE). Both coated and uncoated claddings exhibited similar mechanical response, though Cr-coated specimens showed surface cracking perpendicular to the loading direction, attributable to the hardness mismatch between the coating and substrate. Fracture in all cases remained ductile, and no coating spallation was observed following HFIR irradiation. Complementary efforts were directed toward the fabrication of test specimens from commercially irradiated cladding (rod 47I, ~31.1 GWd/t average burnup). Axial sectioning and computer numerical control(CNC)machining successfully produced ATT geometries suitable for benchmarking against HFIR-irradiated specimens. This capability enables direct comparison of cladding behavior between test reactor and commercial reactor environments, thereby supporting the validation of HFIR as a surrogate irradiation platform for accident tolerant fuel (ATF) development. Once HFIR irradiations are completed in FY26,the relevant comparison tests will be completed. Collectively, the FY 2025 PIE campaign has provided mechanical performance data for irradiated advanced claddings The demonstrated capabilities support the framework for mechanical testing and further evaluations in subsequent years. These efforts will represent an important contribution toward the licensing and deployment of Cr-coated zirconium alloy cladding as a near-term ATF solution.

36 MATERIALS SCIENCE↗

Assembly and Pre-Loading Specifications for the Series Production of the Nb 3 Sn MQXFA Quadrupole Magnets for the HL-LHC

We report the High Luminosity LHC (HL-LHC) Project is planning to install 16 cold-masses made with Nb 3 Sn quadrupole magnets in the LHC Interaction Regions to significantly increase its luminosity. Half of these cold masses are fabricated at BNL, FNAL, and LBNL under the US Accelerator Research Program (AUP). Each cold mass includes two identical Nb 3 Sn quadrupole magnets, called MQXFA with a magnetic length of 4.2 m. Currently, the AUP project has completed the fabrication and test of the first 5 pre-series magnets, and is working on the following 16 magnets for the series production. The brittleness and strain sensitivity of the Nb 3 Sn superconducting material requires a careful definition of the allowable maximum stress in the windings during magnet assembly and pre-load, and a tight control of their variation within the whole coil length. Therefore, a series of assembly and pre-load specifications have been defined with the goals of minimizing the risk of conductor degradation and providing the mechanical support required to reach the nominal current during powering. In this paper we present the specifications defined for the MQXFA magnets and applied during the different assembly phases and during the pre-load process of the first 5 pre-series magnets.

43 PARTICLE ACCELERATORS↗

A semi–automatic analytical methodology for characterizing the energy consumption of MRI systems using load duration curves

Background and purpose: Magnetic resonance imaging (MRI) scanners are a major contributor to greenhouse gas emissions from the healthcare sector, and efforts to improve energy efficiency and reduce energy consumption rely on quantification of the characteristics of energy consumption. The purpose of this work was to develop a semi-automatic analytical methodology for the characterization of the energy consumption of MRI systems using only the load duration curve (LDC). LDCs are a fundamental tool used across various fields to analyze and understand the behavior of loads over time. Methods: An electric current transformer sensor and data logger were installed on two 3T MRI scanners from two vendors, termed M1 (outpatient scanner) and M2 (inpatient/emergency scanner). Data was collected for 1 month (7/11/2023 to 8/11/2023). Active power was calculated, assuming a balanced three-phase system, using the average current measured across all three phases, a 480 V reference voltage for both machines, and vendor-provided power factors. An LDC was constructed for each system by sorting the active power values in descending order and computing the cumulative time (in units of percentage) for each data point. The first derivative of the LDC was then computed (LDC’), smoothed by convolution with a window function (sLDC’), and used to detect transitions between different system modes including (in descending power levels): scan, prepared-to-scan, idle, low-power, and off. The final, segmented LDC was used to measure time (% total time), total energy (kWh), and mean power (kW) for each system mode on both scanners. The method was validated by comparing mean power values, computed using the segmented 1-month LDC, for each nonproductive system mode (i.e., prepared-to-scan, idle, lower-power, and off) against power levels measured after a deliberate system shutdown was performed for each scanner (1 day worth of data). Results: The validation revealed differences in mean power values <1.4% for all nonproductive modes and both scanners. In the scan system mode, the mean power values ranged from 29.8 to 37.2 kW and the total energy consumed for 1 month ranged from 11 106 to 14 466 kWh depending on the scanner. Over the course of 1 month, the portion of time the scanners were in nonproductive modes ranged from 76% to 80% across scanners and the nonproductive energy consumption ranged from 8010 to 6722 kWh depending on the scanner. The M1 (outpatient) scanner consumed 99.9 and 183.9 kWh/day in idle mode for weekdays and weekends, respectively, because the scanner spent 23% more time proportionally in idle mode on the weekends. Conclusions: A semi-automatic method for quantifying energy consumption characteristics of MRI scanners was introduced and validated. This method is relatively simple to implement as it requires only power data from the scanners and avoids the technical challenges associated with extracting and processing scanner log files. Finally, the methodology enables quantitative evaluation of the power, time, and energy characteristics of MRI scanners in scan and nonproductive system modes, providing baseline data and the capability of identifying potential opportunities for enhancing the energy efficiency of MRI scanners.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A numerical study of reversible plasticity using continuum dislocation mechanics

In this contribution, an elasto-viscoplastic fast Fourier transform-based (EVPFFT) numerical implementation of the Mesoscale Field Dislocation Mechanics (MFDM) formulation, called MFDM-EVPFFT, is applied to study the reversible plastic behavior of periodic two-phase crystalline composites with an elasto-viscoplastic plastic matrix and a purely elastic second phase. Additionally, periodic laminate microstructures of this kind with different periods (i.e. sizes) are considered to examine the size dependence of the Bauschinger effect and hardening during cyclic loading. Comparisons with classic composite effects obtained with conventional crystal plasticity are discussed. Specifically, the MFDM-EVPFFT results shed light on the hardening mechanisms due to piling-up/unpiling-up of geometrically-necessary dislocations (GND) during reverse loading.

36 MATERIALS SCIENCE↗

Structural Design and Modeling of MARVEL Primary Coolant System Using the ASME Section III, Division 5, Code

This paper presents the structural design and supporting analysis for the Microreactor Applications Research Validation and Evaluation (MARVEL) primary coolant system (PCS) using the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code Section III, Division 5, rules. MARVEL is a liquid metal–cooled microreactor intended to provide experimental capabilities for the rapid testing and development of microreactor technologies. The PCS utilizes high-temperature sodium-potassium liquid metal as the primary coolant and operates at a design temperature of 570°C, necessitating the consideration of creep-related failure mechanisms. The base metal for the PCS is 316H stainless steel, and the weldments are made with a 16-8-2 filler. The design approach incorporates the current base code rules along with ASME code cases N-924, N-861, and N-862 to address primary load, ratcheting, and creep-fatigue evaluations, respectively. The reactor’s operation involves complex thermal and mechanical interactions due to natural convective flow and differential thermal expansion between components. In conclusion, this paper discusses the structural engineering challenges encountered, such as managing thermal stresses in the distribution plenum and guard vessel, and outlines the strategies implemented to meet the code requirements, including design modifications and operational constraints.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗

TEAMER: Biofouling Analysis for Wave Energy Piston Design - Load Cell Data

Biofouling and corrosion are a major concern for all ocean-deployed components, especially when mechanical motion is involved. Triton has developed the concept of a biofouling mitigation seal as part of the piston sealing assembly for the Triton Wave Energy Converter (TSI-WEC). This mitigation seal has the purpose of preventing the formation of a biofilm on the inside of the piston cylinder. It is hypothesized that the prevention of a biofilm will reduce the amount of macro-biofouling that can occur in the piston assembly. The mitigation seal can also reduce the wear on the main dynamic seal, helping to maintain smooth operation and water-tightness. The cylinder is made from a thermoset composite epoxy, which is resistant to corrosion. However, no studies have researched the material's performance with biofouling. Triton placed two prototype Power Take-Off (PTO) assemblies in a PNNL biofouling tank, one with a biofouling mitigation seal and one without, allowing for an evaluation of seal effectiveness at the prevention of biofouling. In actual WEC operation, wave action would react against the piston, which would drive the linear actuator and electric generator, providing electrical power. In the test setup, this was reversed; a linear actuator was powered to drive the piston in a consistent motion within the cylinder. There are two assemblies: one has a biofouling mitigation seal, the other (control) does not. The following data encompasses a 4 month test period, with load cells being used to monitor piston friction force. Results from this testing will be used to improve seal design and material selection, mitigating risk of premature failure during open water testing and evaluation. This project is part of the TEAMER RFTS 3 (request for technical support) program.

16 TIDAL AND WAVE POWER↗

30 cm Drop Modeling

This report documents modeling and analysis of the 30 cm cask drop scenario for spent nuclear fuel (SNF) casks. The 30 cm drop is an important hypothetical dynamic loading scenario for SNF because 10 CFR 71.71 defines testing requirements for normal conditions of transport (NCT) of SNF packages, and free drop (of 30 cm for the heaviest packages) tests are required. Other NCT conditions specified in 10 CFR 71.71 include vibration normally incident to transport, and DOE’s multimodal transportation test (MMTT) and related analysis concluded that mechanical shock and vibration incident to transportation is negligible for SNF structural integrity because the loads are so small. The SNF cask 30 cm drop has a lot more potential to cause damage to SNF because the kinetic energy and momentum of the system is stopped abruptly during the cask impact event. Test data and model predictions both agree that the SNF assemblies are subjected to significantly higher loads in a cask drop event than normal transportation shock and vibration, but currently there is no reason to expect gross structural damage of the fuel assembly or widespread SNF cladding failure as a result of a 30 cm cask drop. The ultimate goal of this work is to close the knowledge gap of stress profiles on SNF rods and fuel assembly components during NCT, and this report describes the progress toward reaching that goal as well as a plan for finishing this work over the next year or two. Key remaining tasks include refining the fuel assembly model, calculating a selection of specific fuel assembly cases to populate the damage model, and finally running the damage model to predict the full general response range of 30 cm drop loads on SNF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Extended Barton–Bandis model for rock joints under cyclic loading: Formulation and implicit algorithm

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.

58 GEOSCIENCES↗