On the equivalence between traction- and stress-based approaches for the modeling of localized failure in solids

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 (J₂), Drucker–Prager, Mohr–Coulomb and many other frequently employed criteria recovered as its particular cases.     

Self-consistent modeling of large plastic deformation, texture and morphology evolution in semi-crystalline polymers

A self-consistent model for semi-crystalline polymers is proposed to study their constitutive behavior, texture and morphology evolution during large plastic deformation. The material is considered as an aggregate of composite inclusions, each representing a stack of crystalline lamellae with their adjacent amorphous layers. The deformation within the inclusions is volume-averaged over the phases. The interlamellar shear is modeled as an additional slip system with a slip direction depending on the inclusion's stress. Hardening of the amorphous phase due to molecular orientation and, eventually, coarse slip, is introduced via Arruda-Boyce hardening law for the corresponding plastic resistance. The morphology evolution is accounted for through the change of shape of the inclusions under the applied deformation gradient. The overall behavior is obtained via a viscoplastic tangent self-consistent scheme. The model is applied to high density polyethylene (HDPE). The stress-strain response, texture and morphology changes are simulated under different modes of straining and compared to experimental data as well as to the predictions of other models.     

Three-dimensional modeling and numerical analysis of rate-dependent irrecoverable deformation in shape memory alloys

Shape memory alloys (SMAs) provide an attractive solid-state actuation alternative to engineers in various fields due to their ability to exhibit recoverable deformations while under substantial loads. Many constitutive models describing this repeatable phenomenon have been proposed, where some models also capture the effects of rate-independent irrecoverable deformations (i.e., plasticity) in SMAs. In this work, we consider a topic not addressed to date: the generation and evolution of irrecoverable viscoplastic strains in an SMA material. Such strains appear in metals subjected to sufficiently high temperatures. The need to account for these effects in SMAs arises when considering one of two situations: the exposure of a conventional SMA material (e.g., NiTi) to high temperatures for a non-negligible amount of time, as occurs during shape-setting, or the utilization of new high temperature shape memory alloys (HTSMAs), where the elevated transformation temperatures induce transformation and viscoplastic behaviors simultaneously. A new three-dimensional constitutive model based on established SMA and viscoplastic modeling techniques is derived that accounts for these behaviors. The numerical implementation of the model is described in detail. Several finite element analysis (FEA) examples are provided, demonstrating the utility of the new model and its implementation in assessing the effects of viscoplastic behaviors in shape memory alloys.     

A generic approach towards finite growth with examples of athlete's heart, cardiac dilation, and cardiac wall thickening

The objective of this work is to establish a generic continuum-based computational concept for finite growth of living biological tissues. The underlying idea is the introduction of an incompatible growth configuration which naturally introduces a multiplicative decomposition of the deformation gradient into an elastic and a growth part. The two major challenges of finite growth are the kinematic characterization of the growth tensor and the identification of mechanical driving forces for its evolution. Motivated by morphological changes in cell geometry, we illustrate a micromechanically motivated ansatz for the growth tensor for cardiac tissue that can capture both strain-driven ventricular dilation and stress-driven wall thickening. Guided by clinical observations, we explore three distinct pathophysiological cases: athlete's heart, cardiac dilation, and cardiac wall thickening. We demonstrate the computational solution of finite growth within a fully implicit incremental iterative Newton-Raphson based finite element solution scheme. The features of the proposed approach are illustrated and compared for the three different growth pathologies in terms of a generic bi-ventricular heart model.     

Micromechanics, macromechanics and constitutive modeling of the elasto-viscoplastic deformation of rubber-toughened glassy polymers

Under certain conditions, such as sufficiently low temperatures, high loading rates and/or highly triaxial stress states, glassy polymers display an unfavorable characteristic—brittleness. A technique used for reducing the brittleness (increasing the fracture toughness) of these materials is rubber toughening. While there is significant qualitative understanding of the mechanical behavior of rubber-toughened polymers, quantitative modeling tools for the large-strain deformation of rubber-toughened glassy polymers are largely lacking.In this paper, we develop a suite of numerical tools to investigate the mechanical behavior of rubber-toughened glassy polymers, with emphasis on rubber-toughened polycarbonate. The rubber particles are modeled as voids in view of their deformation-induced cavitation early during deformation. A three-dimensional micromechanical model of the heterogeneous microstructure is developed to study the effects of initial rubber particle (void) volume fraction on the underlying elasto-viscoplastic deformation mechanisms in the material, and how these mechanisms influence the macroscopic response of the material. A continuum-level constitutive model is developed for the large-strain elasto-viscoplastic deformation of porous glassy polymers, and it is calibrated against micromechanical modeling results for porous polycarbonate. The constitutive model can be used to study various boundary value problems involving rubber-toughened (porous) glassy polymers. As an example, the case of an axisymmetric notched bar is simulated for the case of polycarbonate with varying levels of initial porosity. The quality of the constitutive model calibration is assessed using a multi-scale modeling approach.