Abstract

This work investigates systematically strain localization and failure mechanics for elastoplastic damage solids. Two complementary methodologies, i.e., traction-based discontinuities localized in an elastic solid and strain localization of a stress-based inelastic softening solid, are addressed. In the former it is assumed a priori that the discontinuity (band) forms with a continuous stress field and along the known orientation. A traction-based failure criterion is introduced to characterize the discontinuity (band) and the orientation is determined from Mohr’s maximization postulate. If the (apparent) displacement jumps are retained as independent variables, the strong/regularized discontinuity approaches follow, requiring constitutive models for both the bulk and discontinuity (band). Elimination of the displacement jumps at the material point level results in the embedded/smeared discontinuity approaches in which an overall inelastic constitutive model fulfilling the static constraint suffices. The second methodology is then adopted to check whether the assumed strain localization can occur and identify its consequences on the resulting approaches. The kinematic constraint guaranteeing stress boundedness/continuity upon strain localization is established for general inelastic softening solids. Application to a unified elastoplastic damage model naturally yields all the ingredients of a localized model for the discontinuity (band), justifying the first methodology. Two dual but not necessarily equivalent approaches, i.e., the traction-based elastoplastic damage model and the stress-based projected discontinuity model, are identified. The former is equivalent to the embedded/smeared discontinuity approaches, whereas in the later the discontinuity orientation and associated failure criterion, not given a priori, are determined consistently from the kinematic constraint. The bi-directional connections and equivalence conditions between the traction- and stress-based approaches are classified. Closed-form 2D results under plane stress condition are also given, with the classical Rankine, Mohr-Coulomb, von Mises and Drucker-Prager criteria analyzed as the illustrative examples. A generic failure criterion of either elliptic, parabolic or hyperbolic type, is then considered in a unified manner, resulting in many failure criteria frequently employed in practice.

Abstract

Damage-induced strain softening is of vital importance for the modeling of localized failure in frictional-cohesive materials. This paper addresses strain localization of damaging solids and the resulting consistent frictional-cohesive crack models. As a supplement to the framework recently established for stress-based continuum material models in rate form (Wu and Cervera 2015, 2016), several classical strain-based damage models, expressed usually in total and secant format, are considered. Upon strain localization of such damaging solids, Maxwell’s kinematics of a strong (or regularized) discontinuity has to be reproduced by the inelastic damage strains, which are defined by a bounded characteristic tensor and an unbounded scalar related to the damage variable. This kinematic constraint yields a set of nonlinear equations from which the discontinuity orientation and damage-type localized cohesive relations can be derived. It is found that for the “Simó and Ju 1987” isotropic damage model, the localization angles and the resulting cohesive model heavily depend on lateral deformations usually ignored in classical crack models for quasi-brittle solids. To remedy this inconsistency, a modified damage model is proposed. Its strain localization analysis naturally results in a consistent frictional-cohesive crack model of damage type, which can be regularized as a classical smeared crack model. The analytical results are numerically verified by the recently-proposed mixed stabilized finite element method, regarding a singly-perforated plate under uniaxial tension. Remarkably, for all of the damage models discussed in this work, the numerically-obtained localization angles agree almost exactly with the closed-form results. This agreement, on the one hand, consolidates the strain localization analysis based on Maxwell’s kinematics and, on the other hand, illustrates versatility of the mixed stabilized finite element method.

Abstract

This work investigates systematically traction- and stress-based approaches for the modeling of strong and regularized discontinuities induced by localized failure in solids. Two complementary methodologies, i.e., discontinuities localized in an elastic solid and strain localization of an inelastic softening solid, are addressed. In the former it is assumed a priori that the discontinuity forms with a continuous stress field and along the known orientation. A traction-based failure criterion is introduced to characterize the discontinuity and the orientation is determined from Mohr's maximization postulate. If the displacement jumps are retained as independent variables, the strong/regularized discontinuity approaches follow, requiring constitutive models for both the bulk and discontinuity. Elimination of the displacement jumps at the material point level results in the embedded/smeared discontinuity approaches in which an overall inelastic constitutive model fulfilling the static constraint suffices. The second methodology is then adopted to check whether the assumed strain localization can occur and identify its consequences on the resulting approaches. The kinematic constraint guaranteeing stress boundedness and continuity upon strain localization is established for general inelastic softening solids. Application to a unified stress-based elastoplastic damage model naturally yields all the ingredients of a localized model for the discontinuity (band), justifying the first methodology. Two dual but not necessarily equivalent approaches, i.e., the traction-based elastoplastic damage model and the stress-based projected discontinuity model, are identified. The former is equivalent to the embedded and smeared discontinuity approaches, whereas in the later the discontinuity orientation and associated failure criterion are determined consistently from the kinematic constraint rather than given a priori. The bi-directional connections and equivalence conditions between the traction- and stress-based approaches are classified. Closed-form results under plane stress condition are also given. A generic failure criterion of either elliptic, parabolic or hyperbolic type is analyzed in a unified manner, with the classical von Mises (J2), Drucker–Prager, Mohr–Coulomb and many other frequently employed criteria recovered as its particular cases.

Abstract

Aiming for the modeling of localized failure in quasi-brittle solids, this paper addresses a thermodynamically consistent plastic-damage framework and the corresponding strain localization analysis. A unified elastoplastic damage model is first presented based on two alternative kinematic decompositions assuming infinitesimal deformations, with the evolution laws of involved internal variables characterized by a dissipative flow tensor. For the strong (or regularized) discontinuity to form in such inelastic quasi-brittle solids and to evolve eventually into a fully softened one, a novel strain localization analysis is then suggested. A kinematic constraint more demanding than the classical discontinuous bifurcation condition is derived by accounting for the traction continuity and the loading/unloading states consistent with the kinematics of a strong (or regularized) discontinuity. More specifically, the strain jumps characterized by Maxwell’s kinematic condition have to be completely inelastic (energy dissipative). Reproduction of this kinematics implies vanishing of the aforesaid dissipative flow tensorial components in the directions orthogonal to the discontinuity orientation. This property allows naturally developing a localized plastic-damage model for the discontinuity (band), with its orientation and the traction-based failure criterion consistently determined a posteriori from the given stress-based counterpart. The general results are then particularized to the 2D conditions of plane stress and plane strain. It is found that in the case of plane stress, strain localization into a strong (or regularized) discontinuity can occur at the onset of strain softening. Contrariwise, owing to an extra kinematic constraint, in the condition of plane strain some continuous inelastic deformations and substantial re-orientation of principal strain directions in general have to take place in the softening regime prior to strain localization. The classical Rankine, Mohr–Coulomb, von Mises (J2) and Drucker–Prager criteria are analyzed as illustrative examples. In particular, both the closed-form solutions for the discontinuity angles validated by numerical simulations and the corresponding traction-based failure criteria are obtained.

Abstract

A general methodology to develop hyper-elastic membrane models applicable to crystalline films one-atom thick is presented. In this method, an extension of the Born rule based on the exponential map is proposed. The exponential map accounts for the fact that the lattice vectors of the crystal lie along the chords of the curved membrane, and consequently a tangent map like the standard Born rule is inadequate. In order to obtain practical methods, the exponential map is locally approximated. The effectiveness of our approach is demonstrated by numerical studies of carbon nanotubes. Deformed configurations as well as equilibrium energies of atomistic simulations are compared with those provided by the continuum membrane resulting from this method discretized by finite elements.