Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic mechanical loading”

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.

At least 613 records · Page 34

LS-DYNA User-Defined Internal Ballistic Modeling

One of the challenges in modeling deflagration of solid propellants with LS-DYNA is its limited capability, which is mostly applicable to airbag systems that use gaseous nitrogen generated by burning sodium azide. To overcome this limitation and enable the modeling of custom propellant grains with specific geometries, perforations, surface inhibitors, impetus, burn rates, and co-volumes, a user-defined burn model must be defined. In this study, a custom internal ballistic analysis code is integrated into LS-DYNA to simulate kinematic systems driven by pyro-mechanical devices such as pyro pushers, cutters, thrusters, separations bolts, ejection seat catapults, etc. Step-by-step guidance is provided on implementing a user-defined loading subroutine and the requirements for compiling a custom LS-DYNA executable along with a simple pyro thruster example.

pyrotechnics↗

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗

Statistical Analysis of Inter-Area Oscillations in the U.S. Eastern Interconnection: A 2017-2023 Perspective

Recent advancements and the accumulation of high-resolution, long-term phasor measurement unit (PMU) data have provided detailed insights into inter-area oscillations in power grids. This study conducts a comprehensive statistical analysis of inter-area oscillations within the United States Eastern Interconnection from 2017 to 2023. Utilizing data captured by the advanced wide-area Frequency Monitoring Network (FNET/GridEye), this investigation examines the occurrence patterns, dominant frequencies, damping ratios, and excitation mechanisms of these oscillations. Our analysis sheds light on the evolving statistical behaviors of inter-area oscillations, offering updated and critical information for grid operators and planners. The insights gained from this study can be instrumental in enhancing the operational resilience of the power network and guiding strategic developments in grid infrastructure to accommodate future challenges. Additionally, the study discusses emerging challenges associated with the modernization of the power grid, including increased renewable penetration, dynamic load variability, and cyber-physical vulnerabilities that complicate oscillation monitoring and control.

Inter-area oscillations↗

Main Shaft Bending Moment Analysis for a Gamesa G97 2-MW Wind Turbine

In a collaborative project between Siemens Gamesa Renewable Energy and the National Renewable Energy Laboratory (NREL), a cost-effective method to measure dynamic main shaft torque and the influence of non-torque loads was applied to the gearbox of the Gamesa G97 2-megawatt (MW) wind turbine at the NREL Flatirons Campus and verified through a field test campaign that took place in 2023. The method involved the installation of fiber-optic sensors on the ring gear of the gearbox. Additionally, traditional strain gauges and a rotary encoder were installed on the main shaft to measure mechanical torque and orthogonal bending moments. We compare the measured total main shaft bending moments from the field test campaign to those predicted by an updated OpenFAST model of the G97 wind turbine for similar environmental conditions. We also assessed the sensitivity of the predicted bending moment to varying levels of turbulence intensity, shear exponent, and yaw error.

17 WIND ENERGY↗

Design of helicopter rotor blades for optimum dynamic characteristics

The possibilities and the limitations of tailoring blade mass and stiffness distributions to give an optimum blade design in terms of weight, inertia, and dynamic characteristics are investigated. Changes in mass or stiffness distribution used to place rotor frequencies at desired locations are determined. Theoretical limits to the amount of frequency shift are established. Realistic constraints on blade properties based on weight, mass moment of inertia size, strength, and stability are formulated. The extent hub loads can be minimized by proper choice of EL distribution is determined. Configurations that are simple enough to yield clear, fundamental insights into the structural mechanisms but which are sufficiently complex to result in a realistic result for an optimum rotor blade are emphasized.

Peters, D. A.↗

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗

Dynamic Reconfiguration of Switchgrass Proteomes in Response to Rust ( Puccinia novopanici ) Infection

Switchgrass (Panicum virgatum L.) can be infected by the rust pathogen (Puccinia novopanici) and results in lowering biomass yields and quality. Label-free quantitative proteomics was conducted on leaf extracts harvested from non-infected and infected plants from a susceptible cultivar (Summer) at 7, 11, and 18 days after inoculation (DAI) to follow the progression of disease and evaluate any plant compensatory mechanisms to infection. Some pustules were evident at 7 DAI, and their numbers increased with time. However, fungal DNA loads did not appreciably change over the course of this experiment in the infected plants. In total, 3830 proteins were identified at 1% false discovery rate, with 3632 mapped to the switchgrass proteome and 198 proteins mapped to different Puccinia proteomes. Across all comparisons, 1825 differentially accumulated switchgrass proteins were identified and subjected to a STRING analysis using Arabidopsis (A. thaliana L.) orthologs to deduce switchgrass cellular pathways impacted by rust infection. Proteins associated with plastid functions and primary metabolism were diminished in infected Summer plants at all harvest dates, whereas proteins associated with immunity, chaperone functions, and phenylpropanoid biosynthesis were significantly enriched. At 18 DAI, 1105 and 151 proteins were significantly enriched or diminished, respectively. Many of the enriched proteins were associated with mitigation of cellular stress and defense.

