Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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,441 records · Page 80

Comparison of NASCAP predictions with experimental data

The NASA charging analyzer program (NASCAP) is a three dimensional, finite element computer code capable of simulating the electrostatic charging of an arbitrary body either in a ground test tank or in the space environment. The code incorporated surface property parameters needed to simulate insulating and conducting materials. These parameters are being updated as required to bring the NASCAP predictions into correspondence with data from ground tests. NASCAP predictions are also being compared with data from the ATS 5 spacecraft. The significance of these results is discussed.

Roche, J. C.↗

Applications of finite element and wave envelope element approximations to turbofan engine noise radiation including flight effects

The problem of acoustic radiation from turbofan engine inlets in flow has not lent itself fully to analysis by numerical means because of the large domains and high frequencies involved. The current work has extended the use of finite elements and wave envelope elements, elements which simulate decay and wavelike behaviour in their interpolation functions, from the no-flow case in which they have been proven, to cases incorporating mean flow. By employing an irrotational mean flow assumption, the acoustics problem has been posed in an axisymmetric formulation in terms of acoustic velocity potential, thus minimizing computer solution storage requirements. The results obtained from the numerical procedures agree well with known analytical solutions, static experimental jet engines inflow data, and also with flight test results.

Parrett, A. V.↗

Quasi-static solution algorithms for kinematically/materially nonlinear thermomechanical problems

This paper develops an algorithmic solution strategy which allows the handling of positive/indefinite stiffness characteristics associated with the pre- and post-buckling of structures subject to complex thermomechanical loading fields. The flexibility of the procedure is such that it can be applied to both finite difference and element-type simulations. Due to the generality of the algorithmic approach developed, both kinematic and thermal/mechanical type material nonlinearity including inelastic effects can be treated. This includes the possibility of handling completely general thermomechanical boundary conditions. To demonstrate the scheme, the results of several benchmark problems is presented.

Padovan, J.↗

Composite transport wing technology development

The design, fabrication, testing, and analysis of stiffened wing cover panels to assess damage tolerance criteria are discussed. The damage tolerance improvements were demonstrated in a test program using full-sized cover panel subcomponents. The panels utilized a hard skin concept with identical laminates of 44-percent 0-degree, 44-percent plus or minus 45-degree, and 12-percent 90-degree plies in the skins and stiffeners. The panel skins were impacted at midbay between the stiffeners, directly over the stiffener, and over the stiffener flange edge. The stiffener blades were impacted laterally. Impact energy levels of 100 ft-lb and 200 ft-lb were used. NASTRAN finite-element analyses were performed to simulate the nonvisible damage that was detected in the panels by nondestructive inspection. A closed-form solution for generalized loading was developed to evaluate the peel stresses in the bonded structure. Two-dimensional delamination growth analysis was developed using the principle of minimum potential energy in terms of closed-form solution for critical strain. An analysis was conducted to determine the residual compressive stress in the panels after impact damage, and the analytical predictions were verified by compression testing of the damaged panels.

Madan, Ram C.↗

Elevated temperature crack growth

A study was performed to examine the applicability of path-independent (P-I) integrals to crack growth problems in hot section components of gas turbine aircraft engines. Alloy 718 was used and the experimental parameters included combined temperature and strain cycling, thermal gradients, elastic-plastic strain levels, and mean strains. A literature review was conducted of proposed P-I integrals, and those capable of analyzing hot section component problems were selected and programmed into the postprocessor of a finite element code. Detailed elastic-plastic finite element analyses were conducted to simulate crack growth and crack closure of the test specimen, and to evaluate the P-I integrals. It was shown that the selected P-I integrals are very effective for predicting crack growth for isothermal conditions.

Kim, K. S.↗

Practical model reduction for flexible structures

The most accepted methods for developing models of flexible structures for the analysis, design and simulation of control systems are finite element methods and extraction from experimental data. Using either of these methods, models are often produced with modes that do not significantly add to the fidelity of the model. This paper presents techniques for eliminating these modes. Particular attention is paid to multiple-input, multiple-output systems. First, control system models developed for flexible structures, using finite element methods and experimental data, are briefly discussed. the shortcomings of using models with unnecessary and/or residual modes are delineated. Then, two techniques for reducing the order of models are presented; the first is applicable to single-input, single-output systems and the second for multiple-input, multiple-output systems. Finally, both are illustrated using model data from the proposed NASA Shuttle-C.

