Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Load Modeling”

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,477 records · Page 82

Multiscale Static Analysis of Notched and Unnotched Laminates Using the Generalized Method of Cells

The generalized method of cells (GMC) is demonstrated to be a viable micromechanics tool for predicting the deformation and failure response of laminated composites, with and without notches, subjected to tensile and compressive static loading. Given the axial [0], transverse [90], and shear [+45/-45] response of a carbon/epoxy (IM7/977-3) system, the unnotched and notched behavior of three multidirectional layups (Layup 1: [0,45,90,-45](sub 2S), Layup 2: [0,60,0](sub 3S), and Layup 3: [30,60,90,-30, -60](sub 2S)) are predicted under both tensile and compressive static loading. Matrix nonlinearity is modeled in two ways. The first assumes all nonlinearity is due to anisotropic progressive damage of the matrix only, which is modeled, using the multiaxial mixed-mode continuum damage model (MMCDM) within GMC. The second utilizes matrix plasticity coupled with brittle final failure based on the maximum principle strain criteria to account for matrix nonlinearity and failure within the Finite Element Analysis--Micromechanics Analysis Code (FEAMAC) software multiscale framework. Both MMCDM and plasticity models incorporate brittle strain- and stress-based failure criteria for the fiber. Upon satisfaction of these criteria, the fiber properties are immediately reduced to a nominal value. The constitutive response for each constituent (fiber and matrix) is characterized using a combination of vendor data and the axial, transverse, and shear responses of unnotched laminates. Then, the capability of the multiscale methodology is assessed by performing blind predictions of the mentioned notched and unnotched composite laminates response under tensile and compressive loading. Tabulated data along with the detailed results (i.e., stress-strain curves as well as damage evolution states at various ratios of strain to failure) for all laminates are presented.

GMC↗

Boride-Carbon Hybrid Technology to Produce Ultra-Wear and Corrosion Resistant Surfaces for Applications in Harsh Conditions (Final Technical Report)

Engineered functional surfaces play an important role to enable new products and manufacturing processes that can endure harsh service conditions such as high impact and contact loads, highly abrasive wear, extreme temperatures, and corrosive environments. Engineered surfaces can also be instrumental to improving the efficient use of energy by reducing frictional losses and extending service life. The main objective of this project was to develop a hybrid surface engineering technology that combines the advantages of a novel ultrafast boriding process with the next generation of superhard carbon coatings. The hypothesis was that this hybrid process will offer an unprecedented combination of wear and corrosion resistance, low frictional losses and affordability for treated parts so that it can be utilized in many applications. During this project, a duplex process was developed that combines the advantages of ultra-fast electrochemical boriding with those of hard tetrahedral amorphous carbon coatings. Both technologies can be combined to form a hybrid technology that is characterized by low friction and wear properties combined with corrosion and fatigue resistance. Good adhesion of both layers to each other was one main goal of this project, that has been achieved with HF1 adhesion through the Rockwell-C adhesion test. In this project, the mechanical properties of the hybrid coating were modeled through a finite-element analysis approach. We can conclude that the FEA model resembles the actual samples and be utilized to predict mechanical behavior under impacts. Based on this model, application-oriented load conditions can be simulated for optimal layer design regarding thickness and mechanical properties. To exemplify, one conclusion that can be drawn from the nanoindentation model is a boride layer thickness of 50 µm is sufficient to effectively support the carbon coating on the identified AISI 1045 low carbon steel substrate material. The duplex treatment yields wear rates as low as 6 x 10 -8 mm 3 N -1 m -1 and a coefficient of friction of 0.14 when tested against a steel counter face in a ball-on-disk test setup. On the other hand, the wear rate of the only-borided AISI 1045 steel was 5 x 10 -5 mm 3 N -1 m -1 , about three orders of magnitude higher than the duplex coating. At the same time, duplex treated samples experience corrosion resistance, which could not be achieved with single-layer carbon coatings. The developed surface treatment withstands a 3-hour exposure to 15% HCl, while the only carbon coated counter sample shows severe delamination of the coating due to pin hole corrosion. The boride layer is chemically stable and pin hole free because it is formed through an electrochemical process under high current densities (700 mA/cm 2 ) and high temperature. Additionally, the hybrid coating led to at least 3x increase in fatigue strength of the steel substrate, which exceeds the target performance of 30% improvement. There are numerous potential applications for the duplex coatings. A representative application is bearing ball coatings for off-shore windmills. Compared to currently employed surface technologies in this field, the initial costs of applying our technology might be higher due to more process steps but the performance benefits lead to an increased life time of treated parts, which will lower the maintenance and replacement costs in the long-term. To validate the technology for this specific application, the team is currently investigating the process of white etching crack initiation of the duplex coating in collaboration with ANL. Overall, this project successfully validated that the boride-carbon hybrid technology can withstand harsh conditions. One possible approach to commercialization under consideration is to transfer the technology to a startup or an existing coatings company.

