Results for “Finite element model” — Page 44

DOE OSTI2025

Microstructure Clones

Background: A material’s microstructure drives its material performance. Contemporary crystal plasticity experiments compare full-field strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which features of the observed deformation are deterministic vs statistical; thus, differences between model and experiment may or may not be significant. Objective: This paper introduces the invention of microstructure clones. Microstructure clones are 2D oligocrystal specimens that have nearly identical microstructures to remedy the aforementioned experimental limitations. Having specimens with nearly identical microstructures will allow for multiple destructive tests of a microstructure (either as repeats or intentionally different experiments), an ability to “see the future” by providing insight into how a specimen will deform, variability quantification, and experimental investigations of response to small microstructural changes. Methods: This work introduces microstructure clones. Repeatability of these clones is demonstrated in tensile bars of pure nickel. Local strain measurements from digital image correlation are compared between clone specimens and compared to results from a crystal plasticity finite element model. Results: Two sets of microstructure clones were tested in this study and displayed very consistent deformation responses within each clone set. Small observed differences in deformation invite investigation into microstructure stochasticity and the effect of small microstructural and loading differences. Conclusions: Microstructure clones represent a significant shift in understanding structure–property relationships. This work reshapes experimental crystal plasticity to allow for experiments that control for specific variables, quantification of microstructural stochasticity (and other sources of stochasticity), and opportunities for replicating experiments.

Crystal Plasticity↗

NASA NTRSDate not supplied

Composite Progressive Delamination Analysis of Thin Ply Buckling

The NASA Advanced Composites Project is a multi-year consortium of government, industry and academia to develop and validate technology that will reduce development and certification of advanced composite materials for aircraft structures. One of its tasks is to simulate and predict composite damage progression accurately through high fidelity finite element models. It included a subtask to verify an approach based on the Virtual Crack Closure Technology (VCCT) as well as Floating Node Method (FNM). This paper discusses predicted damage onset and growth in a thin-ply buckling problem. Sensitivity of the simulations to various initial crack length is studied. The compressive load-displacement curves are plotted and the crack growth critical load is predicted. Overall, the crack progression follows simulation benchmarks well. However, the critical-load overshoot phenomenon is observed in both VCCT and FNM, suggesting further investigation may be worth pursuing.

damage progression↗

DOE OSTI2024

Combined ambient vibration and surface displacement measurements for improved progressive failure monitoring at a toppling rock slab in Utah, USA

Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ~37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

58 GEOSCIENCES↗

DOE OSTI2021

Editors’ Choice—Natural Convection Boundary Layer Thickness at Elevated Chloride Concentrations and Temperatures and the Effects on a Galvanic Couple

The natural convection boundary layer ( δ n c ) and its influence on cathodic current in a galvanic couple under varying electrolytes as a function of concentration (1 − 5.3 M NaCl) and temperature (25 °C−45 °C) were understood. Polarization scans were obtained under quiescent conditions and at defined boundary layer thicknesses using a rotating disk electrode on platinum and stainless steel 304L (SS304L); these were combined to determine δ n c . With increasing chloride concentration and temperature, δ n c decreased. Increased mass transport (Sherwood number) results in a decrease in δ n c , providing a means to predict this important boundary. Using Finite Element Modeling, the cathodic current was calculated for an aluminum alloy/SS304L galvanic couple as a function of water layer ( WL ) thickness and cathode length. Electrolyte domains were delineated, describing (i) dominance of ohmic resistance over mass transport under thin WL , (ii) the transition from thin film to bulk conditions at δ n c , and (iii) dominance of mass transport under thick WL . With increasing chloride concentration, cathodic current decreased due to decreases in mass transport. With increasing temperature, increased cathodic current was related to increases in mass transport and solution conductivity. This study has implications for sample sizing and corrosion prediction under changing environments.

25 ENERGY STORAGE↗

NASA NTRS1983

Propagation of flexural and membrane waves with fluid loaded NASTRAN plate and shell elements

Modeling of flexural and membrane type waves existing in various submerged (or in vacuo) plate and/or shell finite element models that are excited with steady state type harmonic loadings proportioned to e(i omega t) is discussed. Only thin walled plates and shells are treated wherein rotary inertia and shear correction factors are not included. More specifically, the issue of determining the shell or plate mesh size needed to represent the spatial distribution of the plate or shell response is of prime importance towards successfully representing the solution to the problem at hand. To this end, a procedure is presented for establishing guide lines for determining the mesh size based on a simple test model that can be used for a variety of plate and shell configurations such as, cylindrical shells with water loading, cylindrical shells in vacuo, plates with water loading, and plates in vacuo. The procedure for doing these four cases is given, with specific numerical examples present only for the cylindrical shell case.

