Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite deformation”

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 217 records · Page 12

Direct observation on supersonic microprojectile penetration of carbon fiber composites with ultrafast synchrotron X-ray phase contrast imaging

High-speed penetration into carbon fiber composites is of fundamental importance to materials science and impact engineering, but research along this line suffers considerably from the lack of direct experimental observations. Here we investigate such penetration dynamics of a unidirectional carbon fiber reinforced epoxy (UCFRE) composite, with a combination of in situ, ultrafast, synchrotron phase contrast imaging and finite element (FE) analysis. The experiments yield the first direct observation on projectile trajectories and fiber-scale deformation and damage in the UCFRE composites during supersonic microprojectile penetration, for different fiber orientations (0 degrees - 90 degrees from the impact direction) and projectile velocities (600 - 850 m.s -1 ), at unprecedented temporal (~ 100 ps) and spatial (5 μm) resolutions. The maximum penetration depth decreases with increasing fiber orientation angles, as a result of anisotropic damage evolution in the composite sample. Strain localizations are prone to develop along a direction perpendicular to the fiber orientation, while the damage or cavity region, along the fiber direction. FE modeling with a three-dimensional Hashin criterion yields consistent projectile trajectory and cavity morphology with the experimental results. Finally, with increasing fiber orientation angles, damage analyses show a transition in the damage mode from fiber compression to matrix compression damage, in line with the increasing maximum penetration depth.

36 MATERIALS SCIENCE↗

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗

Structural Analysis of Multilayer Tapered Canted-Cosine-Theta Superconducting Magnet for the Interaction Region of Electron Ion Collider

The electron ion collider (EIC) will collide high energy and highly polarized hadron and electron beams with luminosities up to 10 34 cm -2 s -1 . The magnet designs at the interaction region (IR) are challenging due to the close proximity of the hadron and electron beams. A canted cosine theta (CCT) quadrupole magnet was designed due to the space restrictions at the IR. The CCT magnet has a tapered aperture, which allows higher gradients at the interaction point (IP) side of the magnet. The tapered geometry makes the magnet not symmetrical along the axial direction. In conclusion, to validate the design in terms of the mechanical strength and the allowed deformations, a detailed structural mechanical analysis for the tapered CCT quadrupole magnet has been carried out using the finite element method (FEM) considering both the shrinking due to cool down and the electromagnetic force on the conductors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

$T\overline{T}$ deformations and the width of fundamental particles

We provide a simple geometric meaning for deformations of so-called $T\overline{T}$ type in relativistic and non-relativistic systems. Deformations by the cross products of energy and momentum currents in integrable quantum field theories are known to modify the thermodynamic Bethe ansatz equations by a “CDD factor”. In turn, CDD factors may be interpreted as additional, fixed shifts incurred in scattering processes: a finite width added to the fundamental particles (or, if negative, to the free space between them). We suggest that this physical effect is a universal way of understanding $T\overline{T}$ deformations, both in classical and quantum systems. We first show this in non-relativistic systems, with particle conservation and translation invariance, using the deformation formed out of the densities and currents of particles and momentum. This holds at the level of the equations of motion, and for any interaction potential, integrable or not. We then argue, and show by similar techniques in free relativistic particle systems, that $T\overline{T}$ deformations of relativistic systems produce the equivalent phenomenon, accounting for length contractions. We also show that, in both the relativistic and non-relativistic cases, the width of particles is equivalent to a state-dependent change of metric, where the distance function discounts the particles’ widths, or counts the additional free space. This generalises and explains the known field-dependent coordinate change describing $T\overline{T}$ deformations. The results connect such deformations with generalised hydrodynamics, where the relations between scattering shifts, widths of particles and state-dependent changes of metric have been established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗

Study of Microphonic Effects on the C100 Cryomodule for High Energy Electron Beam Accelerators

