Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscoplasticity”

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 271 records · Page 15

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗

Review: Inelastic Constitutive Modeling: Polycrystalline Materials

This article provides a literature review that details the development of inelastic constitutive modeling as it relates to polycrystalline materials. This review distinguishes between inelastic constitutive models that account for nonlinear behavior at the microstructural level, time-independent classic plasticity models, and time-dependent unified models. Particular emphasis is placed on understanding the underlying theoretical framework for unified viscoplasticity models where creep and classical plasticity behavior are considered the result of applied boundary conditions instead of separable rates representing distinct physical mechanisms. This article establishes a clear understanding of the advantages of the unified approach to improve material modeling. This review also discusses recent topics in constitutive modeling that offer new techniques that bridge the gap between the microstructure and the continuum.

36 MATERIALS SCIENCE↗

Deformation and Transformation Textures in the NaMgF3 Neighborite—Post-Perovskite System

The D″ region of the lower mantle, which lies just above the core–mantle boundary, is distinct from the bulk of the lower mantle in that it exhibits complex seismic heterogeneity and seismic anisotropy. Seismic anisotropy in this region is likely to be largely due to the deformation-induced texture (crystallographic preferred orientation) development of the constituent mineral phases. Thus, seismic anisotropy can provide a marker for deformation processes occurring in this dynamic region of the Earth. Post-perovskite-structured (Mg,Fe)SiO3 is believed to be the dominant mineral phase in many regions of the D”. As such, understanding deformation mechanisms and texture development in post-perovskite is important for the interpretation of observed seismic anisotropy. Here, we report on high-pressure diamond anvil cell deformation experiments on NaMgF3 neighborite (perovskite structure) and post-perovskite. During deformation, neighborite develops a 100 texture, as has been previously observed, both in NaMgF3 and MgSiO3 perovskite. Upon transformation to the post-perovskite phase, an initial texture of {130} at high angles to compression is observed, indicating that the {100} planes of perovskite become the ~{130} planes of post-perovskite. Further compression results in the development of a shoulder towards (001) in the inverse pole figure. Plasticity modeling using the elasto-viscoplastic self-consistent code shows this texture evolution to be most consistent with deformation on (001)[100] with some contribution of glide on (100)[010] and (001)<110> in NaMgF3 post-perovskite. The transformation and deformation mechanisms observed in this study in the NaMgF3 system are consistent with the behavior generally observed in other perovskite–post-perovskite systems, including the MgSiO3 system. This shows that NaMgF3 is a good analog for the mantle bridgmanite and MgSiO3 post-perovskite.

Geochemistry & Geophysics↗

Cyclic-stress analysis of notches for supersonic transport conditions

The feasibility of using the finite element method to account for the effects of cyclic load and temperature on local stresses and strains at a notch was demonstrated. The behavior of a notched titanium panel was studied under variable loads and temperatures representative of flight conditions for the lower wing surface of a Supersonic Transport (SST). The analysis was performed with the use of the BOPACE finite-element computer program which provides capability to determine high temperature and large viscoplastic effects caused by cyclic thermal and mechanical loads. The analysis involves the development of the finite-element model as well as determination of the structural behavior of the notched panel. Results are presented for twelve SST flights comprised of five different load-temperature cycles. The results show the approach is feasible, but material response to cyclic loads, temperatures, and hold times requires improved understanding to allow proper modeling of the material.

Source record↗

The 3.3K thrust chamber life prediction

The cause of low-cycle fatigue failure of the 3.3K Thrust Chamber was investigated. This thrust chamber typifies the current trend in rocket nozzle design which calls for high performance coupled with weight and volume limitations as well as the requirement of reusability. The analysis was performed with the BOPACE finite-element computer program which provides capability to determine viscoplastic response of a structure subjected to cyclic thermal and mechanical loading. Results are presented which show the critical region for low-cycle fatigue and the history of strain within that region for each thermo-mechanical loading cycle in the 3.3K thrust chamber. The predicted behavior was used to evaluate the low-cycle fatigue life near the throat plane of the chamber. The results show that BOPACE provides an extremely accurate prediction of structural behavior; the critical region was identified and the life determined from computed strains was within 154 cycles of the observed failure at 1013 cycles.

Armstrong, W. H.↗

Surface temperatures and glassy state investigations in tribology, part 1

