Numerical Modeling of Transient Creep in Polycrystalline Ice

Transient creep, an important deformation mechanism for polycrystalline ice at quasi-static strain rates, is characterized by rate and temperature sensitivity, by isotropic and kinematic strain hardening, as well as by fabric and deformation-induced anisotropy. A physically based constitutive model,...

Full description

Saved in:
Bibliographic Details
Published in:Journal of engineering mechanics 1993-10, Vol.119 (10), p.2011-2035
Main Authors: Shyam Sunder, S, Elvin, Alex, Nanthikesan, S
Format: Article
Language:English
Subjects:
Algorithms
Boundary value problems
Creep
Crystalline materials
Differential equations
Earth sciences
Earth, ocean, space
Engineering and environment geology. Geothermics
Engineering geology
Exact sciences and technology
External geophysics
Finite element method
Integration
Mathematical models
Matrix algebra
Nonlinear equations
Snow. Ice. Glaciers
Strain
Subroutines
TECHNICAL PAPERS
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Transient creep, an important deformation mechanism for polycrystalline ice at quasi-static strain rates, is characterized by rate and temperature sensitivity, by isotropic and kinematic strain hardening, as well as by fabric and deformation-induced anisotropy. A physically based constitutive model, using internal state variables, has been developed by Shyam Sunder and Wu (1989a, b) to describe the multiaxial behavior of ice undergoing transient creep. To solve boundary value problems using this constitutive theory requires the numerical time integration of a coupled set of stiff and highly nonlinear first-order differential equations. A closed-form Newton-Raphson (tangent) formulation, in conjunction with the -method of integration, is developed to solve the constitutive equations. The fully consistent constitutive Jacobian matrix that is used to assemble the finite element tangent stiffness matrix is also established in closed form. This algorithm is implemented as a subroutine in the finite element program ABAQUS and its predictions are verified against experimental data and known solutions. The importance of transient creep is demonstrated by performing simulations of: (1) Arrested subsurface penetration; and (2) in-plane indentation of a floating ice sheet.
ISSN:0733-9399
1943-7889
DOI:10.1061/(ASCE)0733-9399(1993)119:10(2011)