Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Developing a Finite Element Soyuz Seat Environment and Investigating Occupant Injury Metrics for Accelerative Loading Across Spaceflight Vehicle Designs

Landings in the Soyuz vehicle have generated more injuries than models have predicted, and the reason for this remains unknown. Presently, 3 new capsule-type spacecraft—the Orion, the SpaceX Dragon2, and the Boeing CST-100—are in various stages of development. To better understand the injury risk to Soyuz crewmembers, and to mitigate risk of future injuries onboard the new vehicles, we conducted a finite element (FE) modeling effort. At the outset of this study, previously developed FE models were available for the seats in the 3 new vehicles, whereas the Soyuz environment had not been modeled in an extensible manner. Thus, a model of the Soyuz seat was necessary to facilitate comparisons of all seat types.

D. A. Jones↗
DOE OSTI2020

Full circle mechanical dynamic characterization including experimental modal analysis and finite element analysis

During operation, it was observed that a specific mechanical system experienced undesirable vibration and it became necessary to understand and mitigate this phenomenon. This document investigates the tools, methodology, and results of the dynamic characterization of the system. The characterization makes use of the experimental modal analysis (EMA) methods of single input multiple output (SIMO) and single input single output (SISO). The validity of the theory of reciprocity is confirmed to minimize measurement error, cost, and time of repeat testing. Finite element analysis (FEA) is used in choosing transducer and modal impact locations to adequately characterize the system. Single degree of freedom (SDOF) and multiple degree of freedom (MDOF) curve fitting is used to fully characterize the system’s mode shapes and natural frequencies. The EMA characterization results are used to modify and validate the FEA model so that FEA can be used to model potential structural modifications to the system to mitigate the undesirable vibration. Structural modifications are chosen, implemented, and their effectiveness is quantified using EMA. Finally, a qualitative evaluation of the methodology of FEA validation by EMA and tuning of the model to match the experimental results is discussed.

42 ENGINEERING↗
NASA NTRS1987

Comparison of composite rotor blade models: A coupled-beam analysis and an MSC/NASTRAN finite-element model

A methodology was developed for the structural analysis of composite rotor blades. This coupled-beam analysis is relatively simple to use compared with alternative analysis techniques. The beam analysis was developed for thin-wall single-cell rotor structures and includes the effects of elastic coupling. This paper demonstrates the effectiveness of the new composite-beam analysis method through comparison of its results with those of an established baseline analysis technique. The baseline analysis is an MSC/NASTRAN finite-element model built up from anisotropic shell elements. Deformations are compared for three linear static load cases of centrifugal force at design rotor speed, applied torque, and lift for an ideal rotor in hover. A D-spar designed to twist under axial loading is the subject of the analysis. Results indicate the coupled-beam analysis is well within engineering accuracy.

Hodges, Robert V.↗
NASA NTRS1992

Finite element analysis of the stiffness of fabric reinforced composites

The objective of this work is the prediction of all three dimensional elastic moduli of textile fabric reinforced composites. The analysis is general enough for use with complex reinforcing geometries and capable of subsequent improvements. It places no restrictions on fabric microgeometry except that the unit cell be determinate and rectangular. The unit cell is divided into rectangular subcells in which the reinforcing geometries are easier to define and analyze. The analysis, based on inhomogeneous finite elements, is applied to a variety of weave, braid, and knit reinforced composites. Some of these predictions are correlated to test data.

Foye, R. L.↗
NASA NTRS1977

Finite-element solutions for geothermal systems

Vector potential and scalar potential are used to formulate the governing equations for a single-component and single-phase geothermal system. By assuming an initial temperature field, the fluid velocity can be determined which, in turn, is used to calculate the convective heat transfer. The energy equation is then solved by considering convected heat as a distributed source. Using the resulting temperature to compute new source terms, the final results are obtained by iterations of the procedure. Finite-element methods are proposed for modeling of realistic geothermal systems; the advantages of such methods are discussed. The developed methodology is then applied to a sample problem. Favorable agreement is obtained by comparisons with a previous study.

Chen, J. C.↗
NASA NTRS1986

A constitutive law for finite element contact problems with unclassical friction

Techniques for modeling complex, unclassical contact-friction problems arising in solid and structural mechanics are discussed. A constitutive modeling concept is employed whereby analytic relations between increments of contact surface stress (i.e., traction) and contact surface deformation (i.e., relative displacement) are developed. Because of the incremental form of these relations, they are valid for arbitrary load-deformation histories. The motivation for the development of such a constitutive law is that more realistic friction idealizations can be implemented in finite element analysis software in a consistent, straightforward manner. Of particular interest is modeling of two-body (i.e., unlubricated) metal-metal, ceramic-ceramic, and metal-ceramic contact. Interfaces involving ceramics are of engineering importance and are being considered for advanced turbine engines in which higher temperature materials offer potential for higher engine fuel efficiency.

Plesha, M. E.↗
DOE OSTI2023

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1982

