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 667 records · Page 37

Fluid management technology: Liquid slosh dynamics and control

Flight experiments were defined for the Cryogenic On-Orbit Liquid Depot Storage, Acquisition and Transfer Satellite (COLD-SAT) test bed satellite and the Shuttle middeck to help establish the influence of the gravitational environment on liquid slosh dynamics and control. Several analytical and experimental studies were also conducted to support the experiments and to help understand the anticipated results. Both FLOW-3D and NASA-VOF3D computer codes were utilized to simulate low Bond number, small amplitude sloshing, for which the motions are dominated by surface forces; it was found that neither code provided a satisfactory simulation. Thus, a new analysis of low Bond number sloshing was formulated, using an integral minimization technique that will allow the assumptions made about surface physics phenomena to be modified easily when better knowledge becomes available from flight experiments. Several examples were computed by the innovative use of a finite-element structural code. An existing spherical-pendulum analogy of nonlinear, rotary sloshing was also modified for easier use and extended to low-gravity conditions. Laboratory experiments were conducted to determine the requirements for liquid-vapor interface sensors as a method of resolving liquid surface motions in flight experiments. The feasibility of measuring the small slosh forces anticipated in flight experiments was also investigated.

Dodge, Franklin T.↗

Recursive dynamics for flexible multibody systems using spatial operators

Due to their structural flexibility, spacecraft and space manipulators are multibody systems with complex dynamics and possess a large number of degrees of freedom. Here the spatial operator algebra methodology is used to develop a new dynamics formulation and spatially recursive algorithms for such flexible multibody systems. A key feature of the formulation is that the operator description of the flexible system dynamics is identical in form to the corresponding operator description of the dynamics of rigid multibody systems. A significant advantage of this unifying approach is that it allows ideas and techniques for rigid multibody systems to be easily applied to flexible multibody systems. The algorithms use standard finite-element and assumed modes models for the individual body deformation. A Newton-Euler Operator Factorization of the mass matrix of the multibody system is first developed. It forms the basis for recursive algorithms such as for the inverse dynamics, the computation of the mass matrix, and the composite body forward dynamics for the system. Subsequently, an alternative Innovations Operator Factorization of the mass matrix, each of whose factors is invertible, is developed. It leads to an operator expression for the inverse of the mass matrix, and forms the basis for the recursive articulated body forward dynamics algorithm for the flexible multibody system. For simplicity, most of the development here focuses on serial chain multibody systems. However, extensions of the algorithms to general topology flexible multibody systems are described. While the computational cost of the algorithms depends on factors such as the topology and the amount of flexibility in the multibody system, in general, it appears that in contrast to the rigid multibody case, the articulated body forward dynamics algorithm is the more efficient algorithm for flexible multibody systems containing even a small number of flexible bodies. The variety of algorithms described here permits a user to choose the algorithm which is optimal for the multibody system at hand. The availability of a number of algorithms is even more important for real-time applications, where implementation on parallel processors or custom computing hardware is often necessary to maximize speed.

Jain, A.↗

Mixed variational formulations of finite element analysis of elastoacoustic/slosh fluid-structure interaction

A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in a rigid or flexible container by the method of canonical decomposition applied to a modified form of the wave equation in the displacement potential. The general principle is specialized to a mixed two-field principle that contains the fluid displacement potential and pressure as independent fields. This principle contains a free parameter alpha. Semidiscrete finite-element equations of motion based on this principle are displayed and applied to the transient response and free-vibrations of the coupled fluid-structure problem. It is shown that a particular setting of alpha yields a rich set of formulations that can be customized to fit physical and computational requirements. The variational principle is then extended to handle slosh motions in a uniform gravity field, and used to derive semidiscrete equations of motion that account for such effects.

Felippa, Carlos A.↗

Application of the extended boundary condition method to Monte Carlo simulations of scattering of waves by two-dimensional random rough surfaces

