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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1990

Hyperbolic conservation laws and numerical methods

The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.

Leveque, Randall J.↗
NASA NTRS1998

Liquid Motion Experiment Flown on STS-84

During some part or all of each mission, about half of all scientific and commercial spacecraft will spin. For example, commercial spacecraft are made to spin during the transfer maneuver from low Earth orbit to the mission orbit to obtain gyroscopic stiffness. Many spacecraft spin continuously in orbit for the same reason. Other reasons for spinning include controlling the location of liquid propellants within their tanks and distributing solar heat loads. Although spinning has many benefits, it also creates problems because of the unavoidable wobble motion that accompanies spinning. Wobbling makes the spacecraft's flexible components oscillate. The energy dissipated by the internal friction of these components causes the wobbling amplitude to increase continually until, at some point, the attitude control thrusters must be fired to bring the spacecraft's amplitude back to an acceptable level.

Dalton, Penni J.↗
NASA NTRS2013

A Coupled Approach for Structural Damage Detection with Incomplete Measurements

This historical work couples model order reduction, damage detection, dynamic residual/mode shape expansion, and damage extent estimation to overcome the incomplete measurements problem by using an appropriate undamaged structural model. A contribution of this work is the development of a process to estimate the full dynamic residuals using the columns of a spring connectivity matrix obtained by disassembling the structural stiffness matrix. Another contribution is the extension of an eigenvector filtering procedure to produce full-order mode shapes that more closely match the measured active partition of the mode shapes using a set of modified Ritz vectors. The full dynamic residuals and full mode shapes are used as inputs to the minimum rank perturbation theory to provide an estimate of damage location and extent. The issues associated with this process are also discussed as drivers of near-term development activities to understand and improve this approach.

James, George↗
DOE OSTI2022

Implementation of Sacrificial Support Structures for Hybrid Manufacturing of Thin Walls

Thin-walled features can be difficult to produce with traditional machining methods which often rely on excess stock material for stiffness. This challenge is increased in hybrid manufacturing where the feature is already near net shape before machining. Significant workpiece deflection can result in poor geometric and surface finish tolerances on the finished part. A potential solution to this problem is to implement sacrificial support structures to the as-printed geometry. The supports are then machined away during the finishing portion of the hybrid process. In the present work, several different design parameters for these sacrificial supports were evaluated to determine their impact on the quality of representative thin wall geometry samples. The angle, height, and spacing of triangular support structures were varied for each sample and then machined and examined. The addition of these supports relative to an unsupported configuration provided a deflection reduction of around 0.2 mm. Surface roughness was improved by approximately 1.5 µm. Increasing values of support height were found to correspond to reduced wall deflection. Similarly, decreasing values of support angle and support spacing improved geometric accuracy. Efficiency comparisons showed that increases in print time corresponded to rapidly diminishing gains in geometric accuracy but continued to improve surface roughness. Implications for hybrid finishing of additively manufactured thin-walled structures is briefly discussed.

36 MATERIALS SCIENCE↗
NASA NTRS1975

A computerized symbolic integration technique for development of triangular and quadrilateral composite shallow-shell finite elements

Computerized symbolic integration was used in conjunction with group-theoretic techniques to obtain analytic expressions for the stiffness, geometric stiffness, consistent mass, and consistent load matrices of composite shallow shell structural elements. The elements are shear flexible and have variable curvature. A stiffness (displacement) formulation was used with the fundamental unknowns consisting of both the displacement and rotation components of the reference surface of the shell. The triangular elements have six and ten nodes; the quadrilateral elements have four and eight nodes and can have internal degrees of freedom associated with displacement modes which vanish along the edges of the element (bubble modes). The stiffness, geometric stiffness, consistent mass, and consistent load coefficients are expressed as linear combinations of integrals (over the element domain) whose integrands are products of shape functions and their derivatives. The evaluation of the elemental matrices is divided into two separate problems - determination of the coefficients in the linear combination and evaluation of the integrals. The integrals are performed symbolically by using the symbolic-and-algebraic-manipulation language MACSYMA. The efficiency of using symbolic integration in the element development is demonstrated by comparing the number of floating-point arithmetic operations required in this approach with those required by a commonly used numerical quadrature technique.

Anderson, C. M.↗
NASA NTRS1977

Flutter of an elastic plate under tension

The problem of flutter of a tensioned elastic plate is reconsidered on a model of two-dimensional aerodynamic flow over a one-dimensional, infinitely wide plate. The purpose is to clarify the effect of variations in plate bending stiffness on the flutter of a plate under a given tension, and vice versa. On the basis that the flutter mode for the one-dimensional infinitely wide plate is essentially the first plate natural mode for subsonic and low supersonic Mach numbers, and that the dynamic pressure required for flutter is proportional to the square of the first natural mode frequency, it is shown that the membrane paradox disappears for Mach numbers 0 and 1.3, so that reduced bending stiffness reduces the dynamic pressure required for flutter.

Dowell, E. H.↗