Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time to failure”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Current Status and Future Directions of All-Solid-State Batteries with Lithium Metal Anodes, Sulfide Electrolytes, and Layered Ternary Oxide Cathodes

All-solid-state batteries (ASSBs) offer great promise as a next-generation energy storage technology with higher energy density, wider operating temperature range, and improved safety for electric vehicles. ASSBs employing lithium metal anodes (Li), sulfide-based solid-state electrolytes (SSE), and Ni-rich layered transition metal oxide cathodes (LiMO 2 , M = Ni, Mn, Co, Al) are particularly promising due to its superior electrochemical performance compared to other solid-electrolyte systems. However, the battery cycle life at high cathode mass loading and high current is still limited because the failure mechanism is not fully understood. Lithium dendrite growth at the anode or inside a solid electrolyte still represents as a serious risk of cell failure. Interfacial resistance increases attributed to electrolyte decomposition and interfacial void formation at both cathode–electrolyte and anode–electrolyte interfaces lead to gradual capacity fading. In this Review, we present the fundamental challenges and recent scientific understandings of each component in ASSBs. The novel diagnostic tools for these components, especially the interfaces buried under the surface that are often hard for characterization are mainly examined. Lastly, we offer a perspective for future research directions. We hope this Review will provide a timely snapshot of state-of-the-art research progress in ASSBs to accelerate the development of ASSBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

DAQExpert the service to increase CMS data-taking efficiency

The Data Acquisition (DAQ) system of the Compact Muon Solenoid (CMS) experiment at the LHC is a complex system responsible for the data readout, event building and recording of accepted events. Its proper functioning plays a critical role in the data-taking efficiency of the CMS experiment. In order to ensure high availability and recover promptly in the event of hardware or software failure of the subsystems, an expert system, the DAQ Expert, has been developed. It aims at improving the data taking efficiency, reducing the human error in the operations and minimising the on-call expert demand. Introduced in the beginning of 2017, it assists the shift crew and the system experts in recovering from operational faults, streamlining the post mortem analysis and, at the end of Run 2, triggering fully automatic recovery without human intervention. DAQ Expert analyses the real-time monitoring data originating from the DAQ components and the high-level trigger updated every few seconds. It pinpoints data flow problems, and recovers them automatically or after given operator approval. We analyse the CMS downtime in the 2018 run focusing on what was improved with the introduction of automated recovery; present challenges and design of encoding the expert knowledge into automated recovery jobs. Furthermore, we demonstrate the web-based, ReactJS interfaces that ensure an effective cooperation between the human operators in the control room and the automated recovery system. We report on the operational experience with automated recovery.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Independence of Environmental and Mechanical Damages on Silicone Adhesive Stored in a Thermo-Oxidative Environment

Abstract Elastomeric polymer materials hold a special place in today’s industrial sectors with respect to structural application needs. Low quality materials are detrimental to industries like aerospace, and automotive engineering. Assessment of the durability of elastomers for structural applications has been of much interest among the literary circles for decades and new materials keep outperforming the existing ones. Polymeric adhesives are one of the most abundantly used materials in these industries. All polymeric materials get damaged when in contact with aggressive environments in presence of high temperature and oxygen. Commonly referred to as thermo-oxidation, this environment exposes the material to heat and oxidation reactions in presence of oxygen. Resultantly, during service life, the damage to the polymer matrix is primarily caused by two factors: mechanical damage and environmental aging. Environmental aging is an irreversible phenomenon caused by changes in the molecular structure while mechanical damage maintains the shape of the polymer matrix, and the deterioration is mostly due to polymer chain mobility. Environmental aging can be caused by a single environmental agent or by a synergized impact of several environmental elements. Increasing temperature is found to be proportional to decreasing tensile strength and toughness of material. The rate and extent of degradation can be accessed by scrutinizing the changes in constitutive behavior of material through mechanical and chemical properties. Accelerated thermal aging is among the most common modes of process related degradation, leading generally to chain scission, and crosslinking phenomena and reduced resistance to fracture stress, and strain. In this experimental study, our goal is to separate the environmental degradation from mechanical damage. A silicone-based adhesive was aged in thermo-oxidative at (0%RH) aging environment. The damage and decay mechanisms have been used to draw a distinction between environmental degradation and mechanical damage. Material characterization included uniaxial tensile test (failure and cyclic) and scanning electron microscopy (SEM) tests on as-received and aged samples. Aging was conducted at three different temperatures (60°C, 80°C and 95°C) and six different exposure durations (1, 3, 10, 30, 90 and 200 days). This work confirms that environmental damage is superposed on top of the mechanical damage, and thus, they are separable. The total mechanical and environmental damage is a synergized effect of all exposure conditions and parameters involved i.e., aging time, temperature, and oxygen. The chemistry and mechanics of the polymer degradation were found to be in good agreement with each other.

