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

Deformation mechanisms in a precipitation-strengthened ferritic super alloy revealed by in situ neutron dffraction studies at elevated temperatures

The ferritic superalloy Fe–10Ni–6.5Al–10Cr–3.4Mo strengthened by ordered (Ni,Fe)AlB2-type precipitates is a candidate material for ultra-supercritical steam turbine applications above 923 K. Despite earlier success in improving its room-temperature ductility, the creep resistance of this material at high temperatures needs to be further improved, which requires a fundamental understanding of the high-temperature deformation mechanisms at the scales of individual phases and grains. In situ neutron diffraction has been utilized to investigate the lattice strain evolution and the microscopic load-sharing mechanisms during tensile deformation of this ferritic superalloy at elevated temperatures. Finite-element simulations based on the crystal plasticity theory are employed and compared with the experimental results, both qualitatively and quantitatively. Based on these interphase and intergranular load-partitioning studies, it is found that the deformation mechanisms change from dislocation slip to those related to dislocation climb, diffusional flow and possibly grain boundary sliding, below and above 873 K, respectively. Insights into microstructural design for enhancing creep resistance are also discussed.

nuclear (including radiation effects), defects, me↗
DOE OSTI2019

AC Resistance Reduction Using Orthogonal Air Gaps in High-Frequency Inductors: Preprint

This paper presents a relatively simple technique to reduce winding losses due to fringing fields in high-frequency inductors. In high-frequency power electronics, ac inductor winding losses are affected by skin and proximity effects, including uneven current distribution due to fringing magnetic fields around airgaps. It is well known how fringing effects can be mitigated using distributed airgaps, at the expense of non-standard core or winding geometry. In planar magnetics, the orthogonal-airgap approach proposed in this paper combines airgaps in core segments parallel with the windings with conventionally placed airgaps in segments perpendicular to the windings. The approach is developed using a 1D analytical framework and validated by 2D finite-element simulations. Analytical guidelines are presented to optimize the airgaps to achieve minimum ac resistance. As a case study, a planar inductor is designed for an 8~kW SiC-based buck converter operating at 250~kHz. It is shown how the orthogonal airgaps result in more than 45\% reduction in ac resistance and substantially reduced inductor losses compared to the design using standard airgaps. The results are verified by loss measurements on an experimental converter prototype.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Frequency Domain Reflectometry NDE for Aging Cables in Nuclear Power Plants

Cable insulation polymers are among the more susceptible materials to age-related degradation within a nuclear power plant. This is recognized by both regulators and utilities, so all plants have developed cable aging management programs to detect damage before critical component failure in compliance with regulatory guidelines. Although a wide range of tools are available to evaluate cables and cable systems, cable aging management programs vary in how condition monitoring and nondestructive examinations are conducted as utilities search for the most reliable and cost-effective ways to assess cable system condition. Frequency domain reflectometry (FDR) is emerging as one valuable tool to locate and assess damaged portions of a cable system with minimal cost and only requires access in most cases to one of the cable terminal ends. Since laboratory studies to evaluate the use of FDR for inspection of aged cables can be expensive and data interpretation may be confounded by multiple factors which influence results, a model-based approach is desired to parametrically investigate the effect of insulation material damage in a controlled manner. This work describes development of a physics-based FDR model which uses finite element simulations of cable segments in conjunction with cascaded circuit element simulations to efficiently study a cable system. One or more segments of the cable system model have altered physical or electrical properties which represent the degree of damage and the location of the damage in the system. This circuit model is then subjected to a simulated FDR examination. The modeling approach is verified using several experimental cases and by comparing it to a commercial simulator suitable for simulation of some cable configurations. The model is used to examine a broad range of parameters including defect length, defect profile, degree of degradation, number and location of defects, FDR bandwidth, and addition of impedance-matched extensions to minimize the end-shadow effect.

Cable Inspection, Frequency Domain Reflectometry↗
DOE OSTI2017

Effect of Plate Hardening Behavior on the Deformation of Stainless Steel Metal Bellows

