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 397 records · Page 22

A simple rectangular element for 2-dimensional analysis of laminated composites

A simple rectangular finite element was developed for two dimensional analysis of laminated composite materials. The rectangular laminated composite element eliminates the need to add elements to a model simply to account for the material properties of various laminae. This is particularly advantageous for thick laminates with many lamina. Explicit integration in terms of generalized displacements minimizes the algebraic effort required to derive the element stiffness and the thermal load vector. A substitute shape function technique is used to improve the performance of the element in modeling bending type deformation. Results for several example problems are discussed.

Whitcomb, J. D.↗

A simple nonlinear joint model

Hertzian contact theory is applied to a butt joint with specially mismatched bearing surfaces to devise a simple mathematical model of nonlinear axial force-displacement behavior in jointed members. Normalized tangent stiffness-force plots, for several values of a joint imperfection parameter, are presented for the sample case of solid structural members of circular section. The results illustrate the potential problem of high joint compliance at low axial-force levels, as well as the generally desirable stiffening and linearizing effects of preload. A nonlinear oscillator problem based on the static model is also formulated and solved to illustrate the effect of amplitude on natural frequency. As expected, natural frequency is low when amplitude is small. The results call attention to the important roles that tight tolerances and preload are expected to play in the design and fabrication of deployable and erectable truss-type space structures.

Fichter, W. B.↗

A simple rectangular element for two-dimensional analysis of laminated composites

A simple rectangular finite element was developed for two-dimensional analysis of laminated composite materials. The rectangular laminated composite element eliminates the need to add elements to a model simply to account for the material properties of various laminae. This is particularly advantageous for thick laminates with many lamina. Explicit integration in terms of generalized displacements minimizes the algebraic effort required to derive the element stiffness and the thermal load vector. A substitute shape function technique is used to improve the performance of the element in modeling bending type deformation. Results for several example problems are discussed.

Whitcomb, J. D.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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.↗

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. Previously announced in STAR as N83-21388

Yahsi, O. S.↗