Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Enhanced Correlation of SMART Active Flap Rotor Loads

This is a follow-on study to a 2010 correlation effort. Measured data from the SMART rotor test in the NASA Ames 40- by 80- Foot Wind Tunnel are compared with CAMRAD II calculations. As background, during the wind tunnel test, unexpectedly high inboard loads were encountered, and it was hypothesized at that time that due to changes in the flexbeams over the years, the flexbeam properties used in the analysis needed updating. Boeing Mesa, recently updated these properties. This correlation study uses the updated flexbeam properties. Compared to earlier studies, the following two enhancements are implemented: i) the inboard loads (pitchcase and flexbeam loads) correlation is included for the first time (reliable prediction of the inboard loads is a prerequisite for any future anticipated flight-testing); ii) the number of blade modes is increased to better capture the flap dynamics and the pitchcase-flexbeam dynamics. Also, aerodynamically, both the rolled-up wake model and the more complex, multiple trailer wake model are used, with the latter slightly improving the blade chordwise moment correlation. This sensitivity to the wake model indicates that CFD is needed. Three high-speed experimental cases, one uncontrolled free flap case and two commanded flap cases, are considered. The two commanded flap cases include a 2o flap deflection at 5P case and a 0o flap deflection case. For the free flap case, selected modifications to the HH-06 section flap airfoil pitching moment table are implemented. For the commanded 2o flap case, the experimental flap variation is approximately matched by increasing the analytical flap hinge stiffness. This increased flap hinge stiffness is retained for the commanded 0o flap case also, which is treated as a free flap case, but with larger flap hinge stiffness. The change in the mid-span and outboard loads correlation due to the updating of the flexbeam properties is not significant. Increasing the number of blade modes results in an effective, commanded flap hinge stiffness of 4X baseline, not 3X as reported earlier. The inboard loads correlation is reasonable, but needs further study. Overall, the free flap case correlation is reasonable, thus confirming the basic correctness of the current semi-empirical modifications; the correlation for the commanded 2o flap at 5P case and the 0o flap case is also reasonable.

Kottapalli, Sesi↗
NASA NTRS2012

Macro Scale Independently Homogenized Subcells for Modeling Braided Composites

An analytical method has been developed to analyze the impact response of triaxially braided carbon fiber composites, including the penetration velocity and impact damage patterns. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. Currently, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. To determine the stiffness and strength properties required for the constitutive model, a top-down approach for determining the strength properties is merged with a bottom-up approach for determining the stiffness properties. The top-down portion uses global strengths obtained from macro-scale coupon level testing to characterize the material strengths for each subcell. The bottom-up portion uses micro-scale fiber and matrix stiffness properties to characterize the material stiffness for each subcell. Simulations of quasi-static coupon level tests for several representative composites are conducted along with impact simulations.

Blinzler, Brina J.↗
NASA NTRS2012

Utilization of the Generalized Method of Cells to Analyze the Deformation Response of Laminated Ceramic Matrix Composites

In order to practically utilize ceramic matrix composites in aircraft engine components, robust analysis tools are required that can simulate the material response in a computationally efficient manner. The MAC/GMC software developed at NASA Glenn Research Center, based on the Generalized Method of Cells micromechanics method, has the potential to meet this need. Utilizing MAC/GMC, the effective stiffness properties, proportional limit stress and ultimate strength can be predicted based on the properties and response of the individual constituents. In this paper, the effective stiffness and strength properties for a representative laminated ceramic matrix composite with a large diameter fiber are predicted for a variety of fiber orientation angles and laminate orientations. As part of the analytical study, methods to determine the in-situ stiffness and strength properties of the constituents required to appropriately simulate the effective composite response are developed. The stiffness properties of the representative composite have been adequately predicted for all of the fiber orientations and laminate configurations examined in this study. The proportional limit stresses and strains and ultimate stresses and strains were predicted with varying levels of accuracy, depending on the laminate orientation. However, for the cases where the predictions did not have the desired level of accuracy, the specific issues related to the micromechanics theory were identified which could lead to difficulties that were encountered that could be addressed in future work.

Goldberg, Robert K.↗
NASA NTRS2013

Experimental and Analytical Characterization of the Macromechanical Response for Triaxial Braided Composite Materials

