Monte Carlo dynamic analysis for lunar module landing loads
Apollo lunar module structural integrity for lunar landing verified by Monte Carlo dynamic analysis
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.
Apollo lunar module structural integrity for lunar landing verified by Monte Carlo dynamic analysis
Real-time hybrid simulation (RTHS) is an enabling technology that has transformed engineering experimentation and helped researchers expand modeling capabilities. However, breakthroughs are necessary to expand the range of hybrid simulation methods and, thus, enable experiments with loading conditions representing multiple hazards. This paper discusses the development of a new thermomechanical RTHS framework and a systematic approach to determining RTHS control requirements. First, the framework is established using a representative finite element model of a layered structural system subjected to thermal loading. A complete two-layer system model serves as the reference system, and it is then partitioned into a numerical layer and an experimental layer that share interface conditions. Next, a thermal actuator is introduced to impose dynamic thermal loading on the experimental subsystem, serving as a transfer system. Finally, control and performance metrics are defined to evaluate the realization of interface boundary conditions and map this to the RTHS execution. Through an illustrative example considering the influence of temperature on a lunar habitat, we demonstrate how to establish controller requirements for RTHS and demonstrate that this approach can be used to conduct RTHS on structures with thermomechanical loading.
The rapid development of advanced metering infrastructure provides a new data source—building electrical load profiles with high temporal resolution. Electric load profile characterization can generate useful information to enhance building energy modeling and provide metrics to represent patterns and variability of load profiles. Such characterizations can be used to identify changes to building electricity demand due to operations or faulty equipment and controls. In this study, we proposed a two-path approach to analyze high temporal resolution building electrical load profiles: (1) time-domain analysis and (2) frequency-domain analysis. Furthermore, the commonly adopted time-domain analysis can extract and quantify the distribution of key parameters characterizing load shape such as peak-base load ratio and morning rise time, while a frequency-domain analysis can identify major periodic fluctuations and quantify load variability. We implemented and evaluated both paths using whole-year 15-minute interval smart meter data of 188 commercial office building in Northern California. The results from these two paths are consistent with each other and complementary to represent full dynamics of load profiles. The time- and frequency-domain analyses can be used to enhance building energy modeling by: (1) providing more realistic assumptions about building operation schedules, and (2) validating the simulated electric load profiles using the developed variability metrics against the real building load data.
Linear elastic columns dynamic stability under time dependent axial load, investigating almost- sure asymptotic stability
Loads analysis is traditionally performed using dynamically reduced models, which provides the benefit to reduce run time. If the reduced model frequency cutoff is not chosen appropriately, the model will lack the dynamic content required to fully represent the response of the non-reduced model. Standard guidance for reduced model frequency content, provided in NASA-STD-5002, is to solve fixed base modes up to a minimum of 1.5x the model frequency content of interest and to employ static modal truncation methods such as residual vectors, the mode acceleration method, and the residual flexibility method to account for the truncated flexibility of the missing modes. Fixed base modes require the selection of a set of degrees of freedom to be constrained which, if not properly selected, may affect the accuracy of the reduced model by excluding some of the dynamic characteristic of the full model. In this case, the standard NASA guidance would be insufficient, but it may not be readily apparent that a portion of the reduced model response is missing. This error was encountered during an independent verification and validation (IV&V) effort, where it was observed that the resulting dynamic response was lower than the inline analysis. In this specific case, despite following the standard NASA model reduction guidelines in the selection of the frequency cutoff, the inline model still did not fully capture the necessary dynamic content. As part of the IV&V, an alternate reduction methodology was employed using an unconstrained mode acceleration method. The original model initially performed a constrained reduction to twice the frequency content of interest before doing a free-free run, employing the mode acceleration method to account for the truncated modes. In contrast in the IV&V, the unconstrained model reduced directly to the needed frequency content of the free-free run, avoiding any interactions between constraints and dynamic content. To verify the model, the original reduction methodology was used to generate a series of Hurty-Craig-Bampton reductions, each with a higher frequency cutoff than the previous. The results were shown to converge once the frequency cutoff increased past eight times the frequency content of interest. At the same frequency cutoff, the results of the Hurty-Craig-Bampton model converged with the results of the unconstrained mode acceleration model. This comparative study provided confidence that the results of the unconstrained modal acceleration reduced model were correct and that the Hurty-Craig-Bampton needed to increase its frequency cutoff to fully capture the dynamic response.
Loads analysis is traditionally performed using dynamically reduced models, which provides the benefit to reduce run time. If the reduced model frequency cutoff is not chosen appropriately, the model will lack the dynamic content required to fully represent the response of the non-reduced model. Standard guidance for reduced model frequency content, provided in NASA-STD-5002, is to solve fixed base modes up to a minimum of 1.5x the model frequency content of interest and to employ static modal truncation methods such as residual vectors, the mode acceleration method, and the residual flexibility method to account for the truncated flexibility of the missing modes. Fixed base modes require the selection of a set of degrees of freedom to be constrained which, if not properly selected, may affect the accuracy of the reduced model by excluding some of the dynamic characteristic of the full model. In this case, the standard NASA guidance would be insufficient, but it may not be readily apparent that a portion of the reduced model response is missing. This error was encountered during an independent verification and validation (IV&V) effort, where it was observed that the resulting dynamic response was lower than the inline analysis. In this specific case, despite following the standard NASA model reduction guidelines in the selection of the frequency cutoff, the inline model still did not fully capture the necessary dynamic content. As part of the IV&V, an alternate reduction methodology was employed using an unconstrained mode acceleration method. The original model initially performed a constrained reduction to twice the frequency content of interest before doing a free-free run, employing the mode acceleration method to account for the truncated modes. In contrast in the IV&V, the unconstrained model reduced directly to the needed frequency content of the free-free run, avoiding any interactions between constraints and dynamic content. To verify the model, the original reduction methodology was used to generate a series of Hurty-Craig-Bampton reductions, each with a higher frequency cutoff than the previous. The results were shown to converge once the frequency cutoff increased past eight times the frequency content of interest. At the same frequency cutoff, the results of the Hurty-Craig-Bampton model converged with the results of the unconstrained mode acceleration model. This comparative study provided confidence that the results of the unconstrained modal acceleration reduced model were correct and that the Hurty-Craig-Bampton needed to increase its frequency cutoff to fully capture the dynamic response.
In a distributed computing environment, it is important to ensure that the processor workloads are adequately balanced, Among numerous load-balancing algorithms, a unique approach due to Das and Prasad defines a symmetric broadcast network (SBN) that provides a robust communication pattern among the processors in a topology-independent manner. In this paper, we propose and analyze three efficient SBN-based dynamic load-balancing algorithms, and implement them on an SGI Origin2000. A thorough experimental study with Poisson distributed synthetic loads demonstrates that our algorithms are effective in balancing system load. By optimizing completion time and idle time, the proposed algorithms are shown to compare favorably with several existing approaches.
A review of 1980s to early 2010s research supporting and documenting successful missions for on-orbit assembly of large space structures was conducted to bring together flexible space structure dynamics challenges and progress toward solutions. Research in this period focused on issues ranging from ground validation via unique test beds to in-space modal characterization and model adjustment. The phased assembly of the International Space Station (ISS), along with its evolving structural loads and dynamics, provides the central example, with ground and on-orbit tests of other systems providing contributing and contrasting examples of potential value for researchers facing today’s challenges.
For effective integration of building operations into the evolving demand response programs of the power grid, real-time decisions concerning the use of building appliances for grid services must excel on multiple criteria, ranging from the added value to occupants' comfort to the quality of the grid services. In this paper, we present a data-driven stochastic decision-support framework to dynamically rank load control alternatives in a commercial building, addressing the needs of multiple decision criteria (e.g. occupant comfort, grid service quality) under uncertainties in occupancy patterns. We adopt a stochastic multi-criteria decision algorithm recently applied to prioritize residential on/off loads, and extend it to i) consider complex load control decisions (e.g. dimming of lights, changing zone temperature set-points) in a commercial building; and ii) systematically integrate zonal occupancy patterns to better identify short-term (and time-varying) opportunities for grid service participation. We evaluate the performance of the proposed framework for curtailment of air-conditioning, lighting, and plug-loads in a multi-zone commercial office building for a range of design choices. With the help of a prototype system that integrates an interactive \textit{Data Analytics and Visualization} frontend we demonstrate a way for the building operators to monitor and change in real-time the available flexibility in energy consumption and to develop trust in the decision recommendations by interpreting the rationale behind the ranking.
Fluctuating pressure levels were measured on the flap and fuselage of an upper-surface-blown jet-flap airplane configuration in a wind tunnel. The model tested had turbofan engines with a bypass ratio of 3 and a thrust rating of 10 kN. Rectangular nozzles were mounted flush with the upper surface at 35 percent of the wing chord. Test parameters were flap deflection angle, jet impingement angle, angle of attack, free-stream velocity, spanwise location of the engine, and jet dynamic pressure. Load levels were high throughout the jet impingement region, with the highest levels (about 159 dB) occurring on the fuselage and near the knee of the flap. The magnitude of the forward-velocity effect appeared to depend upon the ratio of free-stream and jet velocities. Good agreement was obtained between fluctuating pressure spectra measured at jet dynamic pressures of 7 and 22 kPa when the spectra were scaled by nondimensional functions of dynamic pressure, velocity, and the empirical relationship between dynamic pressure and overall fluctuating pressure level.
On December 9th 2009, the International Space Station (ISS) 2A solar array mast experienced prolonged longeron shadowing during a Soyuz undocking. Analytical reconstruction of induced thermal and dynamic structural loads showed an exceedance of the mast buckling limit. Possible structural damage to the solar array mast could have occurred during this event. A Low fidelity video survey of the 2A mast showed no obvious damage of the mast longerons or battens. The decision was made to conduct an on-orbit dynamic test of the 2A array on December 18th, 2009. The test included thruster pluming on the array while photogrammetry data was recorded. The test was similar to other Dedicated Thruster Firings (DTFs) that were performed to measure structural frequency and damping of a solar array. Results of the DTF indicated lower frequency mast modes than model predictions, thus leading to speculation of mast damage. A detailed nonlinear analysis was performed on the 2A array model to assess possible solutions to modal differences. The setup of the parametric nonlinear trade study included the use of a detailed array model and the reduced mass and stiffness matrices of the entire ISS being applied to the array interface. The study revealed that the array attachment structure is nonlinear and thus was the source of error in the model prediction of mast modes. In addition, a detailed study was performed to determine mast mode sensitivity to mast longeron damage. This sensitivity study was performed to assess if the ISS program has sufficient instrumentation for mast damage detection.
The primary purpose of the Quiet Spike(TradeMark) flight research program was to analyze the aerodynamic, structural, and mechanical proof-of-concept of a large multi-stage telescoping nose spike installed on the National Aeronautics and Space Administration Dryden Flight Research Center (Edwards, California) F-15B airplane. This report describes the preflight stability and control analysis performed to assess the effect of the spike on the stability, controllability, and handling qualities of the airplane; and to develop an envelope expansion approach to maintain safety of flight. The overall flight test objective was to collect flight data to validate the spike structural dynamics and loads model up to Mach 1.8. Other objectives included validating the mechanical feasibility of a morphing fuselage at operational conditions and determining the near-field shock wave characterization. The two main issues relevant to the stability and control objectives were the effects of the spike-influenced aerodynamics on the F-15B airplane flight dynamics, and the air data and angle-of-attack sensors. The analysis covered the sensitivity of the stability margins, and the handling qualities due to aerodynamic variation and the maneuvering limitations of the F-15B Quiet Spike configuration. The results of the analysis and the implications for the flight test program are also presented.
Tests were conducted on the DOE/NASA Mod-O 200-kW horizontal-axis wind turbine in a one-bladed rotor configuration. The objectives of the test were to evaluate the performance, loads, and dynamic characteristics of a one-bladed rotor, and then to compare these parameters with those of an aerodynamically similar two-bladed rotor configuration. Test operations showed that this intermediate-size (15.2-m radius) one-bladed rotor configuration can be operated successfully. Test results show that the one-bladed rotor had cyclic blade loads comparable to those of a two-bladed rotor. A moderate power penalty equivalent to a reduction in windspeed of 1 m/sec occurred with the one-bladed rotor when operated at a rotor speed 50 percent higher than that of the two-bladed rotor.
The SCAPE (SuperConducting Arbitrarily Polarizing Emitter) undulator is under development at the Advanced Photon Source (APS). This new undulator requires a cryostat that will be designed based on expected heat loads. For instance, the expected heating of the beam chamber by electron beam is estimated to be at a level of 182 W – much higher than in planar SCUs. This and other challenges require careful thermal analysis of the LHe-based SCAPE cryostat. A detailed thermal model of the LHe-based SCAPE cryostat has been created in ANSYS. This paper presents calculated cooling capacity and temperatures of the SCAPE cryostat for the static and dynamic heat loads. Index Terms—SCAPE (SuperConducting Arbitrarily Polarizing Emitter), cryocooler, LHe, thermal isolation, thermal conductance.
HIV-1 plasma viral load decays in a biphasic manner during antiretroviral therapy (ART). It was hypothesized that this is due to infection of different cell types, namely CD4+ T cells and macrophages. We studied this possibility directly by modeling the decay of HIV-1 in humanized mice. We utilized previously published data from humanized T-cell only mice (TOM) and myeloid-only mice (MOM) infected with HIV-1 and treated with a potent ART regimen. Viral load decay dynamics were modeled using either a single or a biexponential decay fitted using nonlinear mixed effects techniques. Fits were compared using the corrected Bayesian information criterion (BICc). In TOM, the biphasic model was significantly better than a single-phase decay model (ΔBICc ≈ 16) despite additional parameters. In MOM, the biphasic decay was statistically better, but there was substantial uncertainty because the virus goes below detection very fast. The first-phase half-life was consistent between groups (1.2 days in MOM and 1.3 days in TOM) and similar to the half-life estimated in human infection. The second-phase decay in these mice was minimal likely due to low initial viral loads. Additional analyses with mice containing both CD4+ T cells and macrophages or X4-tropic virus-infected MOM mice confirmed the biphasic pattern, demonstrating the robustness of this result. The biphasic decline in HIV-1 occurs, even with only CD4+ T cells, refuting the hypothesis that distinct cell populations (CD4+ T cells and macrophages) drive each decay phase. These findings support an alternative model in which the observed dynamics arise from intrinsic properties of the viral infection lifecycle rather than from cellular compartmentalization.
In an isolated grid, wind turbines are typically deployed to provide energy to maximize energy production, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) this paper demonstrates, through desktop simulations, advanced wind turbine controls that can be employed to support higher contributions of wind in isolated grids, and to demonstrate ways that wind can play a role in supporting stability of an isolated grid. This isolated grid used in these desktop simulations is comprised of a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a simulated diesel generator (2 MW) and various types of loads (critical, dynamic). We developed a model of the subsystems in MATLAB/Simulink and validated them with available data from real-world components on NREL's Flatirons Campus. These validated models are then configured for various case studies. We compare the output of the desktop simulation with a baseline case with the diesel generator. Active and reactive power control of the wind turbine can help improve frequency and voltage responses in the isolated grid, respectively. By utilizing a small integrated battery energy storage system in the DC-link of the wind turbine, we also demonstrate that wind turbines can help blackstart a critical load comparable to its rated power and support other renewables (e.g. solar PV) come online and pick up an additional load. This report illustrates some of these reliability and resilience services a wind turbine can provide in an isolated grid.
Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.
Development of polyimide two-stage hydraulic actuator rod seals for application in high-performance aircraft was accomplished. The significant portion of the effort was concentrated on optimization of the chevron and K-section second-stage seal geometries to satisfy the requirements for operation at 450 K (350 F) with dynamic pressure loads varying between 200 psig steady-state and 1500 psig impulse cycling. Particular significance was placed on reducing seal gland dimension by efficiently utilizing the fatigue allowables of polyimide materials. Other objectives included investigation of pressure balancing techniques for first-stage polyimide rod seals for 4000 psig 450 K(350 F) environment and fabrication of a modular retainer for the two-stage combination. Seals were fabricated in 0.0254 m (1.0in.) and 0.0635 m (2.5in.) sizes and tested for structural integrity, frictional resistance, and endurance life. Test results showed that carefully designed second stages using polyimides could be made to satisfy the dynamic return pressure requirements of applications in high-performance aircraft. High wear under full system pressure indicated that further research is necessary to obtain an acceptable first-stage design. The modular retainer was successfully tested and showed potential for new actuator applications.