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 73 records · Page 4

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗

Comparison of Test Methods to Determine Failure Parameters for MAT162 Calibration

MAT162, a laminated composite failure material model developed by The Material Sciences Corporation for the commercial finite element software LS-Dyna, is widely used within the aerospace industry to predict damage events under a range of dynamic conditions. The material model involves numerous inputs consisting of both physical material properties and numerical calibration parameters. Due to the large number material card inputs, often there is a lack of uniqueness to MAT162 material cards that limits the predictive capability to only the directly calibrated space. To expand this space, MAT162 requires a prudent and robust calibration process in which significant parameters are calibrated to high confidence damage events observed in experiment. Critical to this success is fully defining the material properties correctly, namely the fiber crush (SFC) and fiber shear (SFS) values, prior to calibrating the numerical parameters. In this paper, the effect of the determination of SFS and SFC on subsequent calibration steps is examined using two different experimental techniques.

Molitor, Matthew↗

Prediction and Verification of Ductile Crack Growth from Simulated Defects in Strength Overmatched Butt Welds

Defects that develop in welds during the fabrication process are frequently manifested as embedded flaws from lack of fusion or lack of penetration. Fracture analyses of welded structures must be able to assess the effect of such defects on the structural integrity of weldments; however, the transferability of R-curves measured in laboratory specimens to defective structural welds has not been fully examined. In the current study, the fracture behavior of an overmatched butt weld containing a simulated buried, lack-of-penetration defect is studied. A specimen designed to simulate pressure vessel butt welds is considered; namely, a center crack panel specimen, of 1.25 inch by 1.25 inch cross section, loaded in tension. The stress-relieved double-V weld has a yield strength 50% higher than that of the plate material, and displays upper shelf fracture behavior at room temperature. Specimens are precracked, loaded monotonically while load-CMOD measurements are made, then stopped and heat tinted to mark the extent of ductile crack growth. These measurements are compared to predictions made using finite element analysis of the specimens using the fracture mechanics code Warp3D, which models void growth using the Gurson-Tvergaard dilitant plasticity formulation within fixed sized computational cells ahead of the crack front. Calibrating data for the finite element analyses, namely cell size and initial material porosities are obtained by matching computational predictions to experimental results from tests of welded compact tension specimens. The R-curves measured in compact tension specimens are compared to those obtained from multi-specimen weld tests, and conclusions as to the transferability of R-curves is discussed.

Nishioka, Owen S.↗

Single Degree-of-Freedom Modeling of SLS Liquid Hydrogen Pre-Valve Flow Guide to Enable Rapid Transient Analysis

A unique single degree-of-freedom approximation technique has been developed to enable rapid application of a temporally-defined multi-spectral semi-narrow-band loading for generation of realistic stress/cycle values compared to a resonant analysis. The technique uses the harmonic analysis at resonance of a high-fidelity finite element model to produce a transfer function, which is then used to calibrate the response of the SDOF model. A standard numerical ordinary differential equation solver is then used to obtain the temporal response, and its histogram is used in a fatigue/fracture model. This technique is related to other SDOF methods used widely in industry, such as Miles' Equation and the Shock Response Spectra, but it is unique in that it produces a realistic time history of the response. The most obvious error in the process, which is the effect of closely-spaced modes, was also assessed using the parallel application of several SDOF models, and the error is shown to be small. The application of this unique and tractable reduced-order methodology has enabled the SLS program to avoid substantial cost and schedule penalties if a redesign or change of material were required. It has also enabled quick analysis of a number of other structures undergoing the same or similar excitation fields, and quick assessment when the excitation and structural configuration has been altered due to design changes in the system.

Brown, Andrew M.↗

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 6061-T651 Aluminum Alloy from Plate

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

36 MATERIALS SCIENCE↗

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗

Calibration of International Space Station (ISS) Node 1 Vibro-Acoustic Model

