Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiff problem”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A novel design for a hybrid space manipulator

Described are the structural design, kinematics, and characteristics of a robot manipulator for space applications and use as an articulate and powerful space shuttle manipulator. Hybrid manipulators are parallel-serial connection robots that give rise to a multitude of highly precise robot manipulators. These manipulators are modular and can be extended by additional modules over large distances. Every module has a hemispherical work space and collective modules give rise to highly dexterous symmetrical work space. Some basic designs and kinematic structures of these robot manipulators are discussed, the associated direct and inverse kinematics formulations are presented, and solutions to the inverse kinematic problem are obtained explicitly and elaborated upon. These robot manipulators are shown to have a strength-to-weight ratio that is many times larger than the value that is currently available with industrial or research manipulators. This is due to the fact that these hybrid manipulators are stress-compensated and have an ultralight weight, yet, they are extremely stiff due to the fact that force distribution in their structure is mostly axial. Actuation is prismatic and can be provided by ball screws for maximum precision.

Shahinpoor, MO↗
NASA NTRS1991

A new uniformly valid asymptotic integration algorithm for elasto-plastic creep and unified viscoplastic theories including continuum damage

A new scheme to integrate a system of stiff differential equations for both the elasto-plastic creep and the unified viscoplastic theories is presented. The method has high stability, allows large time increments, and is implicit and iterative. It is suitable for use with continuum damage theories. The scheme was incorporated into MARC, a commercial finite element code through a user subroutine called HYPELA. Results from numerical problems under complex loading histories are presented for both small and large scale analysis. To demonstrate the scheme's accuracy and efficiency, comparisons to a self-adaptive forward Euler method are made.

Chulya, Abhisak↗
NASA NTRS1994

Analytic Sensitivities for Shape Optimization in Equivalent Plate Structural Wing Models

Equivalent plate modeling techniques based on Ritz analysis with simple polynomials prove to be efficient tools for structural modeling of wings in the preliminary design stage. Accuracy problems are encountered, however, when these models are used to obtain finite difference behavior sensitivities with respect to planform shape. The accuracy problems are associated with the poor numerical conditioning of static and eigenvalue equations. As higher-order polynomials are being used to Improve the analysis itself, the more sensitive is the finite difference derivative to the step size used. This article describes a formulation of wing equivalent plate modeling in which it is simple to obtain analytic, explicit expressions for stiffness and mass matrix elements without the need to perform numerical integration. This formulation leads naturally to analytic expressions for the derivatives of displacements, stresses, and natural frequencies with respect to shape design variables. This article examines the accuracy of finite difference derivatives compared with the analytic derivatives, and shows that In some cases it is impossible to obtain any information of value by finite differences. Analytic sensitivities, in this case, are still sufficiently accurate for design optimization.

Livne, Eli↗
DOE OSTI2021

Galerkin formulations of isogeometric shell analysis: Alleviating locking with Greville quadratures and higher-order elements

Here we propose new quadrature schemes that asymptotically require only four in-plane points for Reissner–Mindlin shell elements and nine in-plane points for Kirchhoff–Love shell elements in B-spline and NURBS-based isogeometric shell analysis, independent of the polynomial degree p of the elements. The quadrature points are Greville abscissae associated with pth-order B-spline basis functions whose continuities depend on the specific Galerkin formulations, and the quadrature weights are calculated by solving a linear moment fitting problem in each parametric direction. The proposed shell element formulations are shown through numerical studies to be rank sufficient and to be free of spurious modes. The studies reveal comparable accuracy, in terms of both displacement and stress, compared with fully integrated spline-based shell elements, while at the same time reducing storage and computational cost associated with forming element stiffness and mass matrices and force vectors. The high accuracy with low computational cost makes the proposed quadratures along with higher-order spline bases, in particular polynomial orders, p = 5 and 6, good choices for alleviating membrane and shear locking in shells.

42 ENGINEERING↗
NASA NTRS1976

Thermal stresses in a spherical pressure vessel having temperature-dependent, transversely isotropic, elastic properties

Rayleigh-Ritz and modified Rayleigh-Ritz procedures are used to construct approximate solutions for the response of a thick-walled sphere to uniform pressure loads and an arbitrary radial temperature distribution. The thermoelastic properties of the sphere are assumed to be transversely isotropic and nonhomogeneous; variations in the elastic stiffness and thermal expansion coefficients are taken to be an arbitrary function of the radial coordinate and temperature. Numerical examples are presented which illustrate the effect of the temperature-dependence upon the thermal stress field. A comparison of the approximate solutions with a finite element analysis indicates that Ritz methods offer a simple, efficient, and relatively accurate approach to the problem.

