SEARCH · Engineering Papers
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.
In the investigation discussed in this paper, the response of FEA cathodes and FEA material candidates to simulated electric propulsion system environments are presented.
This article presents three-dimensional B field solutions for the cylindrical Halbach array in an axial orientation. This arrangement has applications in the design of axial motors and passive axial magnetic bearings and couplers. The analytical model described here assumes ideal magnets with fixed and uniform magnetization. The field component functions are expressed as sums of 2-D definite integrals that are easily computed by a number of mathematical analysis software packages. The analysis is verified with sample calculations and the results are compared to equivalent results from traditional finite-element analysis (FEA). The field solutions are then approximated for use in flux linkage and induced EMF calculations in nearby stator windings by expressing the field variance with angular displacement as pure sinusoidal function whose amplitude depends on radial and axial position. The primary advantage of numerical implementation of the analytical approach presented in the article is that it lends itself more readily to parametric analysis and design tradeoffs than traditional FEA models.
Thermal analyses are performed on the liquid hydrogen (LH2) tank designed for an unmanned aerial vehicle (UAV) powered by solar arrays and a regenerative proton-exchange membrane (PEM) fuel cell. A 14-day cruise mission at a 65,000 ft altitude is considered. Thermal analysis provides the thermal loads on the tank system and the boiling-off rates of LH2. Different approaches are being considered to minimize the boiling-off rates of the LH2. It includes an evacuated multilayer insulation (MLI) versus aerogel insulation on the LH2 tank and aluminum versus stainless steel spacer rings between the inner and outer tank. The resulting boil-off rates of LH2 provided by the one-dimensional model and three-dimensional finite element analysis (FEA) on the tank system are presented and compared to validate the results of the three-dimensional FEA. It concludes that heat flux through penetrations by conduction is as significant as that through insulation around the tank. The tank system with MLI insulation and stainless steel spacer rings result in the lowest boiling-off rate of LH2.
A method is presented for formulating stiffness terms and thermal coefficients of fiber reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN™ FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.
This presentation discusses fractured leads on field-programmable gate array (FPGA) during flight vibration. Actions taken to determine root cause and resolution of the failure include finite element analysis (FEA) and vibration testing and scanning electron microscopy (with X-ray microanalysis) and energy dispersive spectrometry (SEM/EDS) failure assessment. Bonding methods for surface mount parts is assessed, including critical analysis and assessment of random fatigue damage. Regarding ceramic quad flat pack (CQFP) lead fracture, after disassembling the attitude control electronics (ACE) configuration, photographs showed six leads cracked on FPGA RTSX72SU-1 CQ208B package located on the RWIC card. An identical package (FPGA RTSX32SU-1 CQ208B) mounted on the RWIC did not results in cracked pins due to vibration. FPGA lead failure theories include workmanship issues in the lead-forming, material defect in the leads of the FPGA packages, and the insecure mounting of the board in the card guides, among other theories. Studies were conducted using simple calculations to determine the response and fatigue life of the package. Shorter packages exhibited more response when loaded by out-of-plane displacement of PCB while taller packages exhibit more response when loaded by in-plane acceleration of PCB. Additionally, under-fill did not contribute to reducing stress in leads due to out-of-plane PCB loading or from component twisting, as much as corner bonding. The combination of corner bond and under-fill is best to address mechanical and thermal S/C environment. Test results of bonded parts showed reduced (dampened) amplitude and slightly shifted peaks at the un-bonded natural frequency and an additional response at the bonded frequency. Stress due to PCBB out-of-plane loading was decreased on in the corners when only a corner bond was used. Future work may address CQFP fatigue assessment, including the investigation of discrepancy in predicted fatigue damage, as well as comparing fatigue life and fatigue damage cycle ration computed using FEA and Miner's rule to results from a fatigue assessment software program.
The zero-gravity surface figure of optics used in spaceborne astronomical instruments must be known to high accuracy, but earthbound metrology is typically corrupted by gravity sag. Generally, inference of the zero-gravity surface figure from a measurement made under normal gravity requires finite-element analysis (FEA), and for accurate results the mount forces must be well characterized. We describe how to infer the zero-gravity surface figure very precisely using the alternative classical technique of averaging pairs of measurements made with the direction of gravity reversed. We show that mount forces as well as gravity must be reversed between the two measurements and discuss how the St. Venant principle determines when a reversed mount force may be considered to be applied at the same place in the two orientations. Our approach requires no finite-element modeling and no detailed knowledge of mount forces other than the fact that they reverse and are applied at the same point in each orientation. If mount schemes are suitably chosen, zero-gravity optical surfaces may be inferred much more simply and more accurately than with FEA.
A mathematical model, termed the Slug Loss Model, has been developed for describing thermal capacitance (slug) calorimeter behavior when heat losses and other decaying processes are not negligible. This model results in the temperature time slope taking the mathematical form of exponential decay. When data is found to fit well to this model, it allows a heat flux value to be calculated that corrects for the losses and may be a better estimate of the cold wall fully catalytic heat flux, as is desired in arc jet testing. The model was applied to the data from a copper slug calorimeter inserted during a particularly severe high heating rate arc jet run to illustrate its use. The Slug Loss Model gave a cold wall heat flux 15% higher than the value of 2,250 W/sq cm obtained from the conventional approach to processing the data (where no correction is made for losses). For comparison, a Finite Element Analysis (FEA) model was created and applied to the same data, where conduction heat losses from the slug were simulated. The heat flux determined by the FEA model was found to be in close agreement with the heat flux determined by the Slug Loss Model.
A unified framework is presented that enables coupled multiscale analysis of composite structures and associated graphical pre- and postprocessing within the Abaqus/CAE environment. The recently developed, free, Finite Element Analysis--Micromechanics Analysis Code (FEAMAC) software couples NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) with Abaqus/Standard and Abaqus/Explicit to perform micromechanics based FEA such that the nonlinear composite material response at each integration point is modeled at each increment by MAC/GMC. The Graphical User Interfaces (FEAMAC-Pre and FEAMAC-Post), developed through collaboration between SIMULIA Erie and the NASA Glenn Research Center, enable users to employ a new FEAMAC module within Abaqus/CAE that provides access to the composite microscale. FEA IAC-Pre is used to define and store constituent material properties, set-up and store composite repeating unit cells, and assign composite materials as sections with all data being stored within the CAE database. Likewise FEAMAC-Post enables multiscale field quantity visualization (contour plots, X-Y plots), with point and click access to the microscale i.e., fiber and matrix fields).
Three means (i.e., typical foil strain gages, fiber optic sensors, and a digital image correlation (DIC) system) were implemented to measure strains on the back and front surfaces of a longitudinally jointed curved test article subjected to edge-wise compression testing, at NASA Goddard Space Flight Center, according to ASTM C364. The Pre-test finite element analysis (FEA) was conducted to assess ultimate failure load and predict strain distribution pattern throughout the test coupon. The predicted strain pattern contours were then utilized as guidelines for installing the strain measurement instrumentations. The strain gages and fiber optic sensors were bonded on the specimen at locations with nearly the same strain values, as close as possible to each other, so that, comparisons between the measured strains by strain gages and fiber optic sensors, as well as the DIC system are justified. The test article was loaded to failure (at approximately 38 kips), at the strain value of approximately 10,000mu epsilon As a part of this study, the validity of the measured strains by fiber optic sensors is examined against the strain gage and DIC data, and also will be compared with FEA predictions.
Three distinct strain measurement methods (i.e., foil resistance strain gages, fiber optic strain sensors, and a three-dimensional digital image photogrammetry that gives full field strain and displacement measurements) were implemented to measure strains on the back and front surfaces of a longitudinally jointed curved test article subjected to edge-wise compression testing, at NASA Goddard Space Flight Center, according to ASTM C364. The pre-test finite element analysis (FEA) was conducted to assess ultimate failure load and predict strain distribution pattern throughout the test coupon. The predicted strain pattern contours were then utilized as guidelines for installing the strain measurement instrumentations. The foil resistance strain gages and fiber optic strain sensors were bonded on the specimen at locations with nearly the same analytically predicted strain values, and as close as possible to each other, so that, comparisons between the measured strains by strain gages and fiber optic sensors, as well as the three-dimensional digital image photogrammetric system are relevant. The test article was loaded to failure (at ~167 kN), at the compressive strain value of ~10,000 micro epsilon. As a part of this study, the validity of the measured strains by fiber optic sensors is examined against the foil resistance strain gages and the three-dimensional digital image photogrammetric data, and comprehensive comparisons are made with FEA predictions.
A lightweight bulldozer blade prototype has been designed and built to be used as an excavation implement in conjunction with the NASA Chariot lunar mobility platform prototype. The combined system was then used in a variety of field tests in order to characterize structural loads, excavation performance and learn about the operational behavior of lunar excavation in geotechnical lunar simulants. The purpose of this effort was to evaluate the feasibility of lunar excavation for site preparation at a planned NASA lunar outpost. Once the feasibility has been determined then the technology will become available as a candidate element in the NASA Lunar Surface Systems Architecture. In addition to NASA experimental testing of the LANCE blade, NASA engineers completed analytical work on the expected draft forces using classical soil mechanics methods. The Colorado School of Mines (CSM) team utilized finite element analysis (FEA) to study the interaction between the cutting edge of the LANCE blade and the surface of soil. FEA was also used to examine various load cases and their effect on the lightweight structure of the LANCE blade. Overall it has been determined that a lunar bulldozer blade is a viable technology for lunar outpost site preparation, but further work is required to characterize the behavior in 1/6th G and actual lunar regolith in a vacuum lunar environment.
A method is presented for formulating stiffness terms and thermal coefficients of fiber-reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN(TradeMark) FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.
The Stratospheric Aerosol and Gas Experiment (SAGE) III - International Space Station (ISS) instrument will be used to study ozone, providing global, long-term measurements of key components of the Earth's atmosphere for the continued health of Earth and its inhabitants. SAGE III is launched into orbit in an inverted configuration on SpaceX;s Falcon 9 launch vehicle. As one of its four supporting elements, a Contamination Monitoring Package (CMP) mounted to the top panel of the Interface Adapter Module (IAM) box experiences high-frequency response due to structural coupling between the two structures during the SpaceX launch. These vibrations, which were initially observed in the IAM Engineering Development Unit (EDU) test and later verified through finite element analysis (FEA) for the SpaceX launch loads, may damage the internal electronic cards and the Thermoelectric Quartz Crystal Microbalance (TQCM) sensors mounted on the CMP. Three-dimensional (3D) vibration isolators were required to be inserted between the CMP and IAM interface in order to attenuate the high frequency vibrations without resulting in any major changes to the existing system. Wire rope isolators were proposed as the isolation system between the CMP and IAM due to the low impact to design. Most 3D isolation systems are designed for compression and roll, therefore little dynamic data was available for using wire rope isolators in an inverted or tension configuration. From the isolator FEA and test results, it is shown that by using the 3D wire rope isolators, the CMP high-frequency responses have been suppressed by several orders of magnitude over a wide excitation frequency range. Consequently, the TQCM sensor responses are well below their qualification environments. It is indicated that these high-frequency responses due to the typical instrument structural coupling can be significantly suppressed by a vibration passive control using the 3D vibration isolator. Thermal and contamination issues were also examined during the isolator selection period for meeting the SAGE III-ISS instrument requirements.
Stochastic-based, discrete-event progressive damage simulations of ceramic-matrix composite and polymer matrix composite material structures have been enabled through the development of a unique multiscale modeling tool. This effort involves coupling three independently developed software programs: (1) the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), (2) the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program (CARES/ Life), and (3) the Abaqus finite element analysis (FEA) program. MAC/GMC contributes multiscale modeling capabilities and micromechanics relations to determine stresses and deformations at the microscale of the composite material repeating unit cell (RUC). CARES/Life contributes statistical multiaxial failure criteria that can be applied to the individual brittle-material constituents of the RUC. Abaqus is used at the global scale to model the overall composite structure. An Abaqus user-defined material (UMAT) interface, referred to here as "FEAMAC/CARES," was developed that enables MAC/GMC and CARES/Life to operate seamlessly with the Abaqus FEA code. For each FEAMAC/CARES simulation trial, the stochastic nature of brittle material strength results in random, discrete damage events, which incrementally progress and lead to ultimate structural failure. This report describes the FEAMAC/CARES methodology and discusses examples that illustrate the performance of the tool. A comprehensive example problem, simulating the progressive damage of laminated ceramic matrix composites under various off-axis loading conditions and including a double notched tensile specimen geometry, is described in a separate report.
Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.
In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.
An analysis framework is presented that makes available multiscale analysis of composite structures using the open-source FEA solver package CalculiX CrunchiX (CCX). At the center of this framework is the coupling of the Finite Element Analysis - Micromechanics Analysis Code (FEAMAC) library from NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) to allow micromechanics analysis. The results show that the proposed coupling can be used with appropriate care for multiscale FEA simulations of composite materials. The largest error reported in this validation was in a four-point bend test specimen with an error of less than 1% difference in the maximum deflection of the beam.
The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.