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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

A General Reversible Hereditary Constitutive Model: Theoretical Developments - Part 1

Using an internal-variable formalism as a starting point, we describe the viscoelastic extension of a previously-developed viscoplasticity formulation of the complete potential structure type. It is mainly motivated by experimental evidence for the presence of rate/time effects in the so-called quasilinear, reversible, material response range. Several possible generalizations are described, in the general format of hereditary-integral representations for non-equilibrium, stress-type, state variables, both for isotropic as well as anisotropic materials. In particular, thorough discussions are given on the important issues of thermodynamic admissibility requirements for such general descriptions, resulting in a set of explicit mathematical constraints on the associated kernel (relaxation and creep compliance) functions. In addition, a number of explicit, integrated forms are derived, under stress and strain control to facilitate the parametric and qualitative response characteristic studies reported here, as well as to help identify critical factors in the actual experimental characterizations from test data that will be reported in Part II.

Saleeb, A. F.↗
NASA NTRS1998

High Strain Rate Deformation Modeling of a Polymer Matrix Composite: Matrix Constitutive Equations - Part 1

Recently applications have exposed polymer matrix composite materials to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under these extreme conditions. In this first paper of a two part report, background information is presented, along with the constitutive equations which will be used to model the rate dependent nonlinear deformation response of the polymer matrix. Strain rate dependent inelastic constitutive models which were originally developed to model the viscoplastic deformation of metals have been adapted to model the nonlinear viscoelastic deformation of polymers. The modified equations were correlated by analyzing the tensile/ compressive response of both 977-2 toughened epoxy matrix and PEEK thermoplastic matrix over a variety of strain rates. For the cases examined, the modified constitutive equations appear to do an adequate job of modeling the polymer deformation response. A second follow-up paper will describe the implementation of the polymer deformation model into a composite micromechanical model, to allow for the modeling of the nonlinear, rate dependent deformation response of polymer matrix composites.

Goldberg, Robert K.↗
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS1999

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗
NASA NTRS2000

A New Local Debonding Model with Application to the Transverse Tensile and Creep Behavior of Continuously Reinforced Titanium Composites

A new, widely applicable model for local interfacial debonding in composite materials is presented. Unlike its direct predecessors, the new model allows debonding to progress via unloading of interfacial stresses even as global loading of the composite continues. Previous debonding models employed for analysis of titanium matrix composites are surpassed by the accuracy, simplicity, and efficiency demonstrated by the new model. The new model was designed to operate seamlessly within NASA Glenn's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), which was employed to simulate the time- and rate-dependent (viscoplastic) transverse tensile and creep behavior of SiC/Ti composites. MAC/GMC's ability to simulate the transverse behavior of titanium matrix composites has been significantly improved by the new debonding model. Further, results indicate the need for a more accurate constitutive representation of the titanium matrix behavior in order to enable predictions of the composite transverse response, without resorting to recalibration of the debonding model parameters.

Bednarcyk, Brett A.↗
NASA NTRS2000

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗
NASA NTRS2001

The Effect of Heat Treatment on Residual Stress and Machining Distortions in Advanced Nickel Base Disk Alloys

This paper describes an extension of NASA's AST and IDPAT Programs which sought to predict the effect of stabilization heat treatments on residual stress and subsequent machining distortions in the advanced disk alloy, ME-209. Simple "pancake" forgings of ME-209 were produced and given four heat treats: 2075F(SUBSOLVUS)/OIL QUENCH/NO AGE; 2075F/OIL QUENCH/1400F@8HR;2075F/OIL QUENCH/1550F@3HR/l400F@8HR; and 2160F(SUPERSOLVUS)/OIL QUENCH/1550F@3HR/ 1400F@8HR. The forgings were then measured to obtain surface profiles in the heat treated condition. A simple machining plan consisting of face cuts from the top surface followed by measurements of the surface profile opposite the cut were made. This data provided warpage maps which were compared with analytical results. The analysis followed the IDPAT methodology and utilized a 2-D axisymmetric, viscoplastic FEA code. The analytical results accurately tracked the experimental data for each of the four heat treatments. The 1550F stabilization heat treatment was found to significantly reduce residual stresses and subsequent machining distortions for fine grain (subsolvus) ME209, while coarse grain (supersolvus) ME209 would require additional time or higher stabilization temperatures to attain the same degree of stress relief.

