Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stiff problems”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1993

Sensitivity analysis of flutter response of a typical section and a wing in transonic flow

A sensitivity analysis of flutter response of a two-degree of freedom airfoil with plunging and pitching degrees of freedom in transonic flow was performed using a state-space representation of the unsteady aerodynamic behavior. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady airloads, and the aeroelastic system so modeled is solved as an eigenvalue problem to determine the stability. The results of eigenanalysis showed good agreement with flutter calculations performed using a time-integration of the aeroelastic equations. The sensitivities of the flutter speed with respect to the mass and stiffness parameters wee computed by both the analytical and finite difference methods, showing excellent agreement.

Kapania, Rakesh K.↗
NASA NTRS1993

Morphological instability in epitaxially strained dislocation-free solid films - Linear stability theory

The morphological instability of a growing epitaxially strained dislocation-free solid film is analyzed. An evolution equation for the film surface is derived in the dilute limit of vacancies based on surface diffusion driven by a stress-dependent chemical potential. From the time-dependent linear stability problem the conditions for which a growing film is unstable are determined. It is found that the instability is driven by the lattice mismatch between the film and the substrate; however, low temperatures as well as elastically stiff substrates are stabilizing influences. The results also reveal that the critical film thickness for instability depends on the growth rate of the film itself. Detailed comparison with experimental observations indicates that the instability described exhibits many of the observed features of the onset of the 'island instability'.

Spencer, B. J.↗
NASA NTRS2005

Probabilistic Design of a Wind Tunnel Model to Match the Response of a Full-Scale Aircraft

approach is presented for carrying out the reliability-based design of a plate-like wing that is part of a wind tunnel model. The goal is to design the wind tunnel model to match the stiffness characteristics of the wing box of a flight vehicle while satisfying strength-based risk/reliability requirements that prevents damage to the wind tunnel model and fixtures. The flight vehicle is a modified F/A-18 aircraft. The design problem is solved using reliability-based optimization techniques. The objective function to be minimized is the difference between the displacements of the wind tunnel model and the corresponding displacements of the flight vehicle. The design variables control the thickness distribution of the wind tunnel model. Displacements of the wind tunnel model change with the thickness distribution, while displacements of the flight vehicle are a set of fixed data. The only constraint imposed is that the probability of failure is less than a specified value. Failure is assumed to occur if the stress caused by aerodynamic pressure loading is greater than the specified strength allowable. Two uncertain quantities are considered: the allowable stress and the thickness distribution of the wind tunnel model. Reliability is calculated using Monte Carlo simulation with response surfaces that provide approximate values of stresses. The response surface equations are, in turn, computed from finite element analyses of the wind tunnel model at specified design points. Because the response surface approximations were fit over a small region centered about the current design, the response surfaces were refit periodically as the design variables changed. Coarse-grained parallelism was used to simultaneously perform multiple finite element analyses. Studies carried out in this paper demonstrate that this scheme of using moving response surfaces and coarse-grained computational parallelism reduce the execution time of the Monte Carlo simulation enough to make the design problem tractable. The results of the reliability-based designs performed in this paper show that large decreases in the probability of stress-based failure can be realized with only small sacrifices in the ability of the wind tunnel model to represent the displacements of the full-scale vehicle.

Mason, Brian H.↗
NASA NTRS2015

Modeling and Simulation of Vehicle Dynamics on the Surface of Phobos

In this paper we analyze the dynamics of a spacecraft in proximity of Phobos by developing the equations of motion of a test mass in the Phobos rotating frame using a model based on circularly-restricted three body problem, and by analyzing the dynamics of a ATHLETE hopper vehicle interacting with the soil under different soil-interaction conditions. The main conclusion of the numerical studies is that the system response is dominated by the stiffness and damping parameters of the leg springs, with the soil characteristics having a much smaller effect. The system simulations identify ranges of parameters for which the vehicle emerges stably (relying only on the passive viscoelastic damper at each leg) or unstably (needing active attitude control) from the hop. The implication is that further experimental and possibly computational modeling work, as well as site characterization (from precursor missions) will be necessary to obtain validated performance models.

soil-structure interaction↗
NASA NTRS1950

Contribution to the problem of buckling of orthotropic plates, with special reference to plywood

Planar stress-strain relations and bending stress-strain relations are presented for elastic orthotropic plates and specialized to plywood. These relations are used to derive the differential equation and energy expression for the buckling of orthotropic rectangular plates whose principal stiffness directions are not parallel to the plate edges. Buckling analyses are made for the case of pure compression and pure shear of a long plate-strip.

WOODS↗
NASA NTRS1995

Static stability of a three-dimensional space truss

In order to deploy large flexible space structures it is necessary to develop support systems that are strong and lightweight. The most recent example of this aerospace design need is vividly evident in the space station solar array assembly. In order to accommodate both weight limitations and strength performance criteria, ABLE Engineering has developed the Folding Articulating Square Truss (FASTMast) support structure. The FASTMast is a space truss/mechanism hybrid that can provide system support while adhering to stringent packaging demands. However, due to its slender nature and anticipated loading, stability characterization is a critical part of the design process. Furthermore, the dire consequences surely to result from a catastrophic instability quickly provide the motivation for careful examination of this problem. The fundamental components of the space station solar array system are the (1) solar array blanket system, (2) FASTMast support structure, and (3) mast canister assembly. The FASTMast once fully deployed from the canister will provide support to the solar array blankets. A unique feature of this structure is that the system responds linearly within a certain range of operating loads and nonlinearly when that range is exceeded. The source of nonlinear behavior in this case is due to a changing stiffness state resulting from an inability of diagonal members to resist applied loads. The principal objective of this study was to establish the failure modes involving instability of the FASTMast structure. Also of great interest during this effort was to establish a reliable analytical approach capable of effectively predicting critical values at which the mast becomes unstable. Due to the dual nature of structural response inherent to this problem, both linear and nonlinear analyses are required to characterize the mast in terms of stability. The approach employed herein is one that can be considered systematic in nature. The analysis begins with one and two-dimensional failure models of the system and its important components. From knowledge gained through preliminary analyses a foundation is developed for three-dimensional analyses of the FASTMast structure. The three-dimensional finite element (FE) analysis presented here involves a FASTMast system one-tenth the size of the actual flight unit. Although this study does not yield failure analysis results that apply directly to the flight article, it does establish a method by which the full-scale mast can be evaluated.

Shaker, John F.↗
DOE OSTI2022

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗
NASA NTRS1972

Elastic plate spallation

The dynamic finite element method was used to investigate elastic stress waves in a plate. Strain displacement and stress strain relations are discussed along with the stiffness and mass matrix. The results of studying point load, and distributed load over small, intermediate, and large radii are reported. The derivation of finite element matrices, and the derivation of lumped and consistent matrices for one dimensional problems with Laplace transfer solutions are included. The computer program JMMSPALL is also included.

Oline, L.↗
NASA NTRS1984

Analytical results for postbuckling behavior of plates in compression and shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗
NASA NTRS1985

Analytical results for post-buckling behaviour of plates in compression and in shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗
NASA NTRS1989

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