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 325 records · Page 18

Finite element analysis of partly wrinkled membranes

The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.

Miller, R. K.↗

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

Energy Stability of Thermocapillary Convection in Models of the Float Zone Process

The energy-stability of thermocapillary convection in models of the float-zone, crystal-growing process was studied. Stability limits, as functions of pertinent parameters, that will identify conditions which will not allow the existence of an undesirable oscillatory flow instability were determined. Such instabilities may occur in the space processing of semiconductor materials. The determination of the stability limits will involve two sets of numerical computations: (1) solution of the nonlinear governing equations together with the appropriate boundary conditions to determine the basic state (in general, velocity, pressure and temperature fields and the displacement of free surfaces and interfaces); and (2) solution of a nonlinear Euler-Lagrange systems for the energy-stability limit. Both computations, while difficult, should be within the scope of available computer capability and available concepts in numerical analysis. Finite-element methods are attractive candidates for the numerical work.

Neitzel, G. P.↗

Constitutive modeling for isotropic materials (HOST)

This report presents the results of the second year of work on a problem which is part of the NASA HOST Program. Its goals are: (1) to develop and validate unified constitutive models for isotropic materials, and (2) to demonstrate their usefulness for structural analyses of hot section components of gas turbine engines. The unified models selected for development and evaluation are that of Bodner-Partom and Walker. For model evaluation purposes, a large constitutive data base is generated for a B1900 + Hf alloy by performing uniaxial tensile, creep, cyclic, stress relation, and thermomechanical fatigue (TMF) tests as well as biaxial (tension/torsion) tests under proportional and nonproportional loading over a wide range of strain rates and temperatures. Systematic approaches for evaluating material constants from a small subset of the data base are developed. Correlations of the uniaxial and biaxial tests data with the theories of Bodner-Partom and Walker are performed to establish the accuracy, range of applicability, and integability of the models. Both models are implemented in the MARC finite element computer code and used for TMF analyses. Benchmark notch round experiments are conducted and the results compared with finite-element analyses using the MARC code and the Walker model.

Lindholm, U. S.↗

Dynamic gas temperature measurement system

A gas temperature measurement system with compensated frequency response of 1 KHz and capability to operate in the exhaust of a gas turbine combustor was developed. Environmental guidelines for this measurement are presented, followed by a preliminary design of the selected measurement method. Transient thermal conduction effects were identified as important; a preliminary finite-element conduction model quantified the errors expected by neglecting conduction. A compensation method was developed to account for effects of conduction and convection. This method was verified in analog electrical simulations, and used to compensate dynamic temperature data from a laboratory combustor and a gas turbine engine. Detailed data compensations are presented. Analysis of error sources in the method were done to derive confidence levels for the compensated data.

Elmore, D. L.↗

A method for determining acoustic-liner admittance in ducts with sheared flow in two-cross-sectional directions

A method is developed for determining the acoustic admittance of a test liner installed in the wall of a grazing flow impedance tube. The mean flow is permitted flow gradients in both cross-sectional directions of the tube. The unknown admittance value is obtained by solving an eigenvalue problem. This eigenvalue problem results from the application of the finite-element method to the partial differential equation and boundary conditions governing the acoustic field. The credibility of the method is established by comparing results with exact solutions obtained for a constant mean-flow profile and with previous results for cases involving shear in only one cross-sectional direction. Excellent comparisons are obtained in both cases. The analysis is used in conjunction with a limited amount of experimental data and shows that the flow must be accurately modeled in order to determine the acoustic-liner properties.

Watson, W. R.↗

Multidimensional Modeling of Atmospheric Effects and Surface Heterogeneities on Remote Sensing

The overall goal of this project is to establish a modeling capability that allows a quantitative determination of atmospheric effects on remote sensing including the effects of surface heterogeneities. This includes an improved understanding of aerosol and haze effects in connection with structural, angular, and spatial surface heterogeneities. One important objective of the research is the possible identification of intrinsic surface or canopy characteristics that might be invariant to atmospheric perturbations so that they could be used for scene identification. Conversely, an equally important objective is to find a correction algorithm for atmospheric effects in satellite-sensed surface reflectances. The technical approach is centered around a systematic model and code development effort based on existing, highly advanced computer codes that were originally developed for nuclear radiation shielding applications. Computational techniques for the numerical solution of the radiative transfer equation are adapted on the basis of the discrete-ordinates finite-element method which proved highly successful for one and two-dimensional radiative transfer problems with fully resolved angular representation of the radiation field.

Gerstl, S. A. W.↗

The numerical viscosity of entropy stable schemes for systems of conservation laws

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones.