Increasingly, carbon composite structures are being used in aerospace applications. Their highstrength, high-stiffness, and low-weight properties make them good candidates for replacing many aerospace structures currently made of aluminum or steel. Recently, many of the aircraft engine manufacturers have developed new commercial jet engines that will use composite fan cases. Instead of using traditional composite layup techniques, these new fan cases will use a triaxially braided pattern, which improves case performance. The impact characteristics of composite materials for jet engine fan case applications have been an important research topic because Federal regulations require that an engine case be able to contain a blade and blade fragments during an engine blade-out event. Once the impact characteristics of these triaxial braided materials become known, computer models can be developed to simulate a jet engine blade-out event, thus reducing cost and time in the development of these composite jet engine cases. The two main problems that have arisen in this area of research are that the properties for these materials have not been fully determined and computationally efficient computer models, which incorporate much of the microscale deformation and failure mechanisms, are not available. The research reported herein addresses some of the deficiencies present in previous research regarding these triaxial braided composite materials. The current research develops new techniques to accurately quantify the material properties of the triaxial braided composite materials. New test methods are developed for the polymer resin composite constituent and representative composite coupons. These methods expand previous research by using novel specimen designs along with using a noncontact measuring system that is also capable of identifying and quantifying many of the microscale failure mechanisms present in the materials. Finally, using the data gathered, a new hybrid micromacromechanical computer model is created to simulate the behavior of these composite material systems under static and ballistic impact loading using the test data acquired. The model also quantifies the way in which the fiber/matrix interface affects material response under static and impact loading. The results show that the test methods are capable of accurately quantifying the polymer resin under a variety of strain rates and temperature for three loading conditions. The resin strength and stiffness data show a clear rate and temperature dependence. The data also show the hydrostatic stress effects and hysteresis, all of which can be used by researchers developing composite constitutive models for the resins. The results for the composite data reveal noticeable differences in strength, failure strain, and stiffness in the different material systems presented. The investigations into the microscale failure mechanisms provide information about the nature of the different material system behaviors. Finally, the developed computer model predicts composite static strength and stiffness to within 10 percent of the gathered test data and also agrees with composite impact data, where available.

Littell, Justin D.↗
NASA NTRS2016

Effects of Zoledronate and Mechanical Loading during Simulated Weightlessness on Bone Structure and Mechanical Properties

Space flight modulates bone remodeling to favor bone resorption. Current countermeasures include an anti-resorptive drug class, bisphosphonates (BP), and high-force loading regimens. Does the combination of anti-resorptives and high-force exercise during weightlessness have negative effects on the mechanical and structural properties of bone? In this study, we implemented an integrated model to mimic mechanical strain of exercise via cyclical loading (CL) in mice treated with the BP Zoledronate (ZOL) combined with hindlimb unloading (HU). Our working hypothesis is that CL combined with ZOL in the HU model induces additive structural and mechanical changes. Thirty-two C57BL6 mice (male,16 weeks old, n8group) were exposed to 3 weeks of either HU or normal ambulation (NA). Cohorts of mice received one subcutaneous injection of ZOL (45gkg), or saline vehicle, prior to experiment. The right tibia was axially loaded in vivo, 60xday to 9N in compression, repeated 3xweek during HU. During the application of compression, secant stiffness (SEC), a linear estimate of slope of the force displacement curve from rest (0.5N) to max load (9.0N), was calculated for each cycle once per week. Ex vivo CT was conducted on all subjects. For ex vivo mechanical properties, non-CL left femurs underwent 3-point bending. In the proximal tibial metaphysis, HU decreased, CL increased, and ZOL increased the cancellous bone volume to total volume ratio by -26, +21, and +33, respectively. Similar trends held for trabecular thickness and number. Ex vivo left femur mechanical properties revealed HU decreased stiffness (-37),and ZOL mitigated the HU stiffness losses (+78). Data on the ex vivo Ultimate Force followed similar trends. After 3 weeks, HU decreased in vivo SEC (-16). The combination of CL+HU appeared additive in bone structure and mechanical properties. However, when HU + CL + ZOL were combined, ZOL had no additional effect (p0.05) on in vivo SEC. Structural data followed this trend with ZOL not modulating trabecular thickness in CL + NAHU mice. In summary, our integrated model simulates the combination of weightlessness, exercise-induced mechanical strain, and anti-resorptive treatment that astronauts experience during space missions. Based on these results, we conclude that, at the structural and stiffness level, zoledronate treatment during simulated spaceflight does not impede the skeletal response to axial compression. In contrast to our hypothesis, our data show that zoledronate confers no additional mechanical or structural benefit beyond those gained from cyclical loading.

Mechanical Loading↗
NASA NTRS2014

Experimental Investigation of a Flexible Wing with a Variable Camber Continuous Trailing Edge Flap Design

