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 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Error analysis for semi-analytic displacement derivatives with respect to shape and sizing variables

Sensitivity analysis is fundamental to the solution of structural optimization problems. Consequently, much research has focused on the efficient computation of static displacement derivatives. As originally developed, these methods relied on analytical representations for the derivatives of the structural stiffness matrix (K) with respect to the design variables (b sub i). To extend these methods for use with complex finite element formulations and facilitate their implementation into structural optimization programs using the general finite element method analysis codes, the semi-analytic method was developed. In this method the matrix the derivative of K/the derivative b sub i is approximated by finite difference. Although it is well known that the accuracy of the semi-analytic method is dependent on the finite difference parameter, recent work has suggested that more fundamental inaccuracies exist in the method when used for shape optimization. Another study has argued qualitatively that these errors are related to nonuniform errors in the stiffness matrix derivatives. The accuracy of the semi-analytic method is investigated. A general framework was developed for the error analysis and then it is shown analytically that the errors in the method are entirely accounted for by errors in delta K/delta b sub i. Furthermore, it is demonstrated that acceptable accuracy in the derivatives can be obtained through careful selection of the finite difference parameter.

Fenyes, Peter A.↗

NASA NTRS1988

Evaluation of a methodology for model identification in the time domain

A model identification methodology for structural dynamics has been applied to simulated vibrational data as a first step in evaluating its accuracy. The evaluation has taken into account a wide variety of factors which affect the accuracy of the procedure. The effects of each of these factors were observed in both the response time histories and the estimates of the parameters of the model by comparing them with the exact values of the system. Each factor was varied independently but combinations of these have also been considered in an effort to simulate real situations. The results of the tests have shown that for the chain model, the procedure yields robust estimates of the stiffness parameters under the conditions studied whenever uniqueness is ensured. When inaccuracies occur in the results, they are intimately related to non-uniqueness conditions inherent in the inverse problem and not to shortcomings in the methodology.

Beck, R. T.↗

NASA NTRS2005

Application of Ultrasonic Guided Waves for Evaluating Aging Wire Insulation

Aging wiring has become a critical issue to the aerospace and aircraft industries due to Shuttle and aircraft incidents. The problem is that over time the insulation on wire becomes brittle and cracks. This exposes the underlying conductive wire to the potential for short circuits and fire. Popular methods of monitoring aging wire problems focuses on applying electrical sensing techniques that are sensitive to the conductor's condition, but not very sensitive to the wire insulation's condition. Measurement of wire insulation stiffness and ultrasonic properties by ultrasonic guided waves is being examined. Experimental measurements showed that the lowest order extensional mode could be sensitive to stiffness changes in the wire insulation. To test this theory conventional wire samples 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. Results showed that extensional mode phase velocity increased for the samples that were exposed to heat for longer duration.

Anastasi, Robert F.↗

NASA NTRS1975

Numerical solution of quadratic matrix equations for free vibration analysis of structures

This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

Gupta, K. K.↗

NASA NTRS1978

Dynamic model verification of large structural systems

The objective of the present methodology is two-fold: (1) to process test data obtained from either modal survey tests, or slow sine-sweep tests, to extract a set of orthogonal modes best matching the test data while being commensurate with the dynamic model, and (2) to modify submatrices of the dynamic model mass and stiffness matrices to adjust the model to best fit the test data. The method has been implemented using a linear statistical sequential estimator for computation on a CDC computer. Demonstration problems involving Space Shuttle quarter-scale vibration test data and dynamic models have been run. This paper will discuss the general methodology and experience to date.

Lee, L. T.↗

NASA NTRS1981

Development of an orthotropic hole element

A finite element was developed which adequately represents the state of stress in the region around a circular hole in orthotropic material experiencing reasonably general loading. This was achieved with a complementary virtual work formulation of the stiffness and stress matrices for a square element with center circular hole. The assumed stress state provides zero shearing stress on the hole boundary, so the element is suitable for problems involving load transfer without friction. The element has been implemented in the NASTRAN computer program, and sample problem results are presented.