Tensile tests of original plate samples from three types of stainless steel metal bellows were performed at room temperature. The constitutive equations for the three hardening curves were obtained and fitted. The analysis results of the microstructure and fracture morphology of the tensile specimens show that the grain size of the plate with a high logarithmic–exponential hardening rate is uneven and the dimple of the shear fracture is elongated into an ellipse. By contrast, the grain size of the plate with a relatively low linear hardening rate is even and the dimple of the fracture is uniformly equiaxial. Finite element simulations of the hydraulic bulging and repeated limit bending deformation of the metal bellows of the three types of materials were also conducted. The repeated limit bending deformation process was tested experimentally. Although the effect of the hardening exponent on the residual stress of the metal bellows after hydraulic bulging is minimal, this exponent considerably influences the repeated limit bending deformation of the metal bellows after subsequent use. The trough hardening phenomenon is serious in the repeated limit bending process. Moreover, when the hardening exponent of the original plate is high, the resistance to bending fracture at the trough area is poor.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effect of Weld Consumable Conditioning on the Diffusible Hydrogen and Subsequent Residual Stress and Flexural Strength of Multipass Welded P91 Steels

P91 steel weld joint was prepared using the shielded metal arc welding process and four different conditions of weld consumable that provide the different levels of diffusible hydrogen in deposited metal. In the present research, the effects of diffusible hydrogen content on the flexural strength, lower critical stress, and tensile strength of P91 steel welds were also determined with respect to different electrode conditions. To investigate the effect of diffusible hydrogen on multipass welding, top and root side flexural tests were performed. The residual stresses (axial stress and transverse stress) were also measured using the blind hole drilling method for different conditions of welding consumable. The peak value of residual stresses was measured at the center of the weld fusion zone. The maximum value of transverse stress was measured to be 355 MPa for case II (6.21 mL/100 g of diffusible hydrogen), while the maximum axial stress was about 218 MPa for case IV (12.43 mL/100 g of diffusible hydrogen). A three-dimensional finite element simulation was also performed to predict the residual stress distribution and thermal profile along the welded joint. The experimentally determined residual stresses correlated well with the numerically estimated residual stresses. The diffusible hydrogen content was not observed to have any significant effect on the residual stresses. The corrected residual stress values were also predicted by considering the plasticity-induced error. However, the flexural performance of the welded joint was affected by the diffusible hydrogen content. The top and root flexural strength was measured to be optimum for the low level of diffusible hydrogen content, and the values decreased with an increase in diffusible hydrogen content.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Design and Application of CSC-Roll for Heavy Reduction of the Bloom Continuous Casting Process

Center quality control is a key issue affecting large-section bloom production and product quality. The central convex roll could avoid the bloom side part and yield a high reduction amount, which benefit the implement of heavy reduction (HR) for bloom continuous casting process without upgrading the mechanical and hydraulic system of withdrawal and straightening machine. In this work, a curving surface convex roll (CSC-Roll), which could significantly improve bloom center quality, was designed by finite element simulation for a 360 mm × 450 mm #45 steel continuous casting bloom. The optimal roll surface shape of the CSC-Roll was determined based on an analysis of the distribution of the reduction amount in the bloom center and the strain of the bloom centerline and inner surface. The calculation results showed that the optimal size of the CSC-Roll was 200 mm in convex platform region length, 30 mm in convex platform region height, and 80 mm in gradient curvature region length. The bloom heat transfer and deformation behavior in the CSC-Roll reduction process were then analyzed. The results showed that in the CSC-Roll HR process, 61.71 pct of the reduction amount could be transmitted to the inner side of the porosity zone. Macrographs of the 360 mm × 450 mm #45 steel continuous casting bloom proved that the bloom center quality was significantly improved after CSC-Roll HR was applied.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Crash simulation of composite and aluminum helicopter fuselages using a finite-element program

This paper describes a mathematical investigation of the crash-impact responses of an all-composite helicopter cockpit section incorporating an energy absorbing concept and one of conventional aluminum construction, using the Grumman DYCAST finite-element nonlinear structural dynamics computer program. The overall objective of this initial investigation was to explore the application of DYCAST to the design of future helicopter structures. The specific purposes are to assess the value of the finite-element code as a crashworthiness design analysis tool, and to compare the responses of the two fuselage sections. The results indicate that the finite-element simulation is a powerful analysis tool, providing sufficient detail to evaluate individual structural components, and that the crashworthiness of an all-composite cockpit section could be enhanced with energy-absorbing concepts

Winter, R.↗
NASA NTRS1997

