Search NASA⌕ Search

SEARCH · Search NASA

Results for “plane strain”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dual-scale folding in cutting of commercially pure aluminum alloys

We examine a, hitherto, little-studied and curious machining chip morphology, with tell-tale signs of folding at two different length scales, that is common in cutting of certain ductile and highly strain-hardening metals like soft aluminum alloys, tantalum and niobium. This chip morphology does not appear in the usual catalogues of common chip types. The mechanics of formation of the “dual-scale folded chip” is studied in model material systems of commercially pure aluminum alloys (AA 1100 and AA 8040), that prominently exhibit this chip morphology. The flow, folding and associated plastic instabilities are investigated using micro/macro structure observations of the chip in a plane-strain cutting framework, with high-speed in situ imaging and image analysis of material flow; and force measurements. The smaller-scale folding is shown to develop in the primary deformation zone while the larger-scale folding occurs as the chip traverses the rake face of the tool. The resulting chip is composed of irregularly-spaced large folds, superimposed onto which are the quasi periodic small folds. The representative wavelengths of the two folds differ on average by an order of magnitude, 0.1 mm vs. 2 mm. The observations reveal a direct coupling between the material flow and chip morphology, and how specific attributes of the dual-scale folded chip arise from the flow mechanism. Plastic buckling is found to play a key role in the folding at both length scales. The small-scale folds are characteristic of a sinuous plastic flow mode, while the large-scale folding is characterized by buckling and stick-slip along the tool rake face, triggered by adhesive pinning of the chip to the tool. Important consequences of the dualscale folding are very large cutting forces, and force oscillations of large amplitude, despite the alloys being very soft, only ~ 25 HV. Here, the dual-scale folding is why many of these alloys are classified as “gummy” to machine. Since the dual-scale folded chip is associated with large cutting forces and poor surface quality, there is much to be gained by disrupting this flow type in practical machining applications. Methods for controlling the folding to improve machining performance with the gummy alloys are briefly discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Silicon Consortium Project: No-Go on Moir Interferometry for Measuring SEI Strain as a Probe for Calendar Life Testing

Silicon is a promising candidate as a next generation anode to replace or complement graphite electrodes due to its high energy density and low lithiation potential. When silicon is lithiated, it experiences over 300% expansion which stresses the silicon as well as its solid electrolyte interphase (SEI) leading to poor performance. The use of nano-sized silicon has helped to mitigate volume expansion and stress in the silicon, yet the silicon SEI is still both mechanically and chemically unstable. Identifying the mechanical failure mechanism of the SEI will help enhance calendar and cycle life performance through improved SEI design. In situ moiré interferometry was investigated to try and track the in-plane strain in the SEI and silicon electrode for this purpose. Moiré can detect on the order of 10 nm changes in displacement and is therefore a useful tool in the measurement of strain. As the sample undergoes small deformations, large changes in the moiré fringe allow for measurements of displacement below the diffraction limit of light. Figure 1a shows how the moiré fringe changes as the sample grating deforms. As the sample contracts or expands, the frequency of the moiré fringe changes, and this change is proportional to the strain in the sample.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Nodal Constraint, Shear Deformation and Continuity Effects Related to the Modeling of Debonding of Laminates, Using Plate Elements

The effects of several critical assumptions and parameters on the computation of strain energy release rates for delamination and debond configurations modeled with plate elements have been quantified. The method of calculation is based on the virtual crack closure technique (VCCT), and models of the upper and lower surface of the delamination or debond that use two-dimensional (2D) plate elements rather than three-dimensional (3D) solid elements. The major advantages of the plate element modeling technique are a smaller model size and simpler configurational modeling. Specific issues that are discussed include: constraint of translational degrees of freedom, rotational degrees of freedom or both in the neighborhood of the debond front, shear deformation assumptions; and continuity of material properties and section stiffness in the vicinity of the debond front. Where appropriate, the plate element analyses are compared with corresponding two-dimensional plane strain analyses.

Glaessgen, E. H.↗
NASA NTRS1998

Effect of Shear Deformation and Continuity on Delamination Modelling with Plate Elements

The effects of several critical assumptions and parameters on the computation of strain energy release rates for delamination and debond configurations modeled with plate elements have been quantified. The method of calculation is based on the virtual crack closure technique (VCCT), and models that model the upper and lower surface of the delamination or debond with two-dimensional (2D) plate elements rather than three-dimensional (3D) solid elements. The major advantages of the plate element modeling technique are a smaller model size and simpler geometric modeling. Specific issues that are discussed include: constraint of translational degrees of freedom, rotational degrees of freedom or both in the neighborhood of the crack tip; element order and assumed shear deformation; and continuity of material properties and section stiffness in the vicinity of the debond front, Where appropriate, the plate element analyses are compared with corresponding two-dimensional plane strain analyses.

Glaessgen, E. H.↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2023

Direct Observation of Magnetic Bubble Lattices and Magnetoelastic Effects in van der Waals Cr 2 Ge 2 Te 6

Abstract Ferromagnetic van der Waals (vdW) materials are of large current interest for the fundamental study of low‐dimensional magnetism and for potential applications in multilayer heterostructures. Cr 2 Ge 2 Te 6 (CGT) is particularly exciting because it is a ferromagnetic semiconductor with tunable electronic and magnetic properties. Controlling the magnetic domain structure of CGT is a requirement for understanding its novel interface physics and for tuning behavior for potential devices. Herein, cryo‐Lorentz transmission electron microscopy is performed in the temperature range of 12–50K to directly image the magnetic domain structures in CGT. A rich phase diagram of domain structures including stripe domains, magnetic bubble lattices of mixed‐chirality, and topologically‐protected lattices of homochiral magnetic bubbles is observed. The types and chiralities of the bubbles can be controlled by topographical changes in the CGT flakes. Additionally, it is observed that in‐plane strain and magnetoelastic coupling can align and organize both bubble lattices and stripe domains. This study provides insights into creating and controlling complex magnetic domain structures for integration into multilayer heterostructures and for future studies of 2D magnetism.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Coupling Thermomechanical Processing and Alloy Design to Improve Textures in Mg-Zn-Ca Sheet Alloys

Abstract The effect of Ca and Zn additions on the microstructure and texture evolution during thermomechanical processing of Mg-Zn-Ca sheet alloys was systematically investigated and quantified. Plane strain compression testing in a Gleeble thermomechanical simulator was used to physically simulate a 10-pass rolling schedule, while allowing for careful control and monitoring of the processing parameters. Textures in the as-deformed ternary alloy samples demonstrate a weak maximum basal intensity and spreading in the transverse direction. Increasing the Zn content to 3.2 wt.% in the ternary alloys resulted in samples that exhibited weak textures in the as-deformed state. Importantly, static recrystallization (SRX) during post-deformation annealing of these alloys promoted a desirable annular texture, with the c -axis tipped from the normal direction and a lower basal texture intensity. The evolution in texture during SRX is associated with as-deformed microstructures with broad grain orientation spreads and a low degree of recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗
DOE OSTI2022

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗
DOE OSTI2020

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗
DOE OSTI2024

Effects of Zn and Ca on the deformation texture evolution of Mg–Zn–Ca alloys at elevated temperatures

This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.

Crystal plasticity finite element↗
DOE OSTI2021

Plasticity and fracture behavior of Inconel 625 manufactured by laser powder bed fusion: Comparison between as-built and stress relieved conditions

In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