Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-element”

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 901 records · Page 50

Methodology for Predicting the Onset of Widespread Fatigue Damage in Lap-Splice Joints

NASA has conducted an Airframe Structural Integrity Program to develop the methodology to predict the onset of widespread fatigue damage to lap-splice joints of fuselage structures. Several stress analysis codes have been developed or enhanced to analyze the lap-splice-joint configuration. Fatigue lives in lap-splice-joint specimens and fatigue-crack growth in a structural fatigue test article agreed well with calculations from small-crack theory and fatigue-crack growth analyses with the FASTRAN code. Residual-strength analyses of laboratory specimens and wide stiffened panels were predicted quite well from the critical crack-tip-opening angle (CTOA) fracture criterion and elastic-plastic finite-element analyses (two- or three-dimensional codes and the STAGS shell code).

Newman, J. C., Jr.↗

Design, Operation, and Modeling of a Vertical APCVD Reactor for Silicon Carbide Film Growth

An atmospheric pressure chemical vapor deposition (APCVD) reactor utilizing a unique vertical geometry which enables 3C-SiC films to be grown on two, 4-inch diameter Si wafers has been constructed. Contrary to expectations, 3C-SiC films grown in this reactor are thickest at the downstream end of the substrates. To better understand the reason for the thickness distribution on the wafers, an axisymmetric finite-element model of the gas flow in the reactor was constructed. The model uses the ANSYS53 Flowtran package and includes compressible and temperature-dependent fluid properties in laminar or turbulent flow. It does not include reaction chemistry or unsteady flow. The ANSYS53 results predict that the cool, inlet fluid falls through the inlet pipe and the warm, diffuser region like a jet. This jet impinges on top of the susceptor and gets diverted to the reactor side walls, where it flows to the bottom of the reactor, turns, and slowly rises along the face of the susceptor. This may explain why the SiC films are thickest at the downstream side of the wafers, as gas containing fresh reactants first passes over this region. Modeling results are presented for both one atmosphere and one half atmosphere reactor pressure.

DeAnna, Russell G.↗

The Effects of Variable Mass and Geometry, Pretwist, Shear Deformation and Rotatory Inertia on the Resonant Frequencies of Intact Long Bones: A Finite Element Model Analysis

The influence of pretwist, nonuniformities in mass and flexural stiffness, rotatory inertia and shear deformation on the natural frequencies of intact bones is evaluated by means of a linear elastic, finite-element model which has been programmed for solution on the digital computer. Theoretical results are compared to the results on the forced vibration of intact canine radii obtained experimentally by Thompson. Surprisingly, inclusion of fairly large pretwist angles (from -14 to 12 deg for one specimen) had little affect on the first three frequencies of transverse vibration in either the cranial or lateral directions. Inclusion of shear deformation reduced the third-mode frequency in the stiffest (lateral) direction by about six percent, otherwise shear deformation played a minor role in determining natural frequencies. Similarly. rotatory inertia had negligible influence up to the third natural frequency. The predominant influence on the first three natural frequencies of transverse vibration could be attributed to the variations in mass and flexural stiffness along the length of the test specimens. Different effective moduli of elasticity are required to yield correct absolute values for the frequencies which correspond to experimental findings. thus implying the presence of some inhomogeneities in material properties around the bone cross-section and/or along its length.

Young, Donald R.↗

Open-Mode Debonding Analysis of Curved Sandwich Panels Subjected to Heating and Cryogenic Cooling on Opposite Faces

Increasing use of curved sandwich panels as aerospace structure components makes it vital to fully understand their thermostructural behavior and identify key factors affecting the open-mode debonding failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich panels with different radii of curvature. The curved sandwich panels are either simply supported or clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the concave side. The finite-element method was used to study the effects of panel curvature and boundary condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich panels. The critical stress point, where potential debonding failure could initiate, was found to be at the midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and the associated sandwich core depth stretching.

Ko, William L.↗

A Very High Resolution Mid Infrared Spectrograph for SOFIA