Finite element analysis for bearingless rotor blade aeroelasticity

A conventional articulated rotor blade has mechanical flap and lag hinges, a lag damper, and a pitch bearing. In connection with an interest in designs of greater mechanical simplicity and increased maintainability, hingeless and bearingless rotors have been developed. A hingeless blade lacks the hinges and is cantilevered at the hub. It does have a pitch bearing for pitch control. A bearingless design eliminates the hinges and the pitch bearing as well. In the present investigation of bearingless rotor blade characteristics, finite element analysis has been successfully applied to determine the solutions of the nonlinear trim equations and the linearized flutter equations of multiple-load-path blades. The employed formulation is based on Hamilton's principle. The spatial dependence of the equations of motion is discretized by dividing the flexbeams, the torque tube, and the outboard into a number of elements.

Sivaneri, N. T.↗
DOE OSTI2022

Time reversal in a phononic crystal using finite-element modeling and an equivalent circuit model

A phononic crystal acts as a dispersive medium with a phase speed that is lower than the bulk wave speed at frequencies below the resonance of a single resonator. Time reversal is used to compensate for the phase shifts caused by individual resonators as the waves enter the medium and enable focusing of acoustic waves among the crystal. An equivalent circuit, which can predict the dispersion and attenuation of the crystal model, is shown and compared to a full-wave finite-element simulation in frequency and time. The phase shift due to a single resonator is also depicted.

Kingsley, Adam D. (ORCID:0000000226611788)↗
NASA NTRS1975

An interactive graphics package for the automatic node renumbering of finite element matrices

An interactive graphics software package which allows users to display the non-zero structure of large sparse symmetric materials was described and methods used to implement it as a portable FORTRAN callable subroutine were summarized. In particular, the system permits the display of the resulting matrix after reordering the rows and columns, with the reordering scheme either defined by the user or automatically generated by the program with the aim of reducing matrix bandwidth and profile. Although the primary application of the package has been to the finite element analysis of structures, it is equally well suited to the many other areas of engineering and science which use sparse matrices.

Boisvert, R. F.↗
NASA NTRS2012

Influence of Finite Element Size in Residual Strength Prediction of Composite Structures

The sensitivity of failure load to the element size used in a progressive failure analysis (PFA) of carbon composite center notched laminates is evaluated. The sensitivity study employs a PFA methodology previously developed by the authors consisting of Hashin-Rotem intra-laminar fiber and matrix failure criteria and a complete stress degradation scheme for damage simulation. The approach is implemented with a user defined subroutine in the ABAQUS/Explicit finite element package. The effect of element size near the notch tips on residual strength predictions was assessed for a brittle failure mode with a parametric study that included three laminates of varying material system, thickness and stacking sequence. The study resulted in the selection of an element size of 0.09 in. X 0.09 in., which was later used for predicting crack paths and failure loads in sandwich panels and monolithic laminated panels. Comparison of predicted crack paths and failure loads for these panels agreed well with experimental observations. Additionally, the element size vs. normalized failure load relationship, determined in the parametric study, was used to evaluate strength-scaling factors for three different element sizes. The failure loads predicted with all three element sizes provided converged failure loads with respect to that corresponding with the 0.09 in. X 0.09 in. element size. Though preliminary in nature, the strength-scaling concept has the potential to greatly reduce the computational time required for PFA and can enable the analysis of large scale structural components where failure is dominated by fiber failure in tension.

Satyanarayana, Arunkumar↗
NASA NTRS1977

Finite element solutions of free surface flows

A procedure is presented for using NASTRAN to determine the flow field about arbitrarily shaped bodies in the presence of a free surface. The fundamental unknown of the problem is the velocity potential which must satisfy Laplace's equation in the fluid region. Boundary conditions on the free surface may involve second order derivatives in space and time. In cases involving infinite domains either a tractable radiation condition is applied at a truncated boundary or a series expansion is used and matched to the local finite elements. Solutions are presented for harmonic, transient, and steady state problems and compared to either exact solutions or other numerical solutions.

Zarda, P. R.↗
NASA NTRSDate not supplied

Modeling Uprighting of the Orion Crew Module using a Finite Element Method with Wave Coupling and Structural Deformability

The Crew Module Uprighting System (CMUS) for the Orion Crew Module (CM) includes five bags that inflate after splashdown to rotate and upright the CM as necessary. Wave conditions can be severe enough to damage the bags and tethers, which could make the recovery operation hazardous. The CMUS therefore has design limits that ultimately constrain where the CM can land based upon sea conditions. This manuscript describes a full fluid-structure interaction computational model being used for the CMUS. It is a finite element model that uses explicit time integration to directly simulate uprighting, a computational approach usually reserved for short duration, highly transient phenomena and novelly applied here. The goal of this current effort is to determine the feasibility of this approach for CMUS uprighting, and to estimate the loads the CMUS bags experience in different wave conditions. The results from this work will be used to supplement other experimental and computational data in support of NASA’s Artemis program.