Gayda, John↗
NASA NTRS2001

HITEMP Material and Structural Optimization Technology Transfer

The feasibility of adding viscoelasticity and the Generalized Method of Cells (GMC) for micromechanical viscoelastic behavior into the commercial HyperSizer structural analysis and optimization code was investigated. The viscoelasticity methodology was developed in four steps. First, a simplified algorithm was devised to test the iterative time stepping method for simple one-dimensional multiple ply structures. Second, GMC code was made into a callable subroutine and incorporated into the one-dimensional code to test the accuracy and usability of the code. Third, the viscoelastic time-stepping and iterative scheme was incorporated into HyperSizer for homogeneous, isotropic viscoelastic materials. Finally, the GMC was included in a version of HyperSizer. MS Windows executable files implementing each of these steps is delivered with this report, as well as source code. The findings of this research are that both viscoelasticity and GMC are feasible and valuable additions to HyperSizer and that the door is open for more advanced nonlinear capability, such as viscoplasticity.

Collier, Craig S.↗
NASA NTRS2001

Behavior of Three Metallic Alloys Under Combined Axial-Shear Stress at 650 C

Three materials, Inconel 718, Haynes 188, and 316 stainless steel, were tested under an axial-torsional stress state at 650 C. The objective of this study was to quantify the evolution of the material while in the viscoplastic domain. Initial and subsequent yield surfaces were experimentally determined to quantify hardening. Subsequent yield surfaces (yield surfaces taken after a preload) had a well-defined front side, in the prestrain direction, but a poorly defined back side, opposite the prestrain direction. Subsequent yield surfaces exhibited isotropic hardening by expansion of the yield surface, kinematic hardening by translation of the yield surface, and distortional hardening by flattening of the yield surface in the direction opposite to the last prestrain. An existing yield function capable of representing isotropic, kinematic, and distortional hardening was used to fit each yield surface. Four variables are used to describe each surface. These variables evolve as the material state changes and have been regressed to the yield surface data.

Colaiuta, Jason F.↗
NASA NTRS2002

Computational Simulation of the High Strain Rate Tensile Response of Polymer Matrix Composites

A research program is underway to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Under these types of loading conditions, the material response can be highly strain rate dependent and nonlinear. State variable constitutive equations based on a viscoplasticity approach have been developed to model the deformation of the polymer matrix. The constitutive equations are then combined with a mechanics of materials based micromechanics model which utilizes fiber substructuring to predict the effective mechanical and thermal response of the composite. To verify the analytical model, tensile stress-strain curves are predicted for a representative composite over strain rates ranging from around 1 x 10(exp -5)/sec to approximately 400/sec. The analytical predictions compare favorably to experimentally obtained values both qualitatively and quantitatively. Effective elastic and thermal constants are predicted for another composite, and compared to finite element results.

Goldberg, Robert K.↗
NASA NTRS2003

Behavior of Three Metallic Alloys under Combined Axial-Shear Stresses at Elevated Temperature

Type 316 stainless steel, Haynes 188, and Inconel 718 samples were subjected to an axial-shear strain controlled loading history while the specimen temperature was held at 650 C to quantify the evolution of material state under a complex biaxial load path when the material is in the viscoplastic domain. Yield surfaces were constructed in the axial-shear stress plane using a sensitive, 30 x 10(exp -6)m/m, equivalent offset strain definition for the yield strain. Subsequent yield surfaces were constructed at various points along the strain path to define the material evolution. These subsequent yield surface translated, expanded, and distorted relative to the initial yield surface. Each of these very different materials exhibited components of isotropic, kinematic and distortional hardening. Furthermore, subsequent yield surfaces for each material have a very well defined front face and a poorly defined, flattened, back side.

Colaiuta, J. F.↗
NASA NTRS2003

Specific Hardening Function Definition and Characterization of a Multimechanism Generalized Potential-based Viscoelastoplasticity Model

