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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2004

A Six-Node Curved Triangular Element and a Four-Node Quadrilateral Element for Analysis of Laminated Composite Aerospace Structures

Mathematical development and some computed results are presented for Mindlin plate and shell elements, suitable for analysis of laminated composite and sandwich structures. These elements use the conventional 3 (plate) or 5 (shell) nodal degrees of freedom, have no communicable mechanisms, have no spurious shear energy (no shear locking), have no spurious membrane energy (no membrane locking) and do not require arbitrary reduction of out-of-plane shear moduli or under-integration. Artificial out-of-plane rotational stiffnesses are added at the element level to avoid convergence problems or singularity due to flat spots in shells. This report discusses a 6-node curved triangular element and a 4-node quadrilateral element. Findings show that in regular rectangular meshes, the Martin-Breiner 6-node triangular curved shell (MB6) is approximately equivalent to the conventional 8-node quadrilateral with integration. The 4-node quadrilateral (MB4) has very good accuracy for a 4-node element, and may be preferred in vibration analysis because of narrower bandwidth. The mathematical developments used in these elements, those discussed in the seven appendices, have been applied to elements with 3, 4, 6, and 10 nodes and can be applied to other nodal configurations.

Martin, C. Wayne↗

NASA NTRS2006

Validating Finite Element Models of Assembled Shell Structures

The validation of finite element models of assembled shell elements is presented. The topics include: 1) Problems with membrane rotations in assembled shell models; 2) Penalty stiffness for membrane rotations; 3) Physical stiffness for membrane rotations using shell elements with 6 dof per node; and 4) Connections avoiding rotations.

membrane rotations↗

NASA NTRS1993

Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations

LSODE, the Livermore Solver for Ordinary Differential Equations, is a package of FORTRAN subroutines designed for the numerical solution of the initial value problem for a system of ordinary differential equations. It is particularly well suited for 'stiff' differential systems, for which the backward differentiation formula method of orders 1 to 5 is provided. The code includes the Adams-Moulton method of orders 1 to 12, so it can be used for nonstiff problems as well. In addition, the user can easily switch methods to increase computational efficiency for problems that change character. For both methods a variety of corrector iteration techniques is included in the code. Also, to minimize computational work, both the step size and method order are varied dynamically. This report presents complete descriptions of the code and integration methods, including their implementation. It also provides a detailed guide to the use of the code, as well as an illustrative example problem.

Radhakrishnan, Krishnan↗

NASA NTRS1986

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗

NASA NTRS1999

An Elasticity-Based Mesh Scheme Applied to the Computation of Unsteady Three-Dimensional Spoiler and Aeroelastic Problems

This paper presents a modification of the spring analogy scheme which uses axial linear spring stiffness with selective spring stiffening/relaxation. An alternate approach to solving the geometric conservation law is taken which eliminates the need for storage of metric Jacobians at previous time steps. Efficiency and verification are illustrated with several unsteady 2-D airfoil Euler computations. The method is next applied to the computation of the turbulent flow about a 2-D airfoil and wing with two and three- dimensional moving spoiler surfaces, and the results compared with Benchmark Active Controls Technology (BACT) experimental data. The aeroelastic response at low dynamic pressure of an airfoil to a single large scale oscillation of a spoiler surface is computed. This study confirms that it is possible to achieve accurate solutions with a very large time step for aeroelastic problems using the fluid solver and aeroelastic integrator as discussed in this paper.

Bartels, Robert E.↗

NASA NTRS1984

Precise control of flexible manipulators

The design and experimental testing of end point position controllers for a very flexible one link lightweight manipulator are summarized. The latest upgraded version of the experimental set up, and the basic differences between conventional joint angle feedback and end point position feedback are described. A general procedure for application of modern control methods to the problem is outlined. The relationship between weighting parameters and the bandwidth and control stiffness of the resulting end point position closed loop system is shown. It is found that joint rate angle feedback in addition to the primary end point position sensor is essential for adequate disturbance rejection capability of the closed loop system. The use of a low order multivariable compensator design computer code; called Sandy is documented. A solution to the problem of control mode switching between position sensor sets is outlined. The proof of concept for endpoint position feedback for a one link flexible manipulator was demonstrated. The bandwidth obtained with the experimental end point position controller is about twice as fast as the beam's first natural cantilevered frequency, and comes within a factor of four of the absolute physical speed limit imposed by the wave propagation time of the beam.

