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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Flutter and buckling of general laminated plates

An anisotropic buckling and flutter analysis is developed with allowance for both bending-extensional coupling and bending-twisting coupling within the framework of linear small deflection theory for simply supported general laminated plates. The extended Galerkin method is used to obtain approximate solutions to the coupled governing equations. The effects of various anisotropic stiffness parameters on the static and dynamic stability of laminated plates are evaluated, with particular emphasis on assessing the range of applicability of classical orthotropic plate theory. It is shown that bending-extensional coupling and bending-twisting stiffness terms have a destabilizing effect on buckling and flutter, the effect being more pronounced for a small number of layers. For symmetric plates, the number of layers required for orthotropic plate theory to be applicable is generally less for the buckling problem than for flutter. For square plates, aligning the fibers with the direction of airflow over the plate surface results in the highest flutter dynamic pressure.

Sawyer, J. W.↗

NASA NTRS1979

Matrix perturbation for analytical model improvement

A matrix perturbation method is proposed to calculate the Jacobian matrix and to compute the new eigendata for the parameter estimation procedure. The advantages of the method are the applicability to large complex structures without knowing the analytical expressions for the mass and stiffness matrices, and a cost effective approach for the re-computation of the eigendata. This method also allows the use of other measurements such as modal forces, kinetic energy distribution, and strain energy distributions in the estimation procedure. A realistic sample problem is presented to demonstrate the effectiveness of the proposed method.

Chen, J. C.↗

NASA NTRS1982

Nonlinear and tangent stiffness of imperfect beam columns

A curved member under axial load is analyzed using beam column theory to determine nonlinear response and the tangent stiffness associated with small displacements from the nonlinear state. Such a result is suitable for incorporation into a general nonlinear analysis using a corotational coordinate system to describe the rigid body type motion of individual members. The method is applied to buckling problems. Several examples are given to show the accuracy of the method.

Anderson, M. S.↗

NASA NTRS1985

Nonlinear rotordynamics analysis

The special nonlinearities of the Jeffcott equations in rotordynamics are examined. The immediate application of this analysis is directed toward understanding the excessive vibrations recorded in the LOX pump of the SSME during hot firing ground testing. Deadband, side force and rubbing are three possible sources of inducing nonlinearity in the Jeffcott equations. The present analysis initially reduces these problems to the same mathematical description. A special frequency, named the nonlinear natural frequency is defined and used to develop the solutions of the nonlinear Jeffcott equations as asympotic expansions. This nonlinear natural frequency which is the ratio of the cross-stiffness and the damping, plays a major role in determining response frequencies. Numerical solutions are included for comparison with the analysis. Also, nonlinear frequency-response tables are made for a typical range of values.

Day, W. B.↗

NASA NTRS1986

Exact Bernoulli-Euler static stiffness matrix for a range of tapered beam-columns

The static stiffnesses for torsional and flexural deformation of a tapered beam under axial loading are determined analytically, extending the Bernoulli-Euler/Bessel-function approach of Banerjee and Williams (1985). The derivation of the expressions is explained in detail, and numerical results for sample problems are presented in tables and shown to be in excellent agreement with those obtained using stepped beams of 400 or 500 elements.

Banerjee, J. R.↗

NASA NTRS1994

Transfer function modeling of damping mechanisms in distributed parameter models

This work formulates a method for the modeling of material damping characteristics in distributed parameter models which may be easily applied to models such as rod, plate, and beam equations. The general linear boundary value vibration equation is modified to incorporate hysteresis effects represented by complex stiffness using the transfer function approach proposed by Golla and Hughes. The governing characteristic equations are decoupled through separation of variables yielding solutions similar to those of undamped classical theory, allowing solution of the steady state as well as transient response. Example problems and solutions are provided demonstrating the similarity of the solutions to those of the classical theories and transient responses of nonviscous systems.

Slater, J. C.↗

NASA NTRS1978

NDT of DC-10 graphite-epoxy rudder

