Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element (FE) simulation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗
DOE OSTI2023

The Influence of Carbon Fiber Composite Specimen Design Parameters on Artificial Lightning Strike Current Dissipation and Material Thermal Damage

Previous artificial lightning strike direct effect research has examined a broad range of specimen design parameters. No works have studied how such specimen design parameters and electrical boundary conditions impact the dissipation of electric current flow through individual plies. This article assesses the influence of carbon fiber composite specimen design parameters (design parameters = specimen size, shape, and stacking sequence) and electrical boundary conditions on the dissipation of current and the spread of damage resulting from Joule heating. Thermal-electric finite element (FE) modelling is used and laboratory scale (<1 m long) and aircraft scale (>1 m long) models are generated in which laminated ply current dissipation is predicted, considering a fixed artificial lightning current waveform. The simulation results establish a positive correlation between the current exiting the specimen from a given ply and the amount of thermal damage in that ply. The results also establish that the distance to ground, from the strike location to the zero potential boundary conditions (ground), is the controlling factor which dictates the electric current dissipation in each ply. Significantly, this distance to ground is dependent on each of the specimen shape, dimensions, stacking sequence, and location of ground boundary conditions. Therefore, it is not possible to decouple current dissipation and damage from specimen design and boundary condition setup. However, it is possible to define a specimen size for a given specimen shape, stacking sequence, and waveform which limit the influence of specimen dimensions on the resulting current distribution and damage. For a rectangular specimen design which appears in literature multiple times, as 100 × 150 mm and with a stacking sequence of [45/0/-45/90] 4s , a specimen design of greater than 300 × 200 mm is required to limit the influence of specimen dimensions on current distribution and damage.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Conceptual design of a high reactive-power ferroelectric fast reactive tuner

We present a novel design of a ferroelectric fast reactive tuner (FE-FRT) capable of modulating mega-VAR reactive power on a submicrosecond timescale. The high reactive power capability of our design extends the range of applications of reactive tuners to numerous applications. We present a detailed analytical model of the performance of a megawatt-class reactive power device and benchmark it against finite-element method eigenmode and frequency domain electromagnetic simulations. We introduce new features, including an annulus design for the ferroelectric capacitors and capacitive window coupling to the cavity. We consider thermal design issues and nonlinear effects in the ferroelectric. The model covers several configurations, allowing control of the frequency of superconducting and normal-conducting cavities in a variety of applications and frequencies. We calculate that the FE-FRT designed should be capable of handling around 0.45 MVAR of reactive power with around 3 kW of resistive losses, providing a frequency tuning range of 8 kHz in an example of 400 MHz cavity geometry. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Novel Tube Design for Superheater Heat Exchanger Enabled Via Additive Manufacturing

Superheater tubes are critical boiler components that operate at relatively higher temperatures and pressure. Amongst the primary concerns for these tubes is the deposition of ash particles on the tube surface, leading to the reduced thickness of the tube due to material corrosion, consequently causing early creep failure of the component. In this research, a novel tube design has been proposed which resembles a teardrop or ogive shape to reduce the drag and concurrently improve the creep life of the superheater tubes. To administer the practicality of novel tubes, metal additive manufacturing (AM), for instance, laser-powder bed fusion (L-PBF), has been proposed. These unconventional designs were assessed and compared with the baseline circular tube design for mechanical design requirements (hoop stress and creep life) and the particle and flue gas flow characteristics around the differently shaped tubes. A thermomechanical finite element (FE) analysis was performed for hoop stress calculations. This study also emphasizes on effect of circumferential thermal variation on hoop stress distribution in tubes. Therefore, a detailed two-dimensional (2D) thermal simulation has been performed to report the circumferential thermal variation on the tube. A computational fluid dynamics (CFD) analysis coupled with particle tracing was performed for gas flow visualization and particle tracing around the proposed shapes and baseline circular-shaped tube design. The Schlieren optic setup was built and leveraged for qualitative validation of the proposed design. The complete design methodology established in the paper shows teardrop-shaped tubes better in terms of drag and creep life in contrast to the circular-shaped tube.

42 ENGINEERING↗
DOE OSTI2021

