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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

The Simulation of Precipitation Evolutions and Mechanical Properties in Friction Stir Welding with Post-Weld Heat Treatments

A finite element model of friction stir welding capable of re-meshing is used to simulate the temperature variations. Re-meshing of the finite element model is used to maintain a fine mesh resolving the gradients of the solution. The Kampmann–Wagner numerical model for precipitation is then used to study the relation between friction stir welds with post-weld heat treatment (PWHT) and the changes in mechanical properties. Results indicate that the PWHT holding time and PWHT holding temperature need to be optimally designed to obtain FSW with better mechanical properties. Higher precipitate number with lower precipitate sizes gives higher strength in the stirring zone after PWHT. The coarsening of precipitates in HAZ are the main reason to hinder the improvement of mechanical property when PWHT is used.

36 MATERIALS SCIENCE↗
NASA NTRS1987

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗
NASA NTRS1988

The LDCM actuator for vibration suppression

A linear dc motor (LDCM) has been proposed as an actuator for the COFS I mast and the COFS program ground test Mini-Mast. The basic principles of operation of the LDCM as an actuator for vibration suppression in large flexible structures are reviewed. Because of force and stroke limitations, control loops are required to stabilize the actuator, which results in a non-standard actuator-plant configuration. A simulation model that includes LDCM actuator control loops and a finite element model of the Mast is described, with simulation results showing the excitation capability of the actuator.

Ide, Eric N.↗
DOE OSTI2018

Simulation of internal stress waves generated by laser-induced damage in multilayer dielectric gratings

Multilayer dielectric (MLD) gratings used in ultrahigh-intensity laser systems often exhibit a laser-induced damage performance below that of their constituent materials. Reduced performance may arise from fabrication- and/or design-related issues. Finite element models were developed to simulate stress waves in MLD grating structures generated by laser-induced damage events. Here, these models specifically investigate the influence of geometric and material parameters on how stress waves can lead to degradation of material structural integrity that can have adverse effects on its optical performance under subsequent laser irradiation: closer impedance matching of the layer materials reduces maximum interface stresses by ~20% to 30%; increasing sole thickness from 50 nm to 500 nm reduces maximum interface stresses by ~50%.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

TH-AB-BRA-09: Stability Analysis of a Novel Dose Calculation Algorithm for MRI Guided Radiotherapy

Purpose: To determine the iterative deterministic solution stability of the Linear Boltzmann Transport Equation (LBTE) in the presence of magnetic fields. Methods: The LBTE with magnetic fields under investigation is derived using a discrete ordinates approach. The stability analysis is performed using analytical and numerical methods. Analytically, the spectral Fourier analysis is used to obtain the convergence rate of the source iteration procedures based on finding the largest eigenvalue of the iterative operator. This eigenvalue is a function of relevant physical parameters, such as magnetic field strength and material properties, and provides essential information about the domain of applicability required for clinically optimal parameter selection and maximum speed of convergence. The analytical results are reinforced by numerical simulations performed using the same discrete ordinates method in angle, and a discontinuous finite element spatial approach. Results: The spectral radius for the source iteration technique of the time independent transport equation with isotropic and anisotropic scattering centers inside infinite 3D medium is equal to the ratio of differential and total cross sections. The result is confirmed numerically by solving LBTE and is in full agreement with previously published results. The addition of magnetic field reveals that the convergence becomes dependent on the strength of magnetic field, the energy group discretization, and the order of anisotropic expansion. Conclusion: The source iteration technique for solving the LBTE with magnetic fields with the discrete ordinates method leads to divergent solutions in the limiting cases of small energy discretizations and high magnetic field strengths. Future investigations into non-stationary Krylov subspace techniques as an iterative solver will be performed as this has been shown to produce greater stability than source iteration. Furthermore, a stability analysis of a discontinuous finite element space-angle approach (which has been shown to provide the greatest stability) will also be investigated. Dr. B Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization)

60 APPLIED LIFE SCIENCES↗
NASA NTRS2014

Simulating Damage Due to a Lightning Strike Event: Effects of Temperature Dependent Properties on Interlaminar Damage

A multidirectional, carbon fiber-epoxy, composite panel is subjected to a simulated lightning strike, within a finite element method framework, and the effect of material properties on the failure (delamination) response is investigated through a detailed numerical study. The numerical model of the composite panel consists of individual homogenized plies with user-defined, cohesive interface elements between them. Lightning strikes are simulated as an assumed combination of excessive heat and high pressure loadings. It is observed that the initiation and propagation of lightning-induced delamination is a significant function of the temperature dependency of interfacial fracture toughness. This dependency must be defined properly in order to achieve reliable predictions of the present lightning-induced delamination in the composite panel.

