Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peridynamics”

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 109 records · Page 6

Data-driven nonlocal model for fragmentation in the crushing of solids

A technique is proposed for reproducing particle size distributions in three-dimensional simulations of the crushing and comminution of solid materials. The method is designed to produce realistic distributions over a wide range of loading conditions, especially for small fragments. In contrast to most existing methods, the new model does not explicitly treat the small-scale process of fracture. Instead, it uses measured fragment distributions from laboratory tests as the basic material property that is incorporated into the algorithm, providing a data-driven approach. The algorithm is implemented within a nonlocal peridynamic solver, which simulates the underlying continuum mechanics and contact interactions between fragments after they are formed. Finally, the technique is illustrated in reproducing fragmentation data from drop weight testing on sandstone samples.

58 GEOSCIENCES↗

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗

A Review of Local-to-Nonlocal Coupling Methods in Nonlocal Diffusion and Nonlocal Mechanics

Local-to-nonlocal (LtN) coupling refers to a class of methods aimed at combining nonlocal and local modeling descriptions of a given system into a unified coupled representation. This allows to consolidate the accuracy of nonlocal models with the computational expediency of their local counterparts, while often simultaneously removing nonlocal modeling issues such as surface effects. The number and variety of proposed LtN coupling approaches have significantly grown in recent years, yet the field of LtN coupling continues to grow and still has open challenges. This review provides an overview of the state of the art of LtN coupling in the context of nonlocal diffusion and nonlocal mechanics, specifically peridynamics. Furthermore, we present a classification of LtN coupling methods and discuss common features and challenges. The goal of this review is not to provide a preferred way to address LtN coupling but to present a broad perspective of the field, which would serve as guidance for practitioners in the selection of appropriate LtN coupling methods based on the characteristics and needs of the problem under consideration.

97 MATHEMATICS AND COMPUTING↗

Identifying crack tip position and stress intensity factors from displacement data

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, in this work, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

36 MATERIALS SCIENCE↗

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗

Resiliency of Degraded Built Infrastructure

Infrastructure resiliency depends on the ability of infrastructure systems to withstand, adapt, and recover from chronic and extreme stresses. In this white paper, we address the resiliency of infrastructure assets and discuss improving infrastructure stability through development of our understanding of cement and concrete degradation. The resiliency of infrastructure during extreme events relies on the condition, adaptability, and recoverability of built infrastructure (roads, bridges, dams), which serves as the backbone of existing infrastructure systems. Much of the built infrastructure in the US has consistently been rated D+ by the American Society of Civil Engineers (ASCE). Aged infrastructure introduces risk to the system, since unreliable infrastructure increases the likelihood of failures under chronic and extreme stress and are particularly concerning when extreme events occur. To understand and account for this added risk from poor infrastructure quality, more research is needed on (i) how the changing environment alters the aging of new and existing built infrastructure and (ii) how degradation causes unique failure mechanisms. The aging of built infrastructure is based on degradation of the structural materials, such as concrete and steel supports, which causes failure. Current work in cement/concrete degradation is based on (i) the development of high strength and degradation resistance concrete mixtures, (ii) methods of assessing the age and reliability of existing structures, and (3) modeling of structural stability and the microstructural evolution of concrete/cement from degradation mechanisms (sulfide attack, carbonation, decalcification). Sandia National Laboratories (SNL) has made several investments in studying the durability and degradation of cement based materials, including using SNL-developed codes and methodologies (peridynamics, PFLOTRAN) to focus on chemo-mechanical fracture of cement for energy applications. Additionally, a recent collaboration with the University of Colorado Boulder has included fracture of concrete gravity dams, scaling the existing work to applications in full sized infrastructure problems. Ultimately, SNL has the experience in degradation of cementitious materials to extend the current research portfolio and answer concerns about the resilience of aging built infrastructure.

36 MATERIALS SCIENCE↗

Sierra/SolidMechanics 5.0 Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.0 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and $\textit{J}$-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗

Sierra/SolidMechanics 4.58. Capabilities In Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.58 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and /-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations

97 MATHEMATICS AND COMPUTING↗

Local-nonlocal coupling in Emu/PDMS