Alazhary, Sharif↗
DOE OSTI2021

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗
DOE OSTI2023

Neutron Imaging of Al6061 Prepared by Solid-State Friction Stir Additive Manufacturing

Solid-state Friction Stir Additive Manufacturing has recently gained attention as a result of its capacity to fabricate large-scale parts while preserving the mechanical properties of the feedstock material. However, the correlation between the quality of layer-by-layer bonding of the deposited metal and processing parameters has remained unknown. Neutron imaging techniques, with 90% total transmission per cm, are employed for Al6061 parts fabricated by MELD ® Technology as a non-destructive evaluation approach for the first time to investigate the layer-by-layer structure of a stadium-shaped ingot in different sections. The post-processed results show the fabricated parts with an optimized set of processing parameters are void-free. However, the hydrocarbon-based feedstock lubricant segregates between the layers, which consequently may lead to non-uniform weaker mechanical properties along the build direction and stimulate crack initiation during mechanical loading. The tensile test results show 14% lower strain-to-failure values in alleged contaminated areas in transmission imaging results. Additionally, layer bonding is significantly impacted by hot-on-hot and hot-on-cold layer deposition schemes, especially for larger layer thicknesses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AI for Interpreting Nuclear Power Plant Documents for Power Uprates

To reduce the cost and time needed for regulatory compliance, nuclear power plants (NPPs) can utilize artificial intelligence (AI) to assist in interpreting complex and voluminous documents that typically span thousands of pages. Usually, the process of interpreting a plant’s technical specifications (TSs) and associated documents is labor intensive. This study aims to understand what processes state-of-the-art large language models (LLMs) can automate and to identify the pitfalls associated with using LLMs to reduce human labor costs and time. This research uses a recent AI technology called retrieval augmented generation (RAG), which retrieves pages of information from TSs and associated documents to assist with NPP power uprates (cleared to produce more power). LLMs are integral to RAG because they create human-like responses based on the retrieved information, aiding in the interpretation and application processes. A baseline case demonstrates how LLMs can operate successfully for a power uprate application. Then five use cases show five types of potential failures: (1) RAG retrieving the incorrect information, (2) RAG misinterpreting the retrieved information, (3) RAG relying on knowledge not contained in the retrieved information, (4) RAG hallucinating, and (5) RAG refusing to answer. The results of the five use cases suggest that automating the human interpretation of TSs and associated documents with AI should be approached with caution. A subject-matter expert reviewed the AI outputs from the five use cases and concluded that an LLM can produce technical information that is needed to produce power uprate applications in certain instances.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Optimal Gas-Electric Energy System Decarbonization Planning

