SEARCH · Engineering Papers

Results for “rate dependent plasticity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1983

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.↗

NASA NTRS2026

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗

NASA NTRS1987

Effect of grain size on the high temperature properties of B2 aluminides

Measurements of the slow plastic flow behavior of cobalt, iron and nickel B2 crystal structure aluminides were conducted on materials fabricated by metallurical techniques. Due to this processing, the aluminides invariably had small equiaxed grains, ranging in size from about 3 to 60 microns in diameter. Grain size was dependent on the extrusion temperature used for powder consolidation, and it proved to be remarkably stable at elevated temperatures. Mechanical properties of all three aluminides were determined via constant velocity compression testing in air between 1000 and 1400 K at strain rates ranging from approx. 10 to the minus 3 power to 10 to the minus 7 power s (-1).

Whittenberger, J. Daniel↗

NASA NTRS1990

Finite element thermo-viscoplastic analysis of aerospace structures

The time-dependent thermo-viscoplastic response of aerospace structures subjected to intense aerothermal loads is predicted using the finite-element method. The finite-element analysis uses the Bodner-Partom unified viscoplastic constitutive relations to determine rate-dependent nonlinear material behavior. The methodology is verified by comparison with experimental data and other numerical results for a uniaxially-loaded bar. The method is then used (1) to predict the structural response of a rectangular plate subjected to line heating along a centerline, and (2) to predict the thermal-structural response of a convectively-cooled engine cowl leading edge subjected to aerodynamic shock-shock interference heating. Compared to linear elastic analysis, the viscoplastic analysis results in lower peak stresses and regions of plastic deformations.

Pandey, Ajay K.↗

NASA NTRS1990

Finite-element thermo-viscoplastic analysis of aerospace structures

Pandey, Ajay↗

NASA NTRS1985

A survey of unified constitutive theories

The state of the art of time temperature dependent elastic viscoplastic constitutive theories which are based on the unified approach werre assessed. This class of constitutive theories is characterized by the use of kinetic equations and internal variables with appropriate evolutionary equations for treating all aspects of inelastic deformation including plasticity, creep, and stress relaxation. More than 10 such unified theories which are shown to satisfy the uniqueness and stability criteria imposed by Drucker's postulate and Ponter's inequalities are identified. The theories are compared for the types of flow law, kinetic equation, evolutionary equation of the internal variables, and treatment of temperature dependence. The similarities and differences of these theories are outlined in terms of mathematical formulations and illustrated by comparisons of theoretical calculations with experimental results which include monotonic stress-strain curves, cyclic hysteresis loops, creep and stress relaxation rates, and thermomechanical loops. Numerical methods used for integrating these stiff time temperature dependent constitutive equations are reviewed.

Chan, K. S.↗

NASA NTRS1992

Unified Viscoplastic Behavior of Metal Matrix Composites

The need for unified constitutive models was recognized more than a decade ago in the results of phenomenological tests on monolithic metals that exhibited strong creep-plasticity interaction. Recently, metallic alloys have been combined to form high-temperature ductile/ductile composite materials, raising the natural question of whether these metallic composites exhibit the same phenomenological features as their monolithic constituents. This question is addressed in the context of a limited, yet definite (to illustrate creep/plasticity interaction) set of experimental data on the model metal matrix composite (MMC) system W/Kanthal. Furthermore, it is demonstrated that a unified viscoplastic representation, extended for unidirectional composites and correlated to W/Kanthal, can accurately predict the observed longitudinal composite creep/plasticity interaction response and strain rate dependency. Finally, the predicted influence of fiber orientation on the creep response of W/Kanthal is illustrated.

Arnold, S. M.↗

NASA NTRS1983

Mechanics of composite materials: Recent advances; Proceedings of the Symposium, Virginia Polytechnic Institute and State University, Blacksburg, VA, August 16-19, 1982

The present conference on the mechanics of composites discusses microstructure's influence on particulate and short fiber composites' thermoelastic and transport properties, the elastoplastic deformation of composites, constitutive equations for viscoplastic composites, the plasticity and fatigue of metal matrix composites, laminate damping mechanisms, the micromechanical modeling of Kevlar/epoxy composites' time-dependent failure, the variational characterization of waves in composites, and computational methods for eigenvalue problems in composite design. Also discussed are the elastic response of laminates, elastic coupling nonlinear effects in unsymmetrical laminates, elasticity solutions for laminate problems having stress singularities, the mechanics of bimodular composite structures, the optimization of laminated plates and shells, NDE for laminates, the role of matrix cracking in the continuum constitutive behavior of a damaged composite ply, and the energy release rates of various microcracks in short fiber composites.