36 MATERIALS SCIENCE↗

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗

Application of Model-based Prognostics to a Pneumatic Valves Testbed

Pneumatic-actuated valves play an important role in many applications, including cryogenic propellant loading for space operations. Model-based prognostics emphasizes the importance of a model that describes the nominal and faulty behavior of a system, and how faulty behavior progresses in time, causing the end of useful life of the system. We describe the construction of a testbed consisting of a pneumatic valve that allows the injection of faulty behavior and controllable fault progression. The valve opens discretely, and is controlled through a solenoid valve. Controllable leaks of pneumatic gas in the testbed are introduced through proportional valves, allowing the testing and validation of prognostics algorithms for pneumatic valves. A new valve prognostics approach is developed that estimates fault progression and predicts remaining life based only on valve timing measurements. Simulation experiments demonstrate and validate the approach.

prediction↗

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗

Structural Analysis Using Computer Based Methods

The stiffness of a flex hose that will be used in the umbilical arms of the Space Launch Systems mobile launcher needed to be determined in order to properly qualify ground umbilical plate behavior during vehicle separation post T-0. This data is also necessary to properly size and design the motors used to retract the umbilical arms. Therefore an experiment was created to determine the stiffness of the hose. Before the test apparatus for the experiment could be built, the structure had to be analyzed to ensure it would not fail under given loading conditions. The design model was imported into the analysis software and optimized to decrease runtime while still providing accurate rest~lts and allow for seamless meshing. Areas exceeding the allowable stresses in the structure were located and modified before submitting the design for fabrication. In addition, a mock up of a deep space habitat and the support frame was designed and needed to be analyzed for structural integrity under different loading conditions. The load cases were provided by the customer and were applied to the structure after optimizing the geometry. Once again, weak points in the structure were located and recommended design changes were made to the customer and the process was repeated until the load conditions were met without exceeding the allowable stresses. After the stresses met the required factors of safety the designs were released for fabrication.

Dietz, Matthew R.↗

Photogrammetry of a Hypersonic Inflatable Aerodynamic Decelerator

In 2012, two large-scale models of a Hypersonic Inflatable Aerodynamic decelerator were tested in the National Full-Scale Aerodynamic Complex at NASA Ames Research Center. One of the objectives of this test was to measure model deflections under aerodynamic loading that approximated expected flight conditions. The measurements were acquired using stereo photogrammetry. Four pairs of stereo cameras were mounted inside the NFAC test section, each imaging a particular section of the HIAD. The views were then stitched together post-test to create a surface deformation profile. The data from the photogram- metry system will largely be used for comparisons to and refinement of Fluid Structure Interaction models. This paper describes how a commercial photogrammetry system was adapted to make the measurements and presents some preliminary results.

aerodynamic brakes↗

Load Diffusion in Composite Structures

This research has been concerned with load diffusion in composite structures. Fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses. The decay behavior of stresses and other field quantities provides a significant aid towards this process. The results are also amendable to parameter study with a large parameter space and should be useful in structural tailoring studies.

Horgan, Cornelius O.↗

Potential of Data Center Controls in Grid Services

