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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS2016

FEAMAC/CARES Stochastic-Strength-Based Damage Simulation Tool for Ceramic Matrix Composites

Reported here is a coupling of two NASA developed codes: CARES (Ceramics Analysis and Reliability Evaluation of Structures) with the MAC/GMC (Micromechanics Analysis Code/ Generalized Method of Cells) composite material analysis code. The resulting code is called FEAMAC/CARES and is constructed as an Abaqus finite element analysis UMAT (user defined material). Here we describe the FEAMAC/CARES code and an example problem (taken from the open literature) of a laminated CMC in off-axis loading is shown. FEAMAC/CARES performs stochastic-strength-based damage simulation response of a CMC under multiaxial loading using elastic stiffness reduction of the failed elements.

composite structures↗

NASA NTRS2015

Stochastic-Strength-Based Damage Simulation Tool for Ceramic Matrix Composite

ceramic matrix composites↗

DOE OSTI2020

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗

DOE OSTI2021

Multirate linearly-implicit GARK schemes

Many complex applications require the solution of initial-value problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the system, according to their dynamics, in order to achieve increased computational efficiency. The stiff components of the system, fast or slow, are best discretized with implicit base methods in order to ensure numerical stability. To this end, linearly implicit methods are particularly attractive as they solve only linear systems of equations at each step. This paper develops the Multirate GARK-ROS/ROW (MR-GARK-ROS/ROW) framework for linearly-implicit multirate time integration. The order conditions theory considers both exact and approximative Jacobians. The effectiveness of implicit multirate methods depends on the coupling between the slow and fast computations; an array of efficient coupling strategies and the resulting numerical schemes are analyzed. Multirate infinitesimal step linearly-implicit methods, that allow arbitrarily small micro-steps and offer extreme computational flexibility, are constructed. The new unifying framework includes existing multirate Rosenbrock(-W) methods as particular cases, and opens the possibility to develop new classes of highly effective linearly implicit multirate integrators.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Studies on the dynamic characteristics of gas film bearings and dampers

The dynamic characteristics of inherently compensated gas film bearings have been investigated for small excursion ratios. Both circular and rectangular cases have been solved for the stiffness and damping as a function of supply pressure, restrictor coefficient, and squeeze number. The effect of disturbance amplitude has been studied for the inherently compensated strip. Analytical solutions for the simple gas film damper problem have established the effect of disturbance amplitude at low squeeze numbers. These results are applicable to pressurized bearings as limiting case of the restrictor coefficient.

Stiffler, A. K.↗

NASA NTRS1984

Vibration and buckling of general periodic lattice structures

A method is presented for vibration and buckling analysis of arbitrary lattice structures having repetitive geometry in any combination of coordinate directions. The approach is based on exact member theory for representing the stiffness of an individual member subject to axial load, and in the case of vibration, undergoing harmonic oscillation. The method is an extension of previous work that was limited to specific geometries. The resulting eigenvalue problem is of the size associated with the repeating element of the structure. A computer program has been developed incorporating the theory and results are given for vibration of rectangular platforms and a large antenna structure having rotational symmetry. Buckling and vibration results for cable-stiffened rings are also given.

Anderson, M. S.↗

NASA NTRS1985

Oil seal effects and subsynchronous vibrations in high-speed compressors

Oil seals are commonly used in high speed multistage compressors. If the oil seal ring becomes locked up against the fixed portion of the seal, high oil film crosscoupled stiffnesses can result. A method of analysis for determining if the oil seals are locked up or not is discussed. The method is then applied to an oil seal in a compressor with subsynchronous vibration problems.

Allaire, P. E.↗

NASA NTRS1988

Shape sensitivity analysis of flutter response of a laminated wing

A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates' modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities, but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other. The sensitivities of the eigenvalues were used to predict the flutter speed, frequency, and reduced frequency. These approximations were found to be in good agreement with those obtained using a complete reanalysis.

Bergen, Fred D.↗

NASA NTRS1989

Shape sensitivity analysis of flutter response of a laminated wing

Kapania, Rakesh K.↗

NASA NTRS2002

Investigating the Use of Ultrasonic Guided Waves for Aging Wire Insulation Assessment

