Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element Modeling”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Review on failure prediction techniques of composite single lap joint

Adhesive bonding is the most appropriate joining method in construction of composite structures. The use of reliable design and prediction technique will produce better performance of bonded joints. Several papers from recent papers and journals have been reviewed and synthesized to understand the current state of the art in this area. It is done by studying the most relevant analytical solutions for composite adherends with start of reviewing the most fundamental ones involving beam/plate theory. It is then extended to review single lap joint non linearity and failure prediction and finally on the failure prediction on composite single lap joint. The review also encompasses the finite element modelling part as tool to predict the elastic response of composite single lap joint and failure prediction numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

HydroChrono: An Open-Source Hydrodynamics Package for Project Chrono

In this paper we present the development and verification of HydroChrono, a hydrodynamics package for the Project Chrono physics engine. This package includes the implementation of hydrodynamics equations, the added mass for multibody systems, the development of I/O functions as well as a Python API, and comparison against standard reference cases and other existing tools. HydroChrono provides a flexible, fully open-source solution for simulating wave energy converters (WECs), floating offshore wind turbines (FOWTs) platforms, and other hydrodynamic systems. Here we show, via comparisons with existing tools for benchmark verification cases, that HydroChrono accurately models hydrodynamic forces - making it a useful tool for the design and optimization of these systems. Additionally, the integration of HydroChrono with Project Chrono offers access to finite element modeling capabilities and high-fidelity modelling - with Chrono's existing coupling to CFD and SPH codes. This provides numerical modelers with a multifidelity simulation framework for designing and validating these systems. The development of HydroChrono provides a new, open-source solution for simulating hydrodynamic systems. Its compatibility with other simulation tools enables a more streamlined and efficient design process, advancing the field and providing new opportunities for innovation in this area.

BEM↗
NASA NTRS1981

Static shape determination and control for a large space antenna

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and quadratic performance criteria for the comparison of estimates or control forces. Results are illustrated by simulations with a finite element model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic model. These expansions provide approximate solutions in the event that a closed form analytical solution to the system boundary value problem is not available.

Weeks, C.↗
NASA NTRS1993

Influence of local fiber undulations on the global buckling behavior of filament-wound cylinders

The predicted compressive stiffness and buckling strength of filament-wound cylinders using classical lamination theory is significantly higher than those observed experimentally. This discrepancy is partially influenced by the variation of mechanical properties in the region of fiber undulations. These regions are localized geometric defects intrinsic to the filament-winding, weaving, and braiding processes. In the present work, the average mechanical properties of the fiber undulation region are quantified using modified models of woven-fabric composites to account for the 3-dimensional effects. The mechanical properties thus determined can be incorporated as local element properties into global finite-element models. Preliminary results from large-displacement analyses of filament-wound cylinders are relatively more accurate when fiber undulations are accounted for.

Jensen, D. W.↗
NASA NTRS1995

A software tool to design thermal barrier coatings

This paper summarizes work completed for a NASA Phase 1 SBIR program which demonstrated the feasibility of developing a software tool to aid in the design of thermal barrier coating (TBC) systems. Toward this goal, three tasks were undertaken and completed. Task 1 involved the development of a database containing the pertinent thermal and mechanical property data for the top coat, bond coat and substrate materials that comprise a TBC system. Task 2 involved the development of an automated set-up program for generating two dimensional (2D) finite element models of TBC systems. Most importantly, Task 3 involved the generation of a rule base to aid in the design of a TBC system. These rules were based on a factorial design of experiments involving FEM results, and were generated using a Yates analysis. A previous study has indicated the suitability and benefit of applying finite element analysis to perform computer based experiments to decrease but not eliminate physical experiments on TBC's. This program proved feasibility by expanding on these findings by developing a larger knowledge base and developing a procedure to extract rules to aid in TBC design.

Petrus, G.↗
NASA NTRS1995

A software tool to design thermal barrier coatings

This paper summarizes work completed for a NASA Phase 1 SBIR program which demonstrated the feasibility of developing a software tool to aid in the design of thermal barrier coating (TBC) systems. Toward this goal, three tasks were undertaken and completed. Task 1 involved the development of a database containing the pertinent thermal and mechanical property data for the top coat, bond coat and substrate materials that comprise a TBC system. Task 2 involved the development of an automated set-up program for generating two dimensional (2D) finite element models of TBC systems. Most importantly, task 3 involved the generation of a rule base to aid in the design of a TBC system. These rules were based on a factorial design of experiments involving FEM results and were generated using a Yates analysis. A previous study had indicated the suitability and benefit of applying finite element analysis to perform computer based experiments to decrease but not eliminate physical experiments on TBC's. This program proved feasibility by expanding on these findings by developing a larger knowledgebase and developing a procedure to extract rules to aid in TBC design.

Petrus, Gregory↗
NASA NTRS1994

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗
NASA NTRS1992

Fiber shape effects on metal matrix composite behavior

The effects of different fiber shapes on the behavior of metal matrix composites is computationally simulated. A three-dimensional finite element model consisting of a group of nine unidirectional fibers in a three by three unit cell array of a SiC/Ti-15-3 metal matrix composite is used in the analysis. The model is employed to represent five fiber shapes that include a circle, an ellipse, a kidney, and two different cross shapes. The distribution of stresses and the composite material properties, such as moduli, coefficients of thermal expansion, and Poisson's ratios, are obtained from the finite element analysis using the various fiber shapes. Comparisons of these results are used to determine the sensitivity of the composite behavior to the different fiber shapes. In general, fiber dominated properties are not affected by fiber geometry and matrix dominated properties are only moderately affected.

