Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Frequency based localization of structural discrepancies

The intent of modal analysis is to develop a reliable model of a structure by working with the analytical and experimental modal properties of frequency, damping and mode shape. In addition to identifying these modal properties, it would be desirable to determine spatially which parts of the structure are modelled poorly or well. It is shown how the pattern of discrepancies in the analytical and experimental test values for the pole and the driving point zero frequencies of a structure can be linked to discrepancies in the mass or stiffness of the structural elements. The success of the procedure depends on the numerical conditioning of a modal reference matrix. Strategies to insure adequate numerical conditioning require a formulation which avoids geometric and energy storage symmetries of the structure, and ignores structural elements which contribute negligibly small potential or kinetic energy to the excited modes. Physical insight into the numerical conditioning problem is provided by a numerical example and by localization of a mass discrepancy in a real structure based on lab tests.

Shepard, G. D.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS1980

Straightening of a wavy strip: An elastic-plastic contact problem including snap-through

The nonlinear behavior of a wave like deformed metal strip during the levelling process were calculated. Elastic-plastic material behavior as well as nonlinearities due to large deformations were considered. The considered problem lead to a combined stability and contact problem. It is shown that, despite the initially concentrated loading, neglecting the change of loading conditions due to altered contact domains may lead to a significant error in the evaluation of the nonlinear behavior and particularly to an underestimation of the stability limit load. The stability was examined by considering the load deflection path and the behavior of a load-dependent current stiffness parameter in combination with the determinant of the current stiffness matrix.

Fischer, D. F.↗
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
NASA NTRS1993

Flight service evaluation of composite components on the Bell helicopter model 206L

This is the final report on the advanced composite components which were placed in service on the 206L LongRanger helicopters in the continental United States, Canada, and Alaska. This report covers all test data which was gathered, as well as maintenance histories of the parts. The previous reports describe the fabrication, service experiences, and test data through 1986. This report contains information from these references, as well as data gathered after 1986. The status of the 40 sets of components is discussed. Each set consisted of a vertical fin, forward fairing, litter door, and baggage door. Almost 500,000 flight hours were accumulated on the 160 parts, with the high-time part accumulating 14,687 flight hours. Over 60 percent of the parts were destructively tested to measure strength and stiffness retention over the course of the program. The vertical fins had the greatest strength retention followed by the litter doors. The baggage doors had the poorest retention of strength. There was very little difference in property retention between the four primary operating regions: Northwest U.S., Southwest U.S., Gulf of Mexico Coastal Region, and the Northeast U.S. and Eastern Canada Region. The field problems have ranged from two lightning-struck fins to significant delaminations in the baggage doors. There was only one environmentally related field incident, in which the glass windows on the litter doors were found to loosen due to high temperatures experienced in the southwest region.

Wilson, Henry↗
NASA NTRS1997

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗
NASA NTRS1976

On the use of two hardening rules of plasticity in incremental and pseudo force analysis

The tangent stiffness and pseudo force forms of the equations of motion are first derived within the context of a total Lagrangian formulation. After a brief discussion of available incremental theory plasticity models, the small strain formulations and computational procedures of the mechanical sublayer model and combined kinematic-isotropic hardening as used in the general purpose structural analysis program AGGIE I are presented. Several sample problems are then presented along with recommended guidelines for use of the two plasticity models.

Hunsaker, B., Jr.↗
NASA NTRS1993

Computational reacting gas dynamics

In the study of high speed flows at high altitudes, such as that encountered by re-entry spacecrafts, the interaction of chemical reactions and other non-equilibrium processes in the flow field with the gas dynamics is crucial. Generally speaking, problems of this level of complexity must resort to numerical methods for solutions, using sophisticated computational fluid dynamics (CFD) codes. The difficulties introduced by reacting gas dynamics can be classified into three distinct headings: (1) the usually inadequate knowledge of the reaction rate coefficients in the non-equilibrium reaction system; (2) the vastly larger number of unknowns involved in the computation and the expected stiffness of the equations; and (3) the interpretation of the detailed reacting CFD numerical results. The research performed accepts the premise that reacting flows of practical interest in the future will in general be too complex or 'untractable' for traditional analytical developments. The power of modern computers must be exploited. However, instead of focusing solely on the construction of numerical solutions of full-model equations, attention is also directed to the 'derivation' of the simplified model from the given full-model. In other words, the present research aims to utilize computations to do tasks which have traditionally been done by skilled theoreticians: to reduce an originally complex full-model system into an approximate but otherwise equivalent simplified model system. The tacit assumption is that once the appropriate simplified model is derived, the interpretation of the detailed numerical reacting CFD numerical results will become much easier. The approach of the research is called computational singular perturbation (CSP).