Finite Element Aircraft Simulation of Turbulence

A turbulence model has been developed for realtime aircraft simulation that accommodates stochastic turbulence and distributed discrete gusts as a function of the terrain. This model is applicable to conventional aircraft, V/STOL aircraft, and disc rotor model helicopter simulations. Vehicle angular activity in response to turbulence is computed from geometrical and temporal relationships rather than by using the conventional continuum approximations that assume uniform gust immersion and low frequency responses. By using techniques similar to those recently developed for blade-element rotor models, the angular-rate filters of conventional turbulence models are not required. The model produces rotational rates as well as air mass translational velocities in response to both stochastic and deterministic disturbances, where the discrete gusts and turbulence magnitudes may be correlated with significant terrain features or ship models. Assuming isotropy, a two-dimensional vertical turbulence field is created. A novel Gaussian interpolation technique is used to distribute vertical turbulence on the wing span or lateral rotor disc, and this distribution is used to compute roll responses. Air mass velocities are applied at significant centers of pressure in the computation of the aircraft's pitch and roll responses.

McFarland, R. E.↗
NASA NTRS2019

Quadruple Simulations of Thermographic Inspections of Impacted Composites

Thermography has been shown to be a viable technique for inspection of composites. Impact damage in composites typically contains multiple overlapping delaminations at different depths. Understanding the limitations of the thermographic inspection is enhanced by performing simulations of the technique. Most simulations of composite thermographic inspections have focused on simulations of a single delamination at a fixed depth. The quadrupole method has been shown as a viable technique for rapid three-dimensional thermographic simulations of a delamination. This method is expanded to enable rapid simulation of multiple overlapping delaminations at different depths. Quadrupole simulations are compared to finite element simulations of multiple delaminations at different depths. The simulations are also compared to the thermographic measurements on impacted composites where shape and depth of the delaminations are known from x-ray computed tomography data.

Winfree, William P.↗
NASA NTRSDate not supplied

Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. The quadrupole method has been shown as a viable technique for rapidedly simulating the thermal repsonse of layered systems. Often, the effort has focused on a one-dimentional models, in particular for improved analysis of thermal data. For composites, the in-plane heat diffusion often significantly impacts the thermal response of defects of interest, therefore three-dimensional simulations are desirable. This paper discusses the extension of the quadrupole methodology to perform simulations of thermographic responses in three-dimensional configuration. This enables the more realistic simulations of the thermal response of delaminiations in composites. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are known obtained from other sources.

Thermography↗
NASA NTRSDate not supplied

Simulation of Thermographic Responses of Delaminations in Composites with Quadrupole Method

Thermography has been shown to be a viable technique for inspection of composites. The quadrupole method is a valuable technique for rapidly simulating the thermal response of layered systems. Often, the effort has focused on a one-dimensional model, in particular for improved analysis of thermal data. For composites, the in-plane heat diffusion often significantly impacts the thermal response of defects of interest, therefore three-dimensional simulations are desirable. This paper discusses the extension of the quadrupole methodology to perform simulations of thermographic responses in three-dimensional configurations. This enables the more realistic simulations of the thermal response of delaminations in composites. The simulations are compared to finite element simulations of the same inspection configurations.

Thermography↗
DOE OSTI2020

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
NASA NTRS2008

An Approach to Assess Delamination Propagation Simulation Capabilities in Commercial Finite Element Codes

An approach for assessing the delamination propagation simulation capabilities in commercial finite element codes is presented and demonstrated. For this investigation, the Double Cantilever Beam (DCB) specimen and the Single Leg Bending (SLB) specimen were chosen for full three-dimensional finite element simulations. First, benchmark results were created for both specimens. Second, starting from an initially straight front, the delamination was allowed to propagate. The load-displacement relationship and the total strain energy obtained from the propagation analysis results and the benchmark results were compared and good agreements could be achieved by selecting the appropriate input parameters. Selecting the appropriate input parameters, however, was not straightforward and often required an iterative procedure. Qualitatively, the delamination front computed for the DCB specimen did not take the shape of a curved front as expected. However, the analysis of the SLB specimen yielded a curved front as was expected from the distribution of the energy release rate and the failure index across the width of the specimen. Overall, the results are encouraging but further assessment on a structural level is required.

Krueger, Ronald↗