John Puryear↗
DOE OSTI2024

Finite Element Modeling of the Phase Change in Thermally-Grown SiO 2 in SiC Systems for Gas Turbines

The operating lifetimes of SiC-based components in combustion environments are directly linked to the adhesion of the protective environmental barrier coating (EBC) layer. One of the major known failure modes for EBCs is the formation of a thick SiO 2 thermally grown oxide (TGO), which decreases coating adhesion and encourages eventual coating spallation. The effect of the TGO thickness under Yb 2 Si 2 O 7 EBCs on silicon carbide was investigated using finite element models (FEMs) with various interfacial architectures and SiO 2 TGO thicknesses. Further, the FEMs incorporated a user-defined material to simulate the volume contraction of the TGO during the silica phase transformation from β-cristobalite to α-cristobalite upon cooling from the stress-free state at 1350°C to room temperature. Systems with and without a silicon bond coating intermediary layer were assessed. It was shown that the TGO phase transformation stress (1.6–1.7 GPa) dominated the increase in stress in the TGO and EBC layers. Furthermore, it was found that stress increase in the TGO was independent of TGO thickness and interface geometry. These results indicate that stabilization of the TGO to mitigate the phase transformation could dramatically improve the performance of SiC-base components with EBCs.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Finite-element nonlinear transient response computer programs PLATE 1 and CIVM-PLATE 1 for the analysis of panels subjected to impulse or impact loads

Two computer programs are described for predicting the transient large deflection elastic viscoplastic responses of thin single layer, initially flat unstiffened or integrally stiffened, Kirchhoff-Lov ductile metal panels. The PLATE 1 program pertains to structural responses produced by prescribed externally applied transient loading or prescribed initial velocity distributions. The collision imparted velocity method PLATE 1 program concerns structural responses produced by impact of an idealized nondeformable fragment. Finite elements are used to represent the structure in both programs. Strain hardening and strain rate effects of initially isotropic material are considered.

Spilker, R. L.↗
NASA NTRS2002

Finite Element Analysis of Wrinkled Membrane Structures for Sunshield Applications

The deployable sunshield is an example of a gossamer structure envisioned for use on future space telescopes. The basic structure consists of multiple layers of pretensioned, thin-film membranes supported by deployable booms. The prediction and verification of sunshield dynamics has been identified as an area in need of technology development due to the difficulties inherent in predicting nonlinear structural behavior of the membranes and because of the challenges involved. in ground testing of the full-scale structure. This paper describes a finite element analysis of a subscale sunshield that has been subjected to ground testing in support of the Next Generation Space Telescope (NGST) program. The analysis utilizes a nonlinear material model that accounts for wrinkling of the membranes. Results are presented from a nonlinear static preloading analysis and subsequent dynamics analyses to illustrate baseline sunshield structural characteristics. Studies are then described which provide further insight into the effect of membrane. preload on sunshield dynamics and the performance of different membrane modeling techniques. Lastly, a comparison of analytical predictions and ground test results is presented.

Johnston, John D.↗
NASA NTRS2007

Finite Element Development and Specifications of a Patched, Recessed Nomex Core Honeycomb Panel for Increased Sound Transmission Loss

This informal report summarizes the development and the design specifications of a recessed nomex core honeycomb panel in fulfillment of the deliverable in Task Order 13RBE, Revision 10, Subtask 17. The honeycomb panel, with 0.020-inch thick aluminum face sheets, has 0.016-inch thick aluminum patches applied to twenty-five, 6 by 6 inch, quarter inch thick recessed cores. A 10 dB higher transmission loss over the frequency range 250 - 1000 Hz was predicted by a MSC/NASTRAN finite element model when compared with the transmission loss of the base nomex core honeycomb panel. The static displacement, due to a unit force applied at either the core or recessed core area, was of the same order of magnitude as the static displacement of the base honeycomb panel when exposed to the same unit force. The mass of the new honeycomb design is 5.1% more than the base honeycomb panel. A physical model was constructed and is being tested.

Grosveld, Ferdinand W.↗
DOE OSTI2021

Deriving Transmissibility Functions from Finite Elements for Specifications

This work explores deriving transmissibility functions for a missile from a measured location at the base of the fairing to a desired location within the payload. A pressure on the outside of the fairing and the rocket motor’s excitation creates an acceleration at a measured location and a desired location. Typically, the desired location is not measured. In fact, it is typical that the payload may change, but measured acceleration at the base of the fairing is generally similar to previous test flights. Given this knowledge, it is desired to use a finite-element model to create a transmissibility function which relates acceleration from the previous test flight’s measured location at the base of the fairing to acceleration at a location in the new payload. Four methods are explored for deriving this transmissibility, with the goal of finding an appropriate transmissibility when both the pressure and rocket motor excitation are equally present. In conclusion, these methods are assessed using transient results from a simple example problem, and it is found that one of the methods gives good agreement with the transient results for the full range of loads considered.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