Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite deformation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Creep response of the lunar crust in mare regions from an analysis of crater deformation

The settling trends of 318 lunar mare craters are compared with predictions of numerical finite-element models in order to determine the creep response of the upper lunar mare crust. No settling is evident in craters smaller than 5 km in diameter. Settling rates of larger craters increase as function of crater size in a manner suggesting a nonlinear lunar creep response corresponding to the power law epsilon = 8.3 x 10 to the minus 34th sigma squared where epsilon is the strain rate and sigma is the differential stress. However, the observed nonlinearity is probably an apparent nonlinearity resulting from the temperature-induced viscosity decrease with depth due to a lunar crustal temperature gradient of 3 C/km and a creep activation energy of 20 kcal/mole. It is concluded that creep in the lunar medium is essentially Newtonian, and that the effective viscosity of the upper lunar mare is (1.6 plus or minus 0.3) x 10 to the 25th poise.

Kunze, A. W. G.↗
DOE OSTI2026

Design of monolithic piezoelectric bimorph mirrors made from lithium niobate

Recent advances in deformable mirrors based on monolithic piezoelectric substrates, such as lithium niobate, have the potential to improve the image quality of X-ray optical systems on synchrotron and free-electron laser beamlines and microscopes. However, the quantitative relationship between design parameters and the deformed shape has not been readily available in the literature. We present an analytical model, validated through finite element analysis, enabling calculation of tangential and sagittal curvatures based on mirror dimensions, crystallographic orientation, and applied voltage. We demonstrate that through the selection of material orientation, it is possible to achieve different deformed shapes (sphere, cylinder, or hyperbolic paraboloid). This methodology can be generalized to other piezoelectric materials and substrate-integrated actuator systems.