Kalinowski, A. J.↗

DOE OSTI2023

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

DOE OSTI2021

Novel experimental setup for megahertz X-ray diffraction in a diamond anvil cell at the High Energy

Explore the source record for details and available documents.

Liermann, H.P.↗

DOE OSTI2022

Trunk muscle forces and spinal loads while walking in persons with lower limb amputation: Influences of chronic low back pain

Persons with lower limb amputation (LLA) are at high risk for developing chronic low back pain (LBP), often with biomechanical factors considered as likely contributors. Here, trunk and pelvis kinematics, muscle forces, and resultant spinal loads were characterized in persons with LLA, with and without chronic LBP. Thirty-five persons with unilateral LLA – 19 with chronic LBP (“LLA-cLBP”), 16 without LBP (“LLA-nLBP”) – and 15 (uninjured) persons without LBP (“CTR-nLBP”) walked overground (1.3 m/s) while thorax and pelvis kinematics were tracked (and ranges of motion [ROM] computed), and used as inputs for a non-linear finite element model of the spine to estimate global and local muscle forces, and resultant spinal loads. In the frontal and transverse planes, thorax ROM were up to 66.6% smaller in LLA-nLBP versus LLA-cLBP (P < 0.001) and CTR-nLBP (P < 0.001). In the sagittal plane, pelvis ROM was 50.4% smaller in LLA-nLBP versus LLA-cLBP (P = 0.014). LLA-cLBP exhibited 45.5% and 34.2% greater peak local and global muscle forces, respectively, versus CTR-nLBP (P < 0.011). Up to 48.1% greater spinal loads were observed in LLA-cLBP versus CTR-nLBP (P < 0.013); peak compression and local muscle forces were respectively 20.2% and 41.0% larger in LLA-nLBP versus CTR-nLBP (P < 0.005). Despite differences in trunk and pelvis kinematics between LLA-cLBP and LLA-nLBP, trunk muscle forces and spinal loads were similar (P > 0.101) between these groups. Similar loading parameters regardless of LBP presence, while highly dependent on trunk muscle activation strategies, may mitigate further accumulation of mechanical fatigue. Finally, it remains important to understand the temporality of loading with respect to LBP onset following LLA.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRSDate not supplied

Electromagnetic Modeling of Wind tunnel Magnetic Suspension and Balance Systems

An analytical framework for the open-loop behavior of a permanent magnet element levitated within an applied magnetic field are shown, in the context of application to wind tunnel Magnetic Suspension and Balance Systems (MSBS). Various modes of motion are identified, which correlate with observed behavior of the NASA/ODU 6-inch MSBS. This system is being developed as a technology demonstrator for dynamic stability testing of atmospheric entry capsules, with possible application to a supersonic wind tunnel. The analytic equations are populated with field information from experimental measurements and a finite element model. Finally, a system simulation can provide alternative estimates of system dynamic behavior via model linearization, with preliminary comparisons presented. Taken together, the analytical framework, experimental measurement, finite element analysis, and dynamic simulation provide a complete understanding of the characteristics and behavior of the MSBS.

magnetic suspension↗

DOE OSTI2022

Smart material based multilayered microbeam structures for spatial self-deployment and reconfiguration: A residual stress approach

Alleviation of the potentially damaging effects induced by residual stresses was comprehensively investigated in previous research. Here, this paper, however, presents a spatially self-deployable and reconfigurable multilayered microbeam which takes advantage of residual stresses and shape memory effects. Reconfigurable mechanism of a typical four-layered microbeam composed of Pt\Ni 50 Ti 50 \Ni 50 Ti 50 \Pt is introduced, followed by analytical modeling of the maximum distance of the self-deployed gap as functions of variable structural and material parameters, including compressive residual stress in Pt layers and tensile residual stress in Ni 50 Ti 50 layers. Analytical solutions given by the static model agree well with the results obtained via finite element models (FEMs). Fabrication, characterization, and in-situ experiments were carried out to validate the feasibility of deployment of the as-released four-layered microbeam. The maximum distance of the gap was measured to be 41.39 μm at 20 °C, which could be increased to 51.73 μm thanks to controllable reconfiguration driven by shape memory effects. Theoretical analysis of such self-deployment and reconfiguration suggested a tensile residual stress increase by 52 MPa in Ni 50 Ti 50 layers. The multilayered microbeam structure with capabilities of self-deployment and reconfiguration offers great potential for various emerging applications, such as micro robotics, medical drug delivery devices, and intelligent chip scale spacecraft.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Crystal Plasticity Simulation of Magnesium and Its Alloys: A Review of Recent Advances

