Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable load”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Probabilistic simulation of long term behavior in polymer matrix composites

A methodology to compute cumulative probability distribution functions (CDF) of fatigue life for different ratios, r of applied stress to the laminate strength based on first ply failure criteria has been developed and demonstrated. Degradation effects due to long term environmental exposure and mechanical cyclic loads are considered in the simulation process. A unified time-stress dependent multi-factor interaction equation model developed at NASA Lewis Research Center has been used to account for the degradation/aging of material properties due to cyclic loads. Fast probability integration method is used to perform probabilistic simulation of uncertainties. Sensitivity of fatigue life reliability to uncertainties in the primitive random variables are computed and their significance in the reliability based design for maximum life is discussed. The results show that the graphite/epoxy (0/+45/90) deg laminate with ply thickness 0.125 in. has 500,000 cycles life for applied stress to laminate strength ratio of 0.6 and a reliability of 0.999. Also, the fatigue life reliability has been found to be most sensitive to the ply thickness and matrix tensile strength. Tighter quality controls must therefore be enforced on ply thickness and matrix strength in order to achieve high reliability of the structure.

Shah, A. R.↗
NASA NTRS2001

Automated Fiber Placement of PEEK/IM7 Composites with Film Interleaf Layers

The incorporation of thin discrete layers of resin between plies (interleafing) has been shown to improve fatigue and impact properties of structural composite materials. Furthermore, interleafing could be used to increase the barrier properties of composites used as structural materials for cryogenic propellant storage. In this work, robotic heated-head tape placement of PEEK/IM7 composites containing a PEEK polymer film interleaf was investigated. These experiments were carried out at the NASA Langley Research Center automated fiber placement facility. Using the robotic equipment, an optimal fabrication process was developed for the composite without the interleaf. Preliminary interleaf processing trials indicated that a two-stage process was necessary; the film had to be tacked to the partially-placed laminate then fully melted in a separate operation. Screening experiments determined the relative influence of the various robotic process variables on the peel strength of the film-composite interface. Optimization studies were performed in which peel specimens were fabricated at various compaction loads and roller temperatures at each of three film melt processing rates. The resulting data were fitted with quadratic response surfaces. Additional specimens were fabricated at placement parameters predicted by the response surface models to yield high peel strength in an attempt to gage the accuracy of the predicted response and assess the repeatability of the process. The overall results indicate that quality PEEK/lM7 laminates having film interleaves can be successfully and repeatability fabricated by heated head automated fiber placement.

Hulcher, A. Bruce↗
NASA NTRS1997

Growth hormone/IGF-I and/or resistive exercise maintains myonuclear number in hindlimb unweighted muscles

In the present study of rats, we examined the role, during 2 wk of hindlimb suspension, of growth hormone/insulin-like growth factor I (GH/IGF-I) administration and/or brief bouts of resistance exercise in ameliorating the loss of myonuclei in fibers of the soleus muscle that express type I myosin heavy chain. Hindlimb suspension resulted in a significant decrease in mean soleus wet weight that was attenuated either by exercise alone or by exercise plus GH/IGF-I treatment but was not attenuated by hormonal treatment alone. Both mean myonuclear number and mean fiber cross-sectional area (CSA) of fibers expressing type I myosin heavy chain decreased after 2 wk of suspension compared with control (134 vs. 162 myonuclei/mm and 917 vs. 2,076 micron2, respectively). Neither GH/IGF-I treatment nor exercise alone affected myonuclear number or fiber CSA, but the combination of exercise and growth-factor treatment attenuated the decrease in both variables. A significant correlation was found between mean myonuclear number and mean CSA across all groups. Thus GH/IGF-I administration and brief bouts of muscle loading had an interactive effect in attenuating the loss of myonuclei induced by chronic unloading.

NASA Center ARC↗
NASA NTRS2015

Optimization of an Aeroservoelastic Wing with Distributed Multiple Control Surfaces