Cannon, R. H., Jr.↗

NASA NTRS1973

Numerical analysis of stiffened shear webs in the postbuckling range

The postbuckling behavior of shear webs divided into rectangular panels by stiffeners (uprights) was studied numerically, using the STAGS program in which two-dimensional finite differences are used to solve buckling and nonlinear problems. Universal nondimensional parameters, suggested by linear buckling data, are found to hold for the postbuckling range. Results indicate that the postbuckling stiffness of shear webs with isotropic panels is roughly two thirds of the prebuckling stiffness. The postbuckling behavior of shear webs with isotropic and orthotropic material properties is compared. Some practical aspects of solving nonlinear problems of this type by the numerical method employed are examined.

Stein, M.↗

NASA NTRS1993

Application of symbolic and algebraic manipulation software in solving applied mechanics problems

As its name implies, symbolic and algebraic manipulation is an operational tool which not only can retain symbols throughout computations but also can express results in terms of symbols. This report starts with a history of symbolic and algebraic manipulators and a review of the literatures. With the help of selected examples, the capabilities of symbolic and algebraic manipulators are demonstrated. These applications to problems of applied mechanics are then presented. They are the application of automatic formulation to applied mechanics problems, application to a materially nonlinear problem (rigid-plastic ring compression) by finite element method (FEM) and application to plate problems by FEM. The advantages and difficulties, contributions, education, and perspectives of symbolic and algebraic manipulation are discussed. It is well known that there exist some fundamental difficulties in symbolic and algebraic manipulation, such as internal swelling and mathematical limitation. A remedy for these difficulties is proposed, and the three applications mentioned are solved successfully. For example, the closed from solution of stiffness matrix of four-node isoparametrical quadrilateral element for 2-D elasticity problem was not available before. Due to the work presented, the automatic construction of it becomes feasible. In addition, a new advantage of the application of symbolic and algebraic manipulation found is believed to be crucial in improving the efficiency of program execution in the future. This will substantially shorten the response time of a system. It is very significant for certain systems, such as missile and high speed aircraft systems, in which time plays an important role.

Tsai, Wen-Lang↗

NASA NTRS1982

The crack and wedging problem for an orthotropic strip

The plane elasticity problem for an orthotropic strip containing a crack parallel to its boundaries is considered. The problem is formulated under general mixed mode loading conditions. The stress intensity factors depend on two dimensionless orthotropic constants only. For the crack problem the results are given for a single crack and two collinear cracks. The calculated results show that of the two orthotropic constants the influence of the stiffness ratio on the stress intensity factors is much more significant than that of the shear parameter. The problem of loading the strip by a rigid rectangular lengths continuous contact is maintained along the wedge strip interface; at a certain critical wedge length the separation starts at the midsection of the wedge, and the length of the separation zone increases rapidly with increasing wedge length.

Cinar, A.↗

NASA NTRS1983

The crack and wedging problem for an orthotropic strip

Cinar, A.↗

NASA NTRS1973

A penny shaped crack in a filament-reinforced matrix. 2: The crack problem

The elastostatic interaction problem between a penny-shaped crack and a slender inclusion or filament in an elastic matrix was formulated. For a single filament as well as multiple identical filaments located symmetrically around the crack the problem is shown to reduce to a singular integral equation. The solution of the problem is obtained for various geometries and filament-to-matrix stiffness ratios, and the results relating to the angular variation of the stress intensity factor and the maximum filament stress are presented.

Pacella, A. H.↗

NASA NTRS1987

Space Station alpha joint bearing

Perhaps the most critical structural system aboard the Space Station is the Solar Alpha Rotary Joint which helps align the power generation system with the sun. The joint must provide structural support and controlled rotation to the outboard transverse booms as well as power and data transfer across the joint. The Solar Alpha Rotary Joint is composed of two transition sections and an integral, large diameter bearing. Alpha joint bearing design presents a particularly interesting problem because of its large size and need for high reliability, stiffness, and on orbit maintability. The discrete roller bearing developed is a novel refinement to cam follower technology. It offers thermal compensation and ease of on-orbit maintenance that are not found in conventional rolling element bearings. How the bearing design evolved is summarized. Driving requirements are reviewed, alternative concepts assessed, and the selected design is described.

Everman, Michael R.↗

NASA NTRS1988