The rapid proliferation of large data centers brings both challenges and opportunities for grid reliability. The data center resources and their potential flexibility have the potential to contribute resources to grid operations. Through capabilities like energy shifting and resource coordination, data centers can help reduce their net demand on the transmission network, as well as provide additional grid services to support reliable operation on the grid. While transient and long-term grid planning and operations are the scenarios that draw most attention, the quasi-steady state timeseries (QSTS) operation of data centers and grid bring interesting scenarios that can help evaluate the data center controls to aid grid services. This work is focused on modeling data centers for QSTS applications – incorporating the AI data center load profiles and building on the PNNL digital twin model for the thermal management loads to enable simulation studies to reveal the impact of data center controls on grid performance. This includes the integration of a QSTS battery and natural gas generator model to incorporate local resource impacts to the system. The simulation study is performed with a modified IEEE 24-Bus transmission system. Scenarios are focused on evaluating the data center load impacts on the transmission system and leveraging both data center and local generation controls to mitigate those impacts and provide additional grid services. The data center controls revealed the ability to contribute to two main kinds of grid services: preventing congestion on a weak grid by coordinating the data center resources with the collocated BESS and onsite generation; and the ability to help the grid operations during stressed times of operation like during a contingency. Leveraging these and other capabilities has the potential to help data centers become grid responsive assets, aiding in both their integration into the power system and grid reliability.

power grid simulation↗

A Hybrid Reinforcement Learning-MPC Approach for Distribution System Critical Load Restoration

This paper proposes a hybrid control approach for distribution system critical load restoration, combining deep reinforcement learning (RL) and model predictive control (MPC) aiming at maximizing total restored load following an extreme event. RL determines a policy for quantifying operating reserve requirements, thereby hedging against uncertainty, while MPC models grid operations incorporating RL policy actions (i.e., reserve requirements), renewable (wind and solar) power predictions, and load demand forecasts. We formulate the reserve requirement determination problem as a sequential decision-making problem based on the Markov Decision Process (MDP) and design an RL learning environment based on the OpenAI Gym framework and MPC simulation. The RL agent reward and MPC objective function aim to maximize and monotonically increase total restored load and minimize load shedding and renewable power curtailment. The RL algorithm is trained offline using a historical forecast of renewable generation and load demand. The method is tested using a modified IEEE 13-bus distribution test feeder containing wind turbine, photovoltaic, microturbine, and battery. Case studies demonstrated that the proposed method outperforms other policies with static operating reserves.

distribution system↗

Mechanical Comparison of Short Models of Nb$_3$ Sn Low-$\beta$ Quadrupole for the Hi-Lumi LHC

MQXF is the Nb 3 Sn Low-β quadrupole magnet that the HL-LHC project is planning to install in the LHC interaction regions in 2026 to increase the LHC integrated luminosity. The magnet will be fabricated in two different lengths: 4.2 m for MQXFA, built in the US by the Accelerator Upgrade Project (AUP), and 7.15 m for MQXFB, fabricated by CERN. In order to qualify the magnet design and characterize its performance with different conductors, cable geometries and pre-load configurations, five short model magnets, called MQXFS, were fabricated, assembled and tested. We compare the mechanical behavior of short model magnets using experimental data and new numerical models that take into account the measured coil sizes as a function of position.

43 PARTICLE ACCELERATORS↗

Effect of Loading Configuration on Kinematics and Kinetics for Deadlift and Squat Exercises: Case Study in Modeling Exercise Countermeasure Device for Astronauts

This study compares squat and deadlift exercises performed with two different loading configurations: 1) on a novel single-cable resistance exercise countermeasure device (ECD) for spaceflight and 2) with free weights. The results compare joint kinematics and kinetics between different loading configurations for each exercise, and also between the two exercises for each loading configuration. Single-cable versions of the squat (using a harness) and deadlift (using a T-bar) performed on the Hybrid Ultimate Lifting Kit (HULK) ECD have significantly different sagittal plane joint angle kinematics (both peak angle and range of motion) as well as joint kinetics (both peak joint moment and joint impulse) vs. their free weight equivalents at the same load. Differences also exist in hip abduction and rotation. Overall, the single-cable configurations tend to reduce peak joint angles, ranges of motion, peak joint moment and joint impulse vs. free weights. A notable exception is the lumbar joint, which is more heavily loaded for single-cable squats vs. free weight squats. This may have implications for both training benefit and possible risk of injury. Deadlift and squat exercises work the lower body musculature in different ways, with the deadlift emphasizing hip and lumbar extension and the squat emphasizing knee extension. Based on these findings, we would advocate the use of both movements in the exercise prescriptions of astronaut crews on deep-space missions.

