Asymmetric nonlinear dynamic response and buckling of shallow spherical shells
Numerical method and computer program for computing deflection response of thin shells
SEARCH · Search NASA
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.
Numerical method and computer program for computing deflection response of thin shells
Estimates of climatological rainfall are made as part of the validation program for the Tropical Rainfall Measuring Mission. Radar and raingage data collected during AMEX in Darwin, Australia during 1987 and 1988 are analyzed using a nonlinear merging methodology known as disjunctive cokriging (Azimi-Zonooz et al., 1989). The methodology is outlined and the results are used to determine the most advantageous time scale for estimation.
The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.
Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.
The accuracy of chlorophyll estimates by ocean color algorithms is affected by the variability of particulate attenuation; the presence of dissolved organic matter; and the nonlinear inverse relationship between the attenuation coefficient, K, and chlorophyll. Data collected during the Warm Core Rings Program were used to model the downwelling light field and determine the impact of these errors. A possible mechanism for the nonlinearity of K and chlorophyll is suggested; namely, that changing substrate from nitrate-nitrogen to ammonium causes enhanced blue absorption by photosynthetic phytoplankton in oligotrophic surface waters.
A full-scale crash test of the Sikorsky Advanced Composite Airframe Program (ACAP) helicopter was performed in 1999 to generate experimental data for correlation with a crash simulation developed using an explicit nonlinear, transient dynamic finite element code. The airframe was the residual flight test hardware from the ACAP program. For the test, the aircraft was outfitted with two crew and two troop seats, and four anthropomorphic test dummies. While the results of the impact test and crash simulation have been documented fairly extensively in the literature, the focus of this paper is to present the detailed occupant response data obtained from the crash test and to correlate the results with injury prediction models. These injury models include the Dynamic Response Index (DRI), the Head Injury Criteria (HIC), the spinal load requirement defined in FAR Part 27.562(c), and a comparison of the duration and magnitude of the occupant vertical acceleration responses with the Eiband whole-body acceleration tolerance curve.
A full-scale crash test of the Sikorsky Advanced Composite Airframe Program (ACAP) helicopter was performed in 1999 to generate experimental data for correlation with a crash simulation developed using an explicit nonlinear, transient dynamic finite element code. The airframe was the residual flight test hardware from the ACAP program. For the test, the aircraft was outfitted with two crew and two troop seats, and four anthropomorphic test dummies. While the results of the impact test and crash simulation have been documented fairly extensively in the literature, the focus of this paper is to present the detailed occupant response data obtained from the crash test and to correlate the results with injury prediction models. These injury models include the Dynamic Response Index (DRI), the Head Injury Criteria (HIC), the spinal load requirement defined in FAR Part 27.562(c), and a comparison of the duration and magnitude of the occupant vertical acceleration responses with the Eiband whole-body acceleration tolerance curve.
This report documents an effort to develop practical and accurate methods for estimating the fatigue lives of complex aerospace structures subjected to intense random excitations. The emphasis of the current program is to construct analytical schemes for performing fatigue life estimates for structures that exhibit nonlinear vibration behavior and that have numerous resonant modes contributing to the response.
The nonlinear dynamics of propellant sloshing during orbital ascent are usually neglected in the flight control analysis of large boost vehicles under the assumption that the viscous damping of the fluid is sufficient to suppress nonlinear phenomena and confine the fluid to small, planar free surface displacements. In this case, the sloshing dynamics can be modeled using a spring-mass-damper or linearized pendulum mechanical analog. However, large, smooth-wall tanks without significant internal hardware or ring baffles are still susceptible to nonlinear effects. In particular, rotary sloshing can present a risk to flight control as it involves the formation of a stable limit cycle which can lead to undesirable roll coupling. The underlying phenomena of jump resonance does not manifest in linear models, but can be reproduced using a nonlinear spherical pendulum or the Bauer paraboloid model developed during the Apollo/Saturn program. In this paper, a detailed analysis of the rotary sloshing dynamics of these mechanical analogs is presented, and discussed in the context of flight control stability. High-fidelity simulations of a representative boost vehicle are used to verify the semi-analytical predictions of the nonlinear dynamic response.
In order to understand the state of stress and strain in a typical balloon fabricated from thin polyethylene film, experiment data in the literature reviewed. It was determined that the film behaves as a nonlinear viscoelasticity material and should be characterized accordingly. A simple uniaxial, nonlinear viscoelastic model was developed for predicting stress given a certain strain history. The simple model showed good qualitative agreement with results of constant rate, uniaxial accurately predicting stresses for cyclic strain histories typical of balloon flights. A program was outlined which will result in the development of a more complex nonlinear viscoelastic model.
The development is discussed of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors with partial length acoustic liners, injector face baffles, and nonlinear combustion responses. Three primary objectives included: (1) extension of the Green's Function integral-iteration technique previously developed so that nonlinear three dimensional wave propagation effects could be included in stability predictions; (2) development of an analytical technique and a computer program for the prediction of the damping effects caused by injector face baffles; and (3) inclusion of a nonlinear, wave distortation dependent, combustion response model in the global stability analysis and determination of its effect on combustor stability.
The theoretical development of the forced steady state analysis of the structural dynamic response of a turbine engine having nonlinear connecting elements is discussed. Based on modal synthesis, and the principle of harmonic balance, the governing relations are the compatibility of displacements at the nonlinear connecting elements. There are four displacement compatibility equations at each nonlinear connection, which are solved by iteration for the principle harmonic of the excitation frequency. The resulting computer program, TETRA 2, combines the original TETRA transient analysis (with flexible bladed disk) with the steady state capability. A more versatile nonlinear rub or bearing element which contains a hardening (or softening) spring, with or without deadband, is also incorporated.
A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.
The potential role of symbolic manipulation programs in development and solution of the governing equations for rotorcraft dynamics problems is discussed and illustrated. Nonlinear equations of motion for a helicopter rotor blade represented by a rotating beam are developed making use of the computerized symbolic manipulation program MACSYMA. The use of computerized symbolic manipulation allows the analyst to concentrate on more meaningful tasks, such as establishment of physical assumptions, without being sidetracked by the tedious and trivial details of the algebraic manipulations. Furthermore, the resulting equations can be produced, if necessary, in a format suitable for numerical solution. A perturbation-type solution for the resulting dynamical equations is shown to be possible with a combination of symbolic manipulation and standard numerical techniques. This should ultimately lead to a greater physical understanding of the behavior of the solution than is possible with purely numerical techniques. The perturbation analysis of the flapping motion of a rigid rotor blade in forward flight is presented, for illustrative purposes, via computerized symbolic manipulation with a method that bypasses Floquet theory.
Computer programs for elastic analysis of shells of revolution under arbitrary loading with asymmetrical stiffness properties
The digital IPS with wire cable and flex pivot nonlinearity is simulated on the digital computer to determine the effects of varying the sampling period T on the system stability, and to determine a range of optimal values of the parameters of the digital controller. The listing of the computer program is shown as well as the Dahl model for the flex pivot nonlinearity. For the wire cable nonlinearity, two ranges of values were used and the nominal values of the digital controller parameters are included.
Variational equations for maximum range three-dimensional quasisteady glide are given. Nonlinear oscillatory maximum range trajectories obtained with a refined gradient program are approximated by a superposition of quasisteady glide and linearized perturbation equation results. A basic control law is found which is closely followed for maximum cross-range trajectories. The effect of a reradiative heating constraint involving velocity, altitude and angle of attack on a maximum cross-range trajectory for a space shuttle orbiter-type vehicle reentering the earth's atmosphere is investigated numerically.
A retargeting procedure is developed for use as a nonlinear low thrust guidance scheme. The selection of a control program composed of a sequence of inertially fixed thrust-acceleration vectors permits all trajectory computations to be made with closed form expressions, and allows the controls to be represented by constant parameters, thrust-acceleration vectors and thrusting times. By requiring each trajectory to be time optimal, the guidance problem is transformed into a parameter optimization problem which is solved by the conjugate gradient method. The scheme is applied to a low thrust capture mission, and the results of computer simulations are presented.