SEARCH · Engineering Papers
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.
GAMNAS (Geometric and Material Nonlinear Analysis of Structures) is twodimensional finite-element stress-analysis program supporting fracturemechanics studies of debonding and delamination. GAMNAS options include linear, geometricnonlinear, material-nonlinear, and combined geometric- and material-nonlinear analysis.
A technique which employs both linear and nonlinear methods in a multilevel optimization structure to best approximate generalized unsteady aerodynamic forces for arbitrary motion is described. Optimum selection of free parameters is made in a rational function approximation of the aerodynamic forces in the Laplace domain such that a best fit is obtained, in a least squares sense, to tabular data for purely oscillatory motion. The multilevel structure and the corresponding formulation of the objective models are presented which separate the reduction of the fit error into linear and nonlinear problems, thus enabling the use of linear methods where practical. Certain equality and inequality constraints that may be imposed are identified; a brief description of the nongradient, nonlinear optimizer which is used is given; and results which illustrate application of the method are presented.
A nonlinear filter configured to provide phase lead without accompanying gain distortion is analyzed and evaluated. The nonlinear filter is superior to a linear lead/lag compensator in its ability to maintain system stability as open loop crossover frequency is increased. Test subjects subjectively rated the filter as slightly better than a lead/lag compensator in its ability to compensate for delays in a compensatory tracking task. However, the filter does introduce unwanted harmonics. This is particularly noticeable for low frequency pilot inputs. A revised compensation method is proposed which allows such low frequency inputs to bypass the nonlinear filter. A brief analytical and experimental evaluation of the revised filter indicates that further evaluation in more realistic tasks is justified.
The instability of a three dimensional attachment line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time dependent Navier-Stokes equations for the attachment line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment line boundary layer is also investigated.
The fluctuating field of a jet excited by transient mass injection is simulated numerically. The model is developed by expanding the state vector as a mean state plus a fluctuating state. Nonlinear terms are not neglected, and the effect of nonlinearity was studied. A high order numerical method is used to compute the solution. The results show a significant spectral broadening in the flow field due to the nonlinearity. In addition, large scale structures are broken down into smaller scales.
A computational capability for predicting the nonlinear thermomechanical structural response of fiber-reinforced superalloy (FRS) turbine blades is described. This capability is embedded in a special purpose computer code (COBSTRAN) developed at the NASA Lewis Research Center. Special features of this computational capability include accounting for: fiber/matrix reaction, nonlinear and anisotropic material behavior, complex stress distribution due to local and global heterogeneity, and residual stresses due to initial fabrication and/or inelastic behavior during subsequent missions. Numerical results are presented from analyses of a hypothetical FRS turbine blade subjected to a fabrication process and subsequent mission cycle. The results demonstrate the capabilities of this computational tool to; predict local stress/strain response and capture trends of local nonlinear and anisotropic material behavior, relate the effects of this local behavior to the global response of a multilayered fiber-composite turbine blade, and trace material history from fabrication through successive missions.
Nonlinear beam kinematics are developed and applied to the dynamic analysis of a pretwisted, rotating beam element. The common practice of assuming moderate rotations caused by structural deformation in geometric nonlinear analyses of rotating beams was abandoned in the present analysis. The kinematic relations that described the orientation of the cross section during deformation are simplified by systematically ignoring the extensional strain compared to unity in those relations. Open cross section effects such as warping rigidity and dynamics are ignored, but other influences of warp are retained. The beam cross section is not allowed to deform in its own plane. Various means of implementation are discussed, including a finite element formulation. Numerical results obtained for nonlinear static problems show remarkable agreement with experiment.
Nonlinear constitutive relations for high temperature applications were discussed. The state of the art in nonlinear constitutive modeling of high temperature materials was reviewed and the need for future research and development efforts in this area was identified. Considerable research efforts are urgently needed in the development of nonlinear constitutive relations for high temperature applications prompted by recent advances in high temperature materials technology and new demands on material and component performance. Topics discussed include: constitutive modeling, numerical methods, material testing, and structural applications.
The nonlinear vibration response of a double cantilevered beam subjected to pulse loading over a central sector is studied. The initial response is generated in detail to ascertain the energetics of the response. The total energy is used as a gauge of the stability and accuracy of the solution. It is shown that to obtain accurate and stable initial solutions an extremely high spatial and time resolution is required. This requirement was only evident through an examination of the energy of the system. It is proposed, therefore, to use the total energy of the system as a necessary stability and accuracy criterion for the nonlinear response of conservative systems. The results also demonstrate that even for moderate nonlinearities, the effects of membrane forces have a significant influence on the system. It is also shown that while the fundamental response is contained in a first mode envelope, the fluctuations caused by the higher order modes must be resolved.
Nonlinear strain displacement relations for three-dimensional elasticity are determined in orthogonal curvilinear coordinates. To develop a two-dimensional theory, the displacements are expressed by trigonometric series representation through-the-thickness. The nonlinear strain-displacement relations are expanded into series which contain all first and second degree terms. In the series for the displacements only the first few terms are retained. Insertion of the expansions into the three-dimensional virtual work expression leads to nonlinear equations of equilibrium for laminated and thick plates and shells that include the effects of transverse shearing. Equations of equilibrium and buckling equations are derived for flat plates and cylindrical shells. The shell equations reduce to conventional transverse shearing shell equations when the effects of the trigonometric terms are omitted and to classical shell equations when the trigonometric terms are omitted and the shell is assumed to be thin.
It is pointed out that the theory of spiral density waves, invented to explain the spiral structure of disk galaxies, has also been found useful for the study of planetary rings. The linear theory is by now well developed, while the nonlinear theory is less complete. Analytical calculations which include self-gravitation have, so far, obtained results only in the slightly nonlinear regime, or have concentrated on partial effects which are not of primary importance to the physical problem at hand. In the present paper, it is attempted to remedy these shortcomings. The simplest asymptotic ordering which can still yield useful results is adopted. Attention is given to the reduction to a nonlinear integral equation in a single variable, the use of the Wentzel-Kramers-Brillouin-Jeffreys theory, and the replacement of an equation by another which is easier to handle numerically.
A method for designing an automatic flight controller for short and vertical takeoff aircraft is presently being developed at NASA Ames Research Center. This technique involves transformations of nonlinear systems to controllable linear systems and takes into account the nonlinearities of the aircraft. In general, the transformations cannot always be given in closed form. Using partial differential equations, an approximate linear system, called the modified tangent model, was recently introduced. A linear transformation of this tangent model to Brunovsky canonical form can be constructed, and from this an approximation of an exact transformation for the nonlinear system can be found. It is shown that a canonical expansion in Lie brackets about the point x(0) yields the same modified tangent model.
Powerful methods are now available for solving linear parametric inverse problems. However, many inverse problems which arise in geohysics are nonlinear. Fortunately, it is possible to treat most of these with the air of linear perturbation theory and liner inversion. But a convenient method is needed for assessing the importance of nonlinearity in these quasi-linear problems. The present paper provides such a method. Matsu'ura and Jackson (1984) have presented a simple algorithm for evaluating the asymptotic covariance matrix fo estimation errors. In the present investigation, aspects of linear inversion are discussed, taking into account linear parametric inverse problems, nonuniqueness, prior information, confidence limits, conditional and marginal statistics, the relative importance of the prior and observational data, and standardized variables. Attention is also given to nonlinear inversion, and the application of the considered approaches to a number of examples.
An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. With the fitted surface, sensitivity information can be updated at each iteration with less computational effort than that required by either a finite-difference method or integration of the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, and thus provides flexibility to use model equations in any convenient format. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. The degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels and to predict the degree of agreement between CR bounds and search estimates.
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.
The following problem is examined: given a nonlinear control system dot-x(t) = f(x(t)) + the sum to m terms(i=1) u sub i (t)g sub i (x(t)) on R(n) and a point x(0) in R(n), approximate the system near x(0) by a linear system. One approach is to use the usual Taylor series linearization. However, the controllability properties of both the nonlinear and linear systems depend on certain Lie brackets of the vector field under consideration. This suggests that a linear approximation based on Lie bracket matching should be constructed at x(0). In general, the linearizations based on the Taylor method and the Lie bracket approach are different. However, under certain mild assumptions, it is shown that there is a coordinate system for R(n) near x(0) in which these two types of linearizations agree. The importance of this agreement is indicated by examining the time responses of the nonlinear system and its linear approximation and comparing the lower order kernels in Volterra expansions of each.
The concepts of transformation and canonical form have been used in analyzing linear systems. These ideas are extended to nonlinear systems. A coordinate system and a corresponding canonical form are developed for general nonlinear control systems. Their usefulness is demonstrated by showing that every feedback linearizable system becomes a system with only feedback paths in the canonical form. For control design involving a nonlinear system, one approach is to put the system in its canonical form and approximate by that part having only feedback paths.
An analytical method is described for computing the nonlinear thermal and mechanical response of laminated plates. The material model focuses upon the behavior of metal matrix materials by relating the nonlinear composite response to plasticity effects in the matrix. The foundation of the analysis is the unidirectional material model which is used to compute the instantaneous properties of the lamina based upon the properties of the fibers and matrix. The unidirectional model assumes that the fibers properties are constant with temperature and assumes that the matrix can be modelled as a temperature dependent, bilinear, kinematically hardening material. An incremental approach is used to compute average stresses in the fibers and matrix caused by arbitrary mechanical and thermal loads. The layer model is incorporated in an incremental laminated plate theory to compute the nonlinear response of laminated metal matrix composites of general orientation and stacking sequence. The report includes comparisons of the method with other analytical approaches and compares theoretical calculations with measured experimental material behavior. A section is included which describes the limitations of the material model.