Reported here is the ability of utilizing the Energy Finite Element Method (E-FEM) to predict the vibro-acoustic sound fields within the International Space Station (ISS) Node 1 and to compare the results with actual measurements of leak sounds made by a one atmosphere to vacuum leak through a small hole in the pressure wall of the Node 1 STA module during its period of storage at Stennis Space Center (SSC). While the E-FEM method represents a reverberant sound field calculation, of importance to this application is the requirement to also handle the direct field effect of the sound generation. It was also important to be able to compute the sound fields in the ultrasonic frequency range. This report demonstrates the capability of this technology as applied to this type of application.

Zhang, Weiguo↗

Shake Test Results and Dynamic Calibration Efforts for the Large Rotor Test Apparatus

Prior to the full-scale wind tunnel test of the UH-60A Airloads rotor, a shake test was completed on the Large Rotor Test Apparatus. The goal of the shake test was to characterize the oscillatory response of the test rig and provide a dynamic calibration of the balance to accurately measure vibratory hub loads. This paper provides a summary of the shake test results, including balance, shaft bending gauge, and accelerometer measurements. Sensitivity to hub mass and angle of attack were investigated during the shake test. Hub mass was found to have an important impact on the vibratory forces and moments measured at the balance, especially near the UH-60A 4/rev frequency. Comparisons were made between the accelerometer data and an existing finite-element model, showing agreement on mode shapes, but not on natural frequencies. Finally, the results of a simple dynamic calibration are presented, showing the effects of changes in hub mass. The results show that the shake test data can be used to correct in-plane loads measurements up to 10 Hz and normal loads up to 30 Hz.

Large Rotor Test Apparatus↗

Full-Scale Crash Tests and Analyses of Three High-Wing Single

The NASA Emergency Locator Transmitter Survivability and Reliability (ELTSAR) project was initiated in 2014 to assess the crash performance standards for the next generation of ELT systems. Three Cessna 172 aircraft have been acquired to conduct crash testing at NASA Langley Research Center's Landing and Impact Research Facility. Testing is scheduled for the summer of 2015 and will simulate three crash conditions; a flare to stall while emergency landing, and two controlled flight into terrain scenarios. Instrumentation and video coverage, both onboard and external, will also provide valuable data of airframe response. Full-scale finite element analyses will be performed using two separate commercial explicit solvers. Calibration and validation of the models will be based on the airframe response under these varying crash conditions.

Annett, Martin S.↗

Advancements in Constitutive Model Calibration: Leveraging the Power of Full‐Field DIC Measurements and In Situ Load Path Selection for Reliable Parameter Inference

Accurate material characterization and model calibration are essential for computationally supported high-consequence engineering decisions. Historically, characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data are collected for a specific model of interest, (3) use deterministic methods that provide best-fit parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work brings together several recent advancements into an improved workflow called interlaced characterization and calibration (ICC) that advances the state-of-the-art in constitutive model calibration. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) quantifies parameter uncertainty through Bayesian inference and (4) incorporates these advancements into a quasi real-time feedback loop. The ICC framework is demonstrated here on the calibration of a material model using simulated full-field data for an aluminium cruciform specimen being deformed biaxially. The cruciform is actively driven through the myopically preferred load path using Bayesian optimal experimental design, which selects load steps that yield the maximum expected information gain (EIG). Principal component analysis (PCA) is performed on the model predictions of full-field displacements, and fast surrogate models are built to approximate the input-output relationships of the expensive finite element model. Furthermore, the tools developed and demonstrated here show that high-fidelity constitutive models can be efficiently and reliably calibrated with quantified uncertainty, thus supporting credible decision-making and potentially increasing the agility of solid mechanics modelling by enabling utilization of computational simulations at earlier stages of the design cycle.

Bayesian optimal experimental design↗

Calibration of International Space Station (ISS) Node 1 Vibro-Acoustic Model-Report 2

