Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model calibration”

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

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↗
NASA NTRS2019

Structural Analysis Methods for the Roll-Out Solar Array Flight Experiment

The Roll-Out Solar Array (ROSA) flight experiment was launched to the International Space Station (ISS) on June 3rd, 2017. ROSA is an innovative, lightweight solar array with a flexible substrate that makes use of the stored strain energy in its composite structural members to provide deployment without the use of motors. This paper discusses the effort to model the structural dynamics of ROSA using finite element modeling. Two distinct and agnostic approaches were used by separate teams to assess the structural dynamics of the solar array prior to ground vibrational testing and flight testing. Results from each approach are compared to measured dynamics from accelerometers and photogrammetry data gathered on orbit. Advantages and disadvantages of each approach are discussed as are preliminary efforts to calibrate the models to the empirical data for the benefit of future modeling efforts on similar space structures.

Chamberlain, Matthew K.↗
DOE OSTI2020

Simulation of emergent compaction banding fronts caused by frictional boundaries

This study examines the role of boundary friction in promoting heterogeneous compaction in soft rock specimens loaded at high confining pressure outside the domain of compaction localisation. An elastoplastic constitutive model characterised by tunable hardening/softening behaviour is used to conduct the analyses. Finite-element simulations suggest that material instability is a non-necessary condition for the emergence of compaction fronts. Such fronts propagated as a result of a severe deviation in the local responses induced by frictional constraints. Overall, these findings suggest that boundary effects can bias the assessment of the extent of the compaction localisation domain. Experimental countermeasures and informed model calibration procedures are therefore necessary to minimise such bias and enable more accurate predictions of soft rock compaction.

42 ENGINEERING↗
DOE OSTI2023

Results of the Micromorphic Upscaling Workflow for the PSAAP III Year 3 Report

Predicting the mechanical response and failure of heterogeneous materials has proven difficult. Multiscale numerical methods based in higher order continuum theories attempt to bridge the gap between microscale and macroscale structural behavior. Micromorphic continuum theories have shown promise. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling workflows. These include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This work discusses micromorphic upscaling efforts of the University of Colorado Boulder PSAAP III multidisciplinary simulation center (MSC). Verification studies are presented that compare the accuracy of the upscaling workflow with analytical solutions for a trivial stress state and homogeneous material using direct numerical simulations (DNS) conducted in the Ratel finite element method (FEM) and GEOS material point method (MPM) codes. These verification studies consider upscaling using the Micromorphic Filter for a “single filter domain”. Finally, DNS of a heterogeneous composite material is upscaled using a “multiple filter domain” method.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Numerical Analysis of a Radiant Heat Flux Calibration System

A radiant heat flux gage calibration system exists in the Flight Loads Laboratory at NASA's Dryden Flight Research Center. This calibration system must be well understood if the heat flux gages calibrated in it are to provide useful data during radiant heating ground tests or flight tests of high speed aerospace vehicles. A part of the calibration system characterization process is to develop a numerical model of the flat plate heater element and heat flux gage, which will help identify errors due to convection, heater element erosion, and other factors. A 2-dimensional mathematical model of the gage-plate system has been developed to simulate the combined problem involving convection, radiation and mass loss by chemical reaction. A fourth order finite difference scheme is used to solve the steady state governing equations and determine the temperature distribution in the gage and plate, incident heat flux on the gage face, and flat plate erosion. Initial gage heat flux predictions from the model are found to be within 17% of experimental results.

Jiang, Shanjuan↗
NASA NTRS2011

Comparison of Test and Finite Element Analysis for Two Full-Scale Helicopter Crash Tests

Finite element analyses have been performed for two full-scale crash tests of an MD-500 helicopter. The first crash test was conducted to evaluate the performance of a composite deployable energy absorber under combined flight loads. In the second crash test, the energy absorber was removed to establish the baseline loads. The use of an energy absorbing device reduced the impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to analytical results. Details of the full-scale crash tests and development of the system-integrated finite element model are briefly described along with direct comparisons of acceleration magnitudes and durations for the first full-scale crash test. Because load levels were significantly different between tests, models developed for the purposes of predicting the overall system response with external energy absorbers were not adequate under more severe conditions seen in the second crash test. Relative error comparisons were inadequate to guide model calibration. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used for the second full-scale crash test. The calibrated parameter set reduced 2-norm prediction error by 51% but did not improve impact shape orthogonality.

Annett, Martin S.↗
DOE OSTI2023

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗
NASA NTRS2017

A Dynamic Calibration Method for Experimental and Analytical Hub Load Comparison