Smith, C. V.↗

NASA NTRS1980

The solution of variable-step implicit difference equations for dynamic systems analysis

A semi-implicit iterative technique is given for the solution of implicit matrix difference equations which arise in variable-step direct time integration of dynamic systems equations. For all practical purposes, the technique avoids repeated refactoring of the stepsize-dependent solution matrix and achieves the solution convergence for engineering accuracy requirements within a permissible number of iterations. Extension of the present method to nonlinear problems treated by the pseudo-force procedure is straightforward. For nonlinear solution procedures based on the tangent modulus formulation, the technique requires modification to make use of the tangent stiffness matrix. It is noted that the question of whether to include the nonlinear terms within the iteration loop has not yet been resolved.

Park, K. C.↗

NASA NTRS1983

Constrained Self-adaptive Solutions Procedures for Structure Subject to High Temperature Elastic-plastic Creep Effects

This paper will develop a new solution strategy which can handle elastic-plastic-creep problems in an inherently stable manner. This is achieved by introducing a new constrained time stepping algorithm which will enable the solution of creep initiated pre/postbuckling behavior where indefinite tangent stiffnesses are encountered. Due to the generality of the scheme, both monotone and cyclic loading histories can be handled. The presentation will give a thorough overview of current solution schemes and their short comings, the development of constrained time stepping algorithms as well as illustrate the results of several numerical experiments which benchmark the new procedure.

Padovan, J.↗

NASA NTRS1987

An efficient quadrilateral element for plate bending analysis

A simple, shear flexible, quadrilateral plate element is developed based on the Hellinger/Reissner mixed variational principle with independently assumed displacement and stress fields. The crucial point of the selection of appropriate stress parameters is emphasized in the formulation. For this purpose, a set of guidelines is formulated based on the follwoing considerations: (1) suppression of all kinematic deformation modes, (2) the element has a favorable value for the constraint index in the thin plate limit, and (3) element properties are frame-invariant. For computer implementation the components of the element stiffness matrix are evaluated analytically using the symbolic manipulation package MACSYMA. The effectiveness and practical usefulness of the proposed element are demonstrated by the numerical results of a variety of problems involving thin and moderately thick plates under different loading and support conditions.

Saleeb, A. F.↗

NASA NTRS1990

Experimental validation of Order (n) DISCOS

An order (n) algorithm is developed for use in the industry-standard DISCOS (Dynamics Interaction Simulation of Controls and Structures) program for flexible multibody system dynamics simulation. This version should be more computationally efficient than conventional DISCOS for problems where a large number of bodies are involved. Program-dependent changes to DISCOS are described for reducing the number of bodies required to model gear reduction and nonlinear stiffness as internal effects. The gear reduction modification can be extended to general applications, such as gear box, rack and pinion, lever joint, and screw rotation models. Upgrading DISCOS for future use include the addition of robotics capabilities such as event-driven topology changes, surface sliding, pick-and-place, and multiarm hand-off.

Chun, Hon M.↗

NASA NTRS1991

Computational aspects of chemically reacting flows

An implicit, time-dependent approach is taken to solve the inviscid conservative equations in a vertex-based finite-volume formulation. The convective fluxes are estimated by a modified Riemann solver, and higher-order accuracy is obtained by means of the MUSCL. Implicit techniques are adopted to deal with the equation stiffness as a result of the clustered grid and chemical source term. The real-gas Roe matrix and the ADI and LU schemes have been tested on the sample problems: hypersonic flow over blunt configurations, such as the forebody of a double ellipse, and the near wake of a conic aerobrake. Temperature distribution and standoff distance are validated against the shock-fitting, central-differencing results. Equilibrium temperature may be recovered by artificially increasing the chemical rates at the wall. The present method satisfactorily resolves the flow structure, including the strong expansion behind the flat base. The relative merits and performance of the real-gas TVD and implicit schemes are discussed briefly.

Li, C. P.↗

NASA NTRS1992

