Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element analysis”

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 325 records · Page 18

Hoop Tensile Failure Stress of 500 µm-Diameter Ceramic Spheres

The tensile failure stress of a brittle-material sphere is important to the survival of the sphere itself but is particularly relevant to systems that employ such spheres for applications in which the sphere exterior is subjected to tension during system fabrication or operation. That maximum hoop or tangential stress can be measured by diametrally compressing a C-sphere specimen directly machined from such spheres. This combination of specimen geometry and test method was conceived at Oak Ridge National Laboratory around 20 years ago. Although the C-sphere geometry has been tested with larger spheres, the ease and efficacy of its testing are unknown for submillimeter sphere diameters. To consider this, C-sphere specimens were machined from yttria-stabilized zirconia (YSZ) spheres of ∼500 µm diameter. The maximum (hoop or tangential) tensile stress was determined as a function of compressive failure force and specimen dimensions using finite element analysis. Companion effective area and volume as a function of Weibull modulus were computed. Hoop tensile failure stresses were then measured and fractography conducted to examine flaw type and location and to estimate the mirror constant of this YSZ grade. This novel study demonstrates that hoop tensile failure stress can be determined with sub-millimeter brittle-material spheres.

Delia, Daniel [ORNL] (ORCID:0000000209023239)↗

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗

An Analytical Solution for the Initiation and Early Progression of Fretting Wear in Spherical Contacts

Abstract This article derives analytical solutions to calculate the wear volume at the initiation of fretting motion and its early progression over the first few oscillation cycles. The Archard-based model considers a deformable hemisphere that is contact with a deformable flat block. The material pairs investigated are special alloys, the Inconel 617/Incoloy 800H, and Inconel 617/Inconel 617. The analytical study begins with a unidirectional frictional sliding contact, where the local interfacial sliding distance and the nominal sliding distance at the initiation of gross slip are derived. The obtained analytical expressions for unidirectional sliding are then used to derive the corresponding wear volume for the initiation and early progression of gross slip and the wear volume for a general fretting cycle under elastic conditions. These analytical derivations are all verified by the finite element analysis (FEA). The FEA method and the analytical solutions render virtually identical results for both similar and dissimilar material pairs. The effects of plasticity on the wear volume under elastic–plastic conditions are also investigated. It is found that the fretting wear volumes obtained from the FEA simulations, which include plasticity, are close to those obtained from the analytical expressions for purely elastic regimes. All the results are presented in normalized forms, which can be easily generalized and applied to three-dimensional fretting wear of other material pairs.

Engineering↗

Validation and Parametric Investigations of an Internal Permanent Magnet Motor Using a Lumped Parameter Thermal Model

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multiphysics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a lumped parameter thermal network (LPTN) model of an Internal Permanent Magnet synchronous motor within the commercially available MOTOR-CAD modeling environment. The validation is based on temperature comparison with experimental data and with more detailed finite element analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations, and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot-winding and the end-winding - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot-winding and end-winding. As for the LPTN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance.

42 ENGINEERING↗

Design of an Additively Manufactured Recuperator With 800 °C Inlet Temperature for sCO2 Power Cycle Application

Abstract Recuperators with design temperatures at and above 800 °C can further increase the thermal efficiency of supercritical CO2 power cycles by enabling higher turbine exhaust temperatures. Mar-M247 is a well-suited nickel-based superalloy for high temperature service due to its high creep strength that prevents excessive material thickness being required for pressure containment. Additive manufacturing using a high-speed laser-directed energy deposition (L-DED) process presents a promising solution, with build trials demonstrating the ability to produce nonconventional flow channels for enhanced heat transfer. A design process is presented that includes aerothermal and mechanical evaluation to maximize performance within the constraints of the manufacturing process. A 2-D heat transfer network and pressure drop code allows prediction of flow distribution and its effect on overall thermal performance. Established literature correlations, along with CFD simulation, inform the prediction of heat transfer coefficients and friction factors for the flowpaths and enhancement features in the heat exchanger core. Mechanical evaluation using finite element analysis (FEA) modeling with the intent of the ASME Boiler & Pressure Vessel Code (BPVC) Section VIII, Div. 2 assesses the operational safety of the design. The detailed design features annular finned passages that take advantage of helical flow paths to distribute the flow from separated headers to shared heat transfer surfaces. Performance predictions for the recuperator at a 50 kW scale provide insights into the feasibility of the additively manufactured (AM) process to produce recuperators on a commercial scale that extend existing operating envelopes.

Engineering↗

Experimental Verification of Pulse Shaping in Elastic Metamaterials Under Impact Excitation