Resistance Exercise↗

Waters Loads on the XJL-1 Hull as Obtained in Langley Impact Basin, TED No. NACA 2413.3

An investigation was conducted in the Langley impact basin of the water loads on a half scale model of the XJL-1 hull whose forebody has a vee bottom with exaggerated chine flare. The impact loads, moments, and pressures were determined for a range of landing conditions. A normal full-scale landing speed of 86 miles per hour was represented with effective flight paths ranging from 0.6deg to 11.6deg. Landings were made with both fixed trim and free-to-trim mounting of the float over a trim range of -15deg to 12deg into smooth water and into waves having equivalent full-scale length. of 120 feet and heights ranging from 1 to 4 feet. All data and results presented in this report are given in terms of equivalent full-scale values. Summary tables and illustrative plots are used in presenting the material. The following maximum values of load and pressure are those which are apropos for effective flight paths less than 6.5deg which was the maximum value obtained in tests with the XJL-1 hull model representing full-scale landings with vertical velocity of 4.5 feet per second into 4-foot waves. The maximum local pressure on the flat portion of the bottom is 130 pounds per square inch which was measured on a 2-inch-diameter circular area near the step. The maximum local pressure obtained in the curved area near the chines is 200 pounds per square inch. This pressure was also measured near the step.

Steiner, Margaret F.↗

A Hybrid Reinforcement Learning-MPC Approach for Distribution System Critical Load Restoration: Preprint

This paper proposes a hybrid control approach for distribution system critical load restoration, combining deep reinforcement learning (RL) and model predictive control (MPC) aiming at maximizing total restored load following an extreme event. RL determines a policy for quantifying operating reserve requirements, thereby hedging against uncertainty, while MPC models grid operations incorporating RL policy actions, i.e., the reserve requirement, renewable (wind and solar) power predictions, and load demand forecasts. We formulate the reserve requirement determination problem as a sequential decision making problem based on the Markov Decision Process (MDP) and design an RL learning environment based on the OpenAI Gym framework and MPC. The RL agent reward and MPC objective function aim to maximize and monotonically increase total restored load and minimize load shedding and renewable power curtailment. The RL algorithm is trained off-line using historical forecast of renewable generation and load demand. The method is tested using a modified IEEE 13-bus distribution test feeder containing wind turbine, photovoltaic, microturbine and battery. Case studies demonstrated that the proposed method outperforms other operating reserve determination methods.

distribution system↗

Automatic control of a helicopter with a hanging load

An autopilot logic is designed here for controlling a helicopter with a hanging load. A 16th order model for the system is decoupled into four subsystems: (1) a second order system for yawing motion, (2) a second order system for vertical motion, (3) a sixth order system for longitudinal motion, and (4) a sixth order system for lateral motion. A measuring scheme, which could be used in remote areas, is developed and filters are designed to estimate the state variables from these measurements. The autopilot can be used to move the load over short distances without retracting the cables. This is done by automatically shifting the autopilot modes from position-hold (hover) to acceleration-hold to velocity-hold (cruise) to deceleration-hold to velocity-hold (near hover) to position-hold (hover). Use of such an autopilot might save considerable turnaround time. The Sikorsky S-61 helicopter is chosen as an example vehicle. The performance of the controlled system is studied in the presence of longitudinal and lateral winds.

Gupta, N. K.↗

Alfalfa Virtual Building Service: Software Engineering Best Practices Applied to Runtime Interaction with Building Energy Models