The results of visual and ultrasonic nondestructive tests (NDT) on a DC-10 rudder constructed of Thornel 300/5208 graphite-epoxy are presented. A viscous couplant was used for both contact pulse-echo ultrasonics and the Fokker bondtester. Attention is given to the stiffness of cohesive bond structures (the rudder is assembled primarily in a basic six-ply pattern), and to possible delamination. A piezoelectric transducer is employed to sense changes in impedance and signal amplitude. In all, four rudder units were evaluated. The main problems found in the units were porosity and cracking. It is noted that rudder unit four was able to be sufficiently cured to permit static tests on all rib, spar, and skin panel interfaces.

Hagemaier, D. J.↗

NASA NTRS1982

Integrating matrix formulations for vibrations of rotating beams including the effects of concentrated masses

By expressing partial differential equations of motion in matrix notation, utilizing the integrating matrix as a spatial operator, and applying the boundary conditions, the resulting ordinary differential equations can be cast into standard eigenvalue form upon assumption of the usual time dependence. As originally developed, the technique was limited to beams having continuous mass and stiffness properties along their lengths. Integrating matrix methods are extended to treat the differential equations governing the flap, lag, or axial vibrations of rotating beams having concentrated masses. Inclusion of concentrated masses is shown to lead to the same kind of standard eigenvalue problem as before, but with slightly modified matrices.

Lakin, W. D.↗

NASA NTRS1987

Asymptotic expansions in nonlinear rotordynamics

This paper is an examination of special nonlinearities of the Jeffcott equations in rotordynamics. The immediate application of this analysis is directed toward understanding the excessive vibrations recorded in the LOX pump of the SSME during hot-firing ground testing. Deadband, side force, and rubbing are three possible sources of inducing nonlinearity in the Jeffcott equations. The present analysis initially reduces these problems to the same mathematical description. A special frequency, named the nonlinear natural frequency, is defined and used to develop the solutions of the nonlinear Jeffcott equations as singular asymptotic expansions. This nonlinear natural frequency, which is the ratio of the cross-stiffness and the damping, plays a major role in determining response frequencies.

Day, William B.↗

NASA NTRS1994

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗

NASA NTRS1975

NASTRAN nonlinear vibration analysis of beam and frame structures

A capability for the nonlinear vibration analysis of beam and frame structures suitable for use with NASTRAN level 15.5 is described. The nonlinearity considered is due to the presence of axial loads induced by longitudinal end restraints and lateral displacements that are large compared to the beam height. A brief discussion is included of the mathematical analysis and the geometrical stiffness matrix for a prismatic beam (BAR) element. Also included are a brief discussion of the equivalent linearization iterative process used to determine the nonlinear frequency, the required modifications to subroutines DBAR and XMPLBD of the NASTRAN code, and the appropriate vibration capability, four example problems are presented. Comparisons with existing experimental and analytical results show that excellent accuracy is achieved with NASTRAN in all cases.

Mei, C.↗

NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

NASA NTRS1996

CometBoards Users Manual Release 1.0

Several nonlinear mathematical programming algorithms for structural design applications are available at present. These include the sequence of unconstrained minimizations technique, the method of feasible directions, and the sequential quadratic programming technique. The optimality criteria technique and the fully utilized design concept are two other structural design methods. A project was undertaken to bring all these design methods under a common computer environment so that a designer can select any one of these tools that may be suitable for his/her application. To facilitate selection of a design algorithm, to validate and check out the computer code, and to ascertain the relative merits of the design tools, modest finite element structural analysis programs based on the concept of stiffness and integrated force methods have been coupled to each design method. The code that contains both these design and analysis tools, by reading input information from analysis and design data files, can cast the design of a structure as a minimum-weight optimization problem. The code can then solve it with a user-specified optimization technique and a user-specified analysis method. This design code is called CometBoards, which is an acronym for Comparative Evaluation Test Bed of Optimization and Analysis Routines for the Design of Structures. This manual describes for the user a step-by-step procedure for setting up the input data files and executing CometBoards to solve a structural design problem. The manual includes the organization of CometBoards; instructions for preparing input data files; the procedure for submitting a problem; illustrative examples; and several demonstration problems. A set of 29 structural design problems have been solved by using all the optimization methods available in CometBoards. A summary of the optimum results obtained for these problems is appended to this users manual. CometBoards, at present, is available for Posix-based Cray and Convex computers, Iris and Sun workstations, and the VM/CMS system.

Guptill, James D.↗

DOE OSTI2023

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1981

Buckling and vibration of periodic lattice structures