This paper considers the aeroelastic optimization of a subsonic transport wingbox under a variety of static and dynamic aeroelastic constraints. Three types of design variables are utilized: structural variables (skin thickness, stiffener details), the quasi-steady deflection scheduling of a series of control surfaces distributed along the trailing edge for maneuver load alleviation and trim attainment, and the design details of an LQR controller, which commands oscillatory hinge moments into those same control surfaces. Optimization problems are solved where a closed loop flutter constraint is forced to satisfy the required flight margin, and mass reduction benefits are realized by relaxing the open loop flutter requirements.

Stanford, Bret K.↗
NASA NTRS2016

Natural and Anthropogenic Aerosol Trends from Satellite and Surface Observations and Model Simulations over the North Atlantic Ocean from 2002 to 2012

Aerosols contribute to Earths radiative budget both directly and indirectly, and large uncertainties remain in quantifying aerosol effects on climate. Variability in aerosol distribution and properties, as might result from changing emissions and transport processes, must be characterized. In this study, variations in aerosol loading across the eastern seaboard of theUnited States and theNorthAtlanticOcean during 2002 to 2012 are analyzed to examine the impacts of anthropogenic emission control measures using monthly mean data from MODIS, AERONET, and IMPROVE observations and Goddard Chemistry Aerosol Radiation and Transport (GOCART) model simulation.MODIS observes a statistically significant negative trend in aerosol optical depth (AOD) over the midlatitudes (-0.030 decade(sup-1)). Correlation analyses with surface AOD from AERONET sites in the upwind region combined with trend analysis from GOCART component AOD confirm that the observed decrease in the midlatitudes is chiefly associated with anthropogenic aerosols that exhibit significant negative trends from the eastern U.S. coast extending over the western North Atlantic. Additional analysis of IMPROVE surface PM(sub 2.5) observations demonstrates statistically significant negative trends in the anthropogenic components with decreasing mass concentrations over the eastern United States. Finally, a seasonal analysis of observational datasets is performed. The negative trend seen by MODIS is strongest during spring (MAM) and summer (JJA) months. This is supported by AERONET seasonal trends and is identified from IMPROVE seasonal trends as resulting from ammonium sulfate decreases during these seasons.

Jongeward, Andrew R.↗
NASA NTRS2024

SMAnalytics Automated Shape Memory Test Analysis Software User's Manual: UCFTC, UPER, and DSC Test Methods, Version 1.0.2

Over the course of 60 years of shape memory alloy research and development, the properties of these alloys have been measured using various testing methods, which are often customized by the organization performing the test. However, commercial adoption of shape memory alloys in aeronautic actuator applications requires reducing property uncertainty through standardizing test methods. Historically, differential scanning calorimetry (DSC) has been used to measure transformation temperatures under stress-free conditions. This testing method is well established, and tests are run under the applicable standard, ASTM F2004-17: Standard Test Method for Transformation Temperature of Nickel-Titanium Alloys by Thermal Analysis. Two other test methods for measuring strains and transformation temperatures under constant load or free recovery after prestraining have only recently been standardized—ASTM E3097: Standard Test Method for Mechanical Uniaxial Constant Force Thermal Cycling of SMAs (UCFTC) and ASTM E3098: Standard Test Method for Mechanical Uniaxial Pre-Strain and Thermal Free Recovery of SMAs (UPFR). These standards represent a critical step forward in producing reliable material property data for the use of these alloys in aeronautics and other commercial areas. However, no standard programs or software packages were previously available for postprocessing UCFTC- and UPFR-type test data to extract property data. Each organization that performed the testing analyzed data using its own custom manual methods, such as using a plot and a straight edge, spreadsheet plotting methods, or custom code routines, to extract the required values for strains, stresses, and transformation temperatures. SMAnalytics represents the first publicly available uniform software package for processing data generated using the DSC, UCFTC, and UPFR test methods, as well as many modifications of these methods, including loading in martensite versus austenite and performing multiple thermal cycles at stress. In doing so, it provides a tool for reducing potential error or variability in the measured values due to the precision of the technique used and person-to-person variability. It is also expected that this automated data-parsing tool can accelerate the data analysis phase of the test campaign, especially for large data files. This User’s Manual describes how the SMAnalytics software works and the method for its use.

Shape Memory Alloy↗
DOE OSTI2024

ASME Design Code Rule Changes for Nuclear Graphite