Given the previous complete-potential structure framework together with the notion of strain- and stress-partitioning in terms of separate contributions of several submechanisms (viscoelastic and viscoplastic) to the thermodynamic functions (stored energy and dissipation) a detailed viscoelastoplastic multimechanism characterization of a specific hardening functional form of the model is presented and discussed. TIMETAL 21S is the material of choice as a comprehensive test matrix, including creep, relaxation, constant strain-rate tension tests, etc. are available at various temperatures. Discussion of these correlations tests, together with comparisons to several other experimental results, are given to assess the performance and predictive capabilities of the present model particularly with regard to the notion of hardening saturation as well as the interaction of multiplicity of dissipative (reversible/irreversible) mechanisms.

Arnold, S. M.↗
NASA NTRS2003

Microsample Characterization of Coatings for Grcop-84 for High Heat Flux Applications

A multidisciplinary Johns Hopkins University-NASA Glenn team is undertaking a collaborative research program to elucidate and model the thermal stability and mechanical integrity of candidate coatings for GRCop-84. GRCop-84 is a high conductivity, high strength copper alloy that was recently developed at NASA Glenn for use in high temperature, high heat flux applications. With potential applications in rocket motor combustion chamber liners, nozzle ramps and other actively cooled structures, this new material offers great potential for decreasing weight and increasing reliability of third generation reusable launch vehicles. Current emphasis has turned toward the development of environmentally resistant and thermal barrier coatings for this alloy. Metallic coatings such as NiCrAlY and Cu-8-30%Cr have shown promise in: prohibiting blanching, reducing dog-house failures, increasing operating temperatures and decreasing cooling requirements. The focus of this research program is to develop a fundamental understanding of the substrate-coating interactions that occur during thermal cycling (inter-diffusion, viscoplasticity, morphological evolution, crack formation, etc.) and to derive a science-based protocol for future coating selection, optimization and reliability assurance. The microsample tensile testing approach adopted for this study allows us to characterize small-scale and highly scale-specific coatings and properties in a way not possible by conventional means. In addition to providing much needed design data, the integration of microsample testing with detailed microstructural observations provides a mechanistic foundation for coating optimization and life prediction modeling.

Hemker, Kevin↗
NASA NTRS2004

Analysis of Factors Affecting the Performance of RLV Thrust Cell Liners

The reusable launch vehicle (RLV) thrust cell liner, or thrust chamber, is a critical component of the Space Shuttle Main Engine (SSME). It is designed to operate in some of the most severe conditions seen in engineering practice. This requirement, in conjunction with experimentally observed 'dog-house' failure modes characterized by bulging and thinning of the cooling channel wall, provides the motivation to study the factors that influence RLV thrust cell liner performance. Factors or parameters believed to be directly related to the observed characteristic deformation modes leading to failure under in-service loading conditions are identified, and subsequently investigated using the cylindrical version of the higher-order theory for functionally graded materials in conjunction with the Robinson's unified viscoplasticity theory and the power-law creep model for modeling the response of the liner s constituents. Configurations are analyzed in which specific modifications in cooling channel wall thickness or constituent materials are made to determine the influence of these parameters on the deformations resulting in the observed failure modes in the outer walls of the cooling channel. The application of thermal barrier coatings and functional grading are also investigated within this context. Comparison of the higher-order theory results based on the Robinson and power-law creep model predictions has demonstrated that, using the available material parameters, the power-law creep model predicts more precisely the experimentally observed deformation leading to the 'dog-house' failure mode for multiple short cycles, while also providing much improved computational efficiency. However, for a single long cycle, both models predict virtually identical deformations. Increasing the power-law creep model coefficients produces appreciable deformations after just one long cycle that would normally be obtained after multiple cycles, thereby enhancing the efficiency of the analysis. This provides a basis for the development of an accelerated modeling procedure to further characterize dog-house deformation modes in RLV thrust cell liners. Additionally, the results presented herein have demonstrated that the mechanism responsible for deformation leading to 'dog-house' failure modes is driven by pressure, creep/relaxation and geometric effects.

Arnold, Steven M.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗
NASA NTRS1998

Experimental Techniques Verified for Determining Yield and Flow Surfaces

Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.

Lerch, Brad A.↗
NASA NTRS2000

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