The research in this report is divided into two categories: (1) lubricant rheological behavior, and (2) thermal behavior of a simulated elastohydrodynamic contact. The studies of the lubricant rheological behavior consists of high pressure, low shear rate viscosity measurements, viscoelastic transition measurements, by volume dilatometry, dielectric transitions at atmospheric pressure and light scattering transitions. Lubricant shear stress-strain behavior in the amorphous glassy state was measured on several fluids. It appears clear from these investigations that many lubricants undergo viscoplastic transitions in typical EHD contacts and that the lubricant has a limiting maximum shear stress it can support which in turn will determine the traction in the contact except in cases of very low slide-roll ratio. Surface temperature measurements were made for a naphthenic mineral oil and a polyphenyl ether. The maximum surface temperature in these experiments was approximately symmetrical about the zero slide-roll ration except for absolute values of slide-roll ratio greater than about 0.9. Additional surface temperature measurements were made in contacts with rough surfaces where the composite surface roughness was approximately equal to the EHD film thickness. A regression analysis was done to obtain a predictive equation for surface temperatures as a function of pressure, sliding speed, and surface roughness. A correction factor for surface roughness effects to the typical flash temperature analysis was found.

Winer, W. O.↗

Effect of misalignment on mechanical behavior of metals in creep

Application of the endochronic theory of viscoplasticity to creep, creep recovery, and stress relaxation at the small strain and short time range produced the following results: (1) The governing constitutive equations for constant-strain-rate stress-strain behavior, creep, creep recovery, and stress relaxation were derived by imposing appropriate constraints on the general constitutive equation of the endochronic theory. (2) A set of material constants was found which correlate strain-hardening, creep, creep recovery, and stress relaxation. (3) The theory predicts with reasonable accuracy the creep and creep recovery behaviors at short time. (4) The initial strain history prior to the creep stage affects the subsequent creep significantly. (5) A critical stress was established for creep recovery. A computer program, written for the misalignment problem is reported.

Wu, H. C.↗

Endochronic theory of transient creep and creep recovery

Short time creep and creep recovery were investigated by means of the endochronic theory of viscoplasticity. It is shown that the constitutive equations for constant-strain-rate stress-strain behavior, creep, creep recovery, and stress relaxation can all ber derived from the general constitutive equation by imposing appropriate constraints. In this unified approach, the effect of strain-hardening is naturally accounted for when describing creep and creep recovery. The theory predicts with reasonable accuracy the creep and creep recovery behaviors for Aluminum 1100-0 at 150 C. It was found that the strain-rate history at prestraining stage affects the subsequent creep. A critical stress was also established for creep recovery. The theory predicts a forward creep for creep recovery stress greater than the critical stress. For creep recovery stress less than the critical stress, the theory then predicts a normal strain recovery.

Wu, H. C.↗

CELFE: Coupled Eulerian-Lagrangian Finite Element program for high velocity impact. Part 1: Theory and formulation

A 3-D finite element program capable of simulating the dynamic behavior in the vicinity of the impact point, together with predicting the dynamic response in the remaining part of the structural component subjected to high velocity impact is discussed. The finite algorithm is formulated in a general moving coordinate system. In the vicinity of the impact point contained by a moving failure front, the relative velocity of the coordinate system will approach the material particle velocity. The dynamic behavior inside the region is described by Eulerian formulation based on a hydroelasto-viscoplastic model. The failure front which can be regarded as the boundary of the impact zone is described by a transition layer. The layer changes the representation from the Eulerian mode to the Lagrangian mode outside the failure front by varying the relative velocity of the coordinate system to zero. The dynamic response in the remaining part of the structure described by the Lagrangian formulation is treated using advanced structural analysis. An interfacing algorithm for coupling CELFE with NASTRAN is constructed to provide computational capabilities for large structures.

Lee, C. H.↗

Shock Sensitivity of energetic materials

Viscoplastic deformation is examined as the principal source of hot energy. Some shock sensitivity data on a proposed model is explained. A hollow sphere model is used to approximate complex porous matrix of energetic materials. Two pieces of shock sensitivity data are qualitatively compared with results of the proposed model. The first is the p2 tau law. The second is the desensitization of energetic materials by a ramp wave applied stress. An approach to improve the model based on experimental observations is outlined.

Kim, K.↗

Finite-element nonlinear transient response computer programs PLATE 1 and CIVM-PLATE 1 for the analysis of panels subjected to impulse or impact loads

Two computer programs are described for predicting the transient large deflection elastic viscoplastic responses of thin single layer, initially flat unstiffened or integrally stiffened, Kirchhoff-Lov ductile metal panels. The PLATE 1 program pertains to structural responses produced by prescribed externally applied transient loading or prescribed initial velocity distributions. The collision imparted velocity method PLATE 1 program concerns structural responses produced by impact of an idealized nondeformable fragment. Finite elements are used to represent the structure in both programs. Strain hardening and strain rate effects of initially isotropic material are considered.