The Continuous Electron Beam Accelerator Facility (CEBAF) at Thomas Jefferson National Laboratory (JLab) is a particle accelerator which can accelerate an electron beam to relativistic speeds and apply the beam onto target samples. The C100 superconducting radio frequency (SRF) cavity is the primary accelerating structure of the C100 cryomodule, one of the many cryomodules which compose the CEBAF linear accelerator. SRF cavities are particularly sensitive to internal and external vibrations that can result in a phenomenon called microphonics which degrade the operational stability of a cryomodule. The purpose of this thesis is to investigate the significance of mechanical disturbances on the electromagnetic resonant frequency of a C100 SRF cavity. Knowledge of the mechanical resonance of the cavities and cryomodule sheds light into how these disturbances are most easily realized as deformation which causes radio frequency (RF) detuning. Three studies were conducted: the development and hammer test calibration of a Finite Element Analysis (FEA) model of a C100 cavity, the development and hammer test calibration of an FEA model of a C100 cavity string, and the hammer test of the C100-10R cryomodule at the Cryomodule Test Facility (CMTF). The cavity FEA model was found to accurately predict two modes found in two real cavities in a simply supported configuration. The cavity string FEA model leveraged the calibrated cavity FEA model but was not found to accurately predict the modal behavior of a real cavity string. Even so, the modal behavior of the cavity string inside the C100-10R cryomodule was captured during a hammer test while it was partially assembled. Finally, the C100-10R cryomodule was placed in the CMTF to study RF detuning. The RF detuning spectra during hammer hits and background noise was captured. The results of the hammer testing indicate two strong peaks at low frequencies (9-10 Hz and 22-23 Hz). These two frequencies were found to be nearly coincident to four instances of mechanical resonance found during the hammer testing done on the partially-assembled C100-10R. Because of this, these two modes are believed to contribute to RF detuning of the cryomodule. This test event also included the testing of the effectiveness of a configuration of BNNT canisters designed to act as dampers. While these tests show promising results, the lurking variables render these tests somewhat inconclusive.

Hull, Caleb James↗

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗

(2 + δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics

The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2 + δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2 + δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. Finally, this model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling and experimental validation of dynamical effects in Bragg coherent x-ray diffractive imaging of finite crystals

Bragg coherent diffractive imaging (BCDI) is a noninvasive microscopy technique that can visualize the shape and internal lattice deviations of crystals with nanoscale spatial resolution and picometer deformation sensitivity. Its strain imaging capability relies on Fourier transform–based iterative phase retrieval algorithms, which are mostly developed under the kinematical approximation. Such approximation prohibits the application of BCDI on larger crystals, which are commonly seen in most emerging functional materials. Understanding the dynamical effect in BCDI, as well as developing a validated method for modeling BCDI at the dynamical diffraction limit, is crucial for applying BCDI to hierarchical systems that contain micron-sized crystals and grains. Thus we report a comparative study on the impact of dynamical diffraction effects by comparing the reconstruction results from two measurements of the same crystal. Forward simulation is implemented to show subtle changes of interference fringes in the diffraction pattern due to the dynamical diffraction, and is compared directly with the experimental data.

36 MATERIALS SCIENCE↗

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗

Multiscale simulation of spatially correlated microstructure via a latent space representation

When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.

Elastoplasticity↗

Distributed fiber-optic sensing in a subscale high-temperature superconducting dipole magnet

High-temperature superconductors, such as REBa2Cu3O7−x (REBCO, RE = rare earth), are becoming pivotal for high-field magnet technology for future circular colliders and compact fusion reactors. The U.S. Magnet Development Program, in collaboration with industry, is developing REBCO magnet technology using round conductors consisting of multiple REBCO tapes. For these multi-tape cables, traditional instrumentation, such as voltage taps and resistive strain gauges, become insufficient to help measure and understand the performance-limiting factors in these model magnets. Distributed fiber-optic sensing (DFOS) is a potential solution to address this challenge. Although DFOS is well established for various applications, measuring temperature and strain in high-temperature superconducting magnets is in its infancy. Here we report the detailed implementation and test results of DFOS based on Rayleigh scattering in a subscale canted cosθ (CCT) dipole magnet using high-temperature superconducting CORC® wires. We co-wound optical fibers in each layer of the CCT magnet and compared different types of commercial fibers and mold-release agents to reduce the power attenuation in the fibers. The DFOS allowed us to measure mechanical deformation and temperature along the conductor during tests at 77 and 4.2 K. The measured strain agreed quantitively with a finite-element mechanical model of the subscale magnet. Our results indicate that DFOS can effectively identify locations of strain and temperature changes, offering unique insight into magnet performance that can advance our understanding and development of the REBCO magnet technology for high-energy physics and fusion applications.

Luo, Linqing↗

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗

Impact of control blade insertion on the deformation behavior of SiC-SiC channel boxes in BWRs

This article describes the analysis of distortion of a silicon carbide fiber-reinforced, silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of the temperature and neutron flux dependent irradiation-induced swelling of SiC, the SiC-SiC composite channel box can undergo distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermal-hydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that large gradients in fast neutron flux (up to 35–40% across a single axial level) will develop across the channel box. Due to these gradients, the channel box will undergo time-dependent bending for all the control blade positions in the assembly. The time-dependent bowing behavior is dominated by the transient swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The bending will cause temporary interference between the channel box and control blade, and the interference is expected to be most severe for the fully inserted control blade position. The developed stresses due to differential swelling in the channel box exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the tensile strength of the material, and therefore, development of a full, through-thickness crack in the channel box is not expected. Finally, further work is recommended to explore and evaluate the mitigation strategies.

42 ENGINEERING↗

Line operators in Chern-Simons-Matter theories and Bosonization in Three Dimensions II: Perturbative analysis and all-loop resummation

We study mesonic line operators in Chern-Simons theories with bosonic or fermionic matter in the fundamental representation. In this paper, we elaborate on the classification and properties of these operators using all loop resummation of large N perturbation theory. We show that these theories possess two conformal line operators in the fundamental representation. One is a stable renormalization group fixed point, while the other is unstable. They satisfy first-order chiral evolution equations, in which a smooth variation of the path is given by a factorized product of two mesonic line operators. The boundary operators on which the lines can end are classified by their conformal dimension and transverse spin, which we compute explicitly at finite ’t Hooft coupling. We match the operators in the bosonic and fermionic theories. Finally, we extend our findings to the mass deformed theories and discover that the duality still holds true.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deformation twinning and detwinning in extruded Mg-4Al: In-situ experiment and crystal plasticity simulation

Deformation twinning and detwinning in extruded Mg-4Al were investigated using in-situ SEM-DIC experiments and crystal plasticity finite element (CPFE) simulation. In this study, the in-situ SEM-DIC method was used to provide a unique set of data including twin/detwin characteristics and twin area fraction in addition to strain maps. A statistical analysis of the activation of twin variants and twin area fraction during both twinning and detwinning was conducted. A strong correlation was found between twin growth/shrinkage and the Schmid Factor (SF) for individual twin variants Higher twin SF during loading and unloading led to higher twin growth and shrinkage, respectively. However, after the applied compressive strain was removed, the pattern of the twin area fractions of the residual twin variants versus their nominal SFs did not follow the trend observed at the maximum compressive strain. Using a systematic methodology and an advanced twin/detwin model, the PRISMS-Plasticity CPFE simulation was calibrated using experimentally determined stress versus strain and twin area fraction versus strain information. A comprehensive evaluation of the CPFE model was conducted to determine its ability to capture the statistics of twin variants activation and twin area fraction. CPFE accurately captured the statistical aspects of both twinning and detwinning. Finally, it also predicted the first dominant twin variant for 47.5% of the grains and at least one of the two dominant twin variants for 80% of the grains at maximum compressive strain.

36 MATERIALS SCIENCE↗