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.
Analytical models were provided for the orbiting grid structure and the joint dominated beam and computational procedures used in determining the eigen value characteristics.
This paper documents an invited presentation given to The Boeing Company, Seattle, Washington, on September 9, 1997. The audience consisted of structural dynamic and flight test engineers from the Boeing Commercial Airplane Group who were interested in discussing research which may be applied to future flight flutter test programs. A method to compute robust flutter margins is described which is a significant departure from traditional methods. This method uses the structured singular value, mu, to compute a flutter margin which directly accounts for modeling errors such that a worst-case flutter margin is computed with respect to those errors. This method may be applied in several ways. A post-flight application uses data sets from multiple test points to compute worst-case flutter margins and a worst-case flight envelope. An on-line implementation computes flutter margins at each test point to track the flutter margins during a flight test. This on-line implementation is the basis for a flutterometer flight test tool that displays the distance to flutter at a given test point. Such a tool was not previously possible using traditional flutter flight test analysis methods. The F/A-18 System Research Aircraft was used to demonstrate these applications using flight data recorded from test points throughout the flight envelope.
In 2012, the Mars Science Laboratory Mission (MSL) will deploy NASA's largest extra-terrestrial parachute, a technology integral to the safe landing of its advanced robotic explorer on the surface. The supersonic parachute system is a mortar deployed 21.5 m disk-gap-band (DGB) parachute, identical in geometric scaling to the Viking era DGB parachutes of the 1970's. The MSL parachute deployment conditions are Mach 2.3 at a dynamic pressure of 750 Pa. The Viking Balloon Launched Decelerator Test (BLDT) successfully demonstrated a maximum of 700 Pa at Mach 2.2 for a 16.1 m DGB parachute in its AV4 flight. All previous Mars deployments have derived their supersonic qualification from the Viking BLDT test series, preventing the need for full scale high altitude supersonic testing. The qualification programs for Mars Pathfinder, Mars Exploration Rover, and Phoenix Scout Missions were all limited to subsonic structural qualification, with supersonic performance and survivability bounded by the BLDT qualification. The MSL parachute, at the edge of the supersonic heritage deployment space and 33% larger than the Viking parachute, accepts a certain degree of risk without addressing the supersonic environment in which it will deploy. In addition, MSL will spend up to 10 seconds above Mach 1.5, an aerodynamic regime that is associated with a known parachute instability characterized by significant canopy projected area fluctuation and dynamic drag variation. This aerodynamic instability, referred to as "area oscillations" by the parachute community has drag performance, inflation stability, and structural implications, introducing risk to mission success if not quantified for the MSL parachute system. To minimize this risk and as an alternative to a prohibitively expensive high altitude test program, a multi-phase qualification program using computation simulation validated by subscale test was developed and implemented for MSL. The first phase consisted of 2% of fullscale supersonic wind tunnel testing of a rigid DGB parachute with entry-vehicle to validate two high fidelity computational fluid dynamics (CFD) tools. The computer codes utilized Large Eddy Simulation and Detached Eddy Simulation numerical approaches to accurately capture the turbulent wake of the entry vehicle and its coupling to the parachute bow-shock. The second phase was the development of fluid structure interaction (FSI) computational tools to predict parachute response to the supersonic flow field. The FSI development included the integration of the CFD from the first phase with a finite element structural model of the parachute membrane and cable elements. In this phase, a 4% of full-scale supersonic flexible parachute test program was conducted to provide validation data to the FSI code and an empirical dataset of the MSL parachute in a flight-like environment. The final phase is FSI simulations of the full-scale MSL parachute in a Mars type deployment. Findings from this program will be presented in terms of code development and validation, empirical findings from the supersonic testing, and drag performance during supersonic operation.
A NASTRAN analysis of the external tank (ET) substructure of the 1/8-scale space shuttle structural dynamics model is described. The NASTRAN hydroelastic procedures were used to form a model of the liquid oxygen portion of the ET. Large computer storage requirements and running times were required unless these procedures were modified. Several possibilities were demonstrated including the substitution of the real for the complex eigenvalue routine and the use of the OMIT capability to reduce the number of fluid coordinates. A NASTRAN model of the complete ET was formed and reduced to 252 degrees of freedom. A review of the eigenvector extracted, using the unsymmetrical inverse power method, indicated that the structural OMIT's resulted in unsatisfactory modal deflections. Modifications to the NASTRAN hydroelastic capability appeared necessary to generate a successful ET model. An approach is described which, by assuming incompressibility, reduces the fluid representation to a symmetric mass matrix which can be added to the structural mass. The problem is solved using faster and more efficient eigenvalue routines. Using this approach the ET NASTRAN model was analyzed for three separate weight conditions.
Distributed parameter modeling of structural dynamics has been limited to simple spacecraft configurations because of the difficulty of handling several distributed parameter systems linked at their boundaries. Although there is other computer software able to generate such models or complex, flexible spacecraft, unfortunately, neither is suitable for parameter estimation. Because of this limitation the computer software PDEMOD is being developed for the express purposes of modeling, control system analysis, parameter estimation and structure optimization. PDEMOD is capable of modeling complex, flexible spacecraft which consist of a three-dimensional network of flexible beams and rigid bodies. Each beam has bending (Bernoulli-Euler or Timoshenko) in two directions, torsion, and elongation degrees of freedom. The rigid bodies can be attached to the beam ends at any angle or body location. PDEMOD is also capable of performing parameter estimation based on matching experimental modal frequencies and static deflection test data. The underlying formulation and the results of using this approach for test data of the Mini-MAST truss will be discussed. The resulting accuracy of the parameter estimates when using such limited data can impact significantly the instrumentation requirements for on-orbit tests.
A launch vehicle can experience large unsteady aerodynamic forces in the transonic regime that, while usually only lasting for tens of seconds during launch, could be devastating if structural components and electronic hardware are not designed to account for them. These aerodynamic loads are difficult to experimentally measure and even harder to computationally estimate. The current method for estimating buffet loads is through the use of a few hundred unsteady pressure transducers and wind tunnel test. Even with a large number of point measurements, the computed integrated load is not an accurate enough representation of the total load caused by buffeting. This paper discusses an attempt at using a dynamic balance to experimentally determine buffet loads on a generic scale hammer head launch vehicle model tested at NASA Ames Research Center's 11' x 11' transonic wind tunnel. To use a dynamic balance, the structural characteristics of the model needed to be identified so that the natural modal response could be and removed from the aerodynamic forces. A finite element model was created on a simplified version of the model to evaluate the natural modes of the balance flexures, assist in model design, and to compare to experimental data. Several modal tests were conducted on the model in two different configurations to check for non-linearity, and to estimate the dynamic characteristics of the model. The experimental results were used in an inverse force determination technique with a psuedo inverse frequency response function. Due to the non linearity, the model not being axisymmetric, and inconsistent data between the two shake tests from different mounting configuration, it was difficult to create a frequency response matrix that satisfied all input and output conditions for wind tunnel configuration to accurately predict unsteady aerodynamic loads.
Grid adaptation has recently received attention in the computational fluid dynamics (CFD) community as a means to capture the salient features of a flowfield by either moving grid points of a structured or by adding cells in an unstructured manner. An approach based on a background cartesian mesh is investigated from which the geometry is 'cut' out of the mesh. Once the mesh is obtained, a solution on this coarse grid is found, that indicates which cells need to be refined. This process of refining/solving continues until the flow is grid refined in terms of a user specified global parameter (such as drag coefficient etc.). The advantages of this approach are twofold: the generation of the base grid is independent of the topology of the bodies or surfaces around/through which the flow is to be computed, and the resulting grid (in uncut regions) is highly isotropic, so that the truncation error is low. The flow solver (which, along with the grid generation is still under development) uses a completely unstructured data base, and is a finite volume, upwinding scheme. Current and future work will address generating Navier-Stokes suitable grids by using locally aligned and normal face/cell refining. The attached plot shows a simple grid about two turbine blades.
This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.
This work quantifies the potential aeroelastic benefits of tailoring a full-scale wing box structure using tailored thickness distributions, material distributions, or both simultaneously. These tailoring schemes are considered for the wing skins, the spars, and the ribs. Material grading utilizes a spatially-continuous blend of two metals: Al and Al+SiC. Thicknesses and material fraction variables are specified at the 4 corners of the wing box, and a bilinear interpolation is used to compute these parameters for the interior of the planform. Pareto fronts detailing the conflict between static aeroelastic stresses and dynamic flutter boundaries are computed with a genetic algorithm. In some cases, a true material grading is found to be superior to a single-material structure.
Because of inherent low damping and high flexibility, large space structures may require some form of active control of their dynamics. Because of the apparent inability to accurately model the dynamics of these structures, methods for parameter adaptive control are now being developed at Langley. This approach uses a digital computer to process discrete sensor data, identify modal parameters, calculate modal control gains, and then convert the modal forces to real forces. Some of the problems considered are: (1) the possibility that there may be modes to control with limited amounts of hardware, and (2) the required accuracy of identified structural parameters.
A Stirling engine was tested without auxiliaries at NASA-Lewis. Three different regenerator configurations were tested with hydrogen. The test objectives were (1) to obtain steady-state and dynamic engine data, including indicated power, for validation of an existing computer model for this engine; and (2) to evaluate structurally the use of silicon carbide regenerators. This paper presents comparisons of the measured brake performance, indicated mean effective pressure, and cyclic pressure variations with those predicted by the code. The measured data tended to be lower than the computer code predictions. The silicon carbide foam regenerators appear to be structurally suitable, but the foam matrix tested severely reduced performance.
A combined stochastic feedforward and feedback control design methodology was developed. The objective of the feedforward control law is to track the commanded trajectory, whereas the feedback control law tries to maintain the plant state near the desired trajectory in the presence of disturbances and uncertainties about the plant. The feedforward control law design is formulated as a stochastic optimization problem and is embedded into the stochastic output feedback problem where the plant contains unstable and uncontrollable modes. An algorithm to compute the optimal feedforward is developed. In this approach, the use of error integral feedback, dynamic compensation, control rate command structures are an integral part of the methodology. An incremental implementation is recommended. Results on the eigenvalues of the implemented versus designed control laws are presented. The stochastic feedforward/feedback control methodology is used to design a digital automatic landing system for the ATOPS Research Vehicle, a Boeing 737-100 aircraft. The system control modes include localizer and glideslope capture and track, and flare to touchdown. Results of a detailed nonlinear simulation of the digital control laws, actuator systems, and aircraft aerodynamics are presented.
A Stirling engine was tested without auxiliaries at Nasa-Lewis. Three different regenerator configurations were tested with hydrogen. The test objectives were: (1) to obtain steady-state and dynamic engine data, including indicated power, for validation of an existing computer model for this engine; and (2) to evaluate structurally the use of silicon carbide regenerators. This paper presents comparisons of the measured brake performance, indicated mean effective pressure, and cyclic pressure variations from those predicted by the code. The silicon carbide foam generators appear to be structurally suitable, but the foam matrix showed severely reduced performance.
Three design methods, Linear Quadratic Gaussian with Loop Transfer Recovery (LQG/LTR), H-infinity, and mu-synthesis, are used to obtain compensators for suppressing the vibrations of a 10-bay vertical truss structure, a component typical of what may be used to build a large space structure. For the design process the plant dynamic characteristics of the structure were determined experimentally using an identification method. The resulting compensators were implemented on a digital computer and tested for their ability to suppress the first bending mode response of the 10-bay vertical truss. Time histories of the measured motion are presented, and modal damping obtained during the experiments are compared with analytical predictions. The advantages and disadvantages of using the various design methods are discussed.
Many complex systems can be modeled as a collection of interacting bodies, where the relative motion of the bodies may be large. The dynamics of such systems are simulated using multibody dynamic formulations. Many of these treat each body as a rigid component, but recently the flexibility of the components has been incorporated. This paper presents a residual flexibility formulation of the multibody dynamics problem. The formulation is very simple and offers great computational efficiency since it treats each body as a free structure in space, interacting with other bodies only through interface forces. Each body's accelerations can be solved independently, as can each set of interface forces. We have applied the technique successfully to several special applications, and the initial implementation in a general mechanisms code has given excellent results in comparison to a direct finite element representation of flexibility.
The benefits of structural dynamics modeling methods in aerospace structures are reviewed. Four major issues in structural dynamics modeling are discussed which encompass most of its subdisciplines: reduced order modeling, constraints in problems with large motions, computational strategies, and fundamental methods. Directions for future research in these areas are addressed.
Computing speed is leaping forward by several orders of magnitude each decade. Engineers and scientists gathered at a NASA Langley symposium to discuss these exciting trends as they apply to parallel computational methods for large-scale structural analysis and design. Among the topics discussed were: large-scale static analysis; dynamic, transient, and thermal analysis; domain decomposition (substructuring); and nonlinear and numerical methods.
The concept of using a system, consisting of a tow aircraft, glider and tow line, which would enable subsonic flight at altitudes above 24 km (78 kft) has previously been investigated. The preliminary results from these studies seem encouraging. Under certain conditions these studies indicate the concept is feasible. However, the previous studies did not accurately take into account the forces acting on the tow line. Therefore in order to investigate the concept further a more detailed analysis was needed. The code that was selected was the SEADYN cable dynamics computer program which was developed at the Naval Facilities Engineering Service Center. The program is a finite element based structural analysis code that was developed over a period of 10 years. The results have been validated by the Navy in both laboratory and at actual sea conditions. This code was used to simulate arbitrarily-configured cable structures subjected to excitations encountered in real-world operations. The Navy's interest was mainly for modeling underwater tow lines, however the code is also usable for tow lines in air when the change in fluid properties is taken into account. For underwater applications the fluid properties are basically constant over the length of the tow line. For the tow aircraft/glider application the change in fluid properties is considerable along the length of the tow line. Therefore the code had to be modified in order to take into account the variation in atmospheric properties that would be encountered in this application. This modification consisted of adding a variable density to the fluid based on the altitude of the node being calculated. This change in the way the code handled the fluid density had no effect on the method of calculation or any other factor related to the codes validation.