Convergence Analysis of Rate-Dependent J2 Plasticity
End of summer internship write-up/memo
Elmer, Alexis Hope [Sandia National Laboratories (↗
SEARCH · Search NASA
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Light thermal structures and materials for high speed flight
Over the last twenty years, unified viscoplastic constitutive models have evolved to meet this need. These constitutive models provide a means for representing a material's response from the elastic through the plastic range including strain-rate dependent plastic flow, creep, and stress relaxation. Rate-dependent plasticity effects are known to be important at elevated temperatures. The purpose of this paper is to describe computational and experimental research programs underway at the Light Thermal Structures Center focused on the investigation of the response of structures and materials to local heating. In the first part of the paper, finite element thermoviscoplastic analysis is highlighted. In the second part of the paper, the thermal-structures experimental program is outlined.
Thornton, Earl A.↗
Mathematical characterization of mechanical behavior of porous frictional granular media
A new definition of loading and unloading along the yield surface of Roscoe and Burland is introduced. This is achieved by noting that the strain-hardening parameter in the plastic potential function is deduced from the yield locus equation of Roscoe and Burland. The analytical results are compared with the experimental results for plate-bearing and cone-penetrometer problems and close agreements are demonstrated. The wheel-soil interaction is studied under dynamic loading. The rate-dependent plasticity or viscoelastoplastic behavior is considered. This is accomplished by the internal (hidden) variables associated with time-dependent viscous properties directly superimposed with inelastic behavior governed by the yield criteria of Roscoe and Burland. Effects of inertia and energy dissipation are properly accounted for. Example problems are presented.
Chung, T. J.↗
Advances and trends in structures and dynamics; Proceedings of the Symposium, Washington, DC, October 22-25, 1984
Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.
Noor, A. K.↗
Time dependent behavior of a graphite/thermoplastic composite and the effects of stress and physical aging
Two complimentary studies were performed to determine the effects of stress and physical aging on the matrix dominated time dependent properties of IM7/8320 composite. The first of these studies, experimental in nature, used isothermal tensile creep/aging test techniques developed for polymers and adapted them for testing of the composite material. From these tests, the time dependent transverse (S22) and shear (S66) compliance's for an orthotropic plate were found from short term creep compliance measurements at constant, sub-T(sub g) temperatures. These compliance terms were shown to be affected by physical aging. Aging time shift factors and shift rates were found to be a function of temperature and applied stress. The second part of the study relied upon isothermal uniaxial tension tests of IM7/8320 to determine the effects of physical aging on the nonlinear material behavior at elevated temperature. An elastic/viscoplastic constitutive model was used to quantify the effects of aging on the rate-independent plastic and rate-dependent viscoplastic response. Sensitivity of the material constants required by the model to aging time were determined for aging times up to 65 hours. Verification of the analytical model indicated that the effects of prior aging on the nonlinear stress/strain/time data of matrix dominated laminates can be predicted.
Gates, Thomas S.↗
Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems
Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.
Impact↗
Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling
Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.
Dynamic simulation↗
Rate dependent constitutive models for fiber reinforced polymer composites
A literature survey was conducted to assess the state-of-the-art in rate dependent constitutive models for continuous fiber reinforced polymer matrix composite (PMC) materials. Several recent models which include formulations for describing plasticity, viscoelasticity, viscoplasticity, and rate-dependent phenomenon such as creep and stress relaxation are outlined and compared. When appropriate, these comparisons include brief descriptions of the mathematical formulations, the test procedures required for generating material constants, and details of available data comparing test results to analytical predictions.
Gates, Thomas S.↗
Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling
Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .
Close-clearance impellers↗
Impact Simulation of Composite Panels for Aerospace Applications
The aerospace research community has been working towards the development of a comprehensive framework for the design of composite structures at the component, subassembly, assembly, and system levels. Such a framework requires accurate, efficient and easy-to-use mathematical and computational models that are in sync with physical reality. The research work in this paper is designed to help an aerospace structural engineer move towards this goal by simulating dynamic impact to composite panels. First, the existing orthotropic visco-elastic visco-plastic material model in LS-DYNA commercial finite element program, MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE (MAT_213) is enhanced with the development and use of rate-dependent damage data. Using state-of-the-art experimental procedures, impact validation test data are created. Second, the versatility of the constitutive model is examined by its ability to simulate a low velocity (27 ft/s, 8 m/s) impact test, and two high velocity impact tests (385–417 ft/s, 117–127 m/s) that straddle the ballistic limit of the projectile and the composite panels. Third, these impact events are modeled using three different finite element types – under-integrated and fully-integrated solid elements, and fully-integrated thin shell elements. Insight into accuracy versus numerical efficiency is gained by comparing the performance of solid versus thin shell elements. The numerical results show that accurate predictions of these impact events can be obtained with reasonable effort and care.
Simulations↗
A rheological model for loose sands with insights from DEM
A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.
Elasticity↗
Rate-dependent mode I interlaminar crack growth mechanisms in graphite/epoxy and graphite/PEEK
In this paper the mode I fracture behavior of graphite/epoxy and graphite/PEEK composites is examined over four decades of crosshead rates (0.25-250 mm/min). Straight-sided double-cantilever-beam specimens consisting of unidirectional laminates were tested at room temperature. For graphite/epoxy the load-deflection response was linear to fracture, and stable slow crack growth initiating at the highest load level was observed for all rates tested. In contrast, mode I crack growth in the graphite/PEEK material was often unstable and showed stick-slip behavior. Subcritical crack growth occurring prior to the onset of fracture was observed at intermediate displacement rates. A mechanism for the fracture behavior of the graphite/PEEK material (based on viscoelastic, plastic, and microcrack coalescence in the process zone) is proposed and related to the observed rate-dependent phenomena.
Gillespie, J. W., Jr.↗
2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities
The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.
36 MATERIALS SCIENCE↗
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.↗
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↗
Experimental Characterization of Non-Associative Plasticity Flow Rule Coefficients and Post-Peak Stress Degradation for the LS-DYNA MAT213 Model
This project is focused on developing an experimental framework for characterizing non-associative plasticity flow rule coefficients through coupon-scale tests for the LS-DYNA MAT213 model. The main objective is to characterize these coefficients based on the multi-scale (i.e., both microscopic and macroscopic) full-field measurement of the evolution of strain and stress fields. This paper focuses on presenting the experimental work on characterizing the full-scale stress-strain curves of T700/LM-PAEK composites under tension, compression, and shear loads. The experimental data set was intended to build a deformation sub-model in the MAT213 model for the material. The strain data were collected using both microscopic and macroscopic digital image correlation techniques. The microscopic technique was particularly useful for fracture cases under small strains. A preliminary simulation result obtained from the MAT213 model is also presented in the paper. The experimental framework herein will be extended to characterize post-peak stress degradation in the composite material and to develop a damage sub-model for the material. This project will contribute to developing a simulation tool based on the MAT213 model for simulating the rate-dependent impact damage in composites under multi-axial loading.
MAT213↗
Slow growth of cracks in a rate sensitive Tresca solid.
Extension 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 equal 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 it 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.
Wnuk, M. P.↗
C-band SAR observations of marginal ice zone rheology in the Labrador Sea
Synthetic-aperture radar (SAR) imagery collected during the LIMEX'87 field program in the Labrador Sea is used to qualitatively interpret the deformational modes in play throughout the experiment. The ice cover exhibited two distinct rheologies separated by a clear line of shear: a nonlinear viscous outer regime, and a quasi-brittle inner regime. A single constitutive relation capable of effectively modeling both is unlikely within a plastic rate-independent formulation. Rate-dependent effects are discussed in relation to favoring brittle fraction in an otherwise ductile material.
Drinkwater, Mark R.↗