The American Society of Mechanical Engineers Boiler Pressure and Vessel Code (ASME BPVC) Section III, Division 5, Article HHA-3000 outlines graphite core component and graphite core assembly design guidelines. Graphite core components are defined as ?components manufactured from graphite that are installed to form a graphite core assembly within the reactor pressure vessel of a high temperature, graphite moderated fission reactor.? (p. 413) Graphites? inherent defect distributions do not allow for deterministic material reliability. Rather, graphite has variable strength distributions which change by grade. Article HHA-3000 outlines two semi-probabilistic methods, the full and simplified assessments, which set design load limit targets for each of three component structural reliability classes. The Design Task Group was officially recognized as a specialized task group within ASME November of 2023, though we?ve been collaborating since 2022. The purpose of the Design Task Group is to correct, clarify, and make HHA-3000 function as intended. The Design Task Group will sunset once we?ve achieved our objectives. The Design Task Group was specifically told to not write new Code. While there may be more precise and more accurate methods to determine reliability targets, the current methods are conservative, relatively simple to implement, and have thus far been considered satisfactory for setting design reliability targets. Much of the ground-work to write proposal files and background documents for records to make the changes needed to achieve our objective have been completed. The Design Task Group has documented much of their work through papers, presentations, and memorandums. Three memorandums in which INL team members had substantial contributions are found in the Appendices: FEA Modeling for the Baseline Program, Evaluating the Effects on Margin of Updating the Threshold and Shape Parameters in the Full Assessment, and Interpretations of the Full and Simplified Assessments in ASME BPVC. Most of the on-going work to achieve the Design Task Group?s objective will be addressing comments on existing records and moving records through the balloting process. The Design Task Group met bi-weekly mostly through the end of FY2023. Since February 2024, the Design Task Group has mostly been completed with solving and documenting the technical issues associated with the assessments. Unless new tasks are identified, the remaining work of the Design Task Group will be political and editorial.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Static aeroelastic analysis and tailoring of missile control fins

A concept for enhancing the design of control fins for supersonic tactical missiles is described. The concept makes use of aeroelastic tailoring to create fin designs (for given planforms) that limit the variations in hinge moments that can occur during maneuvers involving high load factors and high angles of attack. It combines supersonic nonlinear aerodynamic load calculations with finite-element structural modeling, static and dynamic structural analysis, and optimization. The problem definition is illustrated. The fin is at least partly made up of a composite material. The layup is fixed, and the orientations of the material principal axes are allowed to vary; these are the design variables. The objective is the magnitude of the difference between the chordwise location of the center of pressure and its desired location, calculated for a given flight condition. Three types of constraints can be imposed: upper bounds on static displacements for a given set of load conditions, lower bounds on specified natural frequencies, and upper bounds on the critical flutter damping parameter at a given set of flight speeds and altitudes. The idea is to seek designs that reduce variations in hinge moments that would otherwise occur. The block diagram describes the operation of the computer program that accomplishes these tasks. There is an option for a single analysis in addition to the optimization.

Mcintosh, S. C., Jr.↗
NASA NTRS2005

Energy Finite Element Analysis for Computing the High Frequency Vibration of the Aluminum Testbed Cylinder and Correlating the Results to Test Data

The Energy Finite Element Analysis (EFEA) is a finite element based computational method for high frequency vibration and acoustic analysis. The EFEA solves with finite elements governing differential equations for energy variables. These equations are developed from wave equations. Recently, an EFEA method for computing high frequency vibration of structures either in vacuum or in contact with a dense fluid has been presented. The presence of fluid loading has been considered through added mass and radiation damping. The EFEA developments were validated by comparing EFEA results to solutions obtained by very dense conventional finite element models and solutions from classical techniques such as statistical energy analysis (SEA) and the modal decomposition method for bodies of revolution. EFEA results have also been compared favorably with test data for the vibration and the radiated noise generated by a large scale submersible vehicle. The primary variable in EFEA is defined as the time averaged over a period and space averaged over a wavelength energy density. A joint matrix computed from the power transmission coefficients is utilized for coupling the energy density variables across any discontinuities, such as change of plate thickness, plate/stiffener junctions etc. When considering the high frequency vibration of a periodically stiffened plate or cylinder, the flexural wavelength is smaller than the interval length between two periodic stiffeners, therefore the stiffener stiffness can not be smeared by computing an equivalent rigidity for the plate or cylinder. The periodic stiffeners must be regarded as coupling components between periodic units. In this paper, Periodic Structure (PS) theory is utilized for computing the coupling joint matrix and for accounting for the periodicity characteristics.