Slip and extension twinning are the dominant deformation mechanisms in Magnesium (Mg) and its alloys. Crystal plasticity is a powerful tool to study these deformation mechanisms. Different schemes have incorporated crystal plasticity models to capture different properties, which vary from the simple homogenization Taylor model to the full-scale crystal plasticity finite element model. In the current study, a review of works available in the literature that addresses different properties of Mg and its alloys using crystal plasticity modes is presented. In addition to slip and twinning, detwinning is another deformation mechanism that is activated in Mg and its alloys. The different models that capture detwinning will also be addressed here. Finally, the recent experimental frameworks, such as in-situ neutron diffraction, 3D high energy synchrotron X-ray techniques, and digital image correlation under scanning electron microscopy (SEM-DIC), which are incorporated along crystal plasticity models to investigate the properties of Mg and its alloys, are addressed. Future research directions towards improving the deformation response of Mg and its alloys are identified, which can lead to increased deployment of the lightest structural metal in engineering applications.

36 MATERIALS SCIENCE↗

NASA NTRS1978

Simplified modeling for wind turbine modal analysis using NASTRAN

A detailed finite element model of the MOD-0 wind turbine tower was reduced to six beam elements (stick model). The method used to calculate the properties of the beam elements in the stick model was explained and the accuracy of the stick model in predicting natural frequencies and mode shapes was examined. Computer times were compared and several applications where the stick model was used are described. From results obtained from the MOD-0 tower it is concluded that a tower of this type can be modeled as a simple cantilever beam for modal analysis. However, this model should be limited to tower torsional modes and tower bending modes where the mode shape resembles a cantilever beam first bending mode shape.

Sullivan, T. L.↗

NASA NTRS1993

Time-Accurate Unsteady Flow Simulations Supporting the SRM T+68-Second Pressure Spike Anomaly Investigation (STS-54B)

Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The objectives of the present analysis are to: (1) capture the dynamic responses of the motor combustion gas flow to correlate with available low-frequency (less than 12.5 sample/sec) data and (2) observe the high-frequency (up to 50 Hz) characteristics of the response to determine any potentials for dynamic coupling.

Dougherty, N. S.↗

NASA NTRS1996

Structural Evaluation of a Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump Turbine Blade

Thermal and structural finite-element analyses were performed on the first high pressure fuel turbopump turbine blade of the space shuttle main engine (SSME). A two-dimensional (2-D) finite-element model of the blade and firtree disk attachment was analyzed using the general purpose MARC (finite-element) code. The loading history applied is a typical test stand engine cycle mission, which consists of a startup condition with two thermal spikes, a steady state and a shutdown transient. The blade material is a directionally solidified (DS) Mar-M 246 alloy, the blade rotor is forged with waspalloy material. Thermal responses under steady-state and transient conditions were calculated. The stresses and strains under the influence of mechanical and thermal loadings were also determined. The critical regions that exhibited high stresses and severe localized plastic deformation were the blade-rotor gaps.

Abdul-Aziz, Ali↗

NASA NTRS1999

Global-Local Analysis and Optimization of a Composite Civil Tilt-Rotor Wing

This report gives highlights of an investigation on the design and optimization of a thin composite wing box structure for a civil tilt-rotor aircraft. Two different concepts are considered for the cantilever wing: (a) a thin monolithic skin design, and (b) a thick sandwich skin design. Each concept is examined with three different skin ply patterns based on various combinations of 0, +/-45, and 90 degree plies. The global-local technique is used in the analysis and optimization of the six design models. The global analysis is based on a finite element model of the wing-pylon configuration while the local analysis uses a uniformly supported plate representing a wing panel. Design allowables include those on vibration frequencies, panel buckling, and material strength. The design optimization problem is formulated as one of minimizing the structural weight subject to strength, stiffness, and d,vnamic constraints. Six different loading conditions based on three different flight modes are considered in the design optimization. The results of this investigation reveal that of all the loading conditions the one corresponding to the rolling pull-out in the airplane mode is the most stringent. Also the frequency constraints are found to drive the skin thickness limits, rendering the buckling constraints inactive. The optimum skin ply pattern for the monolithic skin concept is found to be (((0/+/-45/90/(0/90)(sub 2))(sub s))(sub s), while for the sandwich skin concept the optimal ply pattern is found to be ((0/+/-45/90)(sub 2s))(sub s).

Rais-Rohani, Masound↗