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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗

NASA NTRS1975

Long fluid filled bags suspended by line forces

A previous analysis of fluid filled storage bags is extended to the case of a long fluid filled cylindrical membrane supported by uniform line loads. Cross-sectional shape, stiffness of the support system and stress resultants in the membrane are determined. The application of the numerical results to problems arising in the design of nonrigid airships is discussed.

Mullins, M. L.↗

NASA NTRS1977

Application of the aerodynamic energy concept to the selection of transfer functions for flutter suppression and gust alleviation using active controls

Changes are introduced in the aerodynamic energy approach which lead to an increase in the effectiveness of both the trailing-edge (TE) and the leading-edge (LE)-TE control systems. Control laws are determined, using realizable transfer functions, which permit the introduction of aerodynamic damping and stiffness terms in accordance with the requirements of any specific system. It is shown that flutter suppression and gust alleviation problems can successfully be treated either by a TE or by a LE-TE control system. The results obtained are applicable to a very wide class of aircraft operating within the subsonic Mach number range.

Nissim, E.↗

NASA NTRS1948

The buckling of a column on equally spaced deflectional and rotational springs

A solution is presented for the problem of the buckling of a column on equally spaced deflectional and rotational springs. Useful charts, which relate deflectional spring stiffness, rotational spring stiffness, and buckling load, are given for columns having two, three, four, and infinite number of spans.

FORMULAS - BUCKLING↗

NASA NTRS1999

A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

Xu, Kun↗

NASA NTRS1984

Parallel triangularization of substructured finite element problems

Much of the computational effort of the finite element process involves the solution of a system of linear equations. The coefficient matrix of this system, known as the global stiffness matrix, is symmetric, positive definite, and generally sparse. An important technique for reducing the time required to solve this system is substructuring or matrix partitioning. Substructuring is based on the idea of dividing a structure into pieces, each of which can then be analyzed relatively indepenently. As a result of this division, each point in the finite element discretization is either interior to a substructure or on a boundary between substructures. Contributions to the global stiffness matrix from connections between boundary points from the K(bb) matrix are reported. The triangularization of a general K(bb) matrix on a parallel machine is specifically discussed.

Leuze, M. R.↗

NASA NTRS1986

Enhancements to Sperry/NASTRAN

Reviewed is the enhancement to NASTRAN program performed by NUK (Nippon Univac Kaisha, Ltd.) added to Level 15.5. Features discussed include intermediate checkpoint-restart in triangular decomposition, I/O improvement, multibanked memory and new plate element. The first three improvements provide the capability to solve significantly large size problems, while the new elements release the analyst from the cumbersome work to constrain the singularities caused by the lack of stiffness of inplane rotation of old plate elements.

Koga, T.↗

NASA NTRS1988

A perturbation technique for parameter identification in distributed structures

Structures are often characterized by parameters, such as mass and stiffness, that are spatially distributed. Parameter identification of distributed structures is subject to many of the difficulties involved in the modeling problem, and the choice of the model can greatly affect the results of the parameter identification process. Analogously to control spillover in the control of distributed-parameter systems, identification spillover is shown to exist as well and its effect is to degrade the parameter estimates. Moreover, as in modeling by the Rayleigh-Ritz method, it is shown that, for a Rayleigh-Ritz type identification algorithm, an inclusion principle exists in the identification of distributed-parameter systems as well, so that the identified natural frequencies approach the actual natural frequencies monotonically from above.

Meirovitch, L.↗

NASA NTRS1988

A proposed method for enhanced eigen-pair extraction using finite element methods: Theory and application