Buildings are active participants in increasingly complex energy systems. Building Energy Modeling (BEM) has a key role to play in planning and de-risking an equitable energy transition, with BEM-backed "virtual buildings" critical path for diverse applications that include workforce training tools, Hardware-in-the-Loop (HIL) experimentation to study equipment performance under a range of conditions, Control-Hardware-in-the-Loop (CHIL) experimentation to de-risk commercial control implementations at equipment through grid orchestration levels, and integration of dynamic load profiles into grid modeling tools for energy system experimentation at the urban scale. Modeling requirements vary across these applications, but many software engineering tasks do not. The Alfalfa Virtual Building Service (AVBS, see https://github.com/NREL/alfalfa/wiki) is an open-source web service that solves these common tasks robustly in one place, providing a foundational platform for power users to bootstrap their own applications. AVBS abstracts the specifics of runtime interaction with OpenStudio, Modelica, and Spawn of EnergyPlus models behind a unified REST API. Additionally, AVBS provides resources for cloud deployment and scaling to 100s of parallel simulations, a growing library of modular Operational Technology (OT) integrations for emulation of real-world interfaces, and scripts to automate the population of communities of virtual buildings from URBANopt, ResStock and ComStock.

building automation↗

Less Atmospheric Radiative Heating By Dust Due to the Synergy of Coarser Size and Aspherical Shape

Mineral dust aerosols cool and warm the atmosphere by scattering and absorbing solar (shortwave: SW) and thermal (longwave: LW) radiation. However, significant uncertainties remain in dust radiative effects, largely due to differences in the dust size distribution and spectral optical properties simulated in Earth system models. Dust models typically underestimate the coarse dust load (more than 2.5 µm in diameter) and assume a spherical shape, which leads to an overestimate of the fine dust load (less than 2.5 µm) after the dust emissions in the models are scaled to match observed dust aerosol optical depth at 550 nm (DAOD550). Here, we improve the simulated dust properties with data sets that leverage measurements of size-resolved dust concentration, asphericity factor, and refractive index in a coupled global chemical transport model with a radiative transfer module. After the adjustment of size-resolved dust concentration and spectral optical properties, the global and annual average of DAOD550 from the simulation increases from 0.023 to 0.029 and falls within the range of a semi-observationally based estimate (0.030 ± 0.005). The reduction of fine dust load after the adjustment leads to a reduction of the SW cooling at the top of the atmosphere (TOA). To improve agreement against a semi-observationally based estimate of the radiative effect efficiency at TOA, we find that a less absorptive SW dust refractive index is required for coarser aspherical dust. Thus, only a minor difference is estimated for the net global dust radiative effect at TOA (−0.08 vs. −0.00 W m−2 on a global scale). Conversely, our sensitivity simulations reveal that the surface warming is substantially enhanced near the strong dust source regions (less cooling to −0.23 from −0.60 W m−2 on a global scale). Thus, less atmospheric radiative heating is estimated near the major source regions (less heating to 0.15 from 0.59 W m−2 on a global scale), because of enhanced LW warming at the surface by the synergy of coarser size and aspherical shape.

simulated dust properties↗

Inelastic deformation of metal matrix composites

A theoretical model capable of predicting the thermomechanical response of continuously reinforced metal matrix composite laminates subjected to multiaxial loading was developed. A micromechanical model is used in conjunction with nonlinear lamination theory to determine inelastic laminae response. Matrix viscoplasticity, residual stresses, and damage to the fiber/matrix interfacial zone are explicitly included in the model. The representative cell of the micromechanical model is considered to be in a state of generalized plane strain, enabling a quasi two-dimensional analysis to be performed. Constant strain finite elements are formulated with elastic-viscoplastic constitutive equations. Interfacial debonding is incorporated into the model through interface elements based on the interfacial debonding theory originally presented by Needleman, and modified by Tvergaard. Nonlinear interfacial constitutive equations relate interfacial tractions to displacement discontinuities at the interface. Theoretical predictions are compared with the results of an experimental program conducted on silicon carbide/titanium (SiC/Ti) unidirectional, (O4), and angle-ply, (+34)(sub s), tubular specimens. Multiaxial loading included increments of axial tension, compression, torque, and internal pressure. Loadings were chosen in an effort to distinguish inelastic deformation due to damage from matrix plasticity and separate time-dependent effects from time-independent effects. Results show that fiber/matrix debonding is nonuniform throughout the composite and is a major factor in the effective response. Also, significant creep behavior occurs at relatively low applied stress levels at room temperature.

Lissenden, C. J.↗