SEARCH · Engineering Papers

Results for “finite element methods”

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 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Active vibration mitigation of distributed parameter, smart-type structures using Pseudo-Feedback Optimal Control (PFOC)

A new, near-optimal feedback control technique is introduced that is shown to provide excellent vibration attenuation for those distributed parameter systems that are often encountered in the areas of aeroservoelasticity and large space systems. The technique relies on a novel solution methodology for the classical optimal control problem. Specifically, the quadratic regulator control problem for a flexible vibrating structure is first cast in a weak functional form that admits an approximate solution. The necessary conditions (first-order) are then solved via a time finite-element method. The procedure produces a low dimensional, algebraic parameterization of the optimal control problem that provides a rigorous basis for a discrete controller with a first-order like hold output. Simulation has shown that the algorithm can successfully control a wide variety of plant forms including multi-input/multi-output systems and systems exhibiting significant nonlinearities. In order to firmly establish the efficacy of the algorithm, a laboratory control experiment was implemented to provide planar (bending) vibration attenuation of a highly flexible beam (with a first clamped-free mode of approximately 0.5 Hz).

Patten, W. N.↗

NASA NTRS1989

A combined finite element and boundary integral formulation for solution via CGFFT of 2-dimensional scattering problems

A new technique is presented for computing the scattering by 2-D structures of arbitrary composition. The proposed solution approach combines the usual finite element method with the boundary integral equation to formulate a discrete system. This is subsequently solved via the conjugate gradient (CG) algorithm. A particular characteristic of the method is the use of rectangular boundaries to enclose the scatterer. Several of the resulting boundary integrals are therefore convolutions and may be evaluated via the fast Fourier transform (FFT) in the implementation of the CG algorithm. The solution approach offers the principle advantage of having O(N) memory demand and employs a 1-D FFT versus a 2-D FFT as required with a traditional implementation of the CGFFT algorithm. The speed of the proposed solution method is compared with that of the traditional CGFFT algorithm, and results for rectangular bodies are given and shown to be in excellent agreement with the moment method.

Collins, Jeffery D.↗

NASA NTRS1989

Multilayer theory for delamination analysis of a composite curved bar subjected to end forces and end moments

A composite test specimen in the shape of a semicircular curved bar subjected to bending offers an excellent stress field for studying the open-mode delamination behavior of laminated composite materials. This is because the open-mode delamination nucleates at the midspan of the curved bar. The classical anisotropic elasticity theory was used to construct a multilayer theory for the calculations of the stress and deformation fields induced in the multilayered composite semicircular curved bar subjected to end forces and end moments. The radial location and intensity of the open-mode delamination stress were calculated and were compared with the results obtained from the anisotropic continuum theory and from the finite element method. The multilayer theory gave more accurate predictions of the location and the intensity of the open-mode delamination stress than those calculated from the anisotropic continuum theory.

Ko, William L.↗

NASA NTRS1990

Open-mode delamination stress concentrations in horseshoe and elliptic composite curved bars subjected to end forces

The multilayer theory of anisotropic elasticity and a finite element method were used to analyze the open-mode delamination stress concentrations in horseshoe and elliptic laminated composite curved bars. Two types of laminations, solid laminations and sandwich laminations, were analyzed. It was found that the open-mode delamination stress concentration could be greatly increased in these two types of curved bars by decreasing their aspect ratios. The open-mode delamination stress concentration generated in the solid laminations was found to be far more severe than that generated in the sandwich laminations. The horseshoe curved bar may be used to determine both the open-mode delamination strength of solidly laminated composites and the open-mode debonding strength of sandwiched laminated composites. However, the elliptic curved bar is only good for determining the open-mode delamination strength of solidly laminated composites.

Ko, William L.↗

NASA NTRS1990

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

NASA NTRS1990

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow; (2) the energy equation for the structural-thermal response; and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two application studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗

NASA NTRS1990

Fiber pushout test: A three-dimensional finite element computational simulation

A fiber pushthrough process was computationally simulated using three-dimensional finite element method. The interface material is replaced by an anisotropic material with greatly reduced shear modulus in order to simulate the fiber pushthrough process using a linear analysis. Such a procedure is easily implemented and is computationally very effective. It can be used to predict fiber pushthrough load for a composite system at any temperature. The average interface shear strength obtained from pushthrough load can easily be separated into its two components: one that comes from frictional stresses and the other that comes from chemical adhesion between fiber and the matrix and mechanical interlocking that develops due to shrinkage of the composite because of phase change during the processing. Step-by-step procedures are described to perform the computational simulation, to establish bounds on interfacial bond strength and to interpret interfacial bond quality.

Mital, Subodh K.↗

NASA NTRS1989

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Rismantab-Sany, J.↗

NASA NTRS1989

The coupling effects of kinematics and flexibility on the Lagrangian dynamic formulation of open chain deformable links

A nonlinear Lagrangian formulation for the spatial kinematic and dynamic analysis of open chain deformable links consisting of cylindrical joints that connect pairs of flexible links is developed. The special cases of revolute or prismatic joint can also be obtained from the kinematic equations. The kinematic equations are described using a 4x4 matrix method. The configuration of each deformable link in the open loop kinematic chain is identified using a coupled set of relative joint variables, constant geometric parameters, and elastic coordinates. The elastic coordinates define the link deformation with respect to a selected joint coordinate system that is consistent with the kinematic constraints on the boundary of the deformable link. These coordinates can be introduced using approximation techniques such as Rayleigh-Ritz method, finite element technique or any other desired approach. The large relative motion between two neighboring links are defined by a set of joint coordinates which describes the large relative translational and rotational motion between two neighboring joint coordinate systems. The origin of these coordinate systems are rigidly attached to the neighboring links at the joint definition points along the axis of motion.