Spilker, R. L.↗

Effects of misalignment on mechanical behavior of metals in creep

Creep tests were conducted by means of a closed loop servocontrolled materials test system. The strain history prior to creep is carefully monitored. Tests were performed for aluminum alloy 6061-O at 150 C and were monitored by a PDP 11/04 minicomputer at a preset constant plastic strain rate prehistory. The results show that the plastic strain rate prior to creep plays a significant role in creep behavior. The endochronic theory of viscoplasticity was applied to describe the observed creep curves. Intrinsic time and strain rate sensitivity function concepts are employed and modified according to the present observation.

Wu, H. C.↗

Investigation of creep by use of closed loop servo-hydraulic test system

Creep tests were conducted by means of a closed loop servo-controlled materials test system. These tests are different from the conventional creep tests in that the strain history prior to creep may be carefully monitored. Tests were performed for aluminum alloy 6061-0 at 150 C and monitored by a PDP 11/04 minicomputer at a preset constant plastic-strain rate prehistory. The results show that the plastic-strain rate prior to creep plays a significant role in creep behavior. The endochronic theory of viscoplasticity was applied to describe the observed creep curves. The concepts of intrinsic time and strain rate sensitivity function are employed and modified according to the present observation.

Wu, H. C.↗

Operator split methods in the numerical solution of the finite deformation elastoplastic dynamic problem

The spatial formulation of the elastoplastic dynamic problem for finite deformations is considered. A thermodynamic argument leads to an additive decomposition of the spatial rate of deformation tensor and allows an operator split of the evolutionary equations of the problem into elastic and plastic parts. This operator split is taken as the basis for the definition of a global product algorithm. In the context of finite element discretization the product algorithm entails, for every time step, the solution of a nonlinear elastodynamic problem followed by the application of plastic algorithms that operate on the stresses and internal variables at the integration points and bring in the plastic constitutive equations. Suitable plastic algorithms are discussed for the cases of perfect and hardening plasticity and viscoplasticity. The proposed formalism does not depend on any notion of smoothness of the yield surface and is applicable to arbitrary convex elastic regions, with or without corners. The stabiity properties of the global product algorithm are shown to be identical to those of the algorithm used for the integration of the nonlinear elastodynamic problem. Numerical examples illustrate the accuracy of the method.

Pinsky, P. M.↗

On a study of the /Delta T/c and C/asterisk/ integrals for fracture analysis under non-steady creep

Applications of a vector quantity, path-independent integral which has an energy interpretation to the characterization of crack-tip fields in the range from fast to slow crack propagation are examined. The crack tip characterization parameter is defined in terms of a conservation integral for an area around the crack tip in a two-dimensional cracked body. The actual physical interpretation of the parameter is shown to be the difference in crack lengths displayed by two identical bodies which have equal load histories. A steady-state value is obtained for the parameter for cases of steady-state creep and is shown to be related to the standard path-independent integral for macroscopic self-similar crack growth under mode I conditions. A finite element model is developed for viscoplastic material models, using an initial strain approach with steps in a size employed in tangent stiffness methods.

Stonesifer, R. B.↗

Evaluation of Three Constitutive Models for the Prediction of Hastelloy X Elevated Temperature Cyclic Response

The approach for the method development assumes that, for a thermally loaded structure, the overall strain history is defined by linear elastic analysis. The local stress history at a fatigue critical location is then determined from a one dimensional material behavior model and the local strain and temperature conditions. Three material models are currently being evaluated to assess their ability to predict relevant high temperature cyclic material response characteristics. They are: a time independent classical plasticity and creep representation, a time dependent viscoplastic model capable of predicting combined creep and plasticity effects, and an approximate elastic analysis approach that uses a series of stress-strain curves and a cyclic hardening model to determine reverse plasticity.

Moreno, V.↗

Research and development program for the development of advanced time-temperature dependent constitutive relationships. Volume 1: Theoretical discussion

The results of a 10-month research and development program for the development of advanced time-temperature constitutive relationships are presented. The program included (1) the effect of rate of change of temperature, (2) the development of a term to include time independent effects, and (3) improvements in computational efficiency. It was shown that rate of change of temperature could have a substantial effect on the predicted material response. A modification to include time-independent effects, applicable to many viscoplastic constitutive theories, was shown to reduce to classical plasticity. The computation time can be reduced by a factor of two if self-adaptive integration is used when compared to an integration using ordinary forward differences. During the course of the investigation, it was demonstrated that the most important single factor affecting the theoretical accuracy was the choice of material parameters.

Cassenti, B. N.↗