As energy utilities implement climate change mitigation policies, system planners require strategies for achieving affordable emissions reductions. Coordinated planning of electric power and natural gas systems will allow synergistic investments to address cross-sector operational constraints, competing uses for net-zero emissions fuels, and shifts in energy demands across energy carriers. In this study, we develop a novel optimization program that finds the cost-minimizing mix of infrastructure expansion or reduction across gas and electric systems to satisfy sector-specific emissions constraints. Alongside energy supply resources, our framework allows for central-planning of end-use equipment stocks to allow switching between gas and electric appliances upon failure or premature replacement. The proposed model is used to simulate case study scenarios for a benchmark 24-pipe gas network coupled to a 24-node power system test network. We find that electrification of greater than 80% of core gas demands is a component of the least-cost solution for modeled energy systems. Despite this substitution, the gas system is maintained to service difficult-to-electrify customers and to deliver net-zero emissions gas to electricity generators in times of peak electricity demand. Restricting electrification of gas appliances increases reliance on power- to-gas technologies and increases annual costs by 15% in 2040. Neglecting constraints on pipeline blending of hydrogen can produce a misleading result that relies on hydrogen blend fractions of greater than 50%. In all cases, we find the average costs of delivered gas increase nearly 5-fold across the decarbonization transition, highlighting the importance of future work investigating cost-allocation strategies for ensuring an equitable energy transition.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Earth Mover’s Distance as a Metric to Evaluate the Extent of Charge Transfer in Excitations Using Discretized Real-Space Densities

This paper presents a novel theoretical measure, μ EMD , based on the earth mover's distance (EMD), for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground- and excited-state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. The EMD metric is compared against other popular theoretical metrics for describing the extent of electron-hole separation in a wide range of excited states (valence, Rydberg, charge transfer, etc.). Further, the results showcase the EMD metric's effectiveness across all excitation types and suggest that it is useful as an additional tool to characterize electronic excitations. The study also reveals that μ EMD can function as a promising diagnostic tool for predicting the failure of pure exchange-correlation functionals. Specifically, we show statistical relationships among the functional-driven errors, the exact exchange content within the functional, and the magnitude of μ EMD values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Smart manufacturing approach to manufacture bulk nanocrystalline aluminum for lightweight applications

In this research, a smart manufacturing approach was used to enhance the mechanical properties of aluminum (Al) for lightweight applications. The smart manufacturing involved cryomilling of Al powders with and without 5 wt.% magnesium (Mg) powders for varying durations followed by a high-pressure cold spray (HPCS) additive manufacturing process to prepare bulk components. The morphological changes, crystallite size, and composition of the cryomilled powders and cold sprayed (CS’ed) components were examined using scanning electron microscopy (SEM), x-ray diffraction (XRD), and transmission electron microscopy (TEM) techniques. The results showed that the crystallite size reduces with an increase in cryomilling time and the addition of Mg dopant. To test the mechanical properties of the bulk CS’ed components, microhardness tests were performed using a Vickers microhardness tester. Uniaxial tensile tests were also carried out to ascertain the material’s tensile properties. The mechanical testing results showed great improvement in the hardness and tensile strength of CS’ed Al–Mg samples as compared to pure Al samples. Subsequently, fractography analysis of the tensile failed samples was carried out to determine the nature of the failure. Here, the research article also discusses the inherent mechanisms for the improvement in mechanical properties of smart manufactured components as a result of Mg doping and cryomilling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid data-driven and model-informed online tool wear detection in milling machines

Precision machining tool wear is responsible for low product throughput and quality. Monitoring the tool wear online is vital to prevent degradation in machining quality. However, direct real-time tool wear measurement is not practical. This paper presents residual-based anomaly detection models, combining a hybrid model comprised of a physics-based model and a data-driven model (a decision tree or a neural network) to predict signals of interest (e.g., power or forces) under nominal conditions, followed by Page’s cumulative sum test for detecting tool wear on-line using the computer numerical control machine measurements. The most informative features are ranked using dynamic programming and its approximation variants from real-time measurements and machine settings, such as the width of cut, depth of cut, feed rate and spindle speed, that serve as inputs to the predictive models. The baseline nominal model is incrementally updated with experimental data via a gradient boosted adaptation model to generate the residuals that account for discrepancies between the actual machine data under normal conditions and the baseline nominal model predictions. The hybrid model is validated against 20 Mazak milling machine experimental tests and one Haas run-to-failure experiment. The proposed anomaly detector is applied to synthetic data from simulations of the physics-based model at different operating conditions, measurement noise levels, and tool wear levels, and the methods were able to achieve an overall 92% accuracy in data with 1% noise. The anomaly detection methods based on hybrid model reduced the false alarms of either the data-driven or physical-based models alone, and are found to be capable of good online detection of tool wear.

42 ENGINEERING↗
DOE OSTI2025