Changizi, Koorosh↗

NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

NASA NTRS1990

Strain energy release rate analysis of delamination in a tapered laminate subjected to tension load

A tapered composite laminate subjected to tension load was analyzed using the finite-element method. The glass/epoxy laminate has a (+ or - 45)sub 3 group of plies dropped in three distinct steps, each 20 ply-thicknesses apart, thus forming a taper angle of 5.71 degrees. Steep gradients of interlaminar normal and shear stress on a potential delamination interface suggest the existence of stress singularities at the points of material and geometric discontinuities created by the internal plydrops. The delamination was assumed to initiate at the thin end of the taper on a -45/+45 interface and the delamination growth was simulated in both directions, i.e., along the taper and into the thin region. The strain-energy-release rate for a delamination growing into the thin laminate consisted predominantly of mode I (opening) component. For a delamination growing along the tapered region, the strain-energy-release rate was initially all mode I, but the proportion of mode I decreased with increase in delamination size until eventually total G was all mode II. The total G for both delamination tips increased with increase in delamination size, indicating that a delamination initiating at the end of the taper will grow unstably along the taper and into the thin laminate simultaneously.

Salpekar, S. A.↗

NASA NTRS1990

Current research in cavitating fluid films

A review of the current research of cavitation in fluid films is presented. Phenomena and experimental observations include gaseous cavitation, vapor cavitation, and gas entrainment. Cavitation in flooded, starved, and dynamically loaded journal bearings, as well as squeeze films are reviewed. Observations of cavitation damage in bearings and the possibility of cavitation between parallel plates with microasperities were discussed. The transcavity fluid transport process, meniscus motion and geometry or form of the film during rupture, and reformation were summarized. Performance effects were related to heat transfer models in the cavitated region and hysteresis influence on rotor dynamics coefficients. A number of cavitation algorithms was presented together with solution procedures using the finite difference and finite element methods. Although Newtonian fluids were assumed in most of the discussions, the effect of non-Newtonian fluids on cavitation was also discussed.

Brewe, D. E.↗

NASA NTRS1989

Finite element simulations of thermosolutal convection in vertical solidification of binary alloys

Dendritic vertical solidification of a binary alloy is modeled using the finite element method to assess the effect of thermosolutal convection in macrosegregation. The mathematical model assumes steady-state solidification with a planar, undeformable surface defined by the dendrite tips and the eutectic isotherm. The dendritic region is assumed to advance at a constant solidification velocity v. The stability of the modeled system has been investigated and nonlinear calculations performed that show finger-like convection when the system is unstable. Results for lead-tin alloys show that when the system is unstable, convection is only significant in the uppermost part of the mush and is entirely driven by convection in the bulk fluid.

Heinrich, J. C.↗

NASA NTRS1990

Thermocapillary convection in a model float-zone

A finite-element method has been used to study thermocapillary convection in a model of relevance to the float-zone, crystal-growth process. The geometry consists of a pair of horizontal isothermal surfaces with a Boussinesq liquid suspended between them. Because of the variation of surface tension with temperature, the temperature gradient along the free surface drives thermocapillary convection. The free surface is allowed to deform and its location is calculated along with the velocity and temperature fields. Cases in which thermocapillarity dominates and others in which it interacts with an unstable axial buoyancy gradient are treated. In order to obtain information of possible interest to potential microgravity applications, the influence of the Grashof number is investigated.

Neitzel, G. P.↗

NASA NTRS1990

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow, (2) the energy equation for the structural-thermal response, and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two applications studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗

NASA NTRS1989

Application of HOST technology to the SSME HPFTP blade

The development of turbine blade analysis tools as they moved from research contracts under the NASA HOST (Hot Section Technology) program to a PC-based system for preliminary design and to use in evaluating the SSME HPFTP (Space Shuttle Main Engine High Pressure Fuel Turbo Pump) blade is discussed. The heat transfer and structural analysis of the blade, the development of a data base for constitutive modeling of the blade materials, including coatings, and the use of advanced nonlinear finite element methods are addressed.

Mcknight, R. L.↗

NASA NTRS1990

Mode localization in multi-span beams

The influence of numerous effects on mode localization in multi-span beams is investigated. Finite-element methods are used to study localization as a function of: Timoshenko beam effects; beam end conditions; span length, mass, and stiffness imperfection; viscous damping; axial force; transverse support and rotational coupling stiffness; and modeling resolution. Three configurations are studied, commencing with two different two-span models, and culminating in a ten-span configuration resembling lattice-type large space structures. Results indicate that, in addition to the ratio of imperfection to coupling stiffness being an important localization parameter, transverse support stiffness and Timoshenko beam effects greatly affect the tendency of a structure to exhibit localized modes.

Lust, S. D.↗

NASA NTRS1990

Aeroelastic optimization of a helicopter rotor using an efficient sensitivity analysis

To reduce oscillatory hub loads in forward flight, a structural optimization analysis of a hingeless helicopter rotor has been developed and applied. The aeroelastic analysis of the rotor is based on a finite element method in space and time, and linked with automated optimization algorithms. For the optimization analysis two types of structural representation are used: a generic stiffness-distribution and a single-cell thin-walled beam. For the first type, the design variables are nonstructural mass and its placement, chordwise center of gravity offset from the elastic axis, and stiffness. For the second type, width, height and thickness of spar are used as design variables. For the behavior constraints, frequency placement, autorotational inertia and aeroelastic stability of the blade are included. The required sensitivity derivatives are obtained using a direct analytical approach. An optimum oscillatory hub load shows a 25-77 percent reduction for the generic blade, and 30-50 percent reduction for the box-beam.

Lim, Joon W.↗