Brown, H. C.↗
DOE OSTI2021

Structural dynamics modeling of spent nuclear fuel during hypothetical package drop events

The response of spent nuclear fuel (SNF) to hypothetical package drop events is of particular interest in the scope of spent fuel storage and transportation because of the mechanical shock encountered in such scenarios. Previous testing and modeling by the U.S. Department of Energy has demonstrated that the shock and vibration environment of normal shipping and handling conditions (excluding package drop events) is relatively benign and does not challenge the integrity of spent nuclear fuel. Cask drop events are worth considering because SNF packages are required to withstand free drops onto unyielding surfaces as part of their licensing basis. The acceleration experienced during drop events can be orders of magnitude higher, and thus more advanced models are needed to encompass potential nonlinear behavior of the fuel, such as spacer grid buckling and rod-to-rod impact. This work describes a number of finite element models developed to calculate the response of spent nuclear fuel to various hypothetical drop events that have been validated by package and fuel assembly drop tests conducted in the last decade. Sensitivity of the model response to factors such as package drop orientation, secondary impacts, and irradiated material properties as well as their potential impacts to fuel cladding integrity, was also investigated. Cask drops are not expected as a regular occurrence during SNF transportation, but this work helps raise the understanding of SNF mechanical loads to the point of consistency with the package design requirements.

Kadooka, Kevin↗
DOE OSTI2019

Structural Dynamic Analysis of Spent Nuclear Fuel

This unlimited distribution milestone report is the deliverable for M2SF-19PN010202014. This report completes the modeling and analysis related to the Equipos Nucleares Sociedad Anónima/U.S. Department of Energy (ENSA/DOE) spent nuclear fuel (SNF) multimodal transportation test (MMTT) campaign performed in 2017. A previous report studied the shock and vibration data generated during the MMTT and determined effective modeling and analysis methodologies and strategies. This report completes the modeling and analysis activity by applying the structural dynamic analysis models and methods to the Atlas railcar system. This analysis demonstrates that the MMTT configuration provided a conservative, bounding shock and vibration environment for SNF in comparison to finite element model results that are configured to represent the Atlas railcar system and all of its SNF package combinations. Additionally, this report establishes a technical basis for recommending that no further analysis is necessary to demonstrate that SNF cladding will remain intact during shock and vibration loading conditions that occur during normal conditions of transportation (NCT). 10 CFR 71.71 defines package tests to demonstrate acceptable package performance under NCT conditions. The focus of this work is to demonstrate that SNF cladding will remain intact within the package while it is subjected to what 10 CFR 71.71 calls: Vibration normally incident to transport. The MMTT recorded negligible shock and vibrations on the SNF cladding during the test campaign, and this analysis supports the conclusion that intact fuel rods carried in any current SNF transportation package will experience similar, negligible, shock and vibration conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1986

Weight estimation of unconventional structures by structural optimization

Automated techniques are presented that are used in structural optimization technology, with emphasis on modifications of finite element models to obtain an optimal material distribution for minimum weight while satisfying the prescribed design requirements. It is anticipated that the future development of computer aided engineering (CAE) system will provide environments where structural analysis, a design optimization, and weight evaluation modules are integrated, sharing a common data base. Structural optimization capabilities obtained by integrating a finite element structural analysis program and a numerical optimization code are developed and applied to two illustrative examples: marine gear housing structural weight minimization and joined wing structures.

Miura, Hirokazu↗
NASA NTRS2020

Structural Analysis of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the agency’s aging wind tunnel infrastructure. Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper a process for generating structural finite element models of these wind tunnel blades is presented. The model developed in this paper was analyzed under loading conditions from tests conducted on the original blades and the results are compared with legacy strain gage data.

Brian H Mason↗
DOE OSTI2023

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2015

Temperature Profile in Fuel and Tie-Tubes for Nuclear Thermal Propulsion Systems

A finite element method to calculate temperature profiles in heterogeneous geometries of tie-tube moderated LEU nuclear thermal propulsion systems and HEU designs with tie-tubes is developed and implemented in MATLAB. This new method is compared to previous methods to demonstrate shortcomings in those methods. Typical methods to analyze peak fuel centerline temperature in hexagonal geometries rely on spatial homogenization to derive an analytical expression. These methods are not applicable to cores with tie-tube elements because conduction to tie-tubes cannot be accurately modeled with the homogenized models. The fuel centerline temperature directly impacts safety and performance so it must be predicted carefully. The temperature profile in tie-tubes is also important when high temperatures are expected in the fuel because conduction to the tie-tubes may cause melting in tie-tubes, which may set maximum allowable performance. Estimations of maximum tie-tube temperature can be found from equivalent tube methods, however this method tends to be approximate and overly conservative. A finite element model of heat conduction on a unit cell can model spatial dependence and non-linear conductivity for fuel and tie-tube systems allowing for higher design fidelity of Nuclear Thermal Propulsion.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1986

Thermomechanical design criteria for ceramic-coated surfaces

Some early history of ceramic applications is presented. Finite element modeling of components to determine service and fabrication loads found inelastic behavior and residual stresses to be significant to component life. Inelastic behavior mitigates peak strains but enhances residual strains. Results of furnace, Mach 0.3 burner, and engine tests are discussed and categorized into design criteria (loading, geometry, fabrication, materials, analysis, and testing). These design rules and finite element analyses are brought to bear on two test cases: turboshaft engine seals, and rocket thrust chambers.

Mullen, R. L.↗