Lam, S. H.↗
NASA NTRS1988

Precision pointing of scientific instruments on space station: The LFGGREC perspective

An application of Lyapunov function-gradient-generated robustness-enhancing control (LFGGREC) is explored. The attention is directed to a reduced-complexity representation of the pointing problem presented by the system composed of the Space Infrared Telescope Facility gimbaled to a space station configuration. Uncertainties include disturbance forces applied in the crew compartment area and control moments applied to adjacent scientific payloads (modeled as disturbance moments). Also included are uncertainties in gimbal friction and in the structural component of the system, as reflected in the inertia matrix, the damping matrix, and the stiffness matrix, and the effect of the ignored vibrational dynamics of the structure. The emphasis is on the adaptation of LFGGREC to this particular configuration and on the robustness analysis.

Blackwell, C. C.↗
NASA NTRS1989

Shape control of high degree-of-freedom variable geometry trusses

Common static trusses are constrained to permit no relative motion between truss elements. A Variable Geometry Truss (VGT), however, is a truss which contains some number of variable length links. The extensible links allow the truss to change shape in a precise, controllable manner. These changes can also be used to control the vibrational response of a truss structure or to perform robotic tasks. Many geometric configurations, both planar and spatial, are possible candidates for VGT manipulators. Here, only two geometries are discussed; the three degree-of-freedom (DOF) spatial octahedral/octahedral truss and the three DOF planar tetrahedral truss. These truss geometries are used as the fundamental element in a repeating chain of trusses. This results in a highly dexterous manipulator with perhaps 30 to 60 degrees of freedom that retains the favorable stiffness properties of a conventional truss. From a fixed base, this type of manipulator could perform shape or vibration control while extending and snaking through complex passageways or moving around obstacles to perform robotic tasks. The approach taken here is to first concentrate on fully understanding the forward and inverse kinematics of the fundamental elements and then utilizing the insight thus gained to solve the more complex problem of the kinematic chains.

Salerno, R. J.↗
NASA NTRS1991

Control issues of microgravity vibration isolation

Active vibration isolation systems contemplated for microgravity space experiments may be designed to reach given performance requirements in a variety of ways. An analogy to passive isolation systems proves to be illustrative but lacks the flexibility as a design tool of a control systems approach and may lead to poor designs. For example, it is shown that a focus on equivalent stiffness in isolation system design leads to a controller that sacrifices robustness for performance. Control theory as applied to vibration isolation is reviewed and passive analogies are discussed. The loop shaping trade-off is introduced and used to design a single-degree-of-freedom fedback controller. An algebraic control design methodology is contrasted to loop shaping and critiqued. Multi-axis vibration isolation and the problems of decoupled single loop control are introduced through a two-degree-of-freedom example problem. It is shown that center of mass uncertainty may result in instability when decoupled single loop control is used. This results from the ill-conditioned nature of the feedback control design. The use of the Linear Quadratic Regulator synthesis procedure for vibration isolation controller design is discussed.

Knospe, C. R.↗
NASA NTRS1984

Stress intensity factors in a reinforced thick-walled cylinder

An elastic thick-walled cylinder containing a radial crack is considered. It is assumed that the cylinder is reinforced by an elastic membrane on its inner surface. The model is intended to simulate pressure vessels with cladding. The formulation of the problem is reduced to a singular integral equation. Various special cases including that of a crack terminating at the cylinder-reinforcement interface are investigated and numerical examples are given. Results indicate that in the case of the crack touching the interface the crack surface displacement derivative is finite and consequently the stress state around the corresponding crack tip is bounded; and generally, for realistic values of the stiffness parameter, the effect of the reinforcement is not very significant.

Tang, R.↗
NASA NTRS2016

Multifunctional Hybrid Carbon Nanotube/Carbon Fiber Polymer Composites

For aircraft primary structures, carbon fiber reinforced polymer (CFRP) composites possess many advantages over conventional aluminum alloys due to their light weight, higher strengthand stiffness-to-weight ratio, and low life-cycle maintenance costs. However, the relatively low electrical and thermal conductivities of CFRP composites fail to provide structural safety in certain operational conditions such as lightning strikes. Despite several attempts to solve these issues with the addition of carbon nanotubes (CNT) into polymer matrices, and/or by interleaving CNT sheets between conventional carbon fiber (CF) composite layers, there are still interfacial problems that exist between CNTs (or CF) and the resin. In this study, hybrid CNT/CF polymer composites were fabricated by interleaving layers of CNT sheets with Hexcel® IM7/8852 prepreg. Resin concentrations from 1 wt% to 50 wt% were used to infuse the CNT sheets prior to composite fabrication. The interlaminar properties of the resulting hybrid composites were characterized by mode I and II fracture toughness testing (double cantilever beam and end-notched flexure test). Fractographical analysis was performed to study the effect of resin concentration. In addition, multi-directional physical properties like thermal conductivity of the orthotropic hybrid polymer composite were evaluated. Interleaving CNT sheets significantly improved the in-plane (axial and perpendicular direction of CF alignment) thermal conductivity of the hybrid composite laminates by 50 - 400%.