This paper presents the results from an ongoing effort to produce improved correlation between analytical hub force and moment prediction and those measured during wind-tunnel testing on the Aeroelastic Rotor Experimental System (ARES), a conventional rotor testbed commonly used at the Langley Transonic Dynamics Tunnel (TDT). A frequency-dependent transformation between loads at the rotor hub and outputs of the testbed balance is produced from frequency response functions measured during vibration testing of the system. The resulting transformation is used as a dynamic calibration of the balance to transform hub loads predicted by comprehensive analysis into predicted balance outputs. In addition to detailing the transformation process, this paper also presents a set of wind-tunnel test cases, with comparisons between the measured balance outputs and transformed predictions from the comprehensive analysis code CAMRAD II. The modal response of the testbed is discussed and compared to a detailed finite-element model. Results reveal that the modal response of the testbed exhibits a number of characteristics that make accurate dynamic balance predictions challenging, even with the use of the balance transformation.

Kreshock, Andrew R.↗
NASA NTRS2019

TPSAS-NF1676L-32920-DND

In June 2019, a full-scale crash test of a Fokker F28 Fellowship aircraft will be conducted at NASA Langley Research Center?s Landing and Impact Research (LandIR) Facility. The F28 is a high-performance twin-turbo fan narrow-body aircraft with seating in a 3+2 configuration. The MK4000 variant, used in this test, is capable of carrying up to 85 passengers on medium range routes. The F28 was first type certified by the Federal Aviation Administration (FAA) in 1969 and now the majority of the F28 fleet has retired from service in the United States. In 2016, the FAA and NASA Langley Research Center (LaRC) signed an interagency agreement for conducting a research program to obtain test data that will support the development of airframe level crash requirements for transport category airplanes [1]. The objectives of the full-scale crash test can be divided into six categories: (1) To compare and contrast responses in identical aircraft undergoing vertical only to combined vertical and horizontal loading conditions, (2) To examine the effects of horizontal loading on aircraft structure during a crash event, (3) To generate data for the use in calibration of computer modelling efforts, (4) To generate data from onboard Anthropomorphic Test Devices (ATDs) for the evaluation of injury, (5) To obtain data from experimental seats, and (6) To obtain data from new and novel ATD designs including Warrior Injury Assessment MANikin (WIAMan) [2], Test device for Human Occupant Response (THOR), and other newly developed child ATDs. The focus of this presentation will be to document finite element model development of the full-scale F28 aircraft and to present preliminary test-analysis predictions. NASA obtained the full-scale F28 aircraft, plus three fuselage sections and two sets of wings with funding through the NASA Aviation Safety Program in 1998. In addition to the hardware, NASA purchased a full NASTRAN model of the airframe that had been developed by the Dutch manufacturer, Fokker. Beginning in 2016, the NASTRAN model was converted to LS-DYNA? [4, 5] format and modified by combining parts, adding missing parts of internal structure, and including ballast for loading weights, etc. It contains: 89,223 nodes; 24,065 beam elements; 55,404 shell elements; 29,044 solid elements; 740 parts; 80 material definitions; and, 46 Constrained Nodal Rigid Bodies (CNRBs). The aircraft model will be executed in LSDYNA to simulate the test article impact onto a 2-ft.-high bed of soil under combined velocity conditions of 70-ft/s forward and 30-ft/s vertical velocity. Pre-test simulation predictions will be generated and correlated with test data.

Karen E Jackson↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling the interfacial failure and resulting mechanical properties of z-pinned additively manufactured composites

The mechanical properties of Fused Filament Fabrication (FFF) parts are limited in the build direction, albeit advantages such as design flexibility and in-house customization. Here, in this paper, a novel z-pinning approach that enhances mechanical properties in the build direction by depositing material across multiple layers within the part was investigated through numerical simulations and validated experimentally. A finite element model for z-pinned composite structures was developed by assigning fiber-orientation dependent material properties obtained using a micromechanics approach to beads and pins. The bead-to-bead and pin-to-pin adhesions within the z-pinned structures were modeled using a cohesive traction separation law. The properties of cohesive elements for carbon fiber-reinforced polylactic acid (CF-PLA) z-pinned composites were calibrated using tensile experiments. The elastic modulus and tensile strength of CF-PLA z-pinned composites in the build direction were predicted with the developed numerical model. The numerical investigation on various geometrical parameters revealed that the largest pin volume increases the stiffness and tensile strength by 40% and thus, has the greatest influence on the mechanical properties of the z-pinned composites. The effect of z-pin geometrical parameters on the mechanical properties was summarized to aid in the design of z-pinned additively manufactured composite structures.

36 MATERIALS SCIENCE↗
NASA NTRS1990

The effect of friction in the hold down post spherical bearings on hold down post loads

The effect of friction at the connection of the Solid Rocket Booster (SRB) aft skirt and the mobile launch platform (MLP) hold down posts was analyzed. A simplified model of the shuttle response during the Space Shuttle Main Engine (SSME) buildup was constructed. The model included the effect of stick-slip friction for the rotation of the skirt about the spherical bearing. Current finite element models assume the joint is completely frictionless in rotation and therefore no moment is transferred between the skirt and the hold down posts. The model was partially verified against test data and preliminary parameter studies were performed. The parameter studies indicated that the coefficient of friction strongly influenced the moment on the hold down posts. The coefficient of friction had little effect on hold down post vertical loads, however. Further calibration of the model is necessary before the effect of friction on the hold down post horizontal loads can be analyzed.