Vlahopoulos, Nickolas↗
DOE OSTI2022

Cyber-Physical Dynamic System (CPDS) Modeling for Frequency Regulation and AGC Services of Distributed Energy Resources

The substantial integration of renewable energy brings significant challenges to balance the system in real time because of the variability and intermittency of renewable power. For the reliable system operation, the frequency regulation service is used to stabilize the system frequency through automatically balancing the system generation and load. On one hand, with the substantially increasing deployment of renewable energy in electricity system, the requirement of frequency regulation (FR) services increases significantly. On the other hand, the current main resource of FR services, the controllable conventional generation, is continuously decreasing in the system generation mix. This means that in the future high renewable penetration power system, additional and alternative reliable FR services providers such as distributed energy storage (DES) resources should be explored. Although the capability of utility-scale energy storage to provide FR services has been demonstrated, the integrated control and dynamic modeling of distributed energy resource (DER) providing frequency regulation grid services has been rarely explored. There are several challenges to adopt DERs to provide reliable grid services as illustrated in FERC 755. First, the distributed installation of DES requires a comprehensive cyber-physical dynamic system (CPDS) modeling to fully consider the impacts of the communication latency variability on its real time FR provision capability. Unlike the conventional generators, there are two-layer communication between DERs and system operators which increase the communication delay. Second, unlike the conventional generators whose dynamic models have been comprehensively studied, the difference among individual DER components' power-dynamic characteristic brings challenges in the accurate dynamic AGC modeling of its power-frequency relationship. Third, the temporal dependent state of charge uncertainty of DERs challenges DERs' power and frequency regulation capacities scheduling in the look-ahead generation scheduling. Therefore, the aggregator should optimize the frequency regulation provision from individual DER in real time. To overcome these challenges, this project proposes a cyber-physical dynamic system (CPDS) model to handle the uncertainty of DERs two-layer communication latency and power dynamics. The variability of DERs' communication delay and dynamic constraints will be comprehensively modeled. The DERs' AGC model with communication delay is designed to validate DESs' frequency regulation services. Like current performance-based frequency regulation evaluation, the delivery of DERs' frequency regulation will be assessed through post-analysis of the actual AGC response with respect to the AGC control signal from system operators. Consequently, the reliability improvement with DERs providing reliable frequency regulation services can be evaluated from a comprehensive perspective considering all the dynamics of communication and power.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1990

Telefunctioning - An approach to telerobotic manipulations

This study introduces three concepts: telefunctioning, a control method for achieving telefunctioning, and an approach to the analysis of human-robot interaction when telefunctioning governs the system behavior. Telefunctioning facilitates maneuvering of loads by creating a perpetual sense of the load dynamics for the operator. Telefunctioning is defined as a robotic manipulation method in which the dynamic behaviors of the slave robot and the master robot are functions of each other; these functions are the designer's choice and depend on the application. In a subclass of telefunctioning called telepresence, all of the relationships between the master and the slave are specified as unity so that all of the master and slave variables (e.g., position and velocity) are dynamically equal. To create telefunctioning, a minimum number of functions relating the robots' variables is identified. A control architecture which guarantees that the defined functions govern the dynamic behavior of the system is then developed. The stability of the closed-loop system (master robot, slave robot, human, and the load being manipulated) is analyzed, and sufficient conditions for stability are derived.

Kazerooni, H.↗
DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Detailed pressure distribution measurements obtained on several configurations of an aspect-ratio-7 variable twist wing

Detailed pressure distribution measurements were made for 11 twist configurations of a unique, multisegmented wing model having an aspect ratio of 7 and a taper ratio of 1. These configurations encompassed span loads ranging from that of an untwisted wing to simple flapped wings both with and without upper-surface spoilers attached. For each of the wing twist configurations, electronic scanning pressure transducers were used to obtain 580 surface pressure measurements over the wing in about 0.1 sec. Integrated pressure distribution measurements compared favorably with force-balance measurements of lift on the model when the model centerbody lift was included. Complete plots and tabulations of the pressure distribution data for each wing twist configuration are provided.