The paper covers two distinct parts: theory and application. The goal of this work was the reduction of model size with an increase in eigenvalue/vector accuracy. This method is ideal for the condensation of large truss- or beam-type structures. The theoretical approach involves the conversion of a continuum transfer matrix beam element into an 'Exact' dynamic stiffness element. This formulation is implemented in a finite element environment. This results in the need to solve a transcendental eigenvalue problem. Once the eigenvalue is determined the eigenvectors can be reconstructed with any desired spatial precision. No discretization limitations are imposed on the reconstruction. The results of such a combined finite element and transfer matrix formulation is a much smaller FEM eigenvalue problem. This formulation has the ability to extract higher eigenvalues as easily and as accurately as lower eigenvalues. Moreover, one can extract many more eigenvalues/vectors from the model than the number of degrees of freedom in the FEM formulation. Typically, the number of eigenvalues accurately extractable via the 'Exact' element method are at least 8 times the number of degrees of freedom. In contrast, the FEM usually extracts one accurate (within 5 percent) eigenvalue for each 3-4 degrees of freedom. The 'Exact' element results in a 20-30 improvement in the number of accurately extractable eigenvalues and eigenvectors.

Jara-Almonte, J.↗

NASA NTRS1997

Modeling of Structural-Acoustic Interaction Using Coupled FE/BE Method and Control of Interior Acoustic Pressure Using Piezoelectric Actuators

A coupled finite element (FE) and boundary element (BE) approach is presented to model full coupled structural/acoustic/piezoelectric systems. The dual reciprocity boundary element method is used so that the natural frequencies and mode shapes of the coupled system can be obtained, and to extend this approach to time dependent problems. The boundary element method is applied to interior acoustic domains, and the results are very accurate when compared with limited exact solutions. Structural-acoustic problems are then analyzed with the coupled finite element/boundary element method, where the finite element method models the structural domain and the boundary element method models the acoustic domain. Results for a system consisting of an isotropic panel and a cubic cavity are in good agreement with exact solutions and experiment data. The response of a composite panel backed cavity is then obtained. The results show that the mass and stiffness of piezoelectric layers have to be considered. The coupled finite element and boundary element equations are transformed into modal coordinates, which is more convenient for transient excitation. Several transient problems are solved based on this formulation. Two control designs, a linear quadratic regulator (LQR) and a feedforward controller, are applied to reduce the acoustic pressure inside the cavity based on the equations in modal coordinates. The results indicate that both controllers can reduce the interior acoustic pressure and the plate deflection.

Mei, Chuh↗

NASA NTRS2002

The Role of Radial Clearance on the Performance of Foil Air Bearings

Load capacity tests were conducted to determine how radial clearance variations affect the load capacity coefficient of foil air bearings. Two Generation III foil air bearings with the same design but possessing different initial radial clearances were tested at room temperature against an as-ground PS304 coated journal operating at 30,000 rpm. Increases in radial clearance were accomplished by reducing the journal's outside diameter via an in-place grinding system. From each load capacity test the bearing load capacity coefficient was calculated from the rule-of-thumb (ROT) model developed for foil air bearings. The test results indicate that, in terms of the load capacity coefficient, radial clearance has a direct impact on the performance of the foil air bearing. Each test bearing exhibited an optimum radial clearance that resulted in a maximum load capacity coefficient. Relative to this optimum value are two separate operating regimes that are governed by different modes of failure. Bearings operating with radial clearances less than the optimum exhibit load capacity coefficients that are a strong function of radial clearance and are prone to a thermal runaway failure mechanism and bearing seizure. Conversely, a bearing operating with a radial clearance twice the optimum suffered only a 20 percent decline in its maximum load capacity coefficient and did not experience any thermal management problems. However, it is unknown to what degree these changes in radial clearance had on other performance parameters, such as the stiffness and damping properties of the bearings.

Radil, Kevin↗

NASA NTRS1993

Transverse cracking and stiffness reduction in composite laminates

A study of transverse cracking mechanism in composite laminates is presented using a singular hybrid finite element model. The model provides the global structural response as well as the precise local crack-tip stress fields. An elasticity basis for the problem is established by employing Lekhnitskii's complex variable potentials and method of eigenfunction expansion. Stress singularities associated with the transverse crack are obtained by decomposing the deformation into the symmetric and antisymmetric modes and proper boundary conditions. A singular hybrid element is thereby formulated based on the variational principle of a modified hybrid functional to incorporate local crack singularities. Axial stiffness reduction due to transverse cracking is studied. The results are shown to be in very good agreement with the existing experimental data. Comparison with simple shear lag analysis is also given. The effects of stress intensity factors and strain energy density on the increase of crack density are analyzed. The results reveal that the parameters approach definite limits when crack densities are saturated, an evidence of the existence of characteristic damage state.