Aging wiring has become a critical issue to DoD, NASA, FAA, and Industry. The problem is that insulation on environmentally aged wire becomes brittle and cracks. This exposes the underlying conductive wire to the potential for short circuits and fire. The difficulty is that techniques to monitor aging wire problems focus on applying electrical sensing techniques that are not very sensitive to the wire insulation. Thus, the development of methods to quantify and monitor aging wire insulation is highly warranted. Measurement of wire insulation stiffness by ultrasonic guided waves is being examined. Initial laboratory tests were performed on a simple model consisting of a solid cylinder and then a solid cylinder with a polymer coating. Experimental measurements showed that the lowest order axisymmetric mode may be sensitive to stiffness changes in the wire insulation. To test this theory, mil-spec wire samples MIL-W-81381, MIL-W-22759/34, and MIL-W-22759/87 (typically found in aircraft) were heat-damaged in an oven, in a range of heating conditions. The samples were 12, 16, and 20 gauge and the heat-damage introduced material changes in the wire-insulation that made the originally flexible insulation brittle and darker in color. Axisymmetric mode phase velocity increased for the samples that were exposed to heat for longer duration. For example, the phase velocity in the 20-gauge MIL-W-22759/34 wire changed from a baseline value of 2790m/s to 3280m/s and 3530m/s for one-hour exposures to 3490C and 3990C, respectively. Although the heat-damage conditions are not the same as environmental aging, we believe that with further development and refinements, the ultrasonic guided waves can be used to inspect wire-insulation for detrimental environmental aging conditions.

Anastasi, Robert F.↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRS1993

Issac, Jason Cherian ses in transonic flow

Flutter analysis of a two degree of freedom airfoil in compressible flow is performed using a state-space representation of the unsteady aerodynamic behavior. Indicial response functions are used to represent the normal force and moment response of the airfoil. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady air loads and the aeroelastic system so modelled is solved as an eigenvalue problem to determine the stability. The aeroelastic equations are also directly integrated with respect to time and the time-domain results compared with the results from the eigenanalysis. A good agreement is obtained. The derivatives of the flutter speed obtained from the eigenanalysis are calculated with respect to the mass and stiffness parameters by both analytical and finite-difference methods for various transonic Mach numbers. The experience gained from the two degree of freedom model is applied to study the sensitivity of the flutter response of a wing with respect to various shape parameters. The parameters being considered are as follows: (1) aspect ratio; (2) surface area of the wing; (3) taper ratio; and (4) sweep. The wing deflections are represented by Chebyshev polynomials. The compressible aerodynamic state-space model used for the airfoil section is extended to represent the unsteady aerodynamic forces on a generally laminated tapered skewed wing. The aeroelastic equations are solved as an eigenvalue problem to determine the flutter speed of the wing. The derivatives of the flutter speed with respect to the shape parameters are calculated by both analytical and finite difference methods.

Issac, Jason Cherion↗

NASA NTRS1925

The Influence of the Form of a Wooden Beam on Its Stiffness and Strength III : Stresses in Wood Members Subjected to Combined Column and Beam Action

The general purpose in this study was to determine the stresses in a wooden member subjected to combined beam and column action. What may be considered the specific purpose, as it relates more directly to the problem of design, was to determine the particular stress that obtains at maximum load which, for combined loading, does not occur simultaneously with maximum stress.

Newlin, J A↗

DOE OSTI2020

TChem v2.0 - A Software Toolkit for the Analysis of Complex Kinetic Models

TChem is an open source software library for solving complex computational chemistry problems and analyzing detailed chemical kinetic models. The software provides support for: complex kinetic models for gas-phase and surface chemistry; thermodynamic properties based on NASA polynomials; species production/consumption rates; stable time integrator for solving stiff time ordinary differential equations; and, reactor models such as homogenous gas-phase ignition (with analytical Jacobian matrices), continuously stirred tank reactor, plug-flow reactor. This toolkit builds upon earlier versions that were written in C and featured tools for gas-phase chemistry only. The current version of the software was completely refactored in C++, uses an object-oriented programming model, and adopts Kokkos as its portability layer to make it ready for the next generation computing architectures i.e., multi/many core computing platforms with GPU accelerators. We have expanded the range of kinetic models to include surface chemistry and have added examples pertaining to Continuously Stirred Tank Reactors (CSTR) and Plug Flow Reactor (PFR) models to complement the homogenous ignition examples present in the earlier versions. To exploit the massive parallelism available from modern computing platforms, the current software interface is designed to evaluate samples in parallel, which enables large scale parametric studies, e.g. for sensitivity analysis and model calibration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