Mitchell, Jerrel R.↗

Interaction of the sea breeze with a river breeze in an area of complex coastal heating

The interaction of the sea-breeze circulation with a river-breeze circulation in an area of complex coastal heating (east coast of Florida) was studied using a 3D finite-element mesoscale model. The model simulations are compared with temperature and wind fields observed on a typical fall day during the Kennedy Space Center Atmospheric Boundary Layer Experiment. The results from numerical experiments designed to isolate the effect of the river breeze indicate that the convergence in the sea-breeze front is suppressed when it passes over the cooler surface of the rivers.

Zhong, Shiyuan↗

Fiber shape effects on metal matrix composite behavior

The effects of different fiber shapes on the behavior of a SiC/Ti-15 metal matrix composite is computationally simulated. A three-dimensional finite element model consisting of a group of nine unidirectional fibers is used in the analysis. The model is employed to represent five different fiber shapes: a circle, an ellipse, a kidney, and two different cross shapes. The distribution of microstresses and the composite material properties, such as moduli, coefficients of thermal expansion, and Poisson's ratios, are obtained from the finite element analysis for the various fiber shapes. Comparisons of these results are used to determine the sensitivity of the composite behavior to the different fiber shapes and assess their potential benefits. No clear benefits result from different fiber shapes though there are some increases/decreases in isolated properties.

Brown, H. C.↗

Fiber shape effects on metal matrix composite behavior

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

Brown, H. C.↗

State of stress, faulting, and eruption characteristics of large volcanoes on Mars

The formation of a large volcano loads the underlying lithospheric plate and can lead to lithospheric flexure and faulting. In turn, lithospheric stresses affect the stress field beneath and within the volcanic edifice and can influence magma transport. Modeling the interaction of these processes is crucial to an understanding of the history of eruption characteristics and tectonic deformation of large volcanoes. We develop models of time-dependent stress and deformation of the Tharsis volcanoes on Mars. A finite element code is used that simulates viscoelastic flow in the mantle and elastic plate flexural behavior. We calculate stresses and displacements due to a volcano-shaped load emplaced on an elastic plate. Models variously incorporate growth of the volcanic load with time and a detachment between volcano and lithosphere. The models illustrate the manner in which time-dependent stresses induced by lithospheric plate flexure beneath the volcanic load may affect eruption histories, and the derived stress fields can be related to tectonic features on and surrounding martian volcanoes.

Mcgovern, Patrick J.↗

Three-dimensional CTOA and constraint effects during stable tearing in a thin-sheet material

A small strain theory, three-dimensional elastic-plastic finite element analysis was used to simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. Both straight and tunneled cracks were modeled. The tunneled crack front shapes as a function of applied stress were obtained from the fracture surface of tested specimens. The stable crack growth behavior was measured at the specimen surface as a function of applied stress. The fracture simulation modeled the crack tunneling and extension as a function of applied stress. The results indicated that the global constraint factor, alpha(sub g), initially dropped during stable crack growth. After peak applied stress was achieved, alpha(sub g) began to increase slightly. The effect of crack front shape on alpha(sub g) was small, but the crack front shape did greatly influence the local constraint and through-thickness crack-tip opening angle (CTOA) behavior. The surface values of CTOA for the tunneled crack front model agreed well with experimental measurements, showing the same initial decrease from high values during the initial 3mm of crack growth at the specimen's surface. At the same time, the interior CTOA values increased from low angles. After the initial stable tearing region, the CTOA was constant through the thickness. The three-dimensional analysis appears to confirm the potential of CTOA as a two-dimensional fracture criterion.

Dawicke, D. S.↗

Influence of Back-Up Bearings and Support Structure Dynamics on the Behavior of Rotors With Active Supports

This report presents a synopsis of the research work. Specific accomplishments are itemized below: (1) Experimental facilities have been developed. This includes a magnetic bearing test rig and an auxiliary bearing test rig. In addition, components have been designed, constructed, and tested for use with a rotordynamics test rig located at NASA Lewis Research Center. (2) A study of the rotordynamics of an auxiliary bearing supported T-501 engine model was performed. (3) An experimental/simulation study of auxiliary bearing rotordynamics has been performed. (4) A rotordynamical model for a magnetic bearing supported rotor system, including auxiliary bearing effects has been developed and simulation studies performed.(5) A finite element model for a foil bearing has been developed and studies of a rotor supported by foil bearings have been performed. (6) Two students affiliated with this project have graduated with M.S. degrees.