Tadmor, E.↗

Methods for heat transfer and temperature field analysis of the insulated diesel phase 2 progress report

This report describes work done during Phase 2 of a 3 year program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. During Phase 2, radiation heat transfer model was developed, which accounts for soot formation and burn up. A methodology was developed for carrying out the multi-dimensional finite-element heat conduction calculations within the framework of thermodynamic cycle codes. Studies were carried out using the integrated methodology to address key issues in low heat rejection engines. A wide ranging design analysis matrix was covered, including a variety of insulation strategies, recovery devices and base engine configurations. A single cylinder Cummins engine was installed at Purdue University, and it was brought to a full operational status. The development of instrumentation was continued, concentrating on radiation heat flux detector, total heat flux probe, and accurate pressure-crank angle data acquisition.

Morel, T.↗

NASTRAN/FLEXSTAB procedure for static aeroelastic analysis

Presented is a procedure for using the FLEXSTAB External Structural Influence Coefficients (ESIC) computer program to produce the structural data necessary for the FLEXSTAB Stability Derivatives and Static Stability (SD&SS) program. The SD&SS program computes trim state, stability derivatives, and pressure and deflection data for a flexible airplane having a plane of symmetry. The procedure used a NASTRAN finite-element structural model as the source of structural data in the form of flexibility matrices. Selection of a set of degrees of freedom, definition of structural nodes and panels, reordering and reformatting of the flexibility matrix, and redistribution of existing point mass data are among the topics discussed. Also discussed are boundary conditions and the NASTRAN substructuring technique.

Schuster, L. S.↗

On the design of large flexible space structures (LFSS)

For a general finite-element model of an LFSS, a strictly passive compensator results in an exponentially stable feedback system when actuators and sensors are calculated. In the general case (no colocation), necessary and sufficient conditions are stated on the parameter Q for stabilizing a certain number of modes. Conditions for robust stability are given, and it is shown that feedback does not destabilize the unmodeled modes under certain conditions.

Bhaya, A.↗

Longitudinal compressive failure modes in fiber composites End attachment effects on IITRI type test specimens

The end-attachment effects on longitudinal compressive strength of IITRI type specimen unidirectional fiber composites are formally assessed using finite-element analysis (FEA) in conjunction with composite mechanics. Sixteen different cases were analyzed to evaluate end-attachment effects (such as degree of misalignment, type of misalignment, progressive end-tab debonding, and specimen thickness) on stress distribution, peak stresses, buckling loads, and buckling mode shapes. The results obtained from the FEA and comparisons with fractured specimens show that eccentricities induce bending-type stresses which peak near the end-tabs and cause flexural type fracture. Also, guidelines are included for placing back-to-back strain gages to measure the presence/absence of possible end-attachment and eccentricity effects.

Chamis, C. C.↗

A fracture mechanics approach for designing adhesively bonded joints

An analytical and experimental investigation was undertaken to determine if the adhesive debond initiation stress could be predicted for arbitrary joint geometries. The analysis was based upon a threshold total strain-energy-release rate (Gth) concept. Two bonded systems were tested: T300/5208 graphite/epoxy adherends bonded with either EC-3445 or FM-300 adhesive. The Gth for each adhesive was determined from cracked-lap-shear (CLS) specimens by initiation tests. Finite-element analyses of various tapered CLS specimen geometries predicted the specimen stress at which the total strain-energy-release rate (GT) equaled Gth at the joint tip. Experiments verified the predictions. The approach described herein predicts the maximum stress at which an adhesive joint can be cycled yet not debond. Furthermore, total strain-energy-release rate appeared to be the driving parameter for cyclic debonding and debond initiation in structural adhesives. In addition, debond initiation and growth were found to occur with virtually no peel stress present.

Johnson, W. S.↗

Characterization of Mode I and mixed-mode delamination growth in T300/5208 graphite/epoxy

The roles played by Mode I and Mode II strain-energy release rates (G-I and G-II, respectively) in inducing delamination growth under static and fatigue loading were investigated, using T300/5208 graphite/epoxy specimens. Double cantilever beam (DCB) specimens and cracked lap shear (CLS) specimens were used for pure Mode I and mixed-mode tests, respectively. Fatigue-induced delamination growth was characterized by constant-amplitude fatigue tests at a minimum to maximum cyclic load ratio of 0.05 and a frequency of 10 Hz. During the tests, the maximum and minimum strain-energy release rates (Gmax, Gmin) and the delamination growth rate (da/dN) were monitored. Static tests on mixed-mode CLS specimens measured the total strain-energy release rate, which was broken into G-I and G-II components using finite-element analysis. A power-law relationship between da/dN and G-Imax, and da/dN and Gmax were obtained from fatigue test results on DCB and CLS specimens, respectively. The power law for a pure Mode II delamination was derived from CLS results by subtracting the contribution due to G-I.