This paper presents experimental results of a flexible wing wind tunnel model with a variable camber continuous trailing edge flap (VCCTEF) design for drag minimization, tested at the University of Washington Aeronautical Labo- ratory (UWAL). The wind tunnel test was designed to explore the relative merit of the VCCTEF concept for improved cruise efficiency through the use of low-cost aeroelastic model test techniques. The flexible wing model is a 10-scaled model of a typical transport wing and is constructed of woven fabric composites and foam core. The wing structural stiffness in bending is tailored to be half of the stiffness of a Boeing 757-era transport wing, while the torsional stiffness is about the same. This stiffness reduction results in a wing tip deflection of about 10 of the wing semi-span. The VCCTEF is a multi-segment flap design having three chordwise camber segments and five spanwise flap sections for a total of 15 individual flap elements. The three chordwise camber segments can be positioned appropriately to create a desired trailing edge camber. Elastomeric material is used to cover the gaps in between the spanwise flap sections, thereby creating a continuous trailing edge. Wind tunnel data indicate a high degree of data correlation and repeata- bility. The VCCTEF can achieve a drag reduction of up to 6.31 and an improvement in the lift-to-drag ratio (LD) of up to 4.85. The paper also presents two methods for estimating the lift coefficient of a rigid wing using a dynamic pressure correction and an aeroelastic deflection correction. Both methods provide good estimates of the rigid-wing lift coefficient.

Aerodynamics↗
NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗
DOE OSTI2021

Composites and additive manufacturing for high-field coil supports for stellarators

The utilization of composites layered on an additive manufactured substrate, for the production of coil supports for modular coils in small or middle size experimental stellarators is assessed. The focus of the study is a monolithic coil support comprising the coils of a half-period of a stellarator, somewhat similar to the ones in UST_2 and ARIES-CS stellarators. However, the concept may be applicable to quasi-monolithic coil supports (coil forms of NCSX type) or individual coil casings (W7-X type).Coil supports for stellarators require high precision, stiffness and strength for large contorted parts. Traditionally, monolithic coil supports are produced by steel casting/forging and final machining. This production method and material gives accurate, stiff and strong coil supports, but the method may be expensive for monolithic supports due to the geometrical complexity and required accuracy of the structure. In relation to those matters, this work investigates whether a monolithic coil support comprising an additive manufactured resin substrate, which is externally (outward from the coils) surrounded with a thick layer of fibre-reinforced resin, may achieve enough strength and stiffness under middle/high magnetic fields. Finite element calculations are produced to obtain the direction of the principal stresses and their values in compression and tension at different areas of the monolithic support, which is relevant for anisotropic materials. The feasibility of directional application of (carbon) fibres on the winding surface of the stellarator outward from the coils is experimentally tested on a scaled-down additively manufactured prototype of a monolithic support. The strength and stiffness of the composite structure appears sufficient for common magnetic fields in experimental stellarators, and the 3D-composite design and manufacturing was technically feasible.

Queral, Vicente↗
DOE OSTI2022

Characterization of anisotropic mechanical properties of polymer composites from a hybrid additive manufacturing-compression molding process using x-ray computer tomography

In this study, anisotropic stiffness tensors were reconstructed based on fiber orientation distributions obtained from X-ray computer tomography (xCT). A preform was manufactured via a big area additive manufacturing (BAAM) system with carbon fiber (CF) filled acrylonitrile butadiene styrene (ABS). The tailored preform from additive manufacturing (AM) was used in the compression molding (CM) process to produce a low-void high-performance thermoplastic composite panel. An xCT technique was employed to detect the fiber orientations in CF/ABS composites manufactured via three different methods: AM from BAAM, extrusion compression molding (ECM), and AM-CM. The anisotropic stiffness tensor was obtained from the composite panel manufactured via the three manufacturing methods (AM, ECM, and AMCM). A micromechanics theory was used to obtain the orthotropic stiffness tensors of the composite panels and compared with the experimental values. The predicted stiffness tensors of AM and AM-CM composite panels were used to study the deformation characteristics of a steering wheel during airbag deployment by performing finite element analysis (FEA). The approach developed in this study can be utilized for evaluating high-performance composites.

Pokkalla, Deepak↗
DOE OSTI2023

Gearbox assembly

A gearbox assembly includes a first gear coupled to a first shaft, a second gear coupled to a second shaft, and a plurality of planet gears that mesh with the first gear and the second gear. The first shaft includes a first shaft stiff portion and a first shaft flexible portion that includes a greater flexibility than a flexibility of the first shaft stiff portion such that the first gear moves in a radial direction. The second shaft includes a second shaft stiff portion and a second shaft flexible portion that includes a greater flexibility than a flexibility of the second shaft stiff portion such that the second gear moves in the radial direction. The first gear and the second gear move to distribute loads substantially equally on the plurality of planet gears.

O'Connor, Ethan Patrick↗
DOE OSTI2026

Influence of Numerical Modeling Approaches on Damped Behavior of Flexible Beams: Preprint