Grounding of space structures

Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. This paper examines the grounding phenomenon in detail. Numerous stiffness matrices developed by others are examined for rigid body rotation capability, and found lacking. A force imbalance inherent in the formulations examined is the likely cause of the grounding problem, suggesting the need for a directed force formulation.

Bosela, P. A.↗

DOE OSTI2022

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1983

A note on the cracked plates reinforced by a line stiffener

The problem of a cracked plate reinforced by a line stiffener is reconsidered. The original solution of this problem was given in the literature. Also, a variation of the problem with debonding between the plate and the stiffener near the cracked region was reported in the literature. However, the special case of the problem in which the crack tip terminates at the stiffener does not appear to have been studied. In practice, the solution may be necessary in order to assess the crack arrest effectiveness of the stiffener. The problem of a stiffened plate with a crack is reformulated, the asymptotic stress state near the crack tip terminating at the stiffener is examined, and numerical results are given for various stiffness constants.

Yahsi, O. S.↗

NASA NTRS1983

A note on the cracked plates reinforced by a line stiffener

Yahsi, O. S.↗

NASA NTRS1986

Computational methods for the identification of spatially varying stiffness and damping in beams

A numerical approximation scheme for the estimation of functional parameters in Euler-Bernoulli models for the transverse vibration of flexible beams with tip bodies is developed. The method permits the identification of spatially varying flexural stiffness and Voigt-Kelvin viscoelastic damping coefficients which appear in the hybrid system of ordinary and partial differential equations and boundary conditions describing the dynamics of such structures. An inverse problem is formulated as a least squares fit to data subject to constraints in the form of a vector system of abstract first order evolution equations. Spline-based finite element approximations are used to finite dimensionalize the problem. Theoretical convergence results are given and numerical studies carried out on both conventional (serial) and vector computers are discussed.

Banks, H. T.↗

NASA NTRS1985

Modal forced vibration analysis of aerodynamically excited turbosystems

Theoretical aspects of a new capability to determine the vibratory response of turbosystems subjected to aerodynamic excitation are presented. Turbosystems such as advanced turbopropellers with highly swept blades, and axial-flow compressors and turbines can be analyzed using this capability. The capability has been developed and implemented in the April 1984 release of the general purpose finite element program NASTRAN. The dynamic response problem is addressed in terms of the normal modal coordinates of these tuned rotating cyclic structures. Both rigid and flexible hubs/disks are considered. Coriolis and centripetal accelerations, as well as differential stiffness effects are included. Generally non-uniform steady inflow fields and uniform flow fields arbitrarily inclined at small angles with respect to the axis of rotation of the turbosystem are considered sources of aerodynamic excitation. The spatial non-uniformities are considered to be small deviations from a principally uniform inflow. Subsonic and supersonic relative inflows are addressed, with provision for linearly interpolating transonic airloads.

Elchuri, V.↗

NASA NTRS1985

NASTRAN supplemental documentation for modal forced vibration analysis of aerodynamically excited turbosystems

This report is a supplemental NASTRAN document for a new capability to determine the vibratory response of turbosystems subjected to aerodynamic excitation. Supplements of NASTRAN Theoretical, User's, Programmer's, and Demonstration Manuals are included. Turbosystems such as advanced turbopropellers with highly swept blades, and axial-flow compressors and turbines can be analyzed using this capability, which has been developed and implemented in the April 1984 release of the general purpose finite element program NASTRAN. The dynamic response problem is addressed in terms of the normal modal coordinates of these tuned rotating cyclic structures. Both rigid and flexible hubs/disks are considered. Coriolis and centripetal accelerations, as well as differential stiffness effects are included. Generally nonuniform steady inflow fields and uniform flow fields arbitrarily inclined at small angles with respect to the axis of rotation of the turbosystem are considered as the sources of aerodynamic excitation. The spatial nonuniformities are considered to be small deviations from a principally uniform inflow. Subsonic relative inflows are addressed, with provision for linearly interpolating transonic airloads.

Elchuri, V.↗