The extended boundary condition method is applied to Monte Carlo simulations of two-dimensional random rough surface scattering. The numerical results are compared with one-dimensional random rough surfaces obtained from the finite-element method. It is found that the mean scattered intensity from two-dimensional rough surfaces differs from that of one dimension for rough surfaces with large slopes.

Tsang, L.↗

Effect of surface deposits on electromagnetic propagation in uniform ducts

A finite-element Galerkin formulation has been used to study the effect of material surface deposits on the reflective characteristics of straight uniform ducts with PEC (perfectly electric conducting) walls. Over a wide frequency range, the effect of both single and multiple dielectric surface deposits on the duct reflection coefficient were examined. The power reflection coefficient was found to be significantly increased by the addition of deposits on the wall.

Baumeister, Kenneth J.↗

Recent developments in the micromechanics of advanced composites

A summary of recent developments and results in the area of inelastic response of advanced composites obtained with the micromechanics model proposed by Aboudi is presented. The recent developments include investigation of linear viscoelastic response of polymer matrix composites and yielding and elastoplastic response of metal matrix composites. The predictive capability and accuracy of the model is demonstrated through comparison of the predicted results with experimental data and results of finite-element analyses.

Pindera, M.-J.↗

Three-dimensional elastic-plastic analysis of shallow cracks in single-edge-crack-tension specimens

Three-dimensional, elastic-plastic, finite-element results are presented for single-edge crack-tension specimens with several shallow crack-length-to-width ratios (0.05 less than or equal to a/W less than or equal to 0.5). Results showed the need to model the initial yield plateau in the stress-strain behavior to accurately model deformation of the A36 steel specimens. The crack-tip-opening-displacement was found to be linearly proportional to the crack-mouth-opening displacement. A new deformation dependent plastic-eta factor equation is presented for calculating the J-integral from test load-displacement records. This equation was shown to be accurate for all crack lengths considered.

Shivakumar, Kunigal N.↗

An error analysis of least-squares finite element method of velocity-pressure-vorticity formulation for Stokes problem

A theoretical proof of the optimal rate of convergence for the least-squares method is developed for the Stokes problem based on the velocity-pressure-vorticity formula. The 2D Stokes problem is analyzed to define the product space and its inner product, and the a priori estimates are derived to give the finite-element approximation. The least-squares method is found to converge at the optimal rate for equal-order interpolation.

Chang, Ching L.↗

Thermal-structural panel buckling tests

A titanium-matrix-composite (TMC) hat-stiffened panel of 61 cm sq area and 3.175 cm thick was nondestructively tested to 649 C to examine its buckling characteristics. Compressive loads were applied to the panel in a 978.6 kN uniaxial load frame system. High-temperature testing was performed using quartz lamp heating. A single-strain-age force/stiffness buckling prediction technique was developed to predict panel buckling loads. For the monolithic panel, these test predictions correlated within 10 percent with a finite-element buckling analyses performed elsewhere. Comparisons between force/stiffness predictions and analyses for the TMC panel are in progress.

Thompson, Randolph C.↗

The lightcraft project