Holbrook, G. T.↗
DOE OSTI2022

Methods for Computing Physically Realistic Estimates of Electric Water Heater Demand Response Resource Suitable for Bulk Power System Planning Models

Demand response is commonly called on to reduce load during system peak times or to respond to contingency events. In future power systems with higher shares of wind and solar generation (which we describe together as variable generation [VG]), demand response could have more opportunities to provide energy shifting or operating reserve services. This report evaluates the ability of residential electric water heaters, both electric resistance water heaters (ERWHs) and heat pump water heaters (HPWHs), to provide such services starting from detailed whole-building energy models that realistically represent New England single family home stock. We use a parsimonious surrogate model to represent operational flexibility in a form suitable for linear and mixed integer programming. This enables relatively fast determination of aggregate contingency reserve resource, price-taking energy shifting outcomes, and in some cases the determination of aggregate models at the megawatt (MW) scale that can be directly included in large-scale grid models. After selecting modeling methods and parameters through various computational experiments, we find interquartile ranges of contingency reserve resource in ISO-NE for about 603,400 ERWHs of 45 MW - 69 MW for Claim10 (50 minute responses provided with 10 minutes of advanced notification) and 65 MW - 102 MW for Claim30 (30 minute responses provided with 30 minutes of advanced notification), and for about 619,000 HPWHs of 48 MW - 88 MW for Claim10 and 52 MW - 90 MW for Claim30. The overall reserve resource is up to 32% of total load for ERWHs providing Claim10 service, 47% for ERWHs providing Claim30 service, 93% for HPWHs providing Claim10 service, and 97% for HPWHs providing Claim30 service. More work is required to determine if HPWHs are inherently more suitable than ERWHs for providing contingency reserve or if these results reflect idiosyncrasies of the single family home stock model used in this study. The value of this contingency resource in a Near-term VG model of ISO-NE is $\$ 0.40$ to $\$1.20$ per water heater-year, and significantly larger, $\$ 3.80$ to $\$ 5.30$ per water heater-year in a Mid-term VG model of ISONE. Aggregating surrogate models to the MW-scale for energy shifting service is more challenging than for contingency service and we only present such results for ERWHs, because we were unable to determine satisfactory ways to deal with HPWHs' time-varying and path dependent operational characteristics. Individual surrogate models suitable for evaluating the energy shifting resource from both ERWHs and HPWHs are created, however, and dispatched against day-ahead prices from the Near-Term VG and Mid-Term VG models of ISO-NE. The individual surrogate models are able to access and potentially shift all 640 GWh of HPWH load and 1,547 GWh of ERWH load we modeled in two different single family home stock models. In contrast, the most effective model of aggregate ERWH shifting resource we created only captured 34.7% of the total ERWH load. Energy shifting affected by price-taking dispatch against modeled day-ahead energy prices produces per water heater year profits of $\$19.44$ - $\$22.93$ for individual HPWHs, $\$39.11$ - $\$40.54$ for individual ERWHs, and up to $\$4.00$ - $\$4.24$ for aggregated ERWHs, with the variations mainly due to grid conditions (more or less VG). When the supply-side response to these changes is accounted for, the per water heater year production cost savings for ISO-NE are $\$7.50$ to $\$17.70$ for the most effective set of endogenously dispatched aggregate ERWHs, $\$15.60$ to $\$15.70$ for individual ERWHs dispatched against the DA prices, and $\$10.70$ to $\$11.20$ for individual HPWHs dispatched against DA prices. Those ranges primarily represent the difference between Near-Term VG and Mid-Term VG grid conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Low Temperature Combustion Exploration with Negative Valve Overlap