A technique called the splice method for coupling local to peridynamic subregions of a body is described. The method relies on ghost nodes, whose values of displacement are interpolated from nearby physical nodes, to make each subregion visible to the other. In each time step, the nodes in each subregion treat the nodes in the other subregion as boundary conditions. Adaptively changing the subregions is possible through the creation and deletion of ghost nodes. Example problems in 2D and 3D illustrate how the method is used to perform multiscale modeling of fracture and impact events within a larger structure.

42 ENGINEERING↗

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗

Understanding Phase and Interfacial Effects of Spall Fracture in Additively Manufactured Ti-5Al-5V-5Mo-3Cr

Additive manufactured Ti-5Al-5V-5Mo-3Cr (Ti-5553) is being considered as an AM repair material for engineering applications because of its superior strength properties compared to other titanium alloys. Here, we describe the failure mechanisms observed through computed tomography, electron backscatter diffraction (EBSD), and scanning electron microscopy (SEM) of spall damage as a result of tensile failure in as-built and annealed Ti-5553. We also investigate the phase stability in native powder, as-built and annealed Ti-5553 through diamond anvil cell (DAC) and ramp compression experiments. We then explore the effect of tensile loading on a sample containing an interface between a Ti-6Al-V4 (Ti-64) baseplate and additively manufactured Ti-5553 layer. Post-mortem materials characterization showed spallation occurred in regions of initial porosity and the interface provides a nucleation site for spall damage below the spall strength of Ti-5553. Preliminary peridynamics modeling of the dynamic experiments is described. Finally, we discuss further development of Stochastic Parallel PARticle Kinteic Simulator (SPPARKS) Monte Carlo (MC) capabilities to include the integration of alpha (α)-phase and microstructural simulations for this multiphase titanium alloy.

36 MATERIALS SCIENCE↗

Small Particle Impact Damage on Different Glass Substrates

Impact experiments using sand particles were performed on four distinct glass substrates. The sand particles were characterized using the X-Ray micro-CT technique; 3-D reconstruction of the particles was followed by further size and shape analyses. High-speed video footage from impact tests was used to calculate the incoming and rebound velocities of the individual sand impact events, as well as particle volume. Further, video analysis was used in conjunction with optical and scanning electron microscopy to relate the incoming velocity and shape of the particles to subsequent fractures, including both radial and lateral cracks. Analysis was performed using peridynamic simulations.

Waxman, R.↗

NASA ACC High Energy Dynamic Impact Methodology and Outcomes

For High Energy Dynamic Impact (HEDI) events, testing to evaluate the structural response of primary aircraft structure for design and certification is both expensive and time consuming. This paper discusses current work seeking to assess, develop, and validate appropriate analytical models that accurately predict physical response, damage, and failure modes for large scale composite structures in HEDI events. Four state-of-the-art Progressive Damage Analysis (PDA) methods were employed for this phased project: LS-DYNA MAT162, LS-DYNA MAT261, Smoothed Particle Galerkin (SPG), and EMU Peridynamics. Multiple material systems were considered, namely T700/5208 textile-infusion triaxial braid, T800/AMD-825 textile-infusion triaxial braid, IM7/8552 uni-directional tape, and SPG 196-PW/8552 plain-weave fabric. Extensive ballistic impact testing was performed to support this activity and measured results were compared to predictive models for assessment using panel delamination, panel displacement, force at the load cells, and threshold velocity as measures. Ultimately, the work under this activity provided significant progress in advancing the state-of-the-art in the use of PDA for HEDI events. Each material model had favorable performance comparing to test in some parameters and needed improvement in others. With the lessons learned from this activity, significant progress was made in the ability to predict panel behavior for a more general case beyond the flat panel in a ballistic impact event. Subsequent Phase II of the NASA ACC HEDI effort will continue to build on the coupon testing, flat panel ballistic impact testing, and analysis performed to-date with application of the PDA methods for intended material selections to test articles with greater complexity of configuration, curvature, and scale. It is not the intention of this paper to present a full set of data, but rather to give an overview of the NASA HEDI effort and show a small representative subset of the test and analysis results.

Hunziker, Kenneth J.↗

NASA ACC High Energy Dynamic Impact Methodology and Outcomes