Reported here is the capability of the Energy Finite Element Method (E-FEM) to predict the vibro-acoustic sound fields within the International Space Station (ISS) Node 1 and to compare the results with simulated leak sounds. A series of electronically generated structural ultrasonic noise sources were created in the pressure wall to emulate leak signals at different locations of the Node 1 STA module during its period of storage at Stennis Space Center (SSC). The exact sound source profiles created within the pressure wall at the source were unknown, but were estimated from the closest sensor measurement. The E-FEM method represents a reverberant sound field calculation, and of importance to this application is the requirement to correctly handle the direct field effect of the sound generation. It was also important to be able to compute the sound energy fields in the ultrasonic frequency range. This report demonstrates the capability of this technology as applied to this type of application.

Zhang, Weiguo↗

Calibration of Airframe and Occupant Models for Two Full-Scale Rotorcraft Crash Tests

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. Accelerations and kinematic data collected from the crash tests were compared to a system integrated finite element model of the test article. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the second full-scale crash test. This combination of heuristic and quantitative methods was used to identify modeling deficiencies, evaluate parameter importance, and propose required model changes. It is shown that the multi-dimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and co-pilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. This approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and test planning guidance. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, thereby reducing overall development costs.

Annett, Martin S.↗

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗

Modeling the viscoplastic behavior of a semicrystalline polymer

In this study, a complex constitutive relation is identified using inverse modeling with the nominal mechanical response as sole experimental input. The methodology is illustrated for a semicrystalline thermoplastic in the presence of strain localization at finite deformations. The experimental database includes cylindrical tensile bars, compression pins and round notched bars loaded at strain rates spanning up to five decades and temperatures below and above T g . The data is organized into a calibration set and a validation set. The response of tensile specimens is determined using finite element analyses and a two-phase constitutive relation for semicrystalline polymers that accounts for temperature- and rate-sensitive plastic flow, pressure-sensitivity, small-strain softening and large-strain orientational hardening of the amorphous phase, along with the evolution of crystallinity. The large number of constitutive parameters is identified using an optimization tool coupled with the finite element solver and the calibration set from experiments. The methodology is shown to be successful in predicting the response of round notched bars and replicating the effects of temperature and strain rate on the severity of necking in tensile bars. The proposed model identification strategy is both simple and effective in comparison with other elaborate methods that attempt to access intrinsic behavior directly from high-fidelity experimental measurements.

36 MATERIALS SCIENCE↗

Report of METL Mutual Inductance Level Sensor Development for Use in Liquid Metals – FY2023

A robust mutual inductance level sensor (MILS) for use in high temperature liquid metals has been developed at Argonne National Laboratory (ANL). The current design utilizes mineral insulated cables with a 300-series stainless steel sheath, magnesium oxide insulation, and a single copper conductor. Two coils are wrapped on a common core made of a 300-series stainless steel tube. One coil is energized with an alternating current (AC) source, and this electromagnetically couples with the second coil to generate an induced voltage. The voltage is measured to determine the mutual induction between the coils, and the mutual inductance measurement can be used to determine the level of a nearby liquid metal volume. The MILS can be located in an isolating thimble that is fully sealed to a vessel containing liquid metal for ease of maintenance and replacement. When the primary coil is energized with the AC source, the coil not only electromagnetic (EM) couples with the second coil, but also the surrounding liquid metal. The EM coupling with the surrounding liquid metal reduces the EM coupling with the second coil, producing an inverse-linear relationship between liquid metal level and secondary coil voltage. This has all been demonstrated at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at ANL. The most current iteration of the MILS system is the MILS-MK-II. This sensor system has been commissioned in a non-sodium test stand where calibrations were performed using a sodium analog. The calibrations proved to be highly linear and repeatable. The MILS-MK-II has been installed in the METL expansion tank where high temperature sodium tests have been performed. The MILS-MK-II has been calibrated against a known standard at temperature of 300°C, and the calibrations have proved to be highly linear and repeatable. The calibrated MILS-MK-II has been in operation in the METL facility for several thousands of hours at temperatures around 300°C. The experimental data has been used to develop and validate electromagnetic finite element models in COMSOL Multiphysics, and now these models can be used to advance the development of the sensor system. Next steps will include temperature compensation to allow for operation at various temperatures. Additionally, efforts to incorporate a self-calibration methodology are underway.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