Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element 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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Properties of Syntactic Foam for Simulation of Mechanical Insults.

Syntactic foam encapsulation protects sensitive components. The energy mitigated by the foam is calculated with numerical simulations. The properties of a syntactic foam consisting of a mixture of an epoxy-rubber adduct and glass microballoons are obtained from published literature and test results. The conditions and outcomes of the tests are discussed. The method for converting published properties and test results to input for finite element models is described. Simulations of the test conditions are performed to validate the inputs.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Numerical studies of nonlinear ultrasonic guided waves in uniform waveguides with arbitrary cross sections

Nonlinear guided waves have been investigated widely in simple geometries, such as plates, pipe and shells, where analytical solutions have been developed. This paper extends the application of nonlinear guided waves to waveguides with arbitrary cross sections. The criteria for the existence of nonlinear guided waves were summarized based on the finite deformation theory and nonlinear material properties. Numerical models were developed for the analysis of nonlinear guided waves in complex geometries, including nonlinear Semi-Analytical Finite Element (SAFE) method to identify internal resonant modes in complex waveguides, and Finite Element (FE) models to simulate the nonlinear wave propagation at resonant frequencies. Two examples, an aluminum plate and a steel rectangular bar, were studied using the proposed numerical model, demonstrating the existence of nonlinear guided waves in such structures and the energy transfer from primary to secondary modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis

Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Infrared imaging diagnostics for intense pulsed electron beam

Infrared imaging diagnostic method for two-dimensional calorimetric diagnostics has been developed for intense pulsed electron beam (IPEB). By using a 100-μm-thick tungsten film as the infrared heat sink for IPEB, the emitting uniformity of the electron source can be analyzed to evaluate the efficiency and stability of the diode system. Two-dimensional axisymmetric finite element method heat transfer simulation, combined with Monte Carlo calculation, was performed for error estimation and optimization of the method. The test of the method was finished with IPEB generated by explosive emission electron diode with pulse duration (FWHM) of 80 ns, electron energy up to 450 keV, and a total beam current of over 1 kA. The results showed that it is possible to measure the cross-sectional energy density distribution of IPEB with energy sensitivity of 0.1 J/cm{sup 2} and spatial resolution of 1 mm. The technical details, such as irradiation protection of bremsstrahlung γ photons and the functional extensibility of the method were discussed in this work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1993

A comparative study of an ABC and an artificial absorber for truncating finite element meshes

The type of mesh termination used in the context of finite element formulations plays a major role on the efficiency and accuracy of the field solution. The performance of an absorbing boundary condition (ABC) and an artificial absorber (a new concept) for terminating the finite element mesh was evaluated. This analysis is done in connection with the problem of scattering by a finite slot array in a thick ground plane. The two approximate mesh truncation schemes are compared with the exact finite element-boundary integral (FEM-BI) method in terms of accuracy and efficiency. It is demonstrated that both approximate truncation schemes yield reasonably accurate results even when the mesh is extended only 0.3 wavelengths away from the array aperture. However, the artificial absorber termination method leads to a substantially more efficient solution. Moreover, it is shown that the FEM-BI method remains quite competitive with the FEM-artificial absorber method when the FFT is used for computing the matrix-vector products in the iterative solution algorithm. These conclusions are indeed surprising and of major importance in electromagnetic simulations based on the finite element method.

Oezdemir, T.↗
NASA NTRS2013

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi-Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post-landing segment of the Flight Software (FSW). An in-depth trade study was completed to determine optimal characteristics of the touchdown detection method resulting in an algorithm monitoring filtered, lever-arm corrected, 200 Hz Inertial Measurement Unit (IMU) vehicle acceleration magnitude data against a tunable threshold using persistence counter logic. Following the design of the algorithm, high fidelity environment and vehicle simulations, coupled with the actual vehicle FSW, were used to tune the acceleration threshold and persistence counter value to result in adequate performance in detecting touchdown and sufficient safety margin against early detection while descending under parachutes. An analytical approach including Kriging and adaptive sampling allowed for a sufficient number of finite element analysis (FEA) impact simulations to be completed using minimal computation time. The combination of a persistence counter of 10 and an acceleration threshold of approximately 57.3 ft/s2 resulted in an impact performance factor of safety (FOS) of 1.0 and a safety FOS of approximately 2.6 for touchdown declaration. An RCS termination acceleration threshold of approximately 53.1 ft/s(exp)2 with a persistence counter of 10 resulted in an increased impact performance FOS of 1.2 at the expense of a lowered under-parachutes safety factor of 2.2. The resulting tuned algorithm was then tested on data from eight Capsule Parachute Assembly System (CPAS) flight tests, showing an experimental minimum safety FOS of 6.1. The formulated touchdown detection algorithm will be flown on the Orion MPCV FSW during the Exploration Flight Test 1 (EFT-1) mission in the second half of 2014.

Daum, Jared↗
NASA NTRS1988

An elastic-plastic finite-element analysis of the J-resistance curve using a CTOD criterion