Tauchert, T. R.↗
NASA NTRS1988

Frequency based localization of structural discrepancies

The intent of modal analysis is to develop a reliable model of a structure by working with the analytical and experimental modal properties of frequency, damping and mode shape. In addition to identifying these modal properties, it would be desirable to determine spatially which parts of the structure are modelled poorly or well. It is shown how the pattern of discrepancies in the analytical and experimental test values for the pole and the driving point zero frequencies of a structure can be linked to discrepancies in the mass or stiffness of the structural elements. The success of the procedure depends on the numerical conditioning of a modal reference matrix. Strategies to insure adequate numerical conditioning require a formulation which avoids geometric and energy storage symmetries of the structure, and ignores structural elements which contribute negligibly small potential or kinetic energy to the excited modes. Physical insight into the numerical conditioning problem is provided by a numerical example and by localization of a mass discrepancy in a real structure based on lab tests.

Shepard, G. D.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS1980

Straightening of a wavy strip: An elastic-plastic contact problem including snap-through

The nonlinear behavior of a wave like deformed metal strip during the levelling process were calculated. Elastic-plastic material behavior as well as nonlinearities due to large deformations were considered. The considered problem lead to a combined stability and contact problem. It is shown that, despite the initially concentrated loading, neglecting the change of loading conditions due to altered contact domains may lead to a significant error in the evaluation of the nonlinear behavior and particularly to an underestimation of the stability limit load. The stability was examined by considering the load deflection path and the behavior of a load-dependent current stiffness parameter in combination with the determinant of the current stiffness matrix.

Fischer, D. F.↗
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
NASA NTRS1993

Flight service evaluation of composite components on the Bell helicopter model 206L

This is the final report on the advanced composite components which were placed in service on the 206L LongRanger helicopters in the continental United States, Canada, and Alaska. This report covers all test data which was gathered, as well as maintenance histories of the parts. The previous reports describe the fabrication, service experiences, and test data through 1986. This report contains information from these references, as well as data gathered after 1986. The status of the 40 sets of components is discussed. Each set consisted of a vertical fin, forward fairing, litter door, and baggage door. Almost 500,000 flight hours were accumulated on the 160 parts, with the high-time part accumulating 14,687 flight hours. Over 60 percent of the parts were destructively tested to measure strength and stiffness retention over the course of the program. The vertical fins had the greatest strength retention followed by the litter doors. The baggage doors had the poorest retention of strength. There was very little difference in property retention between the four primary operating regions: Northwest U.S., Southwest U.S., Gulf of Mexico Coastal Region, and the Northeast U.S. and Eastern Canada Region. The field problems have ranged from two lightning-struck fins to significant delaminations in the baggage doors. There was only one environmentally related field incident, in which the glass windows on the litter doors were found to loosen due to high temperatures experienced in the southwest region.

Wilson, Henry↗
NASA NTRS1997

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗
NASA NTRS1976

On the use of two hardening rules of plasticity in incremental and pseudo force analysis

The tangent stiffness and pseudo force forms of the equations of motion are first derived within the context of a total Lagrangian formulation. After a brief discussion of available incremental theory plasticity models, the small strain formulations and computational procedures of the mechanical sublayer model and combined kinematic-isotropic hardening as used in the general purpose structural analysis program AGGIE I are presented. Several sample problems are then presented along with recommended guidelines for use of the two plasticity models.

Hunsaker, B., Jr.↗
NASA NTRS1993

Computational reacting gas dynamics

In the study of high speed flows at high altitudes, such as that encountered by re-entry spacecrafts, the interaction of chemical reactions and other non-equilibrium processes in the flow field with the gas dynamics is crucial. Generally speaking, problems of this level of complexity must resort to numerical methods for solutions, using sophisticated computational fluid dynamics (CFD) codes. The difficulties introduced by reacting gas dynamics can be classified into three distinct headings: (1) the usually inadequate knowledge of the reaction rate coefficients in the non-equilibrium reaction system; (2) the vastly larger number of unknowns involved in the computation and the expected stiffness of the equations; and (3) the interpretation of the detailed reacting CFD numerical results. The research performed accepts the premise that reacting flows of practical interest in the future will in general be too complex or 'untractable' for traditional analytical developments. The power of modern computers must be exploited. However, instead of focusing solely on the construction of numerical solutions of full-model equations, attention is also directed to the 'derivation' of the simplified model from the given full-model. In other words, the present research aims to utilize computations to do tasks which have traditionally been done by skilled theoreticians: to reduce an originally complex full-model system into an approximate but otherwise equivalent simplified model system. The tacit assumption is that once the appropriate simplified model is derived, the interpretation of the detailed numerical reacting CFD numerical results will become much easier. The approach of the research is called computational singular perturbation (CSP).