Hashin, Z.↗

NASA NTRS2012

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗

NASA NTRS1987

Nonisothermal elasto-visco-plastic response of shell-type structures

A mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic or static thermomechanical loads is disussed. Among the system responses, which are associated with these load conditions, are thermal buckling and creep buckling. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and have been considered in the development of the mathematical model. Furthermore, this was accommodated in the solution procedures. A complete true ab-inito rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformation, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of the metric and curvature in two coordinate systems, a fixed (spatial) one and a convected (material) coordinate system. The relations between the time derivative and the covariant derivatives (gradients) have been developed for curved space and motion, so that the velocity components supply the connection between the equations of motion and the time rate of change of the metric and curvature tensors.

Simitses, G. J.↗

NASA NTRS1990

Experiment K-6-24, K-6-25, K-6-26. Radiation dosimetry and spectrometry

Radiation experiments flown by the University of San Francisco on the Cosmos 1887 spacecraft were designed to measure the depth dependence of both total dose and heavy particle flux, dose and dose equivalent, down to very thin shielding. Three experiments were flown and were located both inside and outside the Cosmos 1887 spacecraft. Tissue absorbed dose rates of 264 to 0.028 rad d(-1) under shielding of 0.013 to 3.4 g/sq cm of (7)LiF were found outside the spacecraft and 0.025 rad d(-1) inside. Heavy particle fluxes of 3.43 to 1.03 x 10 to the minus 3rd power cm -2 sub s -1 sub sr -1 under shielding of 0.195 to 1.33 g/sq cm plastic were found outside the spacecraft and 4.25 times 10 to the minus 4th power cm -2 sub s -1 sub sr -1 inside (LET infinity H2O greater than or equal to 4 keV/micron m). The corresponding heavy particle dose equivalent rates outside the spacecraft were 30.8 to 19.8 mrem d(-1) and 11.4 mrem d(-1) inside. The large dose and particle fluxes found at small shielding thicknesses emphasize the importance of these and future measurements at low shielding, for predicting radiation effects on space materials and experiments where shielding is minimal and on astronauts during EVA. The Cosmos 1887 mission contained a variety of international radiobiological investigations to which the measurements apply. The high inclination orbit (62 degrees) of this mission provided a radiation environment which is seldom available to U.S. investigators. The radiation measurements will be compared with those of other research groups and also with those performed on the Shuttle, and will be used to refine computer models employed to calculate radiation exposures on other spacecraft, including the Space Station.

Benton, E. V.↗

NASA NTRS1983

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed.

Reed, K. W.↗

NASA NTRS1982

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.↗

NASA NTRS2004

Slip Analysis in a Ni-base Superalloy

A Ni-base superalloy single crystal with Gamma/Gamma' structure was tested at room temperature along the , <110> and <111> directions. Consistent with previously reported investigations, this alloy did not obey the Schmid law and the CRSS (critical resolved shear stress) was noticeably lower for the <111>-oriented samples. Furthermore, the strain hardening rate decreased and the degree of deformation localization increased in the order of <111>, and <110> orientations. The appearance and orientation of deformation traces were found to depend on the loading orientation as well as the amount of strain. In general, when Gamma'-particles were sheared, the traces followed the expected octahedral shear planes. It is demonstrated that the wavy deformation traces that do not follow the {111} planes are associated with changes in the gamma-channels width and the falling off of the gamma-particles. In this paper the evolution of deformation bands are discussed in terms of plastic localization at microscopic, mesoscopic and macroscopic levels.

Westbrooke, Eboni F.↗

NASA NTRS1995

CO2 (dry ice) cleaning system