Richardson, James A.↗
NASA NTRSDate not supplied

Finite Element Thermal Model for Ultrasonic Welding of Thermoplastic Composites

Ultrasonic welding, UW, is a fast and energy-efficient technique for joining thermoplastic composites. It involves the use of high-frequency mechanical vibrations and a static welding force to melt and join adherends. Ultrasonic welding is an enabling technology to reduce the cost and complexity of in-space construction because lightweight thermoplastic composite components can be packaged compactly for launch and then efficiently assembled using supervised autonomous robotic technologies on site. However, the temperatures in space present challenges to UW, and it is critical that the efficacy of process parameters selected for manufacturing in space is understood prior to launch. To this end, a three-dimensional finite element model is presented in this technical presentation. The model incorporates equations for effects of viscoelastic heating and heat transfer on the welding process. The proposed model is applied to predict the temperature distribution in single lap shear, SLS, samples composed of AS4/PEEK (TC1200) composite that were welded using a terrestrial machine as part of a comprehensive weldability study. Thermocouple and infrared spot sensor data from the SLS samples provide empirical temperature measurements for calibration and validation of the UW thermal model. Calibration and validation of the model is an important step given the significant uncertainties in material properties such as the loss modulus and necessary assumptions in the physics implementations which allow the model to converge in an acceptable amount of time. The validated thermal model can be used to simulate the process for the space environment. Future validation sample testing is planned in a vacuum chamber. The result will be a model capable of guiding process parameter selection to ensure acceptable weld bonds when manufacturing in space.

Josh Fody↗
NASA NTRS2022

Finite Element Modeling of Rigid-Flex PCB Robotic Structures for Dynamic Environments

Rigid-flex circuit boards are becoming moreprevalent as the limits are pushed on the size, mass, and geometryof electronic systems. A key aspect of designing a rigid-flexprinted circuit boards (PCB) system is an assessment of thedynamic properties of the PCB and predicting system perfor-mance under dynamic loading. Among current modeling method-ologies for rigid-flex PCB, a simplified modeling methodology thatadequately captures the system dynamics does not exist. This arti-cle presents a novel, computationally efficient approach formodeling rigid-flex PCB systems and the calibration of the mate-rial models via modal testing. The resulting simplified model isable to capture system frequencies, mode shapes, and representa-tive force-displacement behavior. The proposed methodology isused to model NASA Jet Propulsion Laboratory’s Pop-Up FlatFolding Explorer Robot (PUFFER) and assess the sensitivity of asystem model to input parameters.

John Bell↗
DOE OSTI2025

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗
DOE OSTI2024

An active learning framework for the rapid assessment of galvanic corrosion

The current present in a galvanic couple can define its resistance or susceptibility to corrosion. However, as the current is dependent upon environmental, material, and geometrical parameters it is experimentally costly to measure. To reduce these costs, Finite Element (FE) simulations can be used to assess the cathodic current but also require experimental inputs to define boundary conditions. Due to these challenges, it is crucial to accelerate predictions and accurately predict the current output for different environments and geometries representative of in-service conditions. Machine learned surrogate models provides a means to accelerate corrosion predictions. However, a one-time cost is incurred in procuring the simulation and experimental dataset necessary to calibrate the surrogate model. Therefore, an active learning protocol is developed through calibration of a low-cost surrogate model for the cathodic current of an exemplar galvanic couple (AA7075-SS304) as a function of environmental and geometric parameters. The surrogate model is calibrated on a dataset of FE simulations, and calculates an acquisition function that identifies specific additional inputs with the maximum potential to improve the current predictions. This is accomplished through a staggered workflow that not only improves and refines prediction, but identifies the points at which the most information is gained, thus enabling expansion to a larger parameter space. The protocols developed and demonstrated in this work provide a powerful tool for screening various forms of corrosion under in-service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Study of Aircraft Structural Response and Occupant Loading During a Ditching Event Utilizing LS-DYNA Simulation

In this study, the capability of finite element model (FEM) analysis to assess aircraft structural response and occupant injury risk during aircraft water ditching conditions was evaluated. A FEM of the Fokker F28 aircraft, previously validated against full-scale crash test data, was simulated in a realistic water impact environment utilizing recorded conditions from the ditching of US Airways Flight 1549 into the Hudson River. The FEM of the F28 aircraft, was developed from a NASTRAN loads model provided by the vehicle manufacturer, calibrated using component fuselage drop-tests, and validated through test-analysis comparisons from a full-scale crash test of the F28 onto soil. Simulation results are presented in several categories: vehicle kinematics, structural accelerations, damage, occupant responses, and injury risk predictions. Results were consistent with the general outcomes of Flight 1549, indicating that the developed vehicle FEM, validated against land impact test data, may provide an effective tool for evaluating water ditching events.

Jacob Putnam↗