Progressively stringent emission regulations and increasing regulatory demands on fuel economy have led to advanced combustion development. Low temperature combustion (LTC), specifically homogenous charge compression ignition (HCCI), is a promising technology for reducing exhaust emissions and improving efficiency. However, its operating range is limited to low load without boosting and EGR, due to low volumetric efficiency and high pressure rise rates. In addition, effectively controlling the combustion phasing is another challenge in realizing the associated combustion gains. In this work, advanced valve control mechanisms known as continuously variable valve duration (CVVD) and continuously variable valve timing (CVVT) were used for both intake and exhaust valvetrains to enable negative valve overlap (NVO) for trapping hot exhaust residuals and to promote multipoint simultaneous ignition. Heat release phasing was controlled by varying the fueling scheme and by adjusting the amount of NVO. Parametric studies on valve timing and duration, fueling strategy, lambda, spark assist, etc., were carried out first. Afterwards the LTC strategy was proposed and then LTC operation was explored at different engine speeds. Various approaches for extending load limits were summarized and discussed. Finally, combustion performance was compared to that of spark ignition combustion, demonstrating the combustion gains of LTC.

02 PETROLEUM↗
DOE OSTI2022

Building Energy Systems as Behind-the-Meter Resources for Grid Services: Intelligent load control and transactive control and coordination

To mitigate the impacts of climate change, significant reductions in emissions from all sectors of the economy are needed. The electricity generation sector has embarked on an ambitious plan to include renewable generation as part of its decarbonization efforts, and many cities and states are mandating all-electric buildings. While renewable resources will reduce emissions, they are not dispatchable, they vary temporally, and their generation is uncertain. Under these conditions, traditional approaches to managing grid reliability, where supply follows demand, will not be efficient and may not be cost-effective. Further, there is a more efficient alternative for balancing the supply–demand imbalance and for absorbing variability and uncertainty of renewable energy using distributed energy resources (DERs) as opposed to reserve generation. Because buildings consume more than 75% of total U.S. annual electricity consumption, behind-the-meter (BTM) DERs have a load flexibility of 77 GW of power and 90 GWh of virtual energy storage capacity nationwide (Kalsi, 2017). Therefore, some portion of the supply–demand imbalance can be met by these DERs at a lower cost compared to business-as-usual solutions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2012

Reducing Extra-Terrestrial Excavation Forces with Percussion

High launch costs and mission requirements drive the need for low mass excavators with mobility platforms, which in turn have little traction and excavation reaction capacity in low gravity environments. This presents the need for precursor and long term future missions with low mass robotic mining technology to perform In-Situ Resource Utilization (ISRU) tasks. This paper discusses a series of experiments that investigate the effectiveness of a percussive digging device to reduce excavation loads and thereby the mass of the excavator itself. The goal of percussive excavation is to fluidize dry regolith in front of the leading edge of the tool by mechanically separating the microscopic interlocking grains resulting in a reduced force needed to shear the soil. There are several variables involved with this technique; this experiment varied: Impact energy, frequency, and excavation speed and held constant: impact direction, depth of cut, angle of tool, and soil bulk density. The test apparatus consisted of an aluminum truss bridge with a central pivoting arm. Attached to the arm was a winch with a load cell in line that recorded the tension in the cable and therefore the excavation load. The arm could be adjusted for excavation depth which was recorded along with the arm angle relative to the bridge. A percussive mechanism and 30" wide pivoting bucket were attached at the end of the arm simulating a basic backhoe with a percussion direction tangent to the direction of . movement. Internally the mechanism used a set of die springs and barrel cam to produce the percussive blow. By changing the springs and the speed of the motor the impact energy and frequency of percussion could be varied independently. Impact energies from 11.2J to 30.5J and frequencies from 0 BPM to 700 BPM were investigated. A reduction in excavation force of as much as 51% was achieved in this experimental investigation. Smaller percussive digging implements, tested by others, have achieved a reduction of as much as 72%. This paper will examine the effects of impact energy, frequency, scaling and their effect on excavation forces in a dry granular material such as lunar regolith. The past several years have shown an increasing interest in mining space resources both for exploration and commercial enterprises. This work studied the benefits and risks of percussive excavation and prelimin~ry results indicate that this technique may become an enabling technology for extra-terrestrial excavation of regolith and ice.

Schuler, Jason↗
DOE OSTI2021

Parallel projection—An improved return mapping algorithm for finite element modeling of shape memory alloys

Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.

42 ENGINEERING↗