Ion-Exchange Membrane-Centric Durability Testing and Degradation Characterization for Industry-Relevant CO 2 Reduction

Electrochemical CO 2 reduction is a promising conversion process for producing value-added fuels and chemicals from electricity and CO 2 as a sustainable carbon feedstock to domestically produce fuels and chemicals from industrial waste. Having reached industrially viable performance metrics with small-scale CO 2 electrolysis cells, the field must now increasingly focus on extending the device durability of large stacks to achieve equivalent metrics for 35,000+ hours to decrease maintenance and capital costs. Reported device lifetimes have increased in recent years, with the longest stability studies for CO, ethylene, and formic acid production being published in 2024–2025 with operation times of 4500, 1000, and 5200 h, respectively. Unfortunately, significant extension of the device durability is still required. Here, we provide an overview of ion-exchange membranes (IEMs) and provide insight into the variety of degradation mechanisms that must be overcome to enable the community to meet durability targets. In an effort to accelerate the extension of device lifetimes, we propose a general approach for characterizing CO 2 electrolysis cell degradation before and after durability testing to better elucidate the mechanisms and failure modes of IEMs in zero-gap cells. Furthermore, we encourage the adoption of operando characterizations in tandem with accelerated stress and durability tests, postulating that their combined applications will be increasingly valuable. We hope that this perspective motivates future durability studies to evaluate degradation across the entire electrolysis cell.

09 BIOMASS FUELS↗
DOE OSTI2025

Evaluation of LLM-Generated Kokkos Code Using Compile-Time and Run-Time Testing

Due to the growing use of large language models (LLMs) by developers and researchers, it has become essential to reliably evaluate their ability to generate code that uses specialized libraries. We explore the use of compile-time and run-time evaluation of LLM-generated Kokkos code through extending the methods used by OpenAI with the HumanEval dataset. Our evaluation framework is based on the first 40 prompts from the Kokkos138 dataset. We start by discussing two different forms of LLM prompting, using entirely plain English or providing pseudocode for added context. These two methods are used to generate Kokkos code with the Llama-3.1-8B-Instruct and CodeQwen1.5-7B-Chat models. We found that both forms of prompting led to high failure rates and difficulties with reliably parsing LLM-generated code, while prompts with pseudocode for context generally led to improved results on more complicated tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Netostat: analyzing dynamic flow patterns in high-speed networks

Understanding flow traffic patterns in networks, such as the Internet or service provider networks, is crucial to improving their design and building them robustly. However, as networks grow and become more complex, it is increasingly cumbersome and challenging to study how the many flow patterns, sizes and the continually changing source-destination pairs in the network evolve with time. Here, we present Netostat, a visualization-based network analysis tool that uses visual representation and a mathematics framework to study and capture flow patterns, using graph theoretical methods such as clustering, similarity and difference measures. Netostat generates an interactive graph of all traffic patterns in the network, to isolate key elements that can provide insights for traffic engineering. We present results for U.S. and European research networks, ESnet and GEANT, demonstrating network state changes, to identify major flow trends, potential points of failure, and bottlenecks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Prediction of Thermal Conditions of DED With FEA Metal Additive Simulation

This paper presents the integration of wire-arc additive manufacturing (WAAM) using Gas Metal Arc Welding (GMAW) into a machine tool to create a retrofit hybrid computer numeric control (CNC) machine tool. GMAW, along with other direct energy deposition systems, has the capacity to deposit material faster than the excess thermal energy can dissipate. This results in the need to allow the part to cool between consecutive layers, which is the most time-consuming part of the additive process. Finite element analysis (FEA) was used in conjunction with monitored build plate surface temperatures during deposition samples to improve adequate dwell time prediction and to develop a cooling system. A deposition was completed where no dwell time was used and the build plate along with the machine table temperatures were monitored. A second deposition was completed where only one bead was deposited and the traverse speed was increased. The GMAW welder was mounted on a 3-axis CNC machine where two square deposition samples were completed. A FEA model was designed and verified using the monitored samples. The model will be used to determine improved depositions speeds and whether forced cooling would allow for an increased deposition rate without structural failure. It was determined the FEA software can be used to accurately model and predict the thermal response of WAAM AM components.