Tomco Equipment Company has participated in the dry ice (solid carbon dioxide, CO2) cleaning industry for over ten years as a pioneer in the manufacturer of high density, dry ice cleaning pellet production equipment. For over four years Tomco high density pelletizers have been available to the dry ice cleaning industry. Approximately one year ago Tomco introduced the DI-250, a new dry ice blast unit making Tomco a single source supplier for sublimable media, particle blast, cleaning systems. This new blast unit is an all pneumatic, single discharge hose device. It meters the insertion of 1/8 inch diameter (or smaller), high density, dry ice pellets into a high pressure, propellant gas stream. The dry ice and propellant streams are controlled and mixed from the blast cabinet. From there the mixture is transported to the nozzle where the pellets are accelerated to an appropriate blasting velocity. When directed to impact upon a target area, these dry ice pellets have sufficient energy to effectively remove most surface coatings through dry, abrasive contact. The meta-stable, dry ice pellets used for CO2 cleaning, while labeled 'high density,' are less dense than alternate, abrasive, particle blast media. In addition, after contacting the target surface, they return to their equilibrium condition: a superheated gas state. Most currently used grit blasting media are silicon dioxide based, which possess a sharp tetrahedral molecular structure. Silicon dioxide crystal structures will always produce smaller sharp-edged replicas of the original crystal upon fracture. Larger, softer dry ice pellets do not share the same sharp-edged crystalline structures as their non-sublimable counterparts when broken. In fact, upon contact with the target surface, dry ice pellets will plastically deform and break apart. As such, dry ice cleaning is less harmful to sensitive substrates, workers and the environment than chemical or abrasive cleaning systems. Dry ice cleaning system components include: a dry ice pellet supply, a non-reactive propellant gas source, a pellet and propellant metering device, and a media transport and acceleration hose and nozzle arrangement. Dry ice cleaning system operating parameters include: choice of propellant gas, its pressure and temperature, dry ice mass flow rate, dry ice pellet size and shape, and acceleration nozzle configuration. These parameters may be modified to fit different applications. The growth of the dry ice cleaning industry will depend upon timely data acquisition of the effects that independent changes in these parameters have on cleaning rates, with respect to different surface coating and substrate combinations. With this data, optimization of cleaning rates for particular applications will be possible. The analysis of the applicable range of modulation of these parameters, within system component mechanical constraints, has just begun.

Barnett, Donald M.↗

NASA NTRS2001

Viscoplastic Model Development to Account for Strength Differential: Application to Aged Inconel 718 at Elevated Temperature

The magnitude of yield and flow stresses in aged Inconel 718 are observed to be different in tension and compression. This phenomenon, called the Strength differential (SD), contradicts the metal plasticity axiom that the second deviatoric stress invariant alone is sufficient for representing yield and flow. Apparently, at least one of the other two stress invariants is also significant. A unified viscoplastic model was developed that is able to account for the SD effect in aged Inconel 718. Building this model involved both theory and experiments. First, a general threshold function was proposed that depends on all three stress invariants and then the flow and evolution laws were developed using a potential-based thermodynamic framework. Judiciously chosen shear and axial tests were conducted to characterize the material. Shear tests involved monotonic loading, relaxation, and creep tests with different loading rates and load levels. The axial tests were tension and compression tests that resulted in sufficiently large inelastic strains. All tests were performed at 650 C. The viscoplastic material parameters were determined by optimizing the fit to the shear tests, during which the first and the third stress invariants remained zero. The threshold surface parameters were then fit to the tension and compression test data. An experimental procedure was established to quantify the effect of each stress invariant on inelastic deformation. This requires conducting tests with nonproportional three-dimensional load paths. Validation of the model was done using biaxial tests on tubular specimens of aged Inconel 718 using proportional and nonproportional axial-torsion loading. These biaxial tests also helped to determine the most appropriate form of the threshold function; that is, how to combine the stress invariants. Of the set of trial threshold functions, the ones that incorporated the third stress invariant give the best predictions. However, inclusion of the first stress invariant does not significantly improve the model predictions. The model shows excellent predictive capability for nonproportional load paths. Additionally, it reduces to the well-known models of Mises Drucker and Drucker-Prager. The requisite experiments involve reasonably simple load paths in the axial-shear stress plane and hence can be performed on a variety of different materials: be they metallic, geological. polymeric, ceramic or granular. The general form of the threshold function allows representation of inelastic deformation in a range of materials.

Iyer, Saiganesh↗

NASA NTRS1972

Slow and fast motion of cracks in inelastic solids. Part 1: Slow growth of cracks in a rate sensitive tresca solid. Part 2: Dynamic crack represented by the Dugdale model

An extension is proposed of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equals the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus reflects the extent of inelastic behavior of a solid. A concept of apparent surface energy is introduced to account for the geometry dependent and the rate dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared. These are: (1) constant fracture energy (Cherepanov), (2) constant opening displacement at instability (Morozov) and (3) final stretch criterion (Wnuk).

Wnuk, M. P.↗

NASA NTRS2011

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