Dynamic analysis of a flexible spacecraft with rotating components. Volume 3: Program code
For abstract, see N76-10204.
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.
For abstract, see N76-10204.
Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.
A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.
An independent twelve degree-of-freedom simulation of the X-43A separation trajectory was created with the Program to Optimize Simulated trajectories (POST II). This simulation modeled the multi-body dynamics of the X-43A and its booster and included the effect of two pyrotechnically actuated pistons used to push the vehicles apart as well as aerodynamic interaction forces and moments between the two vehicles. The simulation was developed to validate trajectory studies conducted with a 14 degree-of-freedom simulation created early in the program using the Automatic Dynamic Analysis of Mechanics Systems (ADAMS) simulation software. The POST simulation was less detailed than the official ADAMS-based simulation used by the Project, but was simpler, more concise and ran faster, while providing similar results. The increase in speed provided by the POST simulation provided the Project with an alternate analysis tool. This tool was ideal for performing separation control logic trade studies that required the running of numerous Monte Carlo trajectories.
We have pursued several investigations using the Polar/TIDE data set. The first was a comparison of TIDE high-altitude observations with similar ion flux signatures in the Los Alamos National Laboratory (LANL) magnetospheric plasma analyzer (MPA) data. There are several geosynchronously orbiting satellites with LANL MPA instruments onboard. When the satellite is immersed in fresh plasma sheet electrons, the spacecraft voltage drops to several hundred volts negative, and cold ions are accelerated in to the detector. In the LANL MPA spectrograms, the accelerated cold ions appear as a bright, narrow line, following the voltage of the spacecraft. This "ion line" is seen regularly on the nightside, but has not received much attention. The Polar TIDE observations indicated a ubiquitous "lobal wind" in the near-- magnetotail, a field-aligned stream of approx. 100 eV ions flowing out of both polar ionospheres. The interesting result of this research is that the MPA ion line is also peaked in the field-aligned direction, even though the potential well should be uniform in all directions. It is believed that this is evidence that the lobal winds not only populate the high-latitude lobes, but fill the lobes all the way in to the near-Earth plasma sheet (the location of geosynchronous orbit). This activity developed into a full-scale survey of the lobal wind observations in the TIDE database. The universality of these observations with respect to local time and solar wind conditions implies that the ionospheric outflow is supplying the near-Earth plasma sheet at all times, regardless of magnetic activity. We have conducted a statistical study of the characteristics of the lobal wind in these two data sets FIDE and LANL MPA), finding much similarity between them. Using these characteristics as input conditions to our inner magnetosphere ion transport model, we have conducted simulations of the flow of these particles inside of geosynchronous orbit to show the impact these particles will have on the near-Earth space environment. As the study has progressed, these results have been presented several times to the TIDE team during teleconferences. It is planned to present this survey at the GEM 2004 Summer Workshop, and hopefully at additional conferences in the near future. We are presently writing a paper on the lobal wind occurrence statistics, which should be submitted in April or May to the Journal of Geophysical Research for publication. An undergraduate student at the University of Alabama in Huntsville, Will Maddox, working at NASA MSFC with Dr. Paul Craven, is presently creating a database of the values TIDE moments for the lobal wind (under Dr. Liemohn's guidance), and this will be the focus of a follow-on paper.
A hybrid approach for integrated analysis and design of nonlinear multibody systems is proposed. Based on the approach, a general-purpose design tool is developed for nonlinear dynamic systems subjected to nonlinear design constraints. For analysis purposes, second-order nonlinear equations of motion of the system are automatically generated using general-purpose multibody formulations. Once the equations are solved, they are written in first-order form to take advantage of the first-order formulations of design sensitivity analysis for dynamic systems. A nonlinear programming technique is then used to optimize the nonlinear systems, such that a nonlinear cost function is minimized and performance constraints are satisfied. The approach proved to be very general and useful for design and analysis of controlled multibody systems. The approach is applied to the design of passive dynamic controllers for slewing control of multibody systems.
This paper presents the results of a program of analysis and tests to determine the dynamic properties of elastomers as a function of strain and ambient temperature. Measurements were also made to determine the temperature distribution in the elastomer samples during the tests. These measured properties were compared with analytical predictions based on a visco-elastic model designed to take into account the self-heating of the materials as a function of strain. The test method used was well-established Base Excitation Resonant Mass Technique. The specimens tested were two cylindrical button compression specimens and a shear specimen. Strain was shown to be an important parameter in determining the dynamic properties of the elastomers. In general, these properties were much more sensitive to strain than to frequency. The self-heating effect was found to account for a portion of the strain sensitivity of these properties.
Most tank slosh test programs are not conducted with flight hardware in a flight-like environment due to cost and schedule constraints. Instead, subscale tanks are typically tested on the ground. If the “real”, flight liquid in the tank is a cryogenic propellant, sometimes a cryogen is used during testing, but more often a simulant liquid is used. Fluid dynamic similarity analysis, or “scaling analysis” as it is sometimes called, for slosh examines nondimensional numbers to determine which fluid dynamic regime(s) a tank is operating in. If the relevant nondimensional numbers of a scale tank test are similar enough to the real tank’s nondimensional numbers, then it is likely that the scale test data will be similar to what would be obtained with the real tank. This type of analysis is useful for guiding experiment design, test liquid selection, and for understanding the physical relevance of prior test data to a specific application. The authors have performed scaling analyses for many slosh test programs conducted in ground, parabolic flight, suborbital, and in-orbit environments. Basics of how to perform the analysis, along with select examples from ground and microgravity environments, are presented. Suggestions for various assumptions and compromises to make are discussed. Limitations and challenges are also discussed.
An analysis and computer program were prepared to realistically simulate the dynamic behavior of a class of satellites consisting of two end bodies separated by a connecting structure. The shape and mass distribution of the flexible end bodies are arbitrary; the connecting structure is flexible but massless and is capable of deployment and retraction. Fluid flowing in a piping system and rigid moving masses, representing a cargo elevator or crew members, have been modeled. Connecting structure characteristics, control systems, and externally applied loads are modeled in easily replaced subroutines. Subroutines currently available include a telescopic beam-type connecting structure as well as attitude, deployment, spin and wobble control. In addition, a unique mass balance control system was developed to sense and balance mass shifts due to the motion of a cargo elevator. The mass of the cargo may vary through a large range. Numerical results are discussed for various types of runs.
The delivery and demonstration of a computer program for the analysis of aeroelastic and dynamic properties is reported. Approaches to flutter and forced vibration of mistuned discs, and transient aerothermoelasticity are described.
Studies made on the use of NASTRAN for nuclear power plant analysis and design are presented. These studies indicate that NASTRAN could be effectively used for static, dynamic and special purpose problems encountered in the design of such plants. Normal mode capability of NASTRAN is extended through a post-processor program to handle seismic analysis. Static and dynamic substructuring is discussed. Extension of NASTRAN to include the needs in the civil engineering industry is discussed.
CO-ST-IN computer program developed to facilitate study of control/structure interactions; dynamic coupling between control systems and flexible structures. Works in conjunction with other programs to enhance analysis of vibrations. Designed to complement and enhance rather than replace software tools for analysis of control systems and dynamics of structures. FORTRAN 77 program.
Personal computers (PCs) are now used extensively for engineering analysis. their capability exceeds that of mainframe computers of only a few years ago. Programs originally written for mainframes have been ported to PCs to make their use easier. One of these programs is ARDS (Analysis of Rotor Dynamic Systems) which was developed at Arizona State University (ASU) by Nelson et al. to quickly and accurately analyze rotor steady state and transient response using the method of component mode synthesis. The original ARDS program was ported to the PC in 1995. Several extensions were made at ASU to increase the capability of mainframe ARDS. These extensions have also been incorporated into the PC version of ARDS. Each mainframe extension had its own user manual generally covering only that extension. Thus to exploit the full capability of ARDS required a large set of user manuals. Moreover, necessary changes and enhancements for PC ARDS were undocumented. The present document is intended to remedy those problems by combining all pertinent information needed for the use of PC ARDS into one volume.
Techniques are presented for the solution of structural dynamic systems on an electronic digital computer using FORMA (FORTRAN Matrix Analysis). FORMA is a library of subroutines coded in FORTRAN 4 for the efficient solution of structural dynamics problems. These subroutines are in the form of building blocks that can be put together to solve a large variety of structural dynamics problems. The obvious advantage of the building block approach is that programming and checkout time are limited to that required for putting the blocks together in the proper order.
Excessive vibration is one of the most prevalent technical obstacles encountered in the development of new rotorcraft. The inability to predict these vibrations is primarily due to deficiencies in analysis and simulation tools. The Langley Rotorcraft Structural Dynamics Program was instituted in 1984 to meet long term industry needs in the area of rotorcraft vibration prediction. As a part of the Langley program, this research endeavors to develop an efficient means of coupling the rotor to the airframe for preliminary design analysis of helicopter airframe vibrations. The main effort was to modify the existing computer program for modeling the dynamic and aerodynamic behavior of rotorcraft called DYSCO (DYnamic System COupler) to calculate the rotor impedance. DYSCO was recently developed for the U.S. Army and has proven to be adaptable for the inclusion of new solution methods. The solution procedure developed to use DYSCO for the calculation of rotor impedance is presented. Verification of the procedure by comparison with a known solution for a simple wind turbine model is about 75 percent completed, and initial results are encouraging. After the wind turbine impedance is confirmed, the verification effort will continue by comparison to solutions of a more sophisticated rotorcraft model. Future work includes determination of the sensitivity of the rotorcraft airframe vibrations to helicopter flight conditions and rotor modeling assumptions. When completed, this research will ascertain the feasibility and efficiency of the impedance matching method of rotor-airframe coupling for use in the analysis of airframe vibrations during the preliminary rotorcraft design process.
This paper presents the results of a program of analysis and test to determine the dynamic properties of elastomer cartridges operating under a rotating load. These measured properties were compared to predictions based on results of unidirectional tests with the same elastomer material. The test method for the dynamic stiffness and damping measurements was essentially the same as the Base Excitation Resonant Mass Method. The primary difference is that the exciting force used for these most recent tests was exerted by rotating unbalance in a rotational test rig rather than a shake table. The specimens tested were: two rectangular cross-section, continuous ring cartridges of different cross-section and three cylindrical button cartridges of different button thickness. Tests were performed for strains from about 0.0001 to about 0.01 (double amplitude). Material properties and prediction equations determined from reciprocating tests were used to make numerical predictions of stiffness, damping, and loss coefficient for the test elements, with encouraging results. Strain was shown to be an important parameter in determining these dynamic properties, particularly damping and loss coefficient.
NASA structural analysis /NASTRAN/ program applied to large radio telescopes, discussing static and dynamic analysis, problem formulation, finite element method, etc
Recent progress in the area of structural dynamics of large space structures is reviewed. Topics include system identification, large angle slewing of flexible structures, definition of scaling limitations in structural models, and recent results on a tension-stabilized antenna concept known as the hoop-column. Increasingly complex laboratory experiments guide most of the activities leading to realistic technological developments. Theoretical progress in system identification based on system realization theory resulting in unification of several methods is reviewed. Experimental results from implementation of a theoretical large-angle slewing control approach are shown. Status and results of the development of a research computer program for analysis of the transient dynamics of large angle motion of flexible structures are presented. Correlation of results from analysis and vibration tests of the hoop-column antenna concepts are summarized.