Yuan, F. G.↗

NASA NTRS1978

NASTRAN implementation of an isoparametric doubly-curved quadrilateral shell element

A quadrilateral shell element, CQUAD4, was added to level 15.5 and subsequently to level 16.0 of NASTRAN. The element exhibited doubly curved surfaces and used biquadratic interpolation functions. Reduced integration techniques were used to improve the performance of the element in thin shell problems. The creation of several new bulk data items is discussed, along with a special module, GPNORM, to process SHLNORM bulk data cards. In addition to the theoretical basis for the element stiffness matrix, consistent mass and load matrices are presented. Several potential sources of degenerate behavior of the element were investigated. Guidelines for proper use of the element were suggested. Performance of the element on several widely published classical examples was demonstrated. The results showed a significant improvement over presently available NASTRAN shell elements for even the coarsest meshes. Potential applications to two classes of practical problems are discussed.

Potvin, A. B.↗

NASA NTRS1976

Comparison of supercritical and conventional wing flutter characteristics

A wind-tunnel study was undertaken to directly compare the measured flutter boundaries of two dynamically similar aeroelastic models which had the same planform, maximum thickness-to-chord ratio, and as nearly identical stiffness and mass distributions as possible, with one wing having a supercritical airfoil and the other a conventional airfoil. The considerations and problems associated with flutter testing supercritical wing models at or near design lift coefficients are discussed, and the measured transonic boundaries of the two wings are compared with boundaries calculated with a subsonic lifting surface theory.

Farmer, M. G.↗

NASA NTRS1986

Application of the method of lines for solution of the Navier-Stokes equations using nonuniform grid distributions

The method of lines is investigated for the numerical solution of the stream function and velocity form of the Navier-Stokes equations on nonuniform grids. Stiffness characteristics of a linear one-dimensional model equation are examined to establish the feasibility of applying the method to the volicity equation in two dimensions. It is observed that the transformed governing equations in computational domain become stiffer with increasing the number of grid points and the differencing technique affects the stiffness characteristics. The results demonstrate that the method of lines under certain circumstances is highly suitable for the numerical solution of physical problems on domains covered with variable grids.

Tiwari, S. N.↗

NASA NTRS1986

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of elastoplastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, the Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

NASA NTRS1991

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of an elasto-plastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

NASA NTRS1991

Controls of interaction dynamics of orbital assembly

Building structures and spacecraft in orbit will require technologies for positioning, docking/berthing, and joining orbital structures. A fundamental problem underlying the operation of docking and berthing is that of controlling the contact dynamics of mechanical structures actuated by active mechanisms such as robotic devices. Control systems must be designed to control these active mechanisms so that both the free space motions and contact motions are stable and satisfy specifications on position accuracy and bounds on contact forces. For the large orbital structures of the future, the problem of interactive dynamics and control is fundamentally different in several ways than it was for spacecraft docking in the past. First, future space structures must be treated as flexible structures - the operations of docking, berthing, and assembly will need to respect the vibrations of the structures. Second, the assembly of these structures will require multiple-point contact, rather than the essentially single-point positioning of conventional spacecraft docking. Third, some assembly operations require the subassemblies to be brought and held in contact so that successful joining can be accomplished. A preliminary study of contact stability and compliance control design has resulted in the development of an analytical method and a design method to analyze stability. The analytical method analyzes the problem of stability when an actively-controlled structure contacts a passive structure. This method makes it possible to accurately estimate the stiffness of the passive structures with which the contact motion will become unstable. The analytic results suggest that passivity is neither achievable in practice, nor necessary as a design concept. A contact control system need only be passive up to a certain frequency; beyond that frequency the system can be stabilized with sufficiently small gains. With this concept the Center developed a design methodology for achieving desired compliant contact motions. This design method is based on H-infinity norm optimization, which makes it possible to consider both driving point mechanical impedance and systems robustness to modeling uncertainty. A laboratory facility was set up to verify experimentally the analytical and design theory.

Su, Renjeng↗