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.
Information on concurrent algorithms for transient finite element analysis is given in viewgraph form. Information is given on concurrent dynamic algorithms, interprocessor communication, the performance of the BAR problem on the 32 Processor Hypercube, computational efficiency and accuracy analysis.
Dynamic finite element analysis of a real gear housing is presented. The analysis was conducted for the housing without the rotating components (gears, shafts, and bearings). Both rigid and flexible mounting conditions for the gear housing are considered in this analysis. The flexible support simulates the realistic mounting condition on a rotorcraft, and the rigid one is analyzed for comparison purposes. The effect of gear housing stiffeners is also evaluated. The results indicate that the first six natural modes of the flexibly mounted gear housing in the 0 to 200 Hz range correspond to the translational and rotational rigid body vibration modes of the housing. Above this range, the housing plate elastic modes begin to occur. In the case of the rigid mount, only the housing plate elastic modes are observed which are verified by modal analysis experiments. Parametric studies show that the housing plate stiffeners and rigid mounts tend to increase most of the natural frequencies, the lower ones being affected the most.
Stress-intensity factors calculated with a three-dimensional, finite-element analysis for shallow and deep semielliptical surface cracks in finite elastic isotropic plates subjected to tension or bending loads are presented. A wide range of configuration parameters was investigated. The ratio of crack depth to plate thickness ranged from 0.2 to 0.8 and the ratio of crack depth to crack length ranged from 0.2 to 2.0. The effects of plate width on stress-intensity variations along the crack front was also investigated. A wide-range equation for stress-intensity factors along the crack front as a function of crack depth, crack length, plate thickness, and plate width was developed for tension and bending loads. The equation was used to predict patterns of surface-crack growth under tension or bending fatigue loads. A modified form of the equation was also used to correlate surface-crack fracture data for a brittle epoxy material within + or - 10 percent for a wide range of crack shapes and crack sizes.
Alternative techniques for finite element models for transient analyses are explored, using as an example a pinned-pinned beam. Attention is given to the relationships between desired accuracy and the available computer resources. Stiffness, mass, damping, and loading components approximations are examined with note taken of optimization methods. Generation of regular grid results indicates a linear relationship between the number of elements or degrees of freedom and the number of eigenvalues for a given accuracy. Comparisons of the finite element models in terms of accuracy/element or accuracy/degree of freedom were found to be inaccurate for determining model efficiency. Higher accuracies were available from lumped mass and damping models than assumptions of consistent loading, which was preferred to statically equivalent loading.
Stochastic simulation refers to incorporating the effects of design tolerances and uncertainties into the design analysis model and then determining their influence on the design. A high-level evaluation of one such stochastic simulation tool, the MSC.Robust Design tool by MSC.Software Corporation, has been conducted. This stochastic simulation tool provides structural analysts with a tool to interrogate their structural design based on their mathematical description of the design problem using finite element analysis methods. This tool leverages the analyst's prior investment in finite element model development of a particular design. The original finite element model is treated as the baseline structural analysis model for the stochastic simulations that are to be performed. A Monte Carlo approach is used by MSC.Robust Design to determine the effects of scatter in design input variables on response output parameters. The tool was not designed to provide a probabilistic assessment, but to assist engineers in understanding cause and effect. It is driven by a graphical-user interface and retains the engineer-in-the-loop strategy for design evaluation and improvement. The application problem for the evaluation is chosen to be a two-dimensional shell finite element model of a Space Shuttle wing leading-edge panel under re-entry aerodynamic loading. MSC.Robust Design adds value to the analysis effort by rapidly being able to identify design input variables whose variability causes the most influence in response output parameters.
General nonlinear finite element codes were adapted for the purpose of analyzing the dynamics of gas turbine engines. In particular, this adaptation required the development of a squeeze-film damper element software package and its implantation into a representative current generation code. The ADINA code was selected because of prior use of it and familiarity with its internal structure and logic. This objective was met and the results indicate that such use of general purpose codes is viable alternative to specialized codes for general dynamics analysis of engines.
Papers given at the colloquium are presented. Topics discussed include NASTRAN status and plans, computer operations, thermal analysis, NASTRAN applications, vibrations, dynamics, and finite element analysis.
It is found that the standard expressions for elastic stress and displacement in the crack-tip region (i.e., the so-called singular solution) cannot be considered to be approximations that are acceptable in a completely general sense. This conclusion is best illustrated by the instance of a biaxially loaded infinite sheet with a flat horizontal central crack, where the effect of load applied parallel to the plane of the crack appears entirely in the second terms of the series representations for local stresses and displacements. An elastoplastic finite-element analysis of the same biaxially loaded finite specimen geometry shows that the global energy release rate, the J-integral, the plastic stress and strain intensity factors (in the sense of Hilton and Hutchinson), and the size of the crack border region plastic yield, all have pronounced biaxial load dependence.
This paper presents an empirical stress-intensity factor equation for a surface crack as a function of parametric angle, crack depth, crack length, plate thickness and plate width for tension and bending loads. The stress-intensity factors used to develop the equation were obtained from a previous three-dimensional, finite-element analysis of semielliptical surface cracks in finite elastic plates subjected to tension or bending loads. A wide range of configuration parameters was included in the equation. The ratios of crack length to plate thickness and the ratios of crack depth to crack length ranged from 0 to 1.0. The effects of plate width on stress-intensity variations along the crack front were also included. The equation was used to predict patterns of surface-crack growth under tension or bending fatigue loads. The equation was also used to correlate surface-crack fracture data for a brittle epoxy material within + or - 10 percent for a wide range of crack shapes and crack sizes.
Details relating to the design and subsequent vibration and shock tests of a spacecraft component for a severe vibration and shock environment are presented. The design process and analysis method involve the use of finite element analysis coupled with the Modal Strain Energy method with risk graphs to determine the adequacy of the design. The vibration levels experienced by box parts are reduced by the application of passive constrained layer viscoelastic treatments that significantly improve component reliability. All significant internal components were fully instrumented in both the random vibration and shock tests, the latter being done on a mechanical impact pyro simulator. Correlation between analysis and test data is good, validating the modeling and analysis techniques.
Thin-walled rectangular section composite beams are built out of graphite-epoxy prepreg tape using an autoclave molding technique. Cantilevered beams are tested under static loads using a specially-built apparatus and bending slope and elastic twist at a station are measured using a laser system. Measured bending slope and bending-induced twist distributions are correlated with predicted results obtained using a simple beam analysis as well as a refined finite element analysis. Natural vibration characteristics at different rotational speeds are determined experimentally by testing beams in a specially built 10 ft diameter vacuum chamber. The rotor is excited by piezo-electric crystals and response is determined from strain gages mounted on the surface. For composite beams with bending-twisting couplings (symmetric layup), effects of section warping and transverse shear on structural response are found to be substantial. These results were substantiated by the analyses.
A procedure is developed analytically to design and produce an optimized stiffened panel with a centrally located hole. The PASCO panel-buckling analysis and sizing code is employed with the Engineering Analysis Language (EAL) code for finite-element analysis to optimize the design variables. Important design considerations are stiffener height and ply thicknesses in the skin and the stiffener laminates. The EAL code calculates the maximum strains around the hole boundary, and the resulting constraints are combined with the buckling constraint developed by the PASCO code. Cooptimization of the two codes is effected until convergence on the optimal solution is achieved. Results of the analytical procedure are experimentally verified by applying a compressive end-load to graphite-fiber laminated plates with center holes. The experimental results support the analytical procedure, and it is found that for the same weight, the soft-skin designs tested are found to have the best tolerance characteristics and resistance to failure.
Design of a flywheel module, designated the G2 module, is described. The G2 flywheel is a 60,000 RPM, 525 W-hr, 1 kW system designed for a laboratory environment; it will be used for component testing and system demonstrations, with the goal of applying flywheels to aerospace energy storage and integrated power and attitude control (IPACS) applications. G2 has a modular design, which allows for new motors, magnetic bearings, touchdown bearings, and rotors to be installed without a complete redesign of the system. This design process involves several engineering disciplines, and requirements are developed for the speed, energy storage, power level, and operating environment. The G2 rotor system consists of a multilayer carbon fiber rim with a titanium hub on which the other components mount, and rotordynamics analysis is conducted to ensure rigid and flexible rotor modes are controllable or outside of the operating speed range. Magnetic bearings are sized using 1-D magnetic circuit analysis and refined using 3-D finite element analysis. The G2 magnetic bearing system was designed by Texas A&M and has redundancy which allows derated operation after the loss of some components, and an existing liquid cooled two pole permanent magnet motor/generator is used. The touchdown bearing system is designed with a squeeze film damper system allowing spin down from full operating speed in case of a magnetic bearing failure. The G2 flywheel will enable module level demonstrations of component technology, and will be a key building block in system level attitude control and IPACS demonstrations.
The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.
The integrated Finite Element Analysis Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC (Micromechanics Analysis Code/Generalized Method of Cells) multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button-sized disk-shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300 degrees Centigrade to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.
Composite material systems generally exhibit a range of behavior on different length scales (from constituent level to macro); therefore, a multiscale framework is beneficial for the design and engineering of these material systems. The complex nature of the observed composite failure during experiments suggests the need for a three-dimensional (3D) multiscale model to attain a reliable prediction. However, the size of a multiscale three-dimensional finite element model can become prohibitively large and computationally costly. Two-dimensional (2D) models are preferred due to computational efficiency, especially if many different configurations have to be analyzed for an in-depth damage tolerance and durability design study. In this study, various 2D and 3D multiscale analyses will be employed to conduct a detailed investigation into the tensile failure of a given multidirectional, notched carbon fiber reinforced polymer laminate. Threedimensional finite element analysis is typically considered more accurate than a 2D finite element model, as compared with experiments. Nevertheless, in the absence of adequate mesh refinement, large differences may be observed between a 2D and 3D analysis, especially for a shear-dominated layup. This observed difference has not been widely addressed in previous literature and is the main focus of this paper.
Variable-Domain Displacement Transfer Functions were formulated for shape predictions of complex wing structures, for which surface strain-sensing stations must be properly distributed to avoid jointed junctures, and must be increased in the high strain gradient region. Each embedded beam (depth-wise cross section of structure along a surface strain-sensing line) was discretized into small variable domains. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. Through discretization, the embedded beam curvature equation can be piece-wisely integrated to obtain the Variable-Domain Displacement Transfer Functions (for each embedded beam), which are expressed in terms of geometrical parameters of the embedded beam and the surface strains along the strain-sensing line. By inputting the surface strain data into the Displacement Transfer Functions, slopes and deflections along each embedded beam can be calculated for mapping out overall structural deformed shapes. A long tapered cantilever tubular beam was chosen for shape prediction analysis. The input surface strains were analytically generated from finite-element analysis. The shape prediction accuracies of the Variable- Domain Displacement Transfer Functions were then determined in light of the finite-element generated slopes and deflections, and were fofound to be comparable to the accuracies of the constant-domain Displacement Transfer Functions
A finite element analysis of the Iosipescu shear tests for unidirectional and cross-ply composites is presented. It is shown that an iterative analysis procedure must be used to model the fixture-specimen kinematics. The correction factors which are needed to compensate for the nonuniformity of stress distribution in calculating shear modulus are shown to be dependent on the material orthotropic ratio and the finite element loading models. Test section strain distributions representative of typical graphite-epoxy specimens are also presented.