Heinrich, Lauren↗
DOE OSTI2024

Multi-modal characterisation of early-stage, subclinical cardiac deterioration in patients with type 2 diabetes

Abstract Background Type 2 diabetes mellitus (T2DM) is a major risk factor for heart failure with preserved ejection fraction and cardiac arrhythmias. Precursors of these complications, such as diabetic cardiomyopathy, remain incompletely understood and underdiagnosed. Detection of early signs of cardiac deterioration in T2DM patients is critical for prevention. Our goal is to quantify T2DM-driven abnormalities in ECG and cardiac imaging biomarkers leading to cardiovascular disease. Methods We quantified ECG and cardiac magnetic resonance imaging biomarkers in two matched cohorts of 1781 UK Biobank participants, with and without T2DM, and no diagnosed cardiovascular disease at the time of assessment. We performed a pair-matched cross-sectional study to compare cardiac biomarkers in both cohorts, and examined the association between T2DM and these biomarkers. We built multivariate multiple linear regression models sequentially adjusted for socio-demographic, lifestyle, and clinical covariates. Results Participants with T2DM had a higher resting heart rate (66 vs. 61 beats per minute,p < 0.001), longer QTc interval (424 vs. 420ms,p < 0.001), reduced T wave amplitude (0.33 vs. 0.37mV,p < 0.001), lower stroke volume (72 vs. 78ml,p < 0.001) and thicker left ventricular wall (6.1 vs. 5.9mm,p < 0.001) despite a decreased Sokolow-Lyon index (19.1 vs. 20.2mm,p < 0.001). T2DM was independently associated with higher heart rate (beta = 3.11, 95% CI = [2.11,4.10],p < 0.001), lower stroke volume (beta = −4.11, 95% CI = [−6.03, −2.19],p < 0.001) and higher left ventricular wall thickness (beta = 0.133, 95% CI = [0.081,0.186],p < 0.001). Trends were consistent in subgroups of different sex, age and body mass index. Fewer significant differences were observed in participants of non-white ethnic background. QRS duration and Sokolow-Lyon index showed a positive association with the development of cardiovascular disease in cohorts with and without T2DM, respectively. A higher left ventricular mass and wall thickness were associated with cardiovascular outcomes in both groups. Conclusion T2DM prior to cardiovascular disease was linked with a higher heart rate, QTc prolongation, T wave amplitude reduction, as well as lower stroke volume and increased left ventricular wall thickness. Increased QRS duration and left ventricular wall thickness and mass were most strongly associated with future cardiovascular disease. Although subclinical, these changes may indicate the presence of autonomic dysfunction and diabetic cardiomyopathy. Graphical Abstract

Cardiovascular System & Cardiology↗
DOE OSTI2021

A Review of Existing and Emerging Methods for Lithium Detection and Characterization in Li-Ion and Li-Metal Batteries

Whether attempting to eliminate parasitic Li metal plating on graphite (and other Li-ion anodes) or enabling stable, uniform Li metal formation in ‘anode-free’ Li battery configurations, the detection and characterization (morphology, microstructure, chemistry) of Li that cannot be reversibly cycled is essential to understand the behavior and degradation of rechargeable batteries. In this review, various approaches used to detect and characterize the formation of Li in batteries are discussed. Each technique has its unique set of advantages and limitations, and works towards solving only part of the full puzzle of battery degradation. Going forward, multimodal characterization holds the most promise towards addressing two pressing concerns in the implementation of the next generation of batteries in the transportation sector (viz. reducing recharging times and increasing the available capacity per recharge without sacrificing cycle life). Such characterizations involve combining several techniques (experimental- and/or modeling-based) in order to exploit their respective advantages and allow a more comprehensive view of cell degradation and the role of Li metal formation in it. Additionally, it is also discussed which individual techniques, or combinations thereof, can be implemented in real-world battery management systems on-board electric vehicles for early detection of potential battery degradation that would lead to failure.

25 ENERGY STORAGE↗
DOE OSTI2021

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