High-power ferro-electric fast reactive tuner

We present a novel design of a FerroElectric Fast Reactive Tuner (FE-FRT) capable of modulating Mega VAR reactive power on a sub-microsecond time scale. We show detailed analytical estimates of the performance of this device and benchmark these estimates against finite element method eigenmode and frequency domain electromagnetic simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Characterization of Fe-6Si Soft Magnetic Alloy Produced by Laser-Directed Energy Deposition Additive Manufacturing

Abstract Commercial electrical steels, Fe-Si alloys with < 4 wt.% Si, are inexpensive and efficient materials for electrical power conversion. Further efficiency improvements require increasing the silicon concentration to 6 wt.%, at which point the material becomes brittle and difficult to form by conventional rolling and sheet fabrication methods. Additive manufacturing stands to overcome challenges with commercial manufacturing techniques by leveraging near-net-shape fabrication. The wide array of process conditions provides additive manufacturing with increased flexibility, enabling control over the microstructure and mechanical properties. This work explores the microstructures and magnetic properties of ring-shaped Fe-Si alloys produced using concentric and cross-hatch tool paths on a laser-directed energy deposition additive manufacturing system. Concentric-built samples exhibit elongated grain structures while cross-hatch-built samples comprise lower aspect ratio grain structures. Thermal finite element analysis simulations model the stress conditions produced by the different scan path geometries. Microhardness measurements probe the mechanical properties as a function of anneal temperature, providing a qualitative understanding of the intergranular defect density. Soft magnetic properties measured under quasistatic and AC conditions show frequency- and microstructure-dependent coercivity and permeability. Finally, analysis of the core loss quantifies how the build strategies and thermal treatments influence efficiency in electrical power conversion applications. Understanding the influences of scan path geometry and thermal treatment provides a pathway towards application of additively manufactured soft magnetic materials.

Adamczyk, Jesse M.↗
DOE OSTI2020

Biaxial steel plated concrete constitutive models for composite structures: Implementation and validation

The Steel-plated Concrete (SC) technique is an alternative construction technique with faster onsite construction speed, reduced construction time, and increased structural performance. Aiming to predict the force transfer mechanism of SC elements, an innovative constitutive model package is proposed by implementing experimental-based biaxial steel plate concrete models into the nonlinear finite element (FE) model “Membrane Model of SC (MM-CS).” First, the formulation and implementation of the MM-SC is illustrated in detail, including the equilibrium and compatibility equations and the implementation of constitutive models. Next, various experimental data of SC members are selected and simulated using the proposed MM-SC model. In this research, two types of structures, SC panels and framed SC shear walls, and three types of loading conditions, including uniaxial compression, pure shear, and combined axial-flexural-shear tests, are analyzed respectively. Good agreements were obtained between the reported results and the FE simulation results in terms of yield capacity and ultimate capacity, proving the reliability of the implemented analytical constitutive material models and the biaxial membrane model formulations.

42 ENGINEERING↗
DOE OSTI2024

A Study on Particle Trajectory Error in Finite-Element Particle-in-Cell Algorithms

Particle-in-cell (PIC) algorithms are widely used for the simulation of kinetic plasmas. PIC algorithms account for the interaction between charged particles in a plasma and the electromagnetic field in ambient space, including self-field interactions. The objective of this article is to study the error in charged particle trajectories present in finite-element (FE)-based PIC algorithms on unstructured meshes. We study how the trajectory error behaves according to the FE mesh resolution and the matrix solver employed in the FE algorithm. The study is performed by considering a trajectory established by a parabolic electric potential and an axial magnetic force acting on the charged particle. Under a proper choice combination of electric and magnetic field parameters, the 2-D particle trajectories comprise closed orbits. Numerical errors cause small perturbations on the trajectories, with cumulative effects. As a result, the resulting orbital trajectories exhibit distortions including spurious apsidal precession. These distortions provide a clear imprint of the numerical error. Here, we also study the numerical error in a quantitative fashion by computing the distance norm between the trajectories obtained by the exact fields and numerical fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hydro-mechanical coupled modeling of hydraulic fracturing using the mesh fragmentation technique

