ALBERT

Description: Extensive research efforts have been made over the years on the phenomenological representations of constitutive material behavior in the inelastic analysis and life assessment of structures composed of advanced monolithic and composite (CMC, MMC, and PMC) materials. Recently, emphasis has been placed on concurrently addressing three important and related areas of constitutive and degradation modeling; i.e. (i) mathematical formulation, (ii) algorithmic developments for the updating (integrating) of external (e.g. stress) and internal state variable, as well as (iii) parameter estimation for the characterization of the specific model. This concurrent perspective has resulted in; i) the formulation of a fully-associative viscoelastoplastic model (GVIPS), (ii) development of an efficient implicit integration and it's associative, symmetric, consistent tangent stiffness matrix algorithm for integration of the underlying rate flow/evolutionary equations, and iii) a robust, stand-alone, Constitutive Material Parameter Estimator (COMPARE) for automatically characterizing the various time-dependent, nonlinear, material models. Furthermore, to provide a robust multi-scale framework for the deformation and life analysis of structures composed of composite materials, NASA Glenn has aggressively pursued the development of a sufficiently general, accurate, and efficient micromechanics approach known as the generalized method of cells (GMC). This work has resulted in the development of MAC/GMC, a stand-alone micromechanics analysis tool that can easily and accurately design/analyze multiphase (composite) materials subjected to complex histories. MAC/GMC admits generalized, physically based, deformation and damage models for each constituent and provides "closed-form" expressions for the macroscopic composite response in terms of the properties, size, shape, distribution, and response of the individual constituents or phases that comprise the material. Consequently, MAC/GMC can be incorporated directly into a structural finite element code like ABAQUS for cost-effective, micromechanics based, large-scale component design and analysis. Our primary objective here is to report on these recent works conducted over the past decade, in the context of their incorporation into ABAQUS through the various user subroutines. Representative results will be shown to demonstrate the features of the developed schemes.

Description: A global/local time incrementing scheme for viscoplastic analysis of structures is presented. The scheme is very efficient and useful for conducting large scale nonlinear finite element analysis involving viscoplastic materials.

Description: Given the mathematical framework and specific viscoelastic model in Part 1 our primary goal in this second part is focused on model characterization and assessment for the specific titanium alloy, TIMETAL 21S. The model is motivated by experimental evidence suggesting the presence of significant rate/time effects in the so-called quasilinear, reversible, material response range. An explanation of the various experiments performed and their corresponding results are also included. Finally, model correlations and predictions are presented for a wide temperature range.

Description: 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.

Description: Since most advanced material systems (for example metallic-, polymer-, and ceramic-based systems) being currently researched and evaluated are for high-temperature airframe and propulsion system applications, the required constitutive models must account for both reversible and irreversible time-dependent deformations. Furthermore, since an integral part of continuum-based computational methodologies (be they microscale- or macroscale-based) is an accurate and computationally efficient constitutive model to describe the deformation behavior of the materials of interest, extensive research efforts have been made over the years on the phenomenological representations of constitutive material behavior in the inelastic analysis of structures. From a more recent and comprehensive perspective, the NASA Glenn Research Center in conjunction with the University of Akron has emphasized concurrently addressing three important and related areas: that is, 1) Mathematical formulation; 2) Algorithmic developments for updating (integrating) the external (e.g., stress) and internal state variables; 3) Parameter estimation for characterizing the model. This concurrent perspective to constitutive modeling has enabled the overcoming of the two major obstacles to fully utilizing these sophisticated time-dependent (hereditary) constitutive models in practical engineering analysis. These obstacles are: 1) Lack of efficient and robust integration algorithms; 2) Difficulties associated with characterizing the large number of required material parameters, particularly when many of these parameters lack obvious or direct physical interpretations.

Description: The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Description: The issue of developing effective and robust schemes to implement general hyperelastic constitutive models is addressed. To this end, special purpose functions are used to symbolically derive, evaluate, and automatically generate the associated FORTRAN code for the explicit forms of the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid for the entire deformation range. The analytical form of these explicit expressions is given here for the case in which the strain-energy potential is taken as a nonseparable polynomial function of the principle stretches.

Description: The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic constitutive models is addressed. To this end, special purpose functions (running under MACSYMA) are developed for the symbolic derivation, evaluation, and automatic FORTRAN code generation of explicit expressions for the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid over the entire deformation range, since the singularities resulting from repeated principal-stretch values have been theoretically removed. The required computational algorithms are outlined, and the resulting FORTRAN computer code is presented.

Description: The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic constitutive models is addressed. To this end, explicit forms for the corresponding material tangent stiffness tensors are developed, and these are valid for the entire deformation range; i.e., with both distinct as well as repeated principal-stretch values. Throughout the analysis the various implications of the underlying property of separability of the strain-energy functions are exploited, thus leading to compact final forms of the tensor expressions. In particular, this facilitated the treatment of complex cases of uncoupled volumetric/deviatoric formulations for incompressible materials. The forms derived are also amenable for use with symbolic-manipulation packages for systematic code generation.

Description: A computerized procedure for symbolic manipulations and FORTRAN code generation of an elasto-plastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.