Flowers, George T.↗

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗

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↗

Multiphysics Simulation of Active Hypersonic Lip Cooling

This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.

Melis, Matthew E.↗

Modeling and Characterization of Geometric Effects on the Performance of Rainbow and Thunder Actuators

Dome formation in Rainbow and Thunder actuators occurs to relieve thermal expansion mismatch stress between the metallic and piezoelectric layers during cooling from device fabrication temperatures. Accompanying this process is the generation of an internal stress profile within the devices and the development of significant tensile stresses within the surface region of the piezoelectric. These tensile stresses affect the domain configuration (ratio of c-to-a domains), and improve the 90 deg. domain wall movement response of the device in this region of the piezoelectric. This results in improved electromechanical performance compared to standard direct extensional and flextensional devices, presumably because of the contributions of stress to the non-linearity of the piezoelectric d-coefficients. 1 Interestingly, this improvement in response seems counterintuitive; a stress perpendicular to the direction of the applied electric field should impede, not contribute to 90' domain switching. Further consideration of the lower region of the piezoelectric that is under compressive stress thus appears warranted. The specified objectives of the research were to: 1. Conduct finite element and equivalent circuit simulation-based investigations to understand the effects of actuator geometry on internal stress distribution and actuator performance (displacement and load-bearing capabilities). 2. Use the results of the modeling studies to predict the processing conditions (geometry and thickness ratio) required for the fabrication of Rainbow ceramics with optimized performance.

Schwartz, Robert W.↗

Damage Tolerance of Sandwich Plates With Debonded Face Sheets

A nonlinear finite element analysis was performed to simulate axial compression of sandwich beams with debonded face sheets. The load - end-shortening diagrams were generated for a variety of specimens used in a previous experimental study. The energy release rate at the crack tip was computed using the J-integral, and plotted as a function of the load. A detailed stress analysis was performed and the critical stresses in the face sheet and the core were computed. The core was also modeled as an isotropic elastic-perfectly plastic material and a nonlinear post buckling analysis was performed. A Graeco-Latin factorial plan was used to study the effects of debond length, face sheet and core thicknesses, and core density on the load carrying capacity of the sandwich composite. It has been found that a linear buckling analysis is inadequate in determining the maximum load a debonded sandwich beam can carry. A nonlinear post-buckling analysis combined with an elastoplastic model of the core is required to predict the compression behavior of debonded sandwich beams.

Sankar, Bhavani V.↗

Thermal Fault Tolerance Analysis of Carbon Fiber Rope Barrier Systems for Use in the Reusable Solid Rocket Motor ( RSRM) Nozzle Joints

Carbon Fiber Rope (CFR) thermal barrier systems are being considered for use in several RSRM (Reusable Solid Rocket Motor) nozzle joints as a replacement for the current assembly gap close-out process/design. This study provides for development and test verification of analysis methods used for flow-thermal modeling of a CFR thermal barrier subject to fault conditions such as rope combustion gas blow-by and CFR splice failure. Global model development is based on a 1-D (one dimensional) transient volume filling approach where the flow conditions are calculated as a function of internal 'pipe' and porous media 'Darcy' flow correlations. Combustion gas flow rates are calculated for the CFR on a per-linear inch basis and solved simultaneously with a detailed thermal-gas dynamic model of a local region of gas blow by (or splice fault). Effects of gas compressibility, friction and heat transfer are accounted for the model. Computational Fluid Dynamic (CFD) solutions of the fault regions are used to characterize the local flow field, quantify the amount of free jet spreading and assist in the determination of impingement film coefficients on the nozzle housings. Gas to wall heat transfer is simulated by a large thermal finite element grid of the local structure. The employed numerical technique loosely couples the FE (Finite Element) solution with the gas dynamics solution of the faulted region. All free constants that appear in the governing equations are calibrated by hot fire sub-scale test. The calibrated model is used to make flight predictions using motor aft end environments and timelines. Model results indicate that CFR barrier systems provide a near 'vented joint' style of pressurization. Hypothetical fault conditions considered in this study (blow by, splice defect) are relatively benign in terms of overall heating to nozzle metal housing structures.

Clayton, J. Louie↗