The purpose of this project was to study the design of a high resolution mid-infrared spectrograph for the SOFIA airborne observatory. The primary motivation for the construction and use of such an instrument is the study of vibrational transitions of molecules in molecular clouds, circumstellar inflows and outflows, and planetary atmospheres, which are blocked from ground-based observations by absorption by the Earth's atmosphere. The instrument design studied is a cross-dispersed grating spectrograph with several modes of operation giving spectral resolving powers of approximately 4000, 20,000, and 100,000. In the two lower resolution modes a single grating is used in a long-slit spectrograph. At high resolution that grating serves as the cross-disperser for a 1-m long echelon, a very coarsely and steeply ruled grating, diamond-machined in aluminum. Several design study tasks were carried out: 1. An optical and mechanical design was developed. It differed from the preliminary design by using a reflecting focal reducer in the fore-optics and by rearranging the echelon and cross-dispersing chambers to improve the efficiency of the long-slit modes. 2. A finite-element analysis of the dewar and echelon support structures was done. It led us to change the echelon support to improve its rigidity. 3. Further information was obtained which convinced us that Hyperfine Inc. was capable of diamond-machining the echelon. 4. Several possible vendors of readout electronics and software were identified. It was decided that it is more cost-effective to purchase these items rather than developing them ourselves. 5. Observing and calibration procedures were studied. we did not identify significant changes needed from procedures used previously. 6. Possible expansions of instrument capabilities, including wavelength coverage, were studied. It appears unlikely that the echelon will perform well at wavelengths shortward of 5 microns, but the spectral coverage could be extended longward of 28 microns with the use of an Si:Sb detector. No patentable inventions were made during this study. Based on this design study, a proposal was submitted to construct the instrument. It was chosen as one of the first-light instruments for SOFIA.

Lacy, John H.↗

Ice velocity at the ice front of the Filchner-Ronne Ice Shelf, Antarctica, as observed with ERS interferometry

ERS images of the two ends of the ice front of the Filchner-Ronne ice shelf (Antarctica) were utilized interferometrically to study the deformation rate of the ice shelf in response to viscous creep. On the western bank of Berkner Island (BI), near Hemmen Ice Rise (HIR), a time series of ERS data acquired in February 1992, in both ascending and descending mode, and with a three-day time interval, were utilized to map the ice velocity in two dimensions. Finite-element ice-shelf flow simulations are compared with the ERS interferograms to interpret the ice motion in terms of the physical constraints on ice-shelf flow. The efforts to fit artificial interferograms generated with model velocity output suggest that the flow regime is strongly influenced by three processes. First, a void-creation process responsible for rifts at coastal margins tends to uncouple the ice shelf from the ice rise and neighboring coast of BI. Secondly, sea ice within the void space appears to act as a binding agent between discrete ice-shelf fragments, allowing rigid-body rotations of these fragments. Third, strain rates appear to be enhanced in a narrow zone adjacent to HIR, implying significant strain softening along the boundary. It is believed that synthetic aperture radar (SAR) interferogram/model intercomparison represents a powerful impetus towards the development of better, more physically realistic ice-shelf flow models.

Rignot, Eric↗

A New Approach to Overcoming Spatial Aliasing in Structural Damage Detection

Aircraft, reusable launch vehicles,unmanned aircraft, and other advanced structures are being built using lightweight composite materials/metals with design safety factors as low as 1.25. These advanced structures operate in uncertain and severe environments and are susceptible to damage such as delamination, fiber/matrix damage, hydrothermal strain in composite materials, and fatigue and cracking in metals. To ensure human safety and load-bearing integrity these structures must be inspected to detect and locate often invisible damage and faults before they become catastrophic. Conventional methods of non-destructive evaluation sometimes miss significant damage and are time consuming and expensive to perform. In contrast, vibrometry or vibration signature techniques are a global method of structural integrity monitoring that potentially can efficiently detect damage on large structures, including damage that is away from sensor locations, and in the interior of structures. However, a barrier problem in damage detection using vibration measurements is the need to measure the vibration response at a large number of points on the structure. Typically, model reduction or expansion procedures such as Guyan reduction or dynamic expansion are attempted to overcome the problem of insufficient measurements or spatial aliasing. These approaches depend on using information from the healthy model, and thus put error into the reduction used to represent the damaged structure. In this paper, a Frequency Response Function technique is used to detect damage to a fixed-free beam. The technique uses measured frequency response functions from the healthy structure as reference data, and then monitors vibration measurements during the life of the structure to detect damage. In an analytical simulation using a finite-element model of a beam, damage was located using only sparse measurements because the technique uses both rotation and translation measurements from the damaged structure.

Schulz, Mark J.↗

Wing-Body Aeroelasticity Using Finite-Difference Fluid/Finite-Element Structural Equations on Parallel Computers

In recent years significant advances have been made for parallel computers in both hardware and software. Now parallel computers have become viable tools in computational mechanics. Many application codes developed on conventional computers have been modified to benefit from parallel computers. Significant speedups in some areas have been achieved by parallel computations. For single-discipline use of both fluid dynamics and structural dynamics, computations have been made on wing-body configurations using parallel computers. However, only a limited amount of work has been completed in combining these two disciplines for multidisciplinary applications. The prime reason is the increased level of complication associated with a multidisciplinary approach. In this work, procedures to compute aeroelasticity on parallel computers using direct coupling of fluid and structural equations will be investigated for wing-body configurations. The parallel computer selected for computations is an Intel iPSC/860 computer which is a distributed-memory, multiple-instruction, multiple data (MIMD) computer with 128 processors. In this study, the computational efficiency issues of parallel integration of both fluid and structural equations will be investigated in detail. The fluid and structural domains will be modeled using finite-difference and finite-element approaches, respectively. Results from the parallel computer will be compared with those from the conventional computers using a single processor. This study will provide an efficient computational tool for the aeroelastic analysis of wing-body structures on MIMD type parallel computers.