Influence of time domain unsteady aerodynamics on coupled flap-lag-torsional aeroelastic stability and response of rotor blades

This paper describes the incorporation of finite-state, time-domain aerodynamics in a flag-lag-torsional aeroelastic stability and response analysis in forward flight. Improvements to a previous formulation are introduced which eliminate spurious singularities. The methodology for solving the aeroelastic stability and response problems with augmented states, in the time domain, is presented using an implicit formulation. Results describing the aeroelastic behavior of soft and stiff in-plane hingeless rotor blades, in forward flight, are presented to illustrate the sensitivity of both the stability and response problems to time domain unsteady aerodynamics.

Friedmann, P. P.↗

NASA NTRS1992

Automating FEA programming

In this paper we briefly describe a combined symbolic and numeric approach for solving mathematical models on parallel computers. An experimental software system, PIER, is being developed in Common Lisp to synthesize computationally intensive and domain formulation dependent phases of finite element analysis (FEA) solution methods. Quantities for domain formulation like shape functions, element stiffness matrices, etc., are automatically derived using symbolic mathematical computations. The problem specific information and derived formulae are then used to generate (parallel) numerical code for FEA solution steps. A constructive approach to specify a numerical program design is taken. The code generator compiles application oriented input specifications into (parallel) FORTRAN77 routines with the help of built-in knowledge of the particular problem, numerical solution methods and the target computer.

Sharma, Naveen↗

NASA NTRS1994

Error estimation and adaptive mesh refinement for parallel analysis of shell structures

The formulation and application of element-level, element-independent error indicators is investigated. This research culminates in the development of an error indicator formulation which is derived based on the projection of element deformation onto the intrinsic element displacement modes. The qualifier 'element-level' means that no information from adjacent elements is used for error estimation. This property is ideally suited for obtaining error values and driving adaptive mesh refinements on parallel computers where access to neighboring elements residing on different processors may incur significant overhead. In addition such estimators are insensitive to the presence of physical interfaces and junctures. An error indicator qualifies as 'element-independent' when only visible quantities such as element stiffness and nodal displacements are used to quantify error. Error evaluation at the element level and element independence for the error indicator are highly desired properties for computing error in production-level finite element codes. Four element-level error indicators have been constructed. Two of the indicators are based on variational formulation of the element stiffness and are element-dependent. Their derivations are retained for developmental purposes. The second two indicators mimic and exceed the first two in performance but require no special formulation of the element stiffness mesh refinement which we demonstrate for two dimensional plane stress problems. The parallelizing of substructures and adaptive mesh refinement is discussed and the final error indicator using two-dimensional plane-stress and three-dimensional shell problems is demonstrated.

Keating, Scott C.↗

NASA NTRS1966

Large deflection and stability analysis by the direct stiffness method

Matrix numerical procedure by digital computer equipment, and application to elastic stability and deflection problems

ELASTIC STABILITY↗

NASA NTRS1993

An inverse method for computation of structural stiffness distributions of aeroelastically optimized wings

An inverse method has been developed to compute the structural stiffness properties of wings given a specified wing loading and aeroelastic twist distribution. The method directly solves for the bending and torsional stiffness distribution of the wing using a modal representation of these properties. An aeroelastic design problem involving the use of a computational aerodynamics method to optimize the aeroelastic twist distribution of a tighter wing operating at maneuver flight conditions is used to demonstrate the application of the method. This exercise verifies the ability of the inverse scheme to accurately compute the structural stiffness distribution required to generate a specific aeroelastic twist under a specified aeroelastic load.

Schuster, David M.↗

NASA NTRS1986

Development of computer program NAS3D using Vector processing for geometric nonlinear analysis of structures

An algorithm for vectorized computation of stiffness matrices of an 8 noded isoparametric hexahedron element for geometric nonlinear analysis was developed. This was used in conjunction with the earlier 2-D program GAMNAS to develop the new program NAS3D for geometric nonlinear analysis. A conventional, modified Newton-Raphson process is used for the nonlinear analysis. New schemes for the computation of stiffness and strain energy release rates is presented. The organization the program is explained and some results on four sample problems are given. The study of CPU times showed that savings by a factor of 11 to 13 were achieved when vectorized computation was used for the stiffness instead of the conventional scalar one. Finally, the scheme of inputting data is explained.

Mangalgiri, P. D.↗