Lam, S. H.↗
NASA NTRS1988

Precision pointing of scientific instruments on space station: The LFGGREC perspective

An application of Lyapunov function-gradient-generated robustness-enhancing control (LFGGREC) is explored. The attention is directed to a reduced-complexity representation of the pointing problem presented by the system composed of the Space Infrared Telescope Facility gimbaled to a space station configuration. Uncertainties include disturbance forces applied in the crew compartment area and control moments applied to adjacent scientific payloads (modeled as disturbance moments). Also included are uncertainties in gimbal friction and in the structural component of the system, as reflected in the inertia matrix, the damping matrix, and the stiffness matrix, and the effect of the ignored vibrational dynamics of the structure. The emphasis is on the adaptation of LFGGREC to this particular configuration and on the robustness analysis.

Blackwell, C. C.↗
DOE OSTI2024

Formulation of Conductive Inks for DIW Printed Battery Electrodes

Ink formulations for DIW printing of Li ion battery electrodes were optimized for high conductivity, surface area and printability. Our focus is particularly on optimizing Li-ion battery electrodes for low temperature applications by 3-D printing the cathodes into high surface area shapes. The sluggish kinetics of ion transport at the solid-electrolyte phase boundary can be a major problem in batteries at low temperatures, so optimizing for ion transport can help combat these issues. Two families of DIW ink were investigated; aqueous conductive polymer suspension with graphene nano platelets (GNP) ink, and “oil in water” graphene oxide (GO) micro emulsion ink which carbonizes to an open cell graphene aerogel. By iterative formulation changes, the conductivity of PEDOT:PSS based ink was increased over tenfold, and stiffness of GO emulsion ink was increased to provide better printability without sacrificing porosity.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Shape control of high degree-of-freedom variable geometry trusses

Common static trusses are constrained to permit no relative motion between truss elements. A Variable Geometry Truss (VGT), however, is a truss which contains some number of variable length links. The extensible links allow the truss to change shape in a precise, controllable manner. These changes can also be used to control the vibrational response of a truss structure or to perform robotic tasks. Many geometric configurations, both planar and spatial, are possible candidates for VGT manipulators. Here, only two geometries are discussed; the three degree-of-freedom (DOF) spatial octahedral/octahedral truss and the three DOF planar tetrahedral truss. These truss geometries are used as the fundamental element in a repeating chain of trusses. This results in a highly dexterous manipulator with perhaps 30 to 60 degrees of freedom that retains the favorable stiffness properties of a conventional truss. From a fixed base, this type of manipulator could perform shape or vibration control while extending and snaking through complex passageways or moving around obstacles to perform robotic tasks. The approach taken here is to first concentrate on fully understanding the forward and inverse kinematics of the fundamental elements and then utilizing the insight thus gained to solve the more complex problem of the kinematic chains.

Salerno, R. J.↗
NASA NTRS1991

Control issues of microgravity vibration isolation

Active vibration isolation systems contemplated for microgravity space experiments may be designed to reach given performance requirements in a variety of ways. An analogy to passive isolation systems proves to be illustrative but lacks the flexibility as a design tool of a control systems approach and may lead to poor designs. For example, it is shown that a focus on equivalent stiffness in isolation system design leads to a controller that sacrifices robustness for performance. Control theory as applied to vibration isolation is reviewed and passive analogies are discussed. The loop shaping trade-off is introduced and used to design a single-degree-of-freedom fedback controller. An algebraic control design methodology is contrasted to loop shaping and critiqued. Multi-axis vibration isolation and the problems of decoupled single loop control are introduced through a two-degree-of-freedom example problem. It is shown that center of mass uncertainty may result in instability when decoupled single loop control is used. This results from the ill-conditioned nature of the feedback control design. The use of the Linear Quadratic Regulator synthesis procedure for vibration isolation controller design is discussed.

Knospe, C. R.↗
NASA NTRS1984

Stress intensity factors in a reinforced thick-walled cylinder

An elastic thick-walled cylinder containing a radial crack is considered. It is assumed that the cylinder is reinforced by an elastic membrane on its inner surface. The model is intended to simulate pressure vessels with cladding. The formulation of the problem is reduced to a singular integral equation. Various special cases including that of a crack terminating at the cylinder-reinforcement interface are investigated and numerical examples are given. Results indicate that in the case of the crack touching the interface the crack surface displacement derivative is finite and consequently the stress state around the corresponding crack tip is bounded; and generally, for realistic values of the stiffness parameter, the effect of the reinforcement is not very significant.

Tang, R.↗