Ramkumar, R. L.↗

An evaluation of fracture analysis methods

The results of an experimental and predictive round robin on the applications of fracture analysis methods are presented. The objective of the round robin was to verify whether fracture analysis methods currently in use can or cannot predict failure loads on complex structural components containing cracks. Fracture results from tests on a number of compact specimens were used to make the predictions. The accuracy of the prediction methods was evaluated in terms of the variation in the ratio of predicted to experimental failure loads, and the predictions methods are ranked in order of minimum standard error. A range of applicability of the different methods was also considered in assessing their usefulness. For 7075-T651 aluminum alloy, the best methods were: the effective K sub R curve; the critical crack-tip opening displacement (CTOD) criterion using a finite element analysis; and the K sub R curve with the Dugdale model. For the 2024-T351 aluminum alloy, the best methods included: the two-parameter fracture criterion (TPFC); the CTOD parameter using finite element analysis; the K-curve with the Dugdale model; the deformation plasticity failure assessment diagram (DPFAD); and the effective K sub R curve with a limit load condition. For 304 stainless steel, the best methods were the limit load analysis; the CTOD criterion using finite-element analysis TPFC and DPFAD. Some sample experimental results are given in an appendix.

Newman, J. C., Jr.↗

Practical considerations in aeroelastic design

The structural design process for large transport aircraft is described. Critical loads must be determined from a large number of load cases within the flight maneuver envelope. The structural design is also constrained by considerations of producibility, reliability, maintainability, durability, and damage tolerance, as well as impact dynamics and multiple constraints due to flutter and aeroelasticity. Aircraft aeroelastic design considerations in three distinct areas of product development (preliminary design, advanced design, and detailed design) are presented and contrasted. The present state of the art is challenged to solve the practical difficulties associated with design, analysis, and redesign within cost and schedule constraints. The current practice consists of largely independent engineering disciplines operating with unorganized data interfaces. The need is then demonstrated for a well-planned computerized aeroelastic structural design optimization system operating with a common interdisciplinary data base. This system must incorporate automated interfaces between modular programs. In each phase of the design process, a common finite-element model for static and dynamic optimization is required to reduce errors due to modeling discrepancies. As the design proceeds from the simple models in preliminary design to the more complex models in advanced and detailed design, a means of retrieving design data from the previous models must be established.

Rommel, B. A.↗

PIAS: A program for an iterative aeroelastic solution

A Program for an Iterative Aeroelastic Solution (PIAS) is discussed. This will be a modular computer program that combines the use of a finite-element structural analysis code with any linear or nonlinear aerodynamic code. At this point in time, PIAS has been designed but the software has not been written. The idea for this development originated with P. J. (Bud) Bobbitt of the NASA Langley Research Center. There was initial interest in an aeroelastic solution for a separation-induced leading-edge vortex. Some examples of the flow patterns for a low aspect ratio wing are shown. The Leading-Edge Vortex Program, which calculates pressure distributions including the effects of a separation-induced leading-edge vortex, uses an iterative solution method. This led to the concept of an iteration cycle on configuration shape external to the aerodynamic code.

Manro, M. E.↗

STAEBL: Structural tailoring of engine blades, phase 2

The Structural Tailoring of Engine Blades (STAEBL) program was initiated at NASA Lewis Research Center in 1980 to introduce optimal structural tailoring into the design process for aircraft gas turbine engine blades. The standard procedure for blade design is highly iterative with the engineer directly providing most of the decisions that control the design process. The goal of the STAEBL program has been to develop an automated approach to generate structurally optimal blade designs. The program has evolved as a three-phase effort with the developmental work being performed contractually by Pratt & Whitney Aircraft. Phase 1 was intended as a proof of concept in which two fan blades were structurally tailored to meet a full set of structural design constraints while minimizing DOC+I (direct operating cost plus interest) for a representative aircraft. This phase was successfully completed and was reported in reference 1 and 2. Phase 2 has recently been completed and is the basis for this discussion. During this phase, three tasks were accomplished: (1) a nonproprietary structural tailoring computer code was developed; (2) a dedicated approximate finite-element analysis was developed; and (3) an approximate large-deflection analysis was developed to assess local foreign object damage. Phase 3 is just beginning and is designed to incorporated aerodynamic analyses directly into the structural tailoring system in order to relax current geometric constraints.

Hirschbein, M. S.↗