Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elemental analysis”

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 505 records · Page 28

Improved accuracy for finite element structural analysis via a new integrated force method

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Patnaik, Surya N.↗

Applications of FEM and BEM in two-dimensional fracture mechanics problems

A comparison of the finite element method (FEM) and boundary element method (BEM) for the solution of two-dimensional plane strain problems in fracture mechanics is presented in this paper. Stress intensity factors (SIF's) were calculated using both methods for elastic plates with either a single-edge crack or an inclined-edge crack. In particular, two currently available programs, ANSYS for finite element analysis and BEASY for boundary element analysis, were used.

Min, J. B.↗

Coating Development for GRCop-84 Liners for Reusable Launch Vehicles Aided by Modeling Studies

The design of the next generation of reusable launch vehicles calls for using GRCop-84 copper alloy liners based on a composition invented at the NASA Glenn Research Center. Despite its considerable advantage over other copper alloys, it is expected that GRCop-84 will suffer from environmental degradation depending on the type of rocket fuels used and on thermomechanical fatigue. Applying protective coatings on GRCop-84 substrates can minimize or eliminate many of these problems and extend the operational life of the combustion liner. This could increase component reliability, shorten depot maintenance turnaround times, and lower operating costs. Therefore, Glenn is actively pursuing the development of advanced coatings technology for GRCop-84 liners. Technology is being developed in four major areas: (1) new metallic coating compositions, (2) application techniques, (3) test methods, and (4) life prediction design methodology using finite element analysis. The role of finite element analysis in guiding the coating effort is discussed in this report. Thermal analyses were performed at Glenn for different combinations of top- and bondcoat compositions to determine the temperature variation across the coated cross section with the thickness of the top coat. These calculations were conducted for simulated LH2/LO2 booster engine conditions assuming that the bond coat had a constant thickness of 50 m. The preceding graphs show the predicted temperatures at the outer surface of the top coat (hot wall), at the top-coat/bond-coat interface, at the bond-coat/GRCop-84 interface, and at the GRCop-84 cold wall as a function of top-coat thickness for Cu- 26(wt%)Cr top coat (top graph), Ni-17(wt%)Cr-6%Al-0.5%Y top coat and Cu-26%Cr bond coat, and NiAl top coat and Ni bond coat. In all cases, the temperature of the top coat at the hot wall increased with increasing top-coat thickness and with corresponding decreases in the temperatures at the two interfaces and the cold wall. These temperatures are not acutely sensitive to the thermal conductivity of the top coat when it exceeds 25 and 50 W/m/K for low and high heat flux engines. This observation is significant for two reasons. First, several different top-coat compositions can be evaluated as potential protective coatings without loss in the heat-transfer efficiency of the coated system. Second, materials with thermal conductivities less than the critical values of 25 or 50 W/m/K are more likely to act as thermal barrier coatings. The deposition of overlay coatings on GRCop-84 substrates results in the development of residual stresses. The presence of these residual stresses influences the probability of coating spallation, the thermal cycling life, and the fatigue properties of the coated substrate during use. Since it is important to understand how these stresses develop during the vacuum-plasma-spraying coating deposition process, the nature and magnitudes of the cool-down residual stresses were calculated and compared with experimentally determined values across the coated cross section of a disk specimen. The calculations were conducted assuming that the specimen cools down to room temperature from vacuum plasma-spraying temperatures of either 250 or 650 C. The effects of coating the substrate with and without grit blasting were also theoretically examined. The final graph compares the predicted and the experimental results for a GRCop-84 disk coated with about a 50- m-thick Ni bond coat and a 75- to 100- m NiAl top coat, where the curves for NASA-2 assume the presence of a prior residual stress generated by grit blasting under conditions similar to the experimental situation. The predicted cool-down in-plane stresses were compressive in both the NiAl top coat and the Ni bond coat. They were also compressive in the substrate to a depth of about 0.25 mm from the Ni/GRCop-84 interface when the vacuum-plasma-spraying temperature was low. However, using a higher plasma spraying temperaturs likely to leave the substrate under a small tensile stress to counter the compressive stresses in the bond and top coats because of the relaxation of residual stresses generated in the substrate during the grit blasting of its surface prior to spraying. These results suggest that the NiAl and Ni coatings are unlikely to spall after spraying as confirmed by the microstructural observations shown in the following photomicrograph of an as-sprayed specimen. Finally, it is noted that the calculated and experimental results are not in complete agreement, which indicates that both the experimental and modeling techniques need further refinement.

Raj, Sai V.↗

A three-node C deg element for analysis of laminated composite sandwich shells

A three-node flat shell element with C deg rotation fields has been developed for analysis of arbitrary composite shells. The element may consist of any number of orthotropic layers, each layer having different material properties and angular orientation. The formulation includes coupling between bending and extension, which is essential for analysis of unsymmetric laminates. Shearing deflections are included, since laminated and sandwich construction frequently results in shear stiffness much smaller than bending stiffness. Formulation of the element is straightforward, and calculation of its stiffness matrix is simple and fast. Convergence of solutions with mesh refinement is uniform for both thin and thick shells and is insensitive to element shape, although not as rapid as some other elements that lack one or more capabilities of the newly developed element. An experimental verification of the shall element is reported in the appendix.

Martin, C. Wayne↗

A combined experimental and finite element study to predict the failure mechanisms in SiC coated carbon/carbon composites at room and elevated temperatures under flexural loading

Response of quasi-isotropic laminates of SiC coated Carbon/Carbon (C/C) composites have been investigated under flexural loading at various temperatures. Variation of load-deflection behavior with temperatures are studied. Increase in flexural strength and stiffness are observed with the rise in temperature. Extensive analyses through Optical Microscope (OM) and Non-Destructive Evaluation (NDE) have been performed to understand the failure mechanisms. Damage zone is found only within the neighborhood of the loading plane. Isoparametric layered shell elements developed on the basis of the first order shear deformation theory have been used to model the thin laminates of C/C under flexural loading. Large deformation behavior has been considered in the finite element analysis to account for the non-linearities encountered during the actual test. Data generated using finite element analysis are presented to corroborate the experimental findings, and a comparison in respect of displacement and stress-strain behavior are given to check the accuracy of the finite element analysis. Reasonable correlation between the experimental and finite element results have been established.

Mahfuz, Hassan↗

Structural Dynamics and Data Analysis

This project consists of two parts, the first will be the post-flight analysis of data from a Delta IV launch vehicle, and the second will be a Finite Element Analysis of a CubeSat. Shock and vibration data was collected on WGS-5 (Wideband Global SATCOM- 5) which was launched on a Delta IV launch vehicle. Using CAM (CAlculation with Matrices) software, the data is to be plotted into Time History, Shock Response Spectrum, and SPL (Sound Pressure Level) curves. In this format the data is to be reviewed and compared to flight instrumentation data from previous flights of the same launch vehicle. This is done to ensure the current mission environments, such as shock, random vibration, and acoustics, are not out of family with existing flight experience. In family means the peaks on the SRS curve for WGS-5 are similar to the peaks from the previous flights and there are no major outliers. The curves from the data will then be compiled into a useful format so that is can be peer reviewed then presented before an engineering review board if required. Also, the reviewed data will be uploaded to the Engineering Review Board Information System (ERBIS) to archive. The second part of this project is conducting Finite Element Analysis of a CubeSat. In 2010, Merritt Island High School partnered with NASA to design, build and launch a CubeSat. The team is now called StangSat in honor of their mascot, the mustang. Over the past few years, the StangSat team has built a satellite and has now been manifested for flight on a SpaceX Falcon 9 launch in 2014. To prepare for the final launch, a test flight was conducted in Mojave, California. StangSat was launched on a Prospector 18D, a high altitude rocket made by Garvey Spacecraft Corporation, along with their sister satellite CP9 built by California Polytechnic University. However, StangSat was damaged during an off nominal landing and this project will give beneficial insights into what loads the CubeSat experienced during the crash. During the year, the MIHS students generated a SolidWorks (CAD software) geometry model of StangSat. This model will be imported into FEMAP (Finite Element Analysis (FEA) Software) and a finite element model wiiJ be created to predict the loads encountered during the crash of this rocket. This analysis will require learning how to import CAD models into the FEM, mesh and add constraints and concentrated masses to represent components inside the CubeSat frame, such as circuit boards, batteries and accelerometers. During the analysis the loads will be varied, in effort to duplicate the damage to the CubeSat. Results will then be peer reviewed and documented.

Luthman, Briana L.↗

Finite-Element Fracture Analysis of Pins and Bolts

Stress intensities calculated in bending and tension. Finite-element stress-analysis method gives stress-intensity estimates for surface flaws on smooth and threaded round bars. Calculations done for purely tensile and purely bending loads. Results, presented in dimensionless form, useful for determining fatigue lives of bolts and pins.

Nord, K. J.↗

Numerical Analysis Of Interlaminar-Fracture Toughness

Finite-element analysis applied in conjunction with strain-energy and micromechanical concepts. Computational procedure involves local, local-crack-closure, and/or the "unique" local-crack-closure method developed at NASA Lewis Research Center, for mathematical modeling of ENF and MMF. Methods based on three-dimensional finite-element analysis in conjunction with concept of strain-energy-release rate and with micromechanics of composite materials. Assists in interpretation of ENF and MMF fracture tests performed to obtain fracture-toughness parameters, by enabling evaluation of states of stress likely to induce interlaminar fractures.

Chamis, C. C.↗

A study of the response of nonlinear springs

The various phases to developing a methodology for studying the response of a spring-reinforced arch subjected to a point load are discussed. The arch is simply supported at its ends with both the spring and the point load assumed to be at midspan. The spring is present to off-set the typical snap through behavior normally associated with arches, and to provide a structure that responds with constant resistance over a finite displacement. The various phases discussed consist of the following: (1) development of the closed-form solution for the shallow arch case; (2) development of a finite difference analysis to study (shallow) arches; and (3) development of a finite element analysis for studying more general shallow and nonshallow arches. The two numerical analyses rely on a continuation scheme to move the solution past limit points, and to move onto bifurcated paths, both characteristics being common to the arch problem. An eigenvalue method is used for a continuation scheme. The finite difference analysis is based on a mixed formulation (force and displacement variables) of the governing equations. The governing equations for the mixed formulation are in first order form, making the finite difference implementation convenient. However, the mixed formulation is not well-suited for the eigenvalue continuation scheme. This provided the motivation for the displacement based finite element analysis. Both the finite difference and the finite element analyses are compared with the closed form shallow arch solution. Agreement is excellent, except for the potential problems with the finite difference analysis and the continuation scheme. Agreement between the finite element analysis and another investigator's numerical analysis for deep arches is also good.

Hyer, M. W.↗

Application of finite element substructuring to composite micromechanics

Finite element substructuring is used to predict unidirectional fiber composite hygral (moisture), thermal, and mechanical properties. COSMIC NASTRAN and MSC/NASTRAN are used to perform the finite element analysis. The results obtained from the finite element model are compared with those obtained from the simplified composite micromechanics equations. A unidirectional composite structure made of boron/HM-epoxy, S-glass/IMHS-epoxy and AS/IMHS-epoxy are studied. The finite element analysis is performed using three dimensional isoparametric brick elements and two distinct models. The first model consists of a single cell (one fiber surrounded by matrix) to form a square. The second model uses the single cell and substructuring to form a nine cell square array. To compare computer time and results with the nine cell superelement model, another nine cell model is constructed using conventional mesh generation techniques. An independent computer program consisting of the simplified micromechanics equation is developed to predict the hygral, thermal, and mechanical properties for this comparison. The results indicate that advanced techniques can be used advantageously for fiber composite micromechanics.

Caruso, J. J.↗

Comparison of elastic and elastic-plastic structural analyses for cooled turbine blade airfoils

Elastic plastic stress strain states in cooled turbine blade airfoils were calculated by three methods for the initial takeoff transient of an advanced technology aircraft engine. The three analytical methods compared were a three dimensional elastic plastic, finite element analysis, a three dimensional, elastic, finite element analysis, and a one dimensional, elastic plastic, beam theory analysis. Structural analyses were performed for eight cases involving different combinations of mechanical and thermal loading on impingement cooled airfoils with and without leading edge film cooling holes. The von Mises effective total strains at maximum takeoff computed from the elastic and elastic plastic finite element analyses agreed with 9 percent for rotating airfoils and 28 percent for stationary airfoils with the elastic results on the conservative side.

Kaufman, A.↗

Compression creep of filamentary composites

Axial and transverse strain fields induced in composite laminates subjected to compressive creep loading were compared for several types of laminate layups. Unidirectional graphite/epoxy as well as multi-directional graphite/epoxy and graphite/PEEK layups were studied. Specimens with and without holes were tested. The specimens were subjected to compressive creep loading for a 10-hour period. In-plane displacements were measured using moire interferometry. A computer based data reduction scheme was developed which reduces the whole-field displacement fields obtained using moire to whole-field strain contour maps. Only slight viscoelastic response was observed in matrix-dominated laminates, except for one test in which catastrophic specimen failure occurred after a 16-hour period. In this case the specimen response was a complex combination of both viscoelastic and fracture mechanisms. No viscoelastic effects were observed for fiber-dominated laminates over the 10-hour creep time used. The experimental results for specimens with holes were compared with results obtained using a finite-element analysis. The comparison between experiment and theory was generally good. Overall strain distributions were very well predicted. The finite element analysis typically predicted slightly higher strain values at the edge of the hole, and slightly lower strain values at positions removed from the hole, than were observed experimentally. It is hypothesized that these discrepancies are due to nonlinear material behavior at the hole edge, which were not accounted for during the finite-element analysis.

Graesser, D. L.↗

Response of Composite Plates with Inclined Elliptical Notches and Subjected to Axial Compression

An analysis method for predicting the inplane stress states in anisotropic finite plates with an elliptical notch is presented. This method can be used to analyze plates with arbitrary notch orientations with respect to the plate material axes. The analysis results have been validated using finite element analysis results for unnotched composite plates and experimental and finite element analysis results for stiffened composite panels with a skin that has orthotropic properties. The good agreement between these results, until the panel exhibits nonlinear response either due to bending or initiation of damage, indicates that the present analysis method can be used to determine accurately the inplane stress states and stress concentrations at and around an elliptical notch.

Ambur, Damodar R.↗

Delaminations in composite plates under impact loads

A method is presented for calculating the locations, shapes, and sizes of delaminations which occur in a fiber reinforced composite plate subjected to non-penetrating (low velocity) impact of a solid object. The plate may be simply supported, clamped, or free along its edges. A failure model of the delamination formation was developed. This model was then coupled with a finite element analysis. The model and the finite element analysis were then implemented by a computer code (IMPACT-ST) which can be used to estimate the damage initiation load and the locations, shapes, and sizes of the delaminations. Tests were performed measuring the geometries of the delaminations in graphite-epoxy, graphite-toughened epoxy, and graphite-PEEK plates impacted by a projectile with a spherical tip having masses ranging from 0.355 lbm to 0.963 lbm and velocities from 50 in/sec to 225 in/sec. The data were compared to the results of the model, and good agreements were found between the measured and the calculated delamination lengths and widths.

Finn, Scott R.↗

Delaminations in composite plates under transverse static or impact loads - A model

A method is presented for calculating the locations, shapes, and sizes of delaminations which occur in a fiber reinforced composite plate subjected to transverse static or dynamic (impact) loads. The plate may be simply supported, clamped, or free along its edges. A model of the delamination formation was developed. This model was then coupled with a finite element analysis. The model and the finite element analysis were implemented by a computer code which can be used to estimate the load at which damage is initiated as well as the locations, shapes, and sizes of the delaminations.

Finn, Scott R.↗

Application of the p-version of the finite-element method to global-local problems

A brief survey is given of some recent developments in finite-element analysis technology which bear upon the three main research areas under consideration in this workshop: (1) analysis methods; (2) software testing and quality assurance; and (3) parallel processing. The variational principle incorporated in a finite-element computer program, together with a particular set of input data, determines the exact solution corresponding to that input data. Most finite-element analysis computer programs are based on the principle of virtual work. In the following, researchers consider only programs based on the principle of virtual work and denote the exact displacement vector field corresponding to some specific set of input data by vector u(EX). The exact solution vector u(EX) is independent of the design of the mesh or the choice of elements. Except for very simple problems, or specially constructed test problems, vector u(EX) is not known. Researchers perform a finite-element analysis (or any other numerical analysis) because they wish to make conclusions concerning the response of a physical system to certain imposed conditions, as if vector u(EX) were known.

Szabo, Barna A.↗

Stress analysis of the cracked lap shear specimens: An ASTM round robin

This ASTM Round Robin was conducted to evaluate the state of the art in stress analysis of adhesively bonded joint specimens. Specifically, the participants were asked to calculate the strain-energy-release rate for two different geometry cracked lap shear (CLS) specimens at four different debond lengths. The various analytical techniques consisted of 2- and 3-dimensional finite element analysis, beam theory, plate theory, and a combination of beam theory and finite element analysis. The results were examined in terms of the total strain-energy-release rate and the mode I to mode II ratio as a function of debond length for each specimen geometry. These results basically clustered into two groups: geometric linear or geometric nonlinear analysis. The geometric nonlinear analysis is required to properly analyze the CLS specimens. The 3-D finite element analysis gave indications of edge closure plus some mode III loading. Each participant described his analytical technique and results. Nine laboratories participated.

Johnson, W. S.↗