Rensselaer Polytechnic Institute has been developing a transatmospheric 'Lightcraft' technology which uses beamed laser energy to propel advanced shuttle craft to orbit. In the past several years, Rensselaer students have analyzed the unique combined-cycle Lightcraft engine, designed a small unmanned Lightcraft Technology Demonstrator, and conceptualized larger manned Lightcraft - to name just a few of the interrelated design projects. The 1990-91 class carried out preliminary and detailed design efforts for a one-person 'Mercury' Lightcraft, using computer-aided design and finite-element structural modeling techniques. In addition, they began construction of a 2.6 m-diameter, full-scale engineering prototype mockup. The mockup will be equipped with three robotic legs that 'kneel' for passenger entry and exit. More importantly, the articulated tripod gear is crucial for accurately pointing at, and tracking the laser relay mirrors, a maneuver that must be performed just prior to liftoff. Also accomplished were further design improvements on a 6-inch-diameter Lightcraft model (for testing in RPI's hypersonic tunnel), and new laser propulsion experiments. The resultant experimental data will be used to calibrate Computational Fluid Dynamic (CFD) codes and analytical laser propulsion models that can simulate vehicle/engine flight conditions along a transatmospheric boost trajectory. These efforts will enable the prediction of distributed aerodynamic and thruster loads over the entire full-scale spacecraft.

Messitt, Don G.↗

Numerical fracture simulation of compact and bend specimens

An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.

Shivakumar, K. N.↗

CAMELOT 2: Personnel transport between Earth and Mars

A summary of the second phase of a design study fo a cycling spacecraft for regular and frequency personnel transport between Earth and a Mars base is presented. The first-phase study, called project CAMELOT, was a mission analysis and configuration design and was completed by the Space Systems Design class at the University of Michigan in the winter term, 1987. In the second phase, the configuration and component design has been developed by the Space Systems Design class in the winter term, 1988, beyond the initial configuration study. Analysis has led to numerous design modifications and improvements, as well as providing more detailed system definition. Major components of the second phase study include innovative detailed designs such as an interface between the rotating and non-rotating portions of the spacecraft, an elevator system that changes orientation to eliminate Coriolis forces, an electromagnetic radiation shield that significantly reduces the mass of the spacecraft, and a detailed assembly sequence and cost analysis. In addition, modifications to the propulsion and the power generation subsystems, taxi-docking facilities, and internal layout design have led to a more efficient, more reliable, and more comfortable vehicle. Finite-element stress and dynamic analysis has verified much of the system and indicates the validity of the design. Computer simulation of many moving components, including torus rotation, elevator and interface operations, and attitude control mechanisms further validate and support the various systems.

Source record↗

Probabilistic system identification in the time domain

The objective of system identification is to determine reliable dynamical models of a structure by systematically using its measured excitation and response. It brings together in an integrated fashion, experimental, analytical, and computational techniques in structural dynamics. Areas of application for system identification include the following: (1) Model Evaluation--assessing assumptions (linearity and equivalent viscous damping) and techniques (finite-element modeling) used to construct theoretical models of a structure; (2) Model Improvement--updating of a theoretical model to enable more accurate response predictions for possible future loads on the structure, or for control of the structure; (3) Empirical Modelling--developing empirical relationships (nonlinear models) or empirical parameter values (modal damping) because the present state of the art does not provide theoretical results; and (4) Damage Detection and Assessment--continual or episodic updating of a structural model through vibration monitoring to detect and locate any structural damage. It can be argued that since the construction or modification of models using test data is subject to inherent uncertainties, the above problems should be properly treated within a Bayesian probabilistic framework. Such a methodology is presented which allows the precision of the estimates of the model parameters to be computed. It also leads to a guiding principle in applications. Namely, when selecting a single model from a given class of models, one should take the most probable model in the class based on the experimental data. Practical applications of this principle are given which are based on the utilization of measured seismic motions in large civil structures. Examples include the application of a computer program MODE-ID to identify modal properties directly from seismic excitation and response time histories from a nine-story steel-frame building at JPL and from a freeway overpass bridge.

Beck, James L.↗

Approximate minimum-time trajectories for two-link flexible manipulators

The method of recursive quadratic programming was used to generate approximate minimum-time tip trajectories for two-link semi-rigid and flexible manipulator movements in the horizontal plane. The manipulator is modeled with an efficient finite-element scheme for an n-link, m-joint system with bending only in the horizontal plane. Constraints on the trajectory include boundary conditions on position and energy for a rest-to-rest maneuver, straight-line tracking between boundary positions, and motor torque limits. Trajectory comparisons utilize a change in the link stiffness to compare a semi-rigid configuration to a flexible one. The level of bending flexibility necessary to excite significant modal behavior is demonstrated. Applied torques for minimum-time maneuvers are shown to be very similar between configurations and retain much of the qualitative character of rigid-body slewing motion.

Eisler, G. R.↗

Structural dynamic and control considerations for Space Station transportation node concepts

This paper discusses the results of two studies which investigated the expected low frequency dynamic characteristics of conceptual Mars and lunar transportation nodes. Both concepts are based on evolution of the assembly-complete Space Station Freedom configuration by the addition of keels, booms, and modules. Finite-element models of the concepts were developed and a set of undamped modes and frequencies was computed below 2 Hz for each concept. The modes were used as basis vectors for modal analyses describing the dynamic response of the structures to typical reboost maneuvers. The reboost maneuvers had fairly complex loading characteristics since the reaction control system (RCS) jets used for the maneuvers were off-modulated to control pitch and yaw attitude. Interaction of the dynamic response of the structure with the closed-loop attitude control system is investigated for both concepts. The elastic response in the solar dynamic region is investigated to evaluate the severity of the dynamic environment, since the solar dynamic systems have stringent sunpointing requirements which must be maintained during the reboost maneuver.

Cooper, Paul A.↗

Elevation control system model for the DSS 13 antenna

In order to meet the requirements for precision pointing of 34-m antennas, adequate control design and simulation software have to be developed along with a detailed description of the supporting analytical tools. This article describes a control system model for the elevation drive of the DSS 13 antenna. The model allows one to simulate elevation dynamics, cross-coupled dynamics in azimuth and elevation, and RF pointing error. A modal state-space model of the antenna structure was obtained from its finite-element model with a free rotating tipping structure. Model reduction techniques were applied separately for the antenna model and rate-loop model, thereby reducing the system order to one-third of the original one while preserving its dynamic properties. Extensive simulation results illustrate properties of the model.

Gawronski, W.↗

The Effects of Variable Mass and Geometry, Pretwist, Shear Deformation and Rotatory Inertia on the Resonant Frequencies of Intact Long Bones: A Finite Element Model Analysis

The influence of pretwist, nonuniformities in mass and flexural stiffness, rotatory inertia and shear deformation on the natural frequencies of intact bones is evaluated by means of a linear elastic, finite-element model which has been programmed for solution on the digital computer. Theoretical results are compared to the results on the forced vibration of intact canine radii obtained experimentally by Thompson. Surprisingly, inclusion of fairly large pretwist angles (from -14 to 12 deg for one specimen) had little affect on the first three frequencies of transverse vibration in either the cranial or lateral directions. Inclusion of shear deformation reduced the third-mode frequency in the stiffest (lateral) direction by about six percent, otherwise shear deformation played a minor role in determining natural frequencies. Similarly. rotatory inertia had negligible influence up to the third natural frequency. The predominant influence on the first three natural frequencies of transverse vibration could be attributed to the variations in mass and flexural stiffness along the length of the test specimens. Different effective moduli of elasticity are required to yield correct absolute values for the frequencies which correspond to experimental findings. thus implying the presence of some inhomogeneities in material properties around the bone cross-section and/or along its length.

Young, Donald R.↗

Subtle alternating electrocardiographic morphology as an indicator of decreased cardiac electrical stability

Observations from finite-element computer models, together with analytic developments based on percolation theory have suggested that subtle fluctuations of ECG morphology might serve as an indicator diminished cardiac electrical stability. With fixed-rate atrial pacing in canines, we have previously observed a pattern of alternation in T wave energy which correlated with cardiac electrical stability. We report here on a series of 20 canine experiments in which cardiac electrical stability (measured via Ventricular Fibrillation Threshold determination) was compared to a non-degenerate, multidimensional measurement of the degree of alternating activity present in the ECG complex morphology. The decrease in cardiac electrical stability brought on by both coronary artery occlusion and systemic hypothermia was consistently accompanied by subtle alternation in ECG morphology, with the absolute degree of alternating activity being significantly (negatively) correlated with cardiac electrical stability.

Ventricular Fibrillation↗