Lattice booms and platforms composed of flexible members or large diameter rings which may be stiffened by cables in order to support membrane-like antennas or reflector surfaces are the main components of some large space structures. The nature of these structures, repetitive geometry with few different members, makes possible relatively simple solutions for buckling and vibration of a certain class of these structures. Each member is represented by a stiffness matrix derived from the exact solution of the beam column equation. This transcendental matrix gives the current member stiffness at any end load or frequency. Using conventional finite element techniques, equilibrium equations can be written involving displacements and rotations of a typical node and its neighbors. The assumptions of a simple trigonometric mode shape is found to satisfy these equations exactly; thus the entire problem is governed by just one 6 x 6 matrix equation involving the amplitude of the displacement and rotation mode shapes. The boundary conditions implied by this solution are simple supported ends for the column type configurations.

Anderson, M. S.↗

NASA NTRS1984

POSTOP: Postbuckled open-stiffener optimum panels, user's manual

The computer program POSTOP developed to serve as an aid in the analysis and sizing of stiffened composite panels that may be loaded in the postbuckling regime, is intended for the preliminary design of metal or composite panels with open-section stiffeners, subjected to multiple combined biaxial compression (or tension), shear and normal pressure load cases. Longitudinal compression, however, is assumed to be the dominant loading. Temperature, initial bow eccentricity and load eccentricity effects are included. The panel geometry is assumed to be repetitive over several bays in the longitudinal (stiffener) direction as well as in the transverse direction. Analytical routines are included to compute panel stiffnesses, strains, local and panel buckling loads, and skin/stiffener interface stresses. The resulting program is applicable to stiffened panels as commonly used in fuselage, wing, or empennage structures. The capabilities and limitations of the code are described. Instructions required to use the program and several example problems are included.

Biggers, S. B.↗

NASA NTRS1985

A Structural Dynamics Approach to the Simulation of Spacecraft Control/Structure Interaction

A relatively simple approach to the analysis of linear spacecraft control/structure interaction problems is presented. The approach uses a commercially available structural system dynamic analysis package for both controller and plant dynamics, thus obviating the need to transfer data between separate programs. The unilateral coupling between components in the control system block diagram is simulated using sparse matrix stiffness and damping elements available in the structural dynamic code. The approach is illustrated with a series of simple tutorial examples of a rigid spacecraft core with flexible appendages.

Young, J. W.↗

NASA NTRS1988

A stereo triangulation system for structural identification: Analytical and experimental results

Identification of large space structures' distributed mass, stiffness, and energy dissipation characteristics poses formidable analytical, numerical, and implementation difficulties. Development of reliable on-orbit structural identification methods is important for implementing active vibration suppression concepts which are under widespread study in the large space structures community. Near the heart of the identification problem lies the necessity of making a large number of spatially distributed measurements of the structure's vibratory response and the associated force/moment inputs with sufficient spatial and frequency resolution. In the present paper, we discuss a method whereby tens of active or passive (retro-reflecting) targets on the structure are tracked simultaneously by the focal planes of two or more video cameras mounted on an adjacent platform. Triangulation (optical ray intersection) of the conjugate image centroids yield inertial trajectories of each target on the structure. Given the triangulated motion of the targets, we apply and extend methodology developed by Creamer, Junkins, and Juang to identify the frequencies, mode shapes, and updated estimates for the mass/stiffness/damping parameterization of the structure. The methodology is semi-automated, for example, the post experiment analysis of the video imagery to determine the inertial trajectories of the targets typically requires less than thirty minutes of real time. Using methodology discussed herein, the frequency response of a large number of points on the structure (where reflective targets are mounted) on the structure can be determined from optical measurements alone. For comparison purposes, we also utilize measurements from accelerometers and a calibrated impulse hammer. While our experimental work remains in a research stage of development, we have successfully tracked and stereo triangulated 20 targets (on a vibrating cantilevered grid structure) at a sample frequency of 200 HZ, and have established conclusively the feasibility and desirability of this approach. We discuss, in summary, recent advances in analog and digital video processing methodology, actuation methods, and bring them to bear on the structural identification problem. We include a brief discussion of our experimental hardware and some recent experimental results which support the practical feasibility of this structural vibration sensing approach.

Junkins, J. L.↗