Here, we present experimental verification of pulse shaping in elastic metamaterials together with a procedure to design, fabricate, and verify metamaterial pulse shapers under impact excitation. The split Hopkinson pressure bar (SHPB) test, a fundamental dynamic test introduced more than 70 years ago, often incorporates pulse shaping as a means to alter a stress wave, providing the primary motivation for the presented study. Elastic metamaterials hold promise for enhancing conventional pulse shaping abilities and improving capabilities of the SHPB test. We first design the pulse shaper by numerically optimizing its response using finite element analysis. The pulse shaper consists of repeated unit cells based on a combination of a phononic crystal and a local resonator. Then, we fabricate and test pulse shaper candidates to validate the procedural efficacy. An iterative element corrects inaccuracies in input force and material properties and allows convergence on an appropriate pulse shaper. We carry out this procedure by designing pulse shapers fabricated from 3D-printed polylactic acid (PLA) to achieve an extended dwell acceleration pulse shape. In experimental impact tests, the procedure results in rise, dwell, and fall behaviors comparable to that predicted, effectively confirming the efficacy of the presented procedure and verifying the performance of metamaterial-based pulse shapers.

36 MATERIALS SCIENCE↗

Grayscale Digital Light Processing Gradient Printing for Stress Concentration Reduction and Material Toughness Enhancement

Avoiding stress concentrations is essential to achieve robust parts since failure tends to originate at such concentrations. With recent advances in multimaterial additive manufacturing, it is possible to alter the stress (or strain) distribution by adjusting the material properties in selected locations. Here, we investigate the use of grayscale digital light processing (g-DLP) 3D printing to create modulus gradients around areas of high stress. These gradients prevent failure by redistributing high stresses (or strains) to the neighboring material. The improved material distributions are calculated using finite element analysis. The much-enhanced properties are demonstrated experimentally for thin plates with circular, triangular, and elliptical holes. In conclusion, this work suggests that multimaterial additive manufacturing techniques like g-DLP printing provide a unique opportunity to create tougher engineering materials and parts.

36 MATERIALS SCIENCE↗

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗

Validation and Parametric Investigations Using a Lumped Thermal Parameter Model of an Internal Permanent Magnet Motor

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multi-physics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a Lumped Parameter Thermal Network (LPTN) model of an Internal Permanent Magnet (IPM) synchronous motor within the commercially available Motor-CAD® modeling environment. The IPM motor’s stator is studied at steady state, and winding losses are generated by a constant DC current. The validation is based on temperature comparison with experimental data and with more detailed Finite Element Analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% — for the highest temperature in the slot windings and the end windings — was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot windings and end windings. As for the LTPN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance. The latter is often ignored in the literature, whereas its impact on temperature rise was found to be more significant than any other contact resistance within the stator.

47 OTHER INSTRUMENTATION↗

Prediction of Thermal Conditions of DED With FEA Metal Additive Simulation

This paper presents the integration of wire-arc additive manufacturing (WAAM) using Gas Metal Arc Welding (GMAW) into a machine tool to create a retrofit hybrid computer numeric control (CNC) machine tool. GMAW, along with other direct energy deposition systems, has the capacity to deposit material faster than the excess thermal energy can dissipate. This results in the need to allow the part to cool between consecutive layers, which is the most time-consuming part of the additive process. Finite element analysis (FEA) was used in conjunction with monitored build plate surface temperatures during deposition samples to improve adequate dwell time prediction and to develop a cooling system. A deposition was completed where no dwell time was used and the build plate along with the machine table temperatures were monitored. A second deposition was completed where only one bead was deposited and the traverse speed was increased. The GMAW welder was mounted on a 3-axis CNC machine where two square deposition samples were completed. A FEA model was designed and verified using the monitored samples. The model will be used to determine improved depositions speeds and whether forced cooling would allow for an increased deposition rate without structural failure. It was determined the FEA software can be used to accurately model and predict the thermal response of WAAM AM components.

Heinrich, Lauren↗

Sequential infiltration of two-photon polymerized 3D photonic crystals for mid-IR spectroscopic applications

Photonic crystals (PhCs) are spatially organized structures with lattice parameters equivalent to the operational wavelength of light. PhCs have been subject to extensive research efforts in the last two decades and are known for controlling light propagation with applications in sensing and time-delayed communication due to the slow-light phenomenon. Despite their exceptional properties, PhCs are difficult to fabricate using planar micromachining techniques due to their periodic structures. Techniques like two-photon stereolithography have been discussed for PhC fabrication in the literature, but the inherent disadvantage of poor refractive index (RI) contrast results in limited application. In this work, we present sequential infiltration synthesis performed on two-photon stereolithographically printed 3D PhCs for infiltration with zinc oxide to increase the RI of 3D PhCs. Further, finite element analysis was performed over a range of RI contrast values to study the change in photonic bandgap (PBG) with RI contrast. The transmission spectra were recorded on 3D PhCs before and after infiltration to demonstrate the change experimentally. An increase in the PBG width and absorbance is seen postinfiltration due to enhanced RI. This work presents the first, to our knowledge, sequentially infiltrated enhanced 3D PhC fabricated with two-photon stereolithography.

3D printing↗

Ultrastrong colloidal crystal metamaterials engineered with DNA