Byun, Chansup↗

Progress in Integrated Fluid/Structure/Controls Computations Using High Fidelity Equations

Aeroelasticity that involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of aircraft which can experience complex flow/structure interactions. Supersonic transports can experience vortex induced aeroelastic oscillations whereas subsonic transports can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Using these high fidelity methods, design quantities such as structural stresses can be directly computed. Computations using these high fidelity equations require state-of-the-art large computational resources. This paper will describe the transition from potential equations to Euler/Navier-Stokes equations for fluids and from modal equations to finite element equations for structures. Complexities associated with grids for moving control surfaces will be discussed. All contents of the paper will be limited to those already presented at domestic conferences within U.S. such as those sponsored by AIAA. It will not include any materials related to HPCCP computational work.

Guruswamy, Guru P.↗

Thermal Stress Analysis of RCG-Tempered TUFI Tile TPS for Hypersonic Vehicles

This paper presents detailed results from linear and nonlinear finite-element thermal stress analyses of a new tile, Thermal Protection Systems (TPS) concept. A very thin coating of Reaction Cured Glass (RCG) is used to "temper" the surface of Toughened Uni-Piece Fibrous Insulation (TUFI) tiles to improve resistance to thermal shock and thermal cycling effects. The coating also serves to reduce catalytic heating and may improve waterproofing. Calculations include trajectory-based aerothermal heating environments for X-34 wing leading edge TPS designs and arc jet environments for TPS test articles. The nonlinear analyses include the high temperature plasticity of RCG to demonstrate the reuseability of the material.

Milos, Frank S.↗

Numerical Methods Using B-Splines

The seminar will discuss (1) The current range of applications for which B-spline schemes may be appropriate (2) The property of high-resolution and the relationship between B-spline and compact schemes (3) Comparison between finite-element, Hermite finite element and B-spline schemes (4) Mesh embedding using B-splines (5) A method for the incompressible Navier-Stokes equations in curvilinear coordinates using divergence-free expansions.

Shariff, Karim↗

Design-Oriented Analysis of Aircraft Fuselage Structures

A design-oriented analysis capability for aircraft fuselage structures that utilizes equivalent plate methodology is described. This new capability is implemented as an addition to the existing wing analysis procedure in the Equivalent Laminated Plate Solution (ELAPS) computer code. The wing and fuselage analyses are combined to model entire airframes. The paper focuses on the fuselage model definition, the associated analytical formulation and the approach used to couple the wing and fuselage analyses. The modeling approach used to minimize the amount of preparation of input data by the user and to facilitate the making of design changes is described. The fuselage analysis is based on ring and shell equations but the procedure is formulated to be analogous to that used for plates in order to take advantage of the existing code in ELAPS. Connector springs are used to couple the wing and fuselage models. Typical fuselage analysis results are presented for two analytical models. Results for a ring-stiffened cylinder model are compared with results from conventional finite-element analyses to assess the accuracy of this new analysis capability. The connection of plate and ring segments is demonstrated using a second model that is representative of the wing structure for a channel-wing aircraft configuration.

Giles, Gary L.↗

Comparison of Response Surface and Kriging Models for Multidisciplinary Design Optimization

In this paper, we compare and contrast the use of second-order response surface models and kriging models for approximating non-random, deterministic computer analyses. After reviewing the response surface method for constructing polynomial approximations, kriging is presented as an alternative approximation method for the design and analysis of computer experiments. Both methods are applied to the multidisciplinary design of an aerospike nozzle which consists of a computational fluid dynamics model and a finite-element model. Error analysis of the response surface and kriging models is performed along with a graphical comparison of the approximations, and four optimization problems m formulated and solved using both sets of approximation models. The second-order response surface models and kriging models-using a constant underlying global model and a Gaussian correlation function-yield comparable results.

Simpson, Timothy W.↗

Ducted-Fan Engine Acoustic Predictions using a Navier-Stokes Code

A Navier-Stokes computer code is used to predict one of the ducted-fan engine acoustic modes that results from rotor-wake/stator-blade interaction. A patched sliding-zone interface is employed to pass information between the moving rotor row and the stationary stator row. The code produces averaged aerodynamic results downstream of the rotor that agree well with a widely used average-passage code. The acoustic mode of interest is generated successfully by the code and is propagated well upstream of the rotor; temporal and spatial numerical resolution are fine enough such that attenuation of the signal is small. Two acoustic codes are used to find the far-field noise. Near-field propagation is computed by using Eversman's wave envelope code, which is based on a finite-element model. Propagation to the far field is accomplished by using the Kirchhoff formula for moving surfaces with the results of the wave envelope code as input data. Comparison of measured and computed far-field noise levels show fair agreement in the range of directivity angles where the peak radiation lobes from the inlet are observed. Although only a single acoustic mode is targeted in this study, the main conclusion is a proof-of-concept: Navier-Stokes codes can be used both to generate and propagate rotor/stator acoustic modes forward through an engine, where the results can be coupled to other far-field noise prediction codes.