59 BASIC BIOLOGICAL SCIENCES↗

Roller bearing geometry design

A theory of kinematic stabilization of rolling cylinders is extended and applied to the design of cylindrical roller bearings. The kinematic stabilization mechanism puts a reverse skew into the rolling elements by changing the roller taper. Twelve basic bearing modification designs are identified amd modeled. Four have single transverse convex curvature in their rollers while eight have rollers which have compound transverse curvature made up of a central cylindrical band surrounded by symmetric bands with slope and transverse curvature. The bearing designs are modeled for restoring torque per unit axial displacement, contact stress capacity, and contact area including dynamic loading, misalignment sensitivity and roller proportion. Design programs are available which size the single transverse curvature roller designs for a series of roller slopes and load separations and which design the compound roller bearings for a series of slopes and transverse radii of curvature. The compound rollers are proportioned to have equal contact stresses and minimum size. Design examples are also given.

Savage, M.↗

Probabilistic structural mechanics research for parallel processing computers

Aerospace structures and spacecraft are a complex assemblage of structural components that are subjected to a variety of complex, cyclic, and transient loading conditions. Significant modeling uncertainties are present in these structures, in addition to the inherent randomness of material properties and loads. To properly account for these uncertainties in evaluating and assessing the reliability of these components and structures, probabilistic structural mechanics (PSM) procedures must be used. Much research has focused on basic theory development and the development of approximate analytic solution methods in random vibrations and structural reliability. Practical application of PSM methods was hampered by their computationally intense nature. Solution of PSM problems requires repeated analyses of structures that are often large, and exhibit nonlinear and/or dynamic response behavior. These methods are all inherently parallel and ideally suited to implementation on parallel processing computers. New hardware architectures and innovative control software and solution methodologies are needed to make solution of large scale PSM problems practical.

Sues, Robert H.↗

Investigation of Seal-to-Floor Effects on Semi-Span Transonic Models

In an effort to achieve the maximum possible Reynolds number (Re) when conducting production testing for flight loads aerodynamic databases, it has been the preferred practice of The Boeing Company / Commercial Airplanes (BCA) -- Loads and Dynamics Group since the early 1990's to test large scale semi-span models in the 11- By 11-Foot Transonic Wind Tunnel (TWT) leg of the Unitary Plan Wind Tunnel (UPWT) at the NASA Ames Research Center (ARC). There are many problems related to testing large scale semi-span models of high aspect ratio flexible transport wings, such as; floor boundary layer effects, wing spanwise wall effects, solid blockage buoyancy effects, floor mechanical interference effects, airflow under the model effects, or tunnel flow gradient effects. For most of these issues, BCA has developed and implemented either standard testing methods or numerical correction schemes and these will not be discussed in this document. Other researchers have reported on semi-span transonic testing correction issues, however most of the reported research has been for low Mach testing. Some of the reports for low Mach testing address the difficult problem of preventing undesirable airflow under a semi-span model while ensuring unrestricted main balance functionality, however, for transonic models this issue has gone unresolved. BCA has been cognizant for sometime that there are marked differences in wing pressure distributions from semi-span transonic model testing than from full model or flight testing. It has been suspected that these differences are at least in part due to airflow under the model. Previous efforts by BCA to address this issue have proven to be ineffective or inconclusive and in one situation resulted in broken hardware. This paper reports on a Boeing-NASA collaborative investigation based on a series of small tests conducted between June 2006 and November 2007 in the 11 by 11 foot Transonic Wind Tunnel at NASA Ames on three large commercial jet transport configurations to assess the effects of sealing a semi-span model to the floor and to investigate efficient sealing and testing techniques. This document will show how sealing the model to the floor has a small but remarkably far reaching spanwise effect on wing pressures, wing local section forces and wing force summations.