computational mechanics↗
DOE OSTI2014

Finite volume and finite element methods applied to 3D laminar and turbulent channel flows

The work deals with numerical simulations of incompressible flow in channels with rectangular cross section. The rectangular cross section itself leads to development of various secondary flow patterns, where accuracy of simulation is influenced by numerical viscosity of the scheme and by turbulence modeling. In this work some developments of stabilized finite element method are presented. Its results are compared with those of an implicit finite volume method also described, in laminar and turbulent flows. It is shown that numerical viscosity can cause errors of same magnitude as different turbulence models. The finite volume method is also applied to 3D turbulent flow around backward facing step and good agreement with 3D experimental results is obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Multidisciplinary tailoring of hot composite structures

A computational simulation procedure is described for multidisciplinary analysis and tailoring of multilayered multimaterial hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loadings including the effect of aggressive environments. The simulation is based on a 3D finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., the laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermomechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loadings lead to different tailored designs, even those opposite of each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.

Singhal, Surendra N.↗
DOE OSTI2016

Computational homogenization at extreme scales

Multi-scale simulations at extreme scales in terms of both physical length scales and computational resources are presented. In this letter, we introduce a hierarchically parallel computational homogenization solver that employs hundreds of thousands of computing cores and resolves in material length scales (from to ). Simulations of this kind are important in understanding the multi-scale essence of many natural and synthetically made materials. Thus, we present a simulation consisting of 53.8 Billion finite elements with 28.1 Billion nonlinear equations that is solved on 393,216 computing cores (786,432 threads). The excellent parallel performance of the computational homogenization solver is demonstrated by a strong scaling test from 4,096 to 262,144 cores. A fully coupled multi-scale damage simulation shows a complex crack profile at the micro-scale and the macroscopic crack tunneling phenomenon. Such large and predictive simulations are an important step towards Virtual Materials Testing and can aid in development of new material formulations with extreme properties. Furthermore, the high computational efficiency of our computational homogenization solver holds great promise for utilizing the next generation of exascale parallel computing platforms that are expected to accelerate computations through orders of magnitude increase in parallelism rather than speed of each processor.

Engineering↗
NASA NTRS1992

Mixed-mode fracture in unidirectional graphite epoxy composite laminates with central notch

Mixed-mode matrix fracture in central notched off-axis unidirectional composite laminates was investigated. A limited number of unidirectional tensile type specimens with a central, horizontal, notch were tested. Crack initiation and propagation were examined under various local stress fields that were controlled by fiber orientations. The tested specimens were simulated using a two dimensional finite element method with constant strain loading. The strain energy release rates along the crack were evaluated via crack closure technique. The variation of critical strain energy rates with off-axis angle was studied. The results from single (one-sided) and double (two-sided) crack simulations were presented and compared.

Binienda, Wieslaw K.↗
DOE OSTI2015

Finite element analysis of an extended end-plate connection using the T-stub approach

Beam-to-column end-plate bolted connections are usually used as moment-resistant connections in steel framed structures. For this joint type, the deformability is governed by the deformation capacity of the column flange and end-plate under tension and elongation of the bolts. All these elements around the beam tension flange form the tension region of the joint, which can be modeled by means of equivalent T-stubs. In this paper a beam-to-column end-plate bolted connection is substituted with a T-stub of appropriate effective length and it is analyzed using the commercially available finite element software ABAQUS. The performance of the model is validated by comparing the behavior of the T-stub from the numerical simulation with the behavior of the connection as a whole. The moment-rotation curve of the T-stub obtained from the numerical simulation is compared with the behavior of the whole extended end-plate connection, obtained by numerical simulation, experimental tests and analytical approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗
NASA NTRS2011

Three-Dimensional Finite Element Ablative Thermal Response and Thermostructural Design of Thermal Protection Systems