Marzari, Francesco [University of Trento (Italy); ↗
NASA NTRS1998

Parallel Three-Dimensional Computation of Fluid Dynamics and Fluid-Structure Interactions of Ram-Air Parachutes

This is a final report as far as our work at University of Minnesota is concerned. The report describes our research progress and accomplishments in development of high performance computing methods and tools for 3D finite element computation of aerodynamic characteristics and fluid-structure interactions (FSI) arising in airdrop systems, namely ram-air parachutes and round parachutes. This class of simulations involves complex geometries, flexible structural components, deforming fluid domains, and unsteady flow patterns. The key components of our simulation toolkit are a stabilized finite element flow solver, a nonlinear structural dynamics solver, an automatic mesh moving scheme, and an interface between the fluid and structural solvers; all of these have been developed within a parallel message-passing paradigm.

Tezduyar, Tayfun E.↗
NASA NTRS2002

Residual Stress Reversal in Highly Strained Shot Peened Structural Elements

The purpose of this research was to further the understanding of a crack initiation problem in a highly strained pressure containment housing. Finite Element Analysis methods were used to model the behavior of shot peened materials undergoing plastic deformation. Analytical results are in agreement with laboratory tensile tests that simulated the actual housing load conditions. These results further validate the original investigation finding that the shot peened residual stress had reversed, changing from compressive to tensile, and demonstrate that analytical finite element methods can be used to predict this behavior.

Mitchell, William S.↗
NASA NTRS2015

Gear Tooth Root Stresses of a Very Heavily Loaded Gear Pair-Case Study: Orbiter Body Flap Actuator Pinion and Ring Gear

The space shuttle orbiter's body flap actuator gearing was assessed as a case study of the stresses for very heavily loaded external-internal gear pairs (meshing pinion and ring gear). For many applications, using the high point of single tooth contact (HPSTC) to locate the position of the tooth force is adequate for assessing the maximum tooth root stress condition. But for aerospace gearing such an approach may be inadequate for assessing the stress condition while also simultaneously minimizing mass. In this work specialized contact analyses and finite element methods were used to study gear tooth stresses of body flap actuator gears. The analytical solutions considered the elastic deformations as an inherent part of the solutions. The ratio for the maximum tooth stresses using the HPSTC approach solutions relative to the contact analysis and finite element solutions were 1.40 for the ring gear and 1.28 for the pinion gear.

finite element method↗
NASA NTRS1974

Nonlinear finite-element analysis of laminated composite shells

Two aspects of the geometrically nonlinear analysis of laminated composite shells are considered in this paper. The first pertains to the simplifications and computational advantages resulting from the use of group-theoretic methods in conjunction with computerized symbolic (algebraic) manipulation in the element development. The second aspect pertains to the combined effects of shear deformation and nonlinearity on the accuracy and convergence of two-dimensional, shear-flexible stiffness finite element models.

Noor, A. K.↗
DOE OSTI2020

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗
DOE OSTI2021

Bayesian calibration of strength model parameters from Taylor impact data

Materials strength plays a key role in determining the mechanical response of engineered structures. As such, accurate strength models are crucial in simulations involving complex loading conditions, particularly when deformation in the plastic regime is deemed important. In this work, a Gaussian process based surrogate for a finite element simulation of a Taylor impact test is developed and used for Bayesian calibration of the Preston–Tonks–Wallace strength model. Here, the surrogate model is shown to closely approximate the salient features of the Taylor cylinder deformation and is validated against simulation output before being used in the strength model calibration routine. The results show that Taylor impact test data can be used in the calibration of constitutive equations through the use of a combination of data science techniques, namely Gaussian processes and Bayesian inference.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗
DOE OSTI2020

Crystallographic orientation and spatially resolved damage in a dispersion-hardened Al alloy

The in-situ neutron diffraction technique, in combination with both the full-field crystal elasto-viscoplastic finite element model and microstructural characterization, was used to study the deformation-induced damage anisotropy in a commercial Al alloy, subjected to uniaxial tensile and cyclic loading. The simulations capture well the crystallographic-orientation-dependent lattice strain behavior. The hard grains, e.g. those orientated with the <111> and <422> orientations parallel with the uniaxial loading direction (LD), feature large Taylor factors and seem more prone to form damage-related band structures. Their effective elastic moduli decrease with the accumulation of damage, which are different from the soft grains orientated with the <200> orientation along the LD. Correlation between the distribution of voids and that of the residual lattice strain developed after failure may exist. The maximum tensile type residual lattice strain observed after failure may be resulted from the band structure formed in the hardest <111> grains. It was revealed that the band structure triggered by the hard particles could be one of sources of damage. In addition, while the specimen was obviously damaged, a fast stress relief was evidenced after unloading from the tension, especially at the beginning of unloading. Our present investigations provide a novel method for exploring the damage mechanisms of polycrystalline materials during plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2024

Thermal and electric field driven rf breakdown precursor formation on metal surfaces

The phenomenon of electric breakdown poses serious challenges to the design of devices that operate in high electric field environments. Experimental evidence points toward breakdown events that are accompanied by elevated temperatures and dark current spikes, which is attributed to high-asperity nanostructure formation that enhances the local electric field and triggers a runaway process. However, the exact mechanistic origin of such nanostructures under typical macroscopic operational conditions of electric field and magnetic-field-mediated heating remains poorly understood. In this work, we simulate the evolution of a copper surface under the combined action of the electric fields and elevated temperatures. Using a mesoscale curvature-driven surface evolution model, we show how a copper surface can undergo a type of dynamical instability that naturally leads to the formation of sharp asperities in realistic experimental conditions. Exploring the combined effect of fields and temperature rise, we identify the critical regimes that allow for the formation of breakdown precursors. The results show that thermoelastic stresses, while not essential, can significantly lower the critical electric field required for runaway surface instability, which is consistent with experimental observations that thermal effects can increase breakdown rates. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
NASA NTRS1988

Efficient assembly of finite-element subsystems with large relative rotations

A finite element approach is presented for the modeling of rotorcraft undergoing elastic deformation in addition to large rigid body motion with respect to inertial space, with particular attention given to the coupling of the rotor and fuselage subsystems subject to large relative rotations. The component synthesis technique used here allows the coupling of rotors to the fuselage for different rotorcraft configurations. The formulation is general and applicable to any rotorcraft vibration, aeroelasticity, and dynamics problem.

Fuh, Jon-Shen↗
NASA NTRS1989

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗
NASA NTRS1994

Transverse shear stresses and their sensitivity coefficients in multilayered composite panels

A computational procedure is presented for the accurate determination of transverse shear stresses and their sensitivity coefficients in flat multilayered composite panels subjected to mechanical and thermal loads. The sensitivity coefficients measure the sensitivity of the transverse shear stresses to variations in the different lamination and material parameters of the panel. The panel is discretized by using either a three-field mixed finite element model based on a two-dimensional first- order shear deformation plate theory or a two-field degenerate solid element with each of the displacement components having a linear variation throughout the thickness of the laminate. The evaluation of transverse shear stresses can be conveniently divided into two phases. The first phase consists of using a superconvergent recovery technique for evaluating the in-plane stresses in the different layers. In the second phase, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of the transverse shear stresses. The effectiveness of the computational procedure is demonstrated by means of numerical examples of multilayered cross-ply panels subjected to transverse loading, uniform temperature change, and uniform temperature gradient through the thickness of the panel. In each case the standard of the comparison is taken to be the exact solution of the three dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗
NASA NTRS1999

The Considere Condition and Rapid Stretching of Linear and Branched Polymer Melts

We analyze the onset of "necking" and subsequent filament failure during the transient uniaxial elongation of viscoelastic fluid samples in extensional rheometers. In the limit of rapid elongation (such that no molecular relaxation occurs), the external work applied is all stored elastically and the Considere criterion originally developed in solid mechanics can be used to quantitatively predict the critical Hencky strain to failure. By comparing the predictions of the Doi-Edwards model for linear homopolymer melts with those of the "Pom-Pom" model for prototypical branched melts we show that the critical strain to failure in rapid elongation of a rubbery material is intimately linked to the molecular topology of the chain, especially the degree of chain branching. The onset of necking instability is monotonically shifted to larger Hencky strains as the number of branches is increased. Numerical computations at finite Deborah numbers also show that there is an optimal range of deformation rates over which homogeneous extensions can be maintained to large strain. We also consider other rapid homogeneous stretching deformations, such as biaxial and planar stretching, and show that the degree of stabilization afforded by inclusion of material with long-chain branching is a sensitive function of the imposed mode of deformation.

McKinley, Gareth H.↗