Modeling global structural damping in trusses using simple continuum models

Truss beams with members having viscous damping are modeled as continuum Timoshenko beams. Procedures for deriving the equivalent beam stiffnesses and damping are presented. The global damping for the continuum beam is explicitly expressed in terms of the damping coefficients of the individual truss members. The continuum beam model is used to study transient vibration problems and the solutions are compared well with the full scale finite element solutions. The gradient method is used for parameter estimations in conjunction with the Timoshenko beam model. It is shown that the Timoshenko beam model can be easily updated with measured data and the updated model can yield very accurate transient solutions.

Sun, C. T.↗

NASA NTRS1989

Advanced applications of BEM to inelastic analysis of solids

Newly evolved boundary element solution algorithms are utilized for solving a range of elastoplastic and thermoplastic problems. The first algorithm is the typical iterative procedure in which the unknown boundary solution and the initial stress (or strain) rates are found together in an incremental iterative fashion. The second one is a variable stiffness type approach in which the incremental boundary solution is obtained in a direct (noniterative) manner. The third procedure differs from the two prior algorithms in that volume integration is not needed to incorporate the nonlinear effects in the analysis. Instead, initial stress rates are introduced in the boundary element system by way of particular integrals. The methodology is used in a general purpose, multiregion system that employs quadratic isoparametric shape functions to model the geometry and field variables of the body and can admit up to 15 substructured areas of different material properties.

Banerjee, P. K.↗

NASA NTRS1992

An evaluation of some strategies for vibration control of flexible rotors

There is evidence that the reliability of magnetic bearings has achieved an acceptable level in applications when high cost can be tolerated. This acceptability would be enhanced if the inherent capability of magnetic bearings as active control elements were fully used. The technological and commercial promise of magnetic bearings will be fulfilled only if attention is focussed on the control problems associated with their use. The open loop adaptive control algorithm provides an efficient method of controlling the vibration of rotors without the need of a prior knowledge of parameter values. It overcomes the disadvantages normally associated with open loop control while avoiding the problem of instability associated with closed loop control algorithms. The algorithm is conceptually satisfying because it uses the capability of magnetic bearings as fully active vibration control elements rather than limiting them to act as adjustable stiffness and damping elements, as is the case when they are used with local position and velocity feedback.

Burrows, C. R.↗

NASA NTRS2002

Thermostructural Analysis of Carbon Cloth Phenolic Material Tested at the Laser Hardened Material Evaluation Laboratory

RSRM nozzle liner components have been analyzed and tested to explore the occurrence of anomalous material performance known as pocketing erosion. Primary physical factors that contribute to pocketing seem to include the geometric permeability, which governs pore pressure magnitudes and hence load, and carbon fiber high temperature tensile strength, which defines a material limiting capability. The study reports on the results of a coupled thermostructural finite element analysis of Carbon Cloth Phenolic (CCP) material tested at the Laser Hardened Material Evaluation Laboratory (the LHMEL facility). Modeled test configurations will be limited to the special case of where temperature gradients are oriented perpendicular to the composite material ply angle. Analyses were conducted using a transient, one-dimensional flow/thermal finite element code that models pore pressure and temperature distributions and in an explicitly coupled formulation, passes this information to a 2-dimensional finite element structural model for determination of the stress/deformation behavior of the orthotropic fiber/matrix CCP. Pore pressures are generated by thermal decomposition of the phenolic resin which evolve as a multi-component gas phase which is partially trapped in the porous microstructure of the composite. The nature of resultant pressures are described by using the Darcy relationships which have been modified to permit a multi-specie mass and momentum balance including water vapor condensation. Solution to the conjugate flow/thermal equations were performed using the SINDA code. Of particular importance to this problem was the implementation of a char and deformation state dependent (geometric) permeability as describing a first order interaction between the flow/thermal and structural models. Material property models are used to characterize the solid phase mechanical stiffness and failure. Structural calculations were performed using the ABAQUS code. Iterations were made between the two codes involving the dependent variables temperature, pressure and across-ply strain level. Model results comparisons are made for three different surface heat rates and dependent variable sensitivities discussed for the various cases.

Clayton, J. Louie↗