Composites structures are widely used in aerospace and wind energy applications for their excellent stiffness and strength-to-weight properties. In these structures, structural damping is critical to predict vibration amplitudes, performance, and reliability. Structural damping is of particular interest for slender wings, rotorcraft blades, and wind turbine blades that can exhibit complex vibration phenomena and are frequently modeled with geometrically exact beam theory (GEBT). Standard approaches of stiffness proportional or modal damping merely assign user defined values and cannot predict damping behavior. This work compares stiffness proportional damping to two more advanced damping approaches: modal strain energy and Prony series. The modal strain energy approach uses a sectional analysis tool to calculate the beam stiffness and postprocess internal stresses from GEBT simulations. The internal stresses are then used to calculate modal damping factors. The Prony series is implemented within GEBT to directly model viscoelastic behavior of the composites. These approaches are compared by modeling the evolution of the damping factors of a realistic flexible wind turbine blade with varying rotational speed. Discrepancies between the approaches suggest areas for future modeling development, but differences in nonlinear damping values are less than current uncertainties about the magnitude of structural damping.

17 WIND ENERGY↗
NASA NTRS2001

Virtual reality robotic telesurgery simulations using MEMICA haptic system

The authors conceived a haptic mechanism called MEMICA (Remote Mechanical Mirroring using Controlled stiffness and Actuators) that can enable the design of high dexterity, rapid response, and large workspace haptic system. The development of a novel MEMICA gloves and virtual reality models are being explored to allow simulation of telesurgery and other applications. The MEMICA gloves are being designed to provide intuitive mirroring of the conditions at a virtual site where a robot simulates the presence of a human operator. The key components of MEMICA are miniature electrically controlled stiffness (ECS) elements and electrically controlled force and stiffness (ECFS) actuators that are based on the use of Electro-Rheological Fluids (ERF. In this paper the design of the MEMICA system and initial experimental results are presented.

Haptic interfaces↗
DOE OSTI2021

Ion temperature clamping in Wendelstein 7-X electron cyclotron heated plasmas

The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of Active Lengthening and Shortening on Small-Angle X-ray Reflections in Skinned Skeletal Muscle Fibres

Our purpose was to use small-angle X-ray diffraction to investigate the structural changes within sarcomeres at steady-state isometric contraction following active lengthening and shortening, compared to purely isometric contractions performed at the same final lengths. We examined force, stiffness, and the 1,0 and 1,1 equatorial and M3 and M6 meridional reflections in skinned rabbit psoas bundles, at steady-state isometric contraction following active lengthening to a sarcomere length of 3.0 µm (15.4% initial bundle length at 7.7% bundle length/s), and active shortening to a sarcomere length of 2.6 µm (15.4% bundle length at 7.7% bundle length/s), and during purely isometric reference contractions at the corresponding sarcomere lengths. Compared to the reference contraction, the isometric contraction after active lengthening was associated with an increase in force (i.e., residual force enhancement) and M3 spacing, no change in stiffness and the intensity ratio I 1,1 /I 1,0 , and decreased lattice spacing and M3 intensity. Compared to the reference contraction, the isometric contraction after active shortening resulted in decreased force, stiffness, I 1,1 /I 1,0 , M3 and M6 spacings, and M3 intensity. This suggests that residual force enhancement is achieved without an increase in the proportion of attached cross-bridges, and that force depression is accompanied by a decrease in the proportion of attached cross-bridges. Furthermore, the steady-state isometric contraction following active lengthening and shortening is accompanied by an increase in cross-bridge dispersion and/or a change in the cross-bridge conformation compared to the reference contractions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Carbon Fiber and Syntactic Foam Hybrid Materials via Core–Shell Material Extrusion Additive Manufacturing

Abstract Biological materials often employ hybrid architectures, such as the core–shell (C–S) motif present in porcupine quills and plant stems, to achieve unique specific properties and performance. Drawing inspiration from these natural materials, a new method to fabricate lightweight and stiff C–S architected filaments is reported. Specifically, a C–S printhead conducive to printing highly loaded fiber‐filled inks, as well as a new low‐density syntactic foam ink, are utilized to 3D‐print C–S architectures consisting of a syntactic epoxy foam core surrounded by a stiff carbon fiber‐reinforced epoxy composite shell. Effective printing of test specimens and structures with controlled geometry, composition, and architecture is demonstrated. The new foam ink exhibits density as low as 0.68 g cm −3 and C–S structures exhibit up to 25% higher specific flexural stiffness ( E 1/3 /ρ) than either constituent alone. Finally, a new mechanical model is presented to predict this performance improvement while accounting for potential eccentricity of the core.

Pack, Robert C.↗
DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