Rumsey, C. L.↗

Subtle alternating electrocardiographic morphology as an indicator of decreased cardiac electrical stability

Observations from finite-element computer models, together with analytic developments based on percolation theory have suggested that subtle fluctuations of ECG morphology might serve as an indicator diminished cardiac electrical stability. With fixed-rate atrial pacing in canines, we have previously observed a pattern of alternation in T wave energy which correlated with cardiac electrical stability. We report here on a series of 20 canine experiments in which cardiac electrical stability (measured via Ventricular Fibrillation Threshold determination) was compared to a non-degenerate, multidimensional measurement of the degree of alternating activity present in the ECG complex morphology. The decrease in cardiac electrical stability brought on by both coronary artery occlusion and systemic hypothermia was consistently accompanied by subtle alternation in ECG morphology, with the absolute degree of alternating activity being significantly (negatively) correlated with cardiac electrical stability.

Ventricular Fibrillation↗

Gel stretch method: a new method to measure constitutive properties of cardiac muscle cells

Diastolic dysfunction is an important cause of congestive heart failure; however, the basic mechanisms causing diastolic congestive heart failure are not fully understood, especially the role of the cardiac muscle cell, or cardiocyte, in this process. Before the role of the cardiocyte in this pathophysiology can be defined, methods for measuring cardiocyte constitutive properties must be developed and validated. Thus this study was designed to evaluate a new method to characterize cardiocyte constitutive properties, the gel stretch method. Cardiocytes were isolated enzymatically from normal feline hearts and embedded in a 2% agarose gel containing HEPES-Krebs buffer and laminin. This gel was cast in a shape that allowed it to be placed in a stretching device. The ends of the gel were held between a movable roller and fixed plates that acted as mandibles. Distance between the right and left mandibles was increased using a stepper motor system. The force applied to the gel was measured by a force transducer. The resultant cardiocyte strain was determined by imaging the cells with a microscope, capturing the images with a CCD camera, and measuring cardiocyte and sarcomere length changes. Cardiocyte stress was characterized with a finite-element method. These measurements of cardiocyte stress and strain were used to determine cardiocyte stiffness. Two variables affecting cardiocyte stiffness were measured, the passive elastic spring and viscous damping. The passive spring was assessed by increasing the force on the gel at 1 g/min, modeling the resultant stress vs. strain relationship as an exponential [sigma = A/k(ekepsilon - 1)]. In normal cardiocytes, A = 23.0 kN/m2 and k = 16. Viscous damping was assessed by examining the loop area between the stress vs. strain relationship during 1 g/min increases and decreases in force. Normal cardiocytes had a finite loop area = 1.39 kN/m2, indicating the presence of viscous damping. Thus the gel stretch method provided accurate measurements of cardiocyte constitutive properties. These measurements have allowed the first quantitative assessment of passive elastic spring properties and viscous damping in normal mammalian cardiocytes.

Non-NASA Center↗

Mathematical modelling of flow distribution in the human cardiovascular system

The paper presents a detailed model of the entire human cardiovascular system which aims to study the changes in flow distribution caused by external stimuli, changes in internal parameters, or other factors. The arterial-venous network is represented by 325 interconnected elastic segments. The mathematical description of each segment is based on equations of hydrodynamics and those of stress/strain relationships in elastic materials. Appropriate input functions provide for the pumping of blood by the heart through the system. The analysis employs the finite-element technique which can accommodate any prescribed boundary conditions. Values of model parameters are from available data on physical and rheological properties of blood and blood vessels. As a representative example, simulation results on changes in flow distribution with changes in the elastic properties of blood vessels are discussed. They indicate that the errors in the calculated overall flow rates are not significant even in the extreme case of arteries and veins behaving as rigid tubes.

NASA Center JSC↗

Ice-shelf Dynamics Near the Front of Filchner-Ronne Ice Shelf, Antarctica, Revealed by SAR Interferometry: Model/Interferogram Comparison

We compare Earth Remote Sensing (ERS) satellite synthetic-aperture radar (SAR) interferograms with artificial interferograms constructed using output of a finite-element ice-shelf flow model to study the dynamics of Filchner-Ronne Ice Shelf (FRIS) near Hemmen Ice Rise (HIR) where the iceberg-calving front itersects Berkener Island (BI).

Ice-shelf Filchner-Ronne Antarctica SAR Interferom↗