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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Application of the Green's function method for 2- and 3-dimensional steady transonic flows

A Time-Domain Green's function method for the nonlinear time-dependent three-dimensional aerodynamic potential equation is presented. The Green's theorem is being used to transform the partial differential equation into an integro-differential-delay equation. Finite-element and finite-difference methods are employed for the spatial and time discretizations to approximate the integral equation by a system of differential-delay equations. Solution may be obtained by solving for this nonlinear simultaneous system of equations in time. This paper discusses the application of the method to the Transonic Small Disturbance Equation and numerical results for lifting and nonlifting airfoils and wings in steady flows are presented.

Tseng, K.↗

NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗

NASA NTRS1982

Natural convection in melt crystal growth - The influence of flow pattern on solute segregation

The results of two lines of research aimed at calculating the structure of the flows driven by buoyancy in small-scale crystal growth systems and at understanding the coupling between these flows, the shape of the solidification interface, and dopant segregation in the crystal are reviewed. First, finite-element methods are combined with computer-aided methods for detecting multiple steady solutions to analyze the structure of the buoyancy-driven axisymmetric flows in a vertical cylinder heated from below. This system exhibits onset of convection, multiple steady flows, and loss of the primary stable flow beyond a critical value of the Rayleigh number. Second, results are presented for calculations of convection, melt/solid interface shape, and dopant segregation within a vertical ampoule with thermal boundary conditions that represent a prototype of the vertical Bridgman growth system.

Brown, R. A.↗

NASA NTRS1984

Optimal quasi-static shape control for large aerospace antennae

In this paper, an on-line control approach which will adjust the steady-state shape of a large antenna arbitarily close to any achievable desired profile is proposed. The method makes use of distributed parameter system theory and allows refocusing using a limited number of control actuators and sensors. The controller gains are calculated by approximating the solution to an infinite-dimensional optimal quasi-static control problem. A very general convergence result for such quasi-static controllers is proved and applied to the antenna controller to show convergence, using any Galerkin (finite-element) approximation method.

Balas, M. J.↗

NASA NTRS1984

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis has been developed to calculate the strain energy released by delaminating plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, GI, and shear sliding, GII, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow first before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, GI, for a near-surface delamination can be as high as 0.5GII, and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

NASA NTRS1984

Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Harriott, G. M.↗

NASA NTRS1984

Hygrothermal effects on stresses and deformations in a bonded fiber-reinforced plastic/aluminum system

The hygrothermal effect on time dependent mechanical behavior of a model of nonsymmetrical bonded fiber-reinforced plastic (FRP)/aluminum system was investigated analytically. It was found that a linear viscoelastic closed form solution is sufficient for the prediction model deformations for short periods of time and small strains, while a numerical finite-element solution is preferable for determining adhesive stresses and deformational behavior of the model at any time. The thermoelastic numerical solution shows a certain decrease with time of adhesive stresses, which does not seem to affect significantly the global deformations of the model. On the other hand, it was found that the presence of moisture within the polymeric adhesive layer causes significant decrease in adhesive stresses, which affects deformational changes with time, far beyond the level predicted by the static solution.

Yaniv, G.↗

NASA NTRS1984

Influence of damage on the thermal response of graphite-epoxy laminates

A constrained-displacement finite-element approach for studying the influence of transverse cracks and delaminations on the thermal response of laminated composites is presented. Typical results are given in the form of percent-retention curves for the coefficient of thermal expansion as a function of crack density. Cross-ply and quasi-isotropic T300/5208 graphite-epoxy laminates are considered. It is shown that transverse cracks can have a significant influence on the coefficient of thermal expansion, but delaminations located symmetrically about the laminate midplane have no influence on thermal expansion.

Adams, D. S.↗

NASA NTRS1984

Engine cyclic durability by analysis and material testing

The problem of calculating turbine engine component durability is addressed. Nonlinear, finite-element structural analyses, cyclic constitutive behavior models, and an advanced creep-fatigue life prediction method called strainrange partitioning were assessed for their applicability to the solution of durability problems in hot-section components of gas turbine engines. Three different component or subcomponent geometries are examined: a stress concentration in a turbine disk; a louver lip of a half-scale combustor linear; and a squealer tip of a first-stage high-pressure turbine blade. Cyclic structural analyses were performed for all three problems. The computed strain-temperature histories at the critical locations of the combustor linear and turbine blade components were imposed on smooth specimens in uniaxial, strain-controlled, thermomechanical fatigue tests of evaluate the structural and life analysis methods.

Kaufman, A.↗

NASA NTRS1984

A study of interply layer effects on the free-edge stress field of angleplied laminates

The general-purpose finite-element program MSC/NASTRAN is used to study the interply layer effects on the free-edge stress field of symmetric angleplied laminates subjected to uniform tensile stress. The free-edge region is modeled as a separate substructure (superelement) which enables easy mesh refinement and provides the flexibility to move the superelement along the edge. The results indicate that the interply layer reduces the stress intensity significantly at the free edge. Another important observation of the study is that the failures observed near free edges of these types of laminates could have been caused by the interlaminar shear stresses.

Murthy, P. L. N.↗

NASA NTRS1984

Joined-wing research airplane feasibility study

The joined wing is a new type of aircraft configuration which employs tandem wings arranged to form diamond shapes in plan view and front view. Wind-tunnel tests and finite-element structural analyses have shown that the joined wing provides the following advantages over a comparable wing-plus-tail system; lighter weight and higher stiffness, higher span-efficiency factor, higher trimmed maximum lift coefficient, lower wave drag, plus built-in direct lift and direct sideforce control capability. To verify these advantages at full scale a manned research airplane is required. A study has therefore been performed of the feasibility of constructing such an airplane, using the fuselage and engines of the existing NAA AD-1 oblique-wing airplane. Cost and schedule constraints favored converting the AD-1 rather than constructing a totally new airframe. By removing the outboard wing panels the configuration can simulate wings joined at 60, 80, or 100 percent of span. For maximum versatility the aircraft has alternative control surfaces (such as ailerons and elevators on the front and/or rear wings), and a removeable canard to explore canard/joined-wing interactions at high-lift conditions. Design, performance, and flying qualities are discussed.

Wolkovitch, J.↗

NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗

NASA NTRS1984

Dynamic analysis of straight and involute tooth forms

The effect of load speed on straight and involute tooth forms is studied using several finite-element models. It is found that for rapidly rotating gears and sprockets, the load speed along the tooth surface can significantly affect the tooth vibration. Indeed, it is found that for sufficiently high load speeds and for sufficiently slender tooth forms, the tooth deflection can, at times, be directed opposite to the load direction. Comparisons are made of various dynamic models of gear and sprocket teeth. It is shown that for stubby tooth forms there is considerable difference between results obtained with finite element models and results obtained with Timoshenko beam models. Finally, it is shown that gear or sprocket vibrations can be induced by the shape of the tooth form itself. This effect becomes increasingly significant at higher speeds.

Lin, H. H.↗

NASA NTRS1984

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.↗

NASA NTRS1983

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.↗

NASA NTRS1984

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.↗

NASA NTRS1984

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.↗

NASA NTRS1984

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.↗