High fidelity analysis methods known as Progressive Damage Analysis (PDA) methods, capable of reliably predicting the onset and progression of damage in composite materials, are being developed for HEDI event simulation. Four state-of-the-art PDA methods are being investigated: LS-DYNA MAT162, LS-DYNA MAT261, Smooth-Particle Galerkin (SPG), and EMU Peridynamics. Validation of these PDA models aims to follow a building block approach, starting with coupon- and component-level development. As part of the overall effort, material characterization testing was performed to bridge gaps between existing experimental data and the material property inputs required to predict ballistic impact behavior at the component-level. The material models were updated with the results of the coupon-level testing and were then used to predict panel behavior and damage in ballistic impact events. To assess the accuracy of the PDA methods, these pre-test predictions were compared to the results of ballistic impact testing.

Hunziker, Kenneth J.↗

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

Engineering the Interface: Advanced Surface Technologies for Lunar Dust Management and Equipment Longevity

Through the Artemis program, NASA intends to develop a sustainable human foothold on the Moon, ultimately paving the way for crewed exploration of Mars. The Moon's hostile environment poses numerous obstacles, including exposure to radiation, temperature extremes, micrometeoroid threats, and particularly the persistent problem of lunar dust. Lunar dust impacts nearly every aspect of surface operations through adhesion and abrasion mechanisms, with contamination from anthropogenic activities (landing, rovers) far outweighing natural phenomena. Multiple adhesion pathways contribute to surface contamination in the lunar environment, including van der Waals forces, electrostatic forces, chemical reaction, and magnetic forces from elemental iron deposits. Sharp asperities from micrometeoroid bombardment and atmospheric absence increase interaction potential and enable mechanical interlocking. Low cohesion between dust particles exacerbates these challenges, as minimal interaction potential between dust and nearby surfaces overcomes particle cohesion, causing contamination. Lunar dust adhesion mitigation technologies can be categorized as either active, requiring external energy, or passive, relying on intrinsic material properties. Ultrasonic and electrodynamic technologies have been developed to the highest technology readiness level for active approaches. Passive strategies primarily focus on surface chemistry and topography modifications. At NASA Langley Research Center, approaches include surface migration agents to reduce surface energy, topographical modification using laser ablation patterning, and tailored surface conductivity to reduce intrinsic adhesion force. Performance has been evaluated using custom-built ultrasonic and centrifuge instruments. Plume-surface interactions from lunar landers can propel micrometer-sized particles at velocities up to 1000 m s-1.8 These particles pose risks to landers, habitats and infrastructure, leading to erosion, degradation, and reduced component lifespan. A panel recovered from Surveyor III was determined to have been severely abraded because of lunar dust displaced from the Apollo 12 lunar module that landed 160 m away. The performance of metallic surfaces has been evaluated via high velocity single particle impact using the laser-induced project impact test (LIPIT) facility at the University of Utah. Peridynamics modeling, a form of continuum mechanics that uses a nonlocal approach enabling greater simulation capabilities of crack initiation and fracture, has also been utilized to gain greater insight into material response during impact events. Lunar dust contamination challenges extend to power generation systems and moving equipment. Cables, rotation stages, and other mechanisms may experience limited range of motion and reduced lifetime due to dust infiltration. NASA Langley Research Center has evaluated traditional aerospace alloys, softgoods, wear resistant ceramics, and several polymer and polymer composite materials. Test methods have included traditional techniques like Taber abrasion testing, as well as designed test configurations developed in the DUSTE (dust, ultraviolet radiation, and space thermal environmental) chamber that reproduce mechanism functions in operational environment. Beyond laboratory experiments, several flight experiments have been conducted. Materials were exposed to the low Earth orbit environment on the Materials International Space Station Experiment (MISSE) and to the lunar surface environment through the Aegis Aerospace Regolith Adherence Characterization (RAC) payload and the Honeybee Robotics PlanetVac payload. Determining lunar dust's impact on surface exploration and habitation requires comprehensive experimental and computational capabilities combined with lessons learned from initial lunar activities. Identifying the greatest environmental challenges and developing mitigation technologies provides the clearest path toward successfully, expeditiously, and efficaciously completing NASA's mission. This presentation will discuss ongoing efforts at NASA Langley Research Center and collaborator contributions to these critical objectives.

Surface Engineering↗

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↗