Sleppy, Mark A.↗

Visual display and alarm system for wind tunnel static and dynamic loads

A wind tunnel balance monitor and alarm system developed at NASA Ames Research Center will produce several beneficial results. The costs of wind tunnel delays because of inadvertent balance damage and the costs of balance repair or replacement can be greatly reduced or eliminated with better real-time information on the balance static and dynamic loading. The wind tunnel itself will have enhanced utility with the elimination of overly cautious limits on test conditions. The microprocessor-based system features automatic scaling and 16 multicolored LED bargraphs to indicate both static and dynamic components of the signals from eight individual channels. Five individually programmable alarm levels are available with relay closures for internal or external visual and audible warning devices and other functions such as automatic activation of external recording devices, model positioning mechanisms, or tunnel shutdown.

Hanly, Richard D.↗

Review of fatigue and fracture research at NASA Langley Research Center

Most dynamic components in helicopters are designed with a safe-life constant-amplitude testing approach that has not changed in many years. In contrast, the fatigue methodology in other industries has advanced significantly in the last two decades. Recent research at the NASA Langley Research Center and the U.S. Army Aerostructures Directorate at Langley are reviewed relative to fatigue and fracture design methodology for metallic components. Most of the Langley research was directed towards the damage tolerance design approach, but some work was done that is applicable to the safe-life approach. In the areas of testing, damage tolerance concepts are concentrating on the small-crack effect in crack growth and measurement of crack opening stresses. Tests were conducted to determine the effects of a machining scratch on the fatigue life of a high strength steel. In the area of analysis, work was concentrated on developing a crack closure model that will predict fatigue life under spectrum loading for several different metal alloys including a high strength steel that is often used in the dynamic components of helicopters. Work is also continuing in developing a three-dimensional, finite-element stress analysis for cracked and uncracked isotropic and anisotropic structures. A numerical technique for solving simultaneous equations called the multigrid method is being pursued to enhance the solution schemes in both the finite-element analysis and the boundary element analysis. Finally, a fracture mechanics project involving an elastic-plastic finite element analysis of J-resistance curve is also being pursued.

Everett, Richard A., Jr.↗

Visual display and alarm system for wind tunnel static and dynamic loads

A wind tunnel balance monitor and alarm system developed at NASA Ames Research Center will produce several beneficial results. The costs of wind tunnel delays because of inadvertent balance damage and the costs of balance repair or replacement can be greatly reduced or eliminated with better real-time information on the balance static and dynamic loading. The wind tunnel itself will have enhanced utility with the elimination of overly cautious limits on test conditions. The microprocessor-based system features automatic scaling and 16 multicolored LED bargraphs to indicate both static and dynamic components of the signals from eight individual channels. Five individually programmable alarm levels are available with relay closures for internal or external visual and audible warning devices and other functions such as automatic activation of external recording devices, model positioning mechanism, or tunnel shutdown.

Hanly, Richard D.↗

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗

Shock-induced transformation of nitinol shape memory alloy: Effect of stress state on transformation

Due to its numerous practical applications and intriguing phase transformation behavior, shape memory alloys (SMAs) have garnered significant research and development interests. In the past, most studies on the mechanical behavior of SMAs have been conducted under uniaxial stress loadings. Limited research on SMAs under shock loading has not provided conclusive results regarding their transformation behavior and transformation stress under such loading. Additionally, there is a lack of comprehensive understanding regarding the effects of different stress states on transformation behavior. The main objectives of this study are to address these issues. To achieve these objectives, a series of shock wave experiments were designed and conducted. Additionally, quasi-static and dynamic uniaxial stress experiments were carried out to establish a baseline for comparison. The results revealed that the transformation stress under dynamic uniaxial strain shock loading was approximately 1.92 GPa in contrast to 0.5 GPa (quasi-static) to 0.8 GPa (dynamic) observed in uniaxial stress loading. The transformation behavior exhibited noticeable rate sensitivity for both types of loading. There appeared to be a critical strain rate above which the austenite phase was driven to a metastable state. This estimated critical axial strain rate along the loading direction was approximately 2 × 10 3 /s–4 × 10 3 /s for uniaxial stress loading and approximately 2 × 10 6 /s for uniaxial strain loading. The apparent high transformation stress for uniaxial strain loading can likely be attributed to a combination of high-pressure confinement and high strain rate. Furthermore, determining their relative contributions remains an open issue.

36 MATERIALS SCIENCE↗