Kang, Jin Ho↗
NASA NTRS1971

Statistical analysis of the modal properties of large structural systems.

A theory is developed to predict eigenvalue and eigenvector uncertainty in large dynamic models. The uncertainty is based on physical property uncertainty and should not be confused with numerical roundoff, although the method can be extended to include the latter. The theory, when implemented on a computer, is used to analyze the uncertainties in frequencies and mode shapes based on uncertainties in mass, stiffness, modulus of elasticity, etc. The method incorporates a linear statistical model which is quite adequate for handling property uncertainties of 10% or more. The model is not limited to small systems but uses certain statistical assumptions as well as selective matrix manipulations to keep the size of all matrix operations to within the number of degrees of freedom of the system. Examples are given for two longitudinal vibration problems, and the results are supported by a Monte Carlo analysis.

Collins, J. D.↗
NASA NTRS1983

The effect of matrix properties and fiber properties on impact failure mechanics

The low-velocity impact problem in graphite/epoxy composite sheets must be solved before large amounts of that material can be used in commercial aircraft. Many of the low-velocity impacts that affect aircraft parts occur during normal ground operations and maintenance. Service equipment and tools have masses above 1 kg, and at velocities of less than 3 m/s can impact structural parts with energies higher than composites can endure without degradation of stiffness or strength. Simple solutions were developed for large-mass, low-velocity impacts which can be modeled as quasi-static events. Static test data and impact data show that the fiber properties control the impact energy which can be absorbed before penetration. Matrix shear strength and peel resistance control the extent of delamination. Comparison of results from tough matrix and brittle matrix composites show that although tough matrices reduce the extent of delamination, they lead to more fiber damage in the contact area.

Elber, W.↗
NASA NTRS1985

An analysis of cross-coupling of a multicomponent jet engine test stand using finite element modeling techniques

A two axis thrust measuring system was analyzed by using a finite a element computer program to determine the sensitivities of the thrust vectoring nozzle system to misalignment of the load cells and applied loads, and the stiffness of the structural members. Three models were evaluated: (1) the basic measuring element and its internal calibration load cells; (2) the basic measuring element and its external load calibration equipment; and (3) the basic measuring element, external calibration load frame and the altitude facility support structure. Alignment of calibration loads was the greatest source of error for multiaxis thrust measuring systems. Uniform increases or decreases in stiffness of the members, which might be caused by the selection of the materials, have little effect on the accuracy of the measurements. It is found that the POLO-FINITE program is a viable tool for designing and analyzing multiaxis thrust measurement systems. The response of the test stand to step inputs that might be encountered with thrust vectoring tests was determined. The dynamic analysis show a potential problem for measuring the dynamic response characteristics of thrust vectoring systems because of the inherently light damping of the test stand.

Schweikhard, W. G.↗
NASA NTRS1989

Exploiting symmetries in the modeling and analysis of tires

A computational procedure is presented for reducing the size of the analysis models of tires having unsymmetric material, geometry and/or loading. The two key elements of the procedure when applied to anisotropic tires are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic parts; and (2) successive application of the finite-element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The proposed technique has high potential for handling practical tire problems with anisotropic materials, unsymmetric imperfections and asymmetric loading. It is also particularly useful for use with three-dimensional finite-element models of tires.

Noor, Ahmed K.↗
NASA NTRS1989

Integrated aerodynamic-structural design of a transport wing

The integrated aerodynamic-structural design of a subsonic transport wing for minimum weight subject to required range is formulated and solved. The problem requires large computational resources, and two methods are used to alleviate the computational burden. First, a modular sensitivity method that permits the usage of black-box disciplinary software packages, is used to reduce the cost of sensitivity derivatives. In particular, it is shown that derivatives of the aeroelastic response and divergence speed can be calculated without the costly computation of derivatives of aerodynamic influence coefficient and structural stiffness matrices. A sequential approximate optimization is used to further reduce computational cost. The optimization procedure is shown to require a relatively small number of analysis and sensitivity calculations.

Grossman, B.↗