A finite element ablation and thermal response program is presented for simulation of three-dimensional transient thermostructural analysis. The three-dimensional governing differential equations and finite element formulation are summarized. A novel probabilistic design methodology for thermal protection systems is presented. The design methodology is an eight step process beginning with a parameter sensitivity study and is followed by a deterministic analysis whereby an optimum design can determined. The design process concludes with a Monte Carlo simulation where the probabilities of exceeding design specifications are estimated. The design methodology is demonstrated by applying the methodology to the carbon phenolic compression pads of the Crew Exploration Vehicle. The maximum allowed values of bondline temperature and tensile stress are used as the design specifications in this study.

Dec, John A.↗
DOE OSTI2018

Prediction of residual stresses in a multipass pipe weld by a novel 3D finite element approach

Due to enormous computation cost, current residual stress simulation of multipass girth welds are mostly performed using two-dimensional (2D) axisymmetric models. The 2D model can only provide limited estimation on the residual stresses by assuming its axisymmetric distribution. In this study, a highly efficient thermal-mechanical finite element code for three-dimensional (3D) model has been developed based on high performance Graphics Processing Unit (GPU) computers. Our code is further accelerated by considering the unique physics associated with welding processes that are characterized by steep temperature gradient and a moving arc heat source. It is capable of modeling large-scale welding problems that cannot be easily handled by the existing commercial simulation tools. To demonstrate the accuracy and efficiency, our code was compared with a commercial software by simulating a 3D multi-pass girth weld model with over 1 million elements. Our code achieved comparable solution accuracy with respect to the commercial one but with over 100 times saving on computational cost. Moreover, the three-dimensional analysis demonstrated more realistic stress distribution that is not axisymmetric in hoop direction.

Huang, Hui↗
NASA NTRS1977

Status of the Electroforming Shield Design (ESD) project

The utilization of a digital computer to augment electrodeposition/electroforming processes in which nonconducting shielding controls local cathodic current distribution is reported. The primary underlying philosophy of the physics of electrodeposition was presented. The technical approach taken to analytically simulate electrolytic tank variables was also included. A FORTRAN computer program has been developed and implemented. The program utilized finite element techniques and electrostatic theory to simulate electropotential fields and ionic transport.

Fletcher, R. E.↗
NASA NTRS2004

Quantifying Wrinkle Features of Thin Membrane Structures

For future micro-systems utilizing membrane based structures, quantified predictions of wrinkling behavior in terms of amplitude, angle and wavelength are needed to optimize the efficiency and integrity of such structures, as well as their associated control systems. For numerical analyses performed in the past, limitations on the accuracy of membrane distortion simulations have often been related to the assumptions made. This work demonstrates that critical assumptions include: effects of gravity, supposed initial or boundary conditions, and the type of element used to model the membrane. In this work, a 0.2 m x 02 m membrane is treated as a structural material with non-negligible bending stiffness. Finite element modeling is used to simulate wrinkling behavior due to a constant applied in-plane shear load. Membrane thickness, gravity effects, and initial imperfections with respect to flatness were varied in numerous nonlinear analysis cases. Significant findings include notable variations in wrinkle modes for thickness in the range of 50 microns to 1000 microns, which also depend on the presence of an applied gravity field. However, it is revealed that relationships between overall strain energy density and thickness for cases with differing initial conditions are independent of assumed initial conditions. In addition, analysis results indicate that the relationship between wrinkle amplitude scale (W/t) and structural scale (L/t) is independent of the nonlinear relationship between thickness and stiffness.

Jacobson, Mindy B.↗
DOE OSTI2019

Simulation of Fast-Rise Transients in a Large-Power Transformer Winding

A high-altitude electromagnetic pulse (HEMP) may pose a significant threat to the electrical power grid. In this context, the present study considers here the potential impact of the HEMP early-time component (E1) on the winding of a typical transmission power line transformer. A time-domain, fully electromagnetic, three-dimensional, finite-element model was used to simulate the propagation of an E1-induced voltage surge into the transformer winding. The analysis is focused on the first few turns of the winding, the ones that are impacted by the largest electric field resulting from the applied pulse. The model follows the transient dynamics, while the pulse is traveling along the conductor, and the electric field is established across the insulation gaps through capacitive coupling. The simulations also considered the effect of an electrically floating, “shield” winding conductor. This is used as a countermeasure to improve the voltage uniformity across the winding during fast-rise transients. In this case, the propagation of the E1-induced pulse is compared with that from a slower rising, lightning standard waveform. The simulation results show that the presence of the shield conductor can effectively reduce the peak electric field between winding turns, in the presence of either E1 or lightning waveforms.

24 POWER TRANSMISSION AND DISTRIBUTION↗