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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
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↗
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.↗
NASA NTRS1978

An algorithm for optimal structural design with frequency constraints

The paper presents a finite element method for minimum weight design of structures with lower-bound constraints on the natural frequencies, and upper and lower bounds on the design variables. The design algorithm is essentially an iterative solution of the Kuhn-Tucker optimality criterion. The three most important features of the algorithm are: (1) a small number of design iterations are needed to reach optimal or near-optimal design, (2) structural elements with a wide variety of size-stiffness may be used, the only significant restriction being the exclusion of curved beam and shell elements, and (3) the algorithm will work for multiple as well as single frequency constraints. The design procedure is illustrated with three simple problems.

Kiusalaas, J.↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1986

Variable-Force Eddy-Current Damper

Variable damping achieved without problems of containing viscous fluids. Eddy-current damping obtained by moving copper or aluminum conductors through magnetic fields. Position of magnet carrier determines amount of field engagement and, therefore, amount of damping. Three advantages of concept: Magnitudes of stiffness and damping continously varied from maximum to zero without bringing rotor or shaft to stop; used in rotating machines not having viscous fluids available such as lubricating oils; produces sizable damping forces in machines that pump liquid hydrogen at - 246 degrees C and liquid oxygen at - 183 degrees C and are compact in size.

Cunningham, R. E.↗
NASA NTRS1989

Transient response of lattice structures based on exact member theory

The computer program BUNVIS-RG, which treats vibration and buckling of lattice structures using exact member stiffness matrices, has been extended to calculate the exact modal mass and stiffness quantities that can be used in a conventional transient response analysis based on modes. The exact nature of the development allows inclusion of local member response without introduction of any interior member nodes. Results are given for several problems in which significant interaction between local and global response occurs.

Anderson, Melvin S.↗
NASA NTRS1976

Optimum design of composite laminates with thermal effects

An analytical approach to determine an optimum laminate for a variety of thermal and mechanical loading combinations is presented. The analysis is performed for a linear elastic material under static mechanical and uniform thermal loadings. The problem is restricted to a unit width and length laminate with angle orientations resulting in an orthotropic, symmetric, and balanced configuration. An objective function defining total strain energy, is formulated and an optimum laminate design determined subject to constraints on stiffness, average coefficient of thermal expansion, and strength. The objective function is formulated in terms of the orientation angles, number of plies, and material properties. The method presented has, in varying degrees, shown that the design of a laminate can be accomplished using strain energy minimization as the primary criteria. The results of various combinations of applied constraints in the optimized design process are presented and discussed.

Markham, L. R.↗
NASA NTRS1977

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗
NASA NTRS1983

Vibration control of large linear quadratic symmetric systems

Some unique properties on a class of the second order lambda matrices were found and applied to determine a damping matrix of the decoupled subsystem in such a way that the damped system would have preassigned eigenvalues without disturbing the stiffness matrix. The resulting system was realized as a time invariant velocity only feedback control system with desired poles. Another approach using optimal control theory was also applied to the decoupled system in such a way that the mode spillover problem could be eliminated. The procedures were tested successfully by numerical examples.

Jeon, G. J.↗
NASA NTRS1984

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1991

Optimal identification using inconsistent modal data

This work examines techniques under the general approach of optimal-update identification which produce optimally adjusted, or updated, property matrices (i.e., mass, stiffness and damping matrices) to more closely match the structure modal response. For practical applications, the techniques must perform when the modal response is inconsistent with other constraints on the desired model. An alternate view of the optimal-update problem is presented that leads to new techniques for addressing inconsistent data. Viewpoints used for previously published techniques are also examined to explore issues in optimal-update identification.

Smith, Suzanne Weaver↗
NASA NTRS1993

A consistent model reduction of measured modal parameters for reduced-order active control

The problem of synthesizing reduced-order linear models of vibrating structures for the design of fixed-order dynamic feedback control is investigated. The present technique builds on a recently developed procedure for constructing an objective set of mass and stiffness matrices from measured modal parameters that are akin to the Craig-Bampton synthesized ones obtained from finite element models. The constructed mass and stiffness matrices are determined directly from the identification of experimental data, however, rather than through correlation or reconciliation of a finite element model. A model truncation criterion is then applied to the identified minimum-order mass and stiffness model to satisfy certain observability/controllability requirements for the reduced model. Numerical examples illustrate the effectiveness of the proposed technique for synthesizing reduced-order controllers from system realizations of experimental data. The dynamic performance of the resulting closed-loop models is assessed using the known full-order structural dynamics and compared with existing model reduction techniques.

Alvin, K. F.↗
NASA NTRS1998

On the Buckling of Imperfect Anisotropic Shells with Elastic Edge Supports Under Combined Loading Part I:: Theory and Numerical Analysis - Pt. 1

A rigorous solution is presented for the case of stiffened anisotropic cylindrical shells with general imperfections under combined loading, where the edge supports are provided by symmetrical or unsymmetrical elastic rings. The circumferential dependence is eliminated by a truncated Fourier series. The resulting nonlinear 2-point boundary value problem is solved numerically via the "Parallel Shooting Method". The changing deformation patterns resulting from the different degrees of interaction between the given initial imperfections and the specified end rings are displayed. Recommendations are made as to the minimum ring stiffnesses required for optimal load carrying configurations.

Arbocz, Johann↗