The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.

58 GEOSCIENCES↗
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
DOE OSTI2022

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗
DOE OSTI2022

Mach Conference, 2022.

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiC/SiC) composite tubes due to their superior mechanical properties. The SiC/SiC tubes are known to remain stable even after prolonged exposure to radiation. However, there are many variables in the manufacturing of SiC/SiC composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. The intention is to identify the distribution of elastic constants of SiC/SiC tubes which would help model the performance of the novel cladding material. To do so, a sensitivity analysis will be carried out to determine the dominant variables that influence the elastic constants of SiC/SiC tubes. A chosen number of combinations of dominant variables will be sampled through the Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in the circumferential and longitudinal direction, Poisson's ratio, and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used as a surrogate to FE analysis. Once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The uncertainty can be quantified by calculating the coefficient of variation (CV) based on MCS-generated data which will be presented at the conference.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Environmentally Assisted Fatigue in Light Water Reactor Environment

This report summarizes the Environmentally Assisted Fatigue (EAF) research conducted at ANL under the US DOE Light Water Reactor Sustainability (LWRS) program. Starting from a rich background in theoretical and experimental EAF, ANL previously developed an approach to evaluate fatigue performance of reactor materials in light water reactor environments with the correction factor F en . The approach was based on a large body of experimental work performed at ANL and elsewhere, and was consistent with American Society of Mechanical Engineers (ASME)’s methodology governing the design and construction of reactor components. In recent years, the program was focused on component fatigue prediction and made several major and fundamental contributions in this area. These accomplishments help meet the needs identified by the industry concerning component level fatigue predictions in complex, transient conditions. The main contribution of the ANL program involved the development of a system-level model for estimating residual strain and life of nuclear reactor coolant system components under connected-system-thermal-mechanical boundary conditions. The goal was to predict the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the finite element (FE) model results, the strain residuals, strain amplitudes and resulting fatigue lives of reactor coolant system (RCS) components were predicted. The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. In addition, the simulated component-level strain profile can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests. Building upon the system-level model, ANL developed a digital twin (DT) framework to predict the structural states and associated fatigue life of components in real-time. This framework is a comprehensive system designed to predict the structural states and fatigue lives of reactor components. It includes multiple models and integrates artificial intelligence (AI), machine learning (ML), and FE based modeling tools to evaluate the structural states and fatigue lives.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Carbide Particle Redistribution in U-10Mo Alloy during Hot and Cold Rolling Processes: Synthetic Microstructure-Based Finite-Element Approach

The influence of initial carbide particle distribution features on particle fracture and the evolution of stringer distribution within an alloy of uranium with 10 weight percent molybdenum during cold rolling was investigated in this study. This work will provide guidance to casting optimization work and identify actionable process changes to eliminate the surface features directly attributed to carbides. In a previous study, various synthetic microstructure-based finite-element (FE) models were generated with consideration of different parameters (i.e., particle shape, particle volume fraction, and particle distribution), and used to examine the effects of hot rolling on the microstructure evolution. In this study, microstructures after ~80% hot rolling were used as the initial microstructures for the subsequent cold rolling simulation. Plane-strain compression loading was employed in the FE models to simulate cold rolling with reductions from ~80% to ~94%. Stringer analyses with the consideration of particle fracture were then performed with the simulated cold rolling results. The results of void evolution analyses show that microstructures with elliptical particles incorporate a much larger (almost double) void volume fraction than those with circular particles, indicating that the elliptical particles are fractured more easily than the circular ones. The results from the fracture mode analyses show that, under the adopted simple compression loading assumption, almost all the particles (i.e., >99%) were fractured by tensile separation. The results of stringer analyses show that, during cold rolling, the stringer volume fractions increase substantially due to the particle fracture, especially for the microstructures with the elliptical particles, and that the stringer volume fractions from stringers comprising many particles tend to increase faster. The maximum stringer lengths for the microstructures with circular particles tend to remain constant or decrease during cold rolling, whereas those for elliptical particles tend to slightly increase after a slight decrease. The number of particles at the surface increases very rapidly due to particle fracture during cold rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