Lattice-based constructs, often made by additive manufacturing, are attractive for many applications. Typically, such constructs are made from microscale or larger elements; however, smaller nanoscale components can lead to more unusual properties, including greater strength, lighter weight, and unprecedented resiliencies. Here, solid and hollow nanoparticles (nanoframes and nanocages; frame size: ~15 nanometers) were assembled into colloidal crystals using DNA, and their mechanical strengths were studied. Nanosolid, nanocage, and nanoframe lattices with identical crystal symmetries exhibit markedly different specific stiffnesses and strengths. Unexpectedly, the nanoframe lattice is approximately six times stronger than the nanosolid lattice. Nanomechanical experiments, electron microscopy, and finite element analysis show that this property results from the buckling, densification, and size-dependent strain hardening of nanoframe lattices. Last, these unusual open architectures show that lattices with structural elements as small as 15 nanometers can retain a high degree of strength, and as such, they represent target components for making and exploring a variety of miniaturized devices.

42 ENGINEERING↗

Three-dimensional nucleation and growth of deformation twins in magnesium

At two-thirds the weight of aluminum, magnesium alloys have the potential to reduce the fuel consumption of transportation vehicles. These advancements depend on our ability to optimize the desirable versus undesirable effects of deformation twins, which are three-dimensional (3D) microstructural domains that form under mechanical stresses. Previously only characterized through surface or thin-film measurements, we here present 3D in situ characterization of deformation twinning inside an embedded grain over mesoscopic fields of view using dark-field x-ray microscopy supported by crystal plasticity finite element analysis. The results revealed the role of triple junctions on twin nucleation and the sequence and irregularity of twin growth and showed that twin-grain junctions, twin-twin junctions, and twin boundaries were the sites of localized dislocation accumulation.

Lee, Sangwon [Univ. of Michigan, Ann Arbor, MI (Un↗

Tetrahedral Trees: A Family of Hierarchical Spatial Indexes for Tetrahedral Meshes

In this work, we address the problem of performing efficient spatial and topological queries on large tetrahedral meshes with arbitrary topology and complex boundaries. Such meshes arise in several application domains, such as 3D Geographic Information Systems (GISs), scientific visualization, and finite element analysis. To this aim, we propose Tetrahedral trees, a family of spatial indexes based on a nested space subdivision (an octree or a kD-tree) and defined by several different subdivision criteria. We provide efficient algorithms for spatial and topological queries on Tetrahedral trees and compare to state-of-the-art approaches. Our results indicate that Tetrahedral trees are an improvement over R*-trees for querying tetrahedral meshes; they are more compact, faster in many queries, and stable at variations of construction thresholds. They also support spatial queries on more general domains than topological data structures, which explicitly encode adjacency information for efficient navigation but have difficulties with domains with a non-trivial geometric or topological shape.

97 MATHEMATICS AND COMPUTING↗

Precautions of Using Three-Electrode Configuration to Measure Electrode Overpotential in Solid Oxide Electrochemical Cells: Insights from Finite Element Modeling

Accurate determination of electrode overpotentials is essential to assess the performance of the electrode and understand the rate-limiting steps involved. Three-electrode configuration with the use of a reference electrode at a fixed potential is a standard way to measure overpotential of a specific electrode in liquid electrochemical systems. However, application of such three-electrode configuration to solid electrochemical cells for overpotential determination is not straightforward and requires extra caution. Here we report a theoretical Finite Element Analysis on the geometrical requirements for which the reference electrode can be applied to anode- or electrolyte-supported solid oxide button cells. The modeling results suggest that the symmetry of the working and counter electrodes is the key factor determining if a reference electrode is suitable to use. For anode-supported fuel cells with asymmetrical working and counter electrodes, reference electrode cannot be used under all circumstances. To use reference electrode for overpotential measurements, electrolyte-supported cells with symmetrical semicircular-shaped electrodes are preferred. A data processing procedure has also been presented to obtain the electrode overpotential from the measured potential using the three-electrode scheme in solid oxide electrochemical cells.

25 ENERGY STORAGE↗

inp2cjw

inp2cjw is a utility program for hierarchical finite element analysis which converts Abaqus formatted (inp) model output from CUBIT to a format (cjw) for HFEM and HFEM-Hydro modeling codes while introducing volume, facet, and edge material properties into the model. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-4551 O

Weiss, ChesterJ↗

Cubit v.16.X

SAND2021-3051 O CUBIT is a full-featured software toolkit for robust generation of two- and three-dimensional finite element meshes (grids) and geometry preparation. Its main goal is to reduce the time required to generate meshes - particularly large hex meshes of complicated, interlocking assemblies. It is a solid-modeler-based preprocessor that meshes volumes and surfaces for finite element analysis. CUBIT also includes state-of-the-art smoothing algorithms and provides an extensive suite of tools for geometry decomposition and mesh generation. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Stimpson, Clint↗

exodus_helper

SAND2022-7288 O The exodus_helper Python package is a set of tools that allows a user to query and create datasets in the Exodus II format for use in finite element analysis. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Alleman, Coleman↗