The fracture process in compact and bend specimens was simulated using a two-dimensional finite-element analysis of the J-resistance (J-R) curve and a CTOD criterion. The J-R curves were calculated from the numerical results for each specimen type using several different methods. In general, the J-R curves obtained for the bend specimens were found to be higher than those for the compact specimens, especially beyond maximum load. However, below the maximum load, the modified deformation theory of plasticity and the contour-integral J(r) method resulted in very similar J-R curves for both specimen types.

Newman, J. C., Jr.↗
DOE OSTI2018

In-plane crashworthiness of bio-inspired hierarchical honeycombs

Biological tissues like bone, wood, and sponge possess hierarchical cellular topologies, which are lightweight and feature an excellent energy absorption capability. Here we present a system of bio-inspired hierarchical honeycomb structures based on hexagonal, Kagome, and triangular tessellations. The hierarchical designs and a reference regular honeycomb configuration are subjected to simulated in-plane impact using the nonlinear finite element code LS-DYNA. The numerical simulation results show that the triangular hierarchical honeycomb provides the best performance compared to the other two hierarchical honeycombs, and features more than twice the energy absorbed by the regular honeycomb under similar loading conditions. We also propose a parametric study correlating the microstructure parameters (hierarchical length ratio r and the number of sub cells N) to the energy absorption capacity of these hierarchical honeycombs. The triangular hierarchical honeycomb with N = 2 and r = 1/8 shows the highest energy absorption capacity among all the investigated cases, and this configuration could be employed as a benchmark for the design of future safety protective systems.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Acoustically Tailored Composite Rotorcraft Fuselage Panels

A rotorcraft roof sandwich panel has been redesigned to optimize sound power transmission loss (TL) and minimize structure-borne sound for frequencies between 1 and 4 kHz where gear meshing noise from the transmission has the most impact on speech intelligibility. The roof section, framed by a grid of ribs, was originally constructed of a single honeycomb core/composite face sheet panel. The original panel has coincidence frequencies near 700 Hz, leading to poor TL across the frequency range of 1 to 4 kHz. To quiet the panel, the cross section was split into two thinner sandwich subpanels separated by an air gap. The air gap was sized to target the fundamental mass-spring-mass resonance of the double panel system to less than 500 Hz. The panels were designed to withstand structural loading from normal rotorcraft operation, as well as 'man-on-the-roof' static loads experienced during maintenance operations. Thin layers of VHB 9469 viscoelastomer from 3M were also included in the face sheet ply layups, increasing panel damping loss factors from about 0.01 to 0.05. Measurements in the NASA SALT facility show the optimized panel provides 6-11 dB of acoustic transmission loss improvement, and 6-15 dB of structure-borne sound reduction at critical rotorcraft transmission tonal frequencies. Analytic panel TL theory simulates the measured performance quite well. Detailed finite element/boundary element modeling of the baseline panel simulates TL slightly more accurately, and also simulates structure-borne sound well.

Hambric, Stephen↗
NASA NTRSDate not supplied

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of the limitations of the technique can be obtained by simulating the inspection process. Time domain finite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

thermography↗
NASA NTRSDate not supplied

Realistic Thermographic Simulation of Impact Damage with Quadrupole Method

Flash thermography has been shown to be an effective method for detection of delaminations in carbon fiber reinforced polymer (CFRP) composite structures. Improved understanding of thelimitations of the technique can be obtained by simulating the inspection process. Time domainfinite difference and finite element methods are well suited for such simulations but can be computationally intensive. An alternate method is to solve the Laplace transform of the heat equation, then invert the Laplace transform to produce a time domain response. Often this is referred to as the thermal quadrupole method and is well suited for simulating flash thermography. The quadrupole method has been used extensively for simulating one-dimensional heat flow in multilayer systems with contact resistances at the interfaces. It is also applicable for three dimensional configurations, in particular for simulation of changes in thermographic response due to delaminations in composites. This presentation discusses three-dimensional simulations with realistic shapes and contact resistances at the interfaces between layers. The results of these simulations are compared to thermography measurements on composites specimens with impact induced delaminations.

Thermography↗
DOE OSTI2020

Factors and Considerations for Modeling Loss of a GaN-based Inverter

The article investigates the impacts of four often-neglected factors on the loss model of a GaN-based full-bridge inverter: parasitic capacitance of the devices, dynamics of junction temperature (Tj) under time-varying power dissipation (Ploss), case temperature estimation, and detailed considerations of the passive components. Procedures to calculate the converter loss considering the above factors are proposed and implemented. A 4.5-kW hard-switching inverter prototype using gallium nitrite (GaN) high-electron-mobility transistors is used to experimentally demonstrate the impact of each factor on the converter loss model. Furthermore, it is found that the accuracy of a converter loss model is mainly affected by the passive components at the light load condition, whereas the thermal and loss models of the active components become the major factors as the output power increases. The results show that after considering the above factors, the converter loss discrepancy between calculation and measurement can be reduced from 30.6 W (28%) to 2.5 W (less than 3%) at heavy load (Po = 4.5 kW), while at the light load condition (Po = 500 W), it is reduced from 3.9 W (28%) to 2.6 W (16%). Furthermore, the difference between simulated and measured case temperature of the GaN devices is within 6 °C.

42 ENGINEERING↗
NASA NTRS1990

Computational technology for flight vehicles; Proceedings of the Symposium, Washington, DC, Nov. 5-7, 1990

Recent advances in computational fluid mechanics are discussed in reviews and reports. Sections are devoted to (1) the modeling of local phenomena and edge effects in solids, (2) stochastic modeling and simulation of fracture toughness, and (3) partitioning strategy and new finite elements. Particular attention is given to global and local finite-element/spectral-boundary-element techniques for failure analysis, simulations of microfracture in metal-matrix composites, fatigue analysis of cracked anisotropic plates under stochastic loading, mathematical modeling for the analysis of nonlinear aircraft dynamics, physical and mathematical modeling of wave propagation in the Ariane 5 VEB structure, partitioning based on symmetry transformations, an FEM approach to adaptive reliability assurance, and time-domain FEMs for the large rotational dynamics of multibody systems.

Noor, Ahmed K.↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting in real time the residual strength of flight structures with discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. A residual strength test of a metallic, integrally-stiffened panel is simulated to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data would, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high-fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
NASA NTRS2004

Test-Analysis Correlation of a Crash Simulation of a Vertical Drop Test of a Commuter-Category Aircraft

A finite element model of an ATR42-300 commuter-class aircraft was developed and a crash simulation was executed. Analytical predictions were correlated with data obtained from a 30-feet per second (9.14-meters per second) vertical drop test of the aircraft. The purpose of the test was to evaluate the structural response of the aircraft when subjected to a severe, but survivable, impact. The aircraft was configured with seats, dummies, luggage, and other ballast. The wings were filled with 8,700 lb. (3,946 kilograms) of water to represent the fuel. The finite element model, which consisted of 57,643 nodes and 62,979 elements, was developed from direct measurements of the airframe geometry. The seats, dummies, luggage, simulated engines and fuel, and other ballast were represented using concentrated masses. The model was executed in LS-DYNA, a commercial finite element code for performing explicit transient dynamic simulations. Analytical predictions of structural deformation and selected time-history responses were correlated with experimental data from the drop test to validate the simulation.

Jackson, Karen E.↗
DOE OSTI2017

Characterization of Flame Cut Heavy Steel: Modeling of Temperature History and Residual Stress Formation

Heavy steel plates are used in demanding applications that require both high strength and hardness. An important step in the production of such components is cutting the plates with a cost-effective thermal cutting method such as flame cutting. Flame cutting is performed with a controlled flame and oxygen jet, which burns the steel and forms a cutting edge. However, the thermal cutting of heavy steel plates causes several problems. A heat-affected zone (HAZ) is generated at the cut edge due to the steep temperature gradient. Consequently, volume changes, hardness variations, and microstructural changes occur in the HAZ. In addition, residual stresses are formed at the cut edge during the process. In the worst case, unsuitable flame cutting practices generate cracks at the cut edge. The flame cutting of thick steel plate was modeled using the commercial finite element software ABAQUS. The results of modeling were verified by X-ray diffraction-based residual stress measurements and microstructural analysis. The model provides several outcomes, such as obtaining more information related to the formation of residual stresses and the temperature history during the flame cutting process. In addition, an extensive series of flame cut samples was designed with the assistance of the model.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Assessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module

Orion Crew Module (CM) touchdown detection is critical to activating the post-landing sequence that safe?s the Reaction Control Jets (RCS), ensures that the vehicle remains upright, and establishes communication with recovery forces. In order to accommodate safe landing of an unmanned vehicle or incapacitated crew, an onboard automated detection system is required. An Orion-specific touchdown detection algorithm was developed and evaluated to differentiate landing events from in-flight events. The proposed method will be used to initiate post-landing cutting of the parachute riser lines, to prevent CM rollover, and to terminate RCS jet firing prior to submersion. The RCS jets continue to fire until touchdown to maintain proper CM orientation with respect to the flight path and to limit impact loads, but have potentially hazardous consequences if submerged while firing. The time available after impact to cut risers and initiate the CM Up-righting System (CMUS) is measured in minutes, whereas the time from touchdown to RCS jet submersion is a function of descent velocity, sea state conditions, and is often less than one second. Evaluation of the detection algorithms was performed for in-flight events (e.g. descent under chutes) using hi-fidelity rigid body analyses in the Decelerator Systems Simulation (DSS), whereas water impacts were simulated using a rigid finite element model of the Orion CM in LS-DYNA. Two touchdown detection algorithms were evaluated with various thresholds: Acceleration magnitude spike detection, and Accumulated velocity changed (over a given time window) spike detection. Data for both detection methods is acquired from an onboard Inertial Measurement Unit (IMU) sensor. The detection algorithms were tested with analytically generated in-flight and landing IMU data simulations. The acceleration spike detection proved to be faster while maintaining desired safety margin. Time to RCS jet submersion was predicted analytically across a series of simulated Orion landing conditions. This paper details the touchdown detection method chosen and the analysis used to support the decision.

Gay, Robert S.↗