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.
Nonlinear behavior of phase-locked loops with rapidly varying phase error is examined by using computer phase-plane analysis. The phase variation is modeled by a sinusoidal function. Threshold loop parameters are presented for both sinusoidal and sawtooth phase comparators.
A comprehensive development of nondimensional parameters and equations for nonlinear and bifurcations analyses of quasi-shallow shells, based on the Donnell-Mushtari-Vlasov theory for thin anisotropic shells, is presented. A complete set of field equations for geometrically imperfect shells is presented in terms general of lines-of-curvature coordinates. A systematic nondimensionalization of these equations is developed, several new nondimensional parameters are defined, and a comprehensive stress-function formulation is presented that includes variational principles for equilibrium and compatibility. Bifurcation analysis is applied to the nondimensional nonlinear field equations and a comprehensive set of bifurcation equations are presented. An extensive collection of tables and figures are presented that show the effects of lamina material properties and stacking sequence on the nondimensional parameters.
An easily computerized analytical method is developed for identifying critical airplane maneuvers in which nonlinear rotational coupling effects may cause sudden jumps in the response to pilot's control inputs. Fifth and ninth degree polynomials for predicting multiple pseudo-steady states of roll-coupled maneuvers are derived. The program calculates the pseudo-steady solutions and their stability. The occurrence of jump-like responses for several airplanes and a variety of maneuvers is shown to correlate well with the appearance of multiple stable solutions for critical control combinations. The analysis is extended to include aerodynamics nonlinear in angle of attack.
An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.
A single-engine, propeller-driven, general-aviation model was incorporated into the nonlinear simulation and into the linear analysis of root loci and frequency response. Full-scale wind tunnel data provided its aerodynamic model, and the thrust model included the airspeed dependent effects of power and propeller efficiency. Also, the parameters of the Jet Transport model were changed to correspond more closely to the Boeing 727. In order to study their effects on steady-state repsonse to vertical wind inputs, altitude and total specific energy (air-relative and inertial) feedback capabilities were added to the nonlinear and linear models. Multiloop system design goals were defined. Attempts were made to develop controllers which achieved these goals.
Bearing-load response for a pin-loaded hole is studied within the context of two-dimensional finite element analyses. Pin-loaded-hole configurations are representative of mechanically connected structures, such as a stiffener fastened to a rib of an isogrid panel, that are idealized as part of a larger structural component. Within this context, the larger structural component may be idealized as a two-dimensional shell finite element model to identify load paths and high stress regions. Finite element modeling and analysis aspects of a pin-loaded hole are considered in the present paper including the use of linear and nonlinear springs to simulate the pin-bearing contact condition. Simulating pin-connected structures within a two-dimensional finite element analysis model using nonlinear spring or gap elements provides an effective way for accurate prediction of the local effective stress state and peak forces.
The present paper takes the basic mathematical approach to analyze the stability of Riccati equations. First, a simple stability criterion like that of the linear constant coefficient system is sought for the nonlinear constant Riccati Equation. Then, a complete analysis on the limiting covariance P follows. The analysis is focussed on the main theorem. Finally, the computational aspects of the Riccati equations are analyzed with examples. The application of this analysis is to provide a simple stability criterion for the Kalman filter computation to test and predict filter convergence and thus modify filter parameters accordingly.
Nonlinear equations of motion are derived for the rim inertial measuring system, which is mounted in a strapped down configuration on a carrier vehicle. The measurement and equations derived from angular rates and linear accelerations were interfaced with the dynamic model of a user defined carried vehicle. Two methods are presented for rate measurement, and the results of the nonlinear computer simulation are presented.
Two computer programs are described for predicting the transient large deflection elastic viscoplastic responses of thin single layer, initially flat unstiffened or integrally stiffened, Kirchhoff-Lov ductile metal panels. The PLATE 1 program pertains to structural responses produced by prescribed externally applied transient loading or prescribed initial velocity distributions. The collision imparted velocity method PLATE 1 program concerns structural responses produced by impact of an idealized nondeformable fragment. Finite elements are used to represent the structure in both programs. Strain hardening and strain rate effects of initially isotropic material are considered.
An analysis is presented which examines the modulation of different instability modes satisfying the triad resonance condition in time and space in a three-dimensional boundary layer flow. Detuning parameters are used for the wave numbers and the frequencies. The nonparallelism of the mean flow is taken into account in the analysis. At the leading-edge region of an infinite swept wing, different resonant triads are investigated that are comprised of travelling crossflow, vertical vorticity and Tollmein-Schlichting modes. The spatial evolution of the resonating triad components are studied.
The present analysis of free liquid drops' inviscid oscillations proceeds through solution of Bernoulli's equation to obtain the free surface shape and of Laplace's equation for the velocity potential field. Results thus obtained encompass drop-shape sequences, pressure distributions, particle paths, and the temporal evolution of kinetic and surface energies; accuracy is verified by the near-constant drop volume and total energy, as well as the diminutiveness of mass and momentum fluxes across drop surfaces. Further insight into the nature of oscillations is provided by Fourier power spectrum analyses of mode interactions and frequency shifts.
A solution procedure for solving nonlinear time-marching problems is presented. The nonsymmetric systems of equations arising from a Newton-type linearization of these time-marching problems are solved using an iterative strategy based on the generalized minimal residual (GMRES) algorithm. Matrix-free techniques leading to reduction in storage are presented. Incorporation of a linesearch algorithm in the Newton-GMRES scheme is discussed. An automatic time-increment control strategy is developed to increase the stability of the time-marching process. High-speed flow computations demonstrate the effectiveness of these algorithms.
We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.
A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.
During conceptual design studies of advanced aircraft, the usual practice is to use linear theory to calculate the aerodynamic characteristics of candidate rigid (nonflexible) geometric external shapes. Recent developments and improvements in computational methods, especially computational fluid dynamics (CFD), provide significantly improved capability to generate detailed analysis data for the use of all disciplines involved in the evaluation of a proposed aircraft design. A multidisciplinary application of such analysis methods to calculate the effects of nonlinear aerodynamics and static aeroelasticity on the mission performance of a fighter aircraft concept is described. The aircraft configuration selected for study was defined in a previous study using linear aerodynamics and rigid geometry. The results from the previous study are used as a basis of comparison for the data generated herein. Aerodynamic characteristics are calculated using two different nonlinear theories, potential flow and rotational (Euler) flow. The aerodynamic calculations are performed in an iterative procedure with an equivalent plate structural analysis method to obtain lift and drag data for a flexible (nonrigid) aircraft. These static aeroelastic data are then used in calculating the combat and mission performance characteristics of the aircraft.
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.
Analysis based on the Jeffcott model is presented to explain 1/2 speed and 1/3 speed whirling motion occurring in rotors which are subject to periodic normal-loose or normal-tight radial stiffness variations. The normal-loose stiffness variation results due to bearing-clearance effects, while normal-tight stiffness variations result from rubbing over a portion of a rotor's orbit. The results demonstrate that 1/2 speed subharmonic motion can be explained as either a linear parametric-excitation phenomenon or as a stable nonlinear subharmonic motion. The 1/3 speed motion is shown to be possible due to the radial stiffness nonlinearity. A linear parametric-excitation analysis demonstrates that during a normal-tight rubbing condition, Coulomb damping significantly widens the potential range of unstable speeds.
A three-dimensional, full potential flow analysis with realistic hover wake geometry is presented for the prediction of hovering rotor airloads. The method of analysis is based on the concept of matching inner and outer domain solutions in three dimensions. The inner domain nonlinear solution is obtained using a finite difference analysis and the outer domain solution is based on prescribed wake methodology. This formulation which includes three-dimensional wake influence, was initially validated using a fixed-wing analysis, and has been extended to hovering rotor flight. Detailed chordwise and spanwise loading results are compared with subsonic and transonic test results from two rotor configurations to illustrate the predictive capabilities of the analysis. The extension of the method to steady-level forward flight is also discussed.