Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

Robust optimization of uncertain multistage inventory systems with inexact data in decision rules

In production-inventory problems customer demand is often subject to uncertainty. Therefore, it is challenging to design production plans that satisfy both demand and a set of constraints on e.g. production capacity and required inventory levels. Adjustable robust optimization (ARO) is a technique to solve these dynamic (multistage) production-inventory problems. In ARO, the decision in each stage is a function of the data on the realizations of the uncertain demand gathered from the previous periods. These data, however, are often inaccurate; there is much evidence in the information management literature that data quality in inventory systems is often poor. Reliance on data “as is” may then lead to poor performance of “data-driven” methods such as ARO. In this paper, we remedy this weakness of ARO by introducing a model that treats past data itself as an uncertain model parameter. We show that computational tractability of the robust counterparts associated with this extension of ARO is still maintained. The benefits of the new model are demonstrated by a numerical test case of a well-studied production-inventory problem. Our approach is also applicable to other ARO models outside the realm of production-inventory planning.     

Microstructural instabilities: dislocation substructures with pronounced mechanical effects

Microstructure evolution is largely dominated by the internal stress fields that appear upon the appearance of inhomogeneous structures in a material. The hardening behaviour of metals physically originates from such a complex microstructure evolution. As deformation proceeds, statistically homogeneous distributions of dislocations in grains become unstable, which constitutes the driving force for the development of a pronounced dislocation substructure. The dislocation structure already appears at early stages of deformation due to the statistical trapping of dislocations. Cell walls contain dislocation dipoles and multipoles with high dislocation densities and enclose cell-interior regions with a considerably smaller dislocation density. The presence and evolution of such a dislocation arrangement in the material influence the mechanical response of the material and is commonly associated with the transient hardening after strain path changes. This contribution introduces a micromechanical continuum model of the dislocation cell structure based on the physics of the dislocation interactions. The approximation of the internal stress field in such a microstructure and the impact on the macroscopic mechanical response are the main items investigated here.     

Modelling of the internal stress in dislocation cell structures

The nonuniform distribution of dislocations in metals gives rise to material anisotropy and internal stresses that determine the mechanical response. This paper proposes a micromechanical model of a dislocation cell structure that accounts for the material inhomogeneity and incorporates the internal stresses in a physically-based manner. A composite model is employed to describe the material with its dislocation cell structure. The internal stress is obtained as a natural result of plastic deformation incompatibility and incorporated in the composite model. Applications of this model enable the prediction of the mechanical behavior of metals under various nonuniform deformations. The implementation of the model is relatively straightforward, allowing easy use in macroscopic engineering computations.     

Gradient Estimation Using Lagrange Interpolation Polynomials

We use Lagrange interpolation polynomials to obtain good gradient estimations. This is e.g. important for nonlinear programming solvers. As an error criterion, we take the mean squared error, which can be split up into a deterministic error and a stochastic error. We analyze these errors using N-times replicated Lagrange interpolation polynomials. We show that the mean squared error is of order if we replicate the Lagrange estimation procedure N times and use 2d evaluations in each replicate. As a result, the order of the mean squared error converges to N ⁻¹ if the number of evaluation points increases to infinity. Moreover, we show that our approach is also useful for deterministic functions in which numerical errors are involved. We provide also an optimal division between the number of gridpoints and replicates in case the number of evaluations is fixed. Further, it is shown that the estimation of the derivatives is more robust when the number of evaluation points is increased. Finally, test results show the practical use of the proposed method. We thank Jack Kleijnen, Gül Gürkan, and Peter Glynn for useful remarks on an earlier version of this paper. We thank Henk Norde for the proof of Lemma 2.2.     

On the role of dislocation conservation in higher-order crystal plasticity

The interactions between individual dislocations contribute significantly to size effects as observed at the plastic deformation of miniaturized structures. When employing a crystal plasticity framework, these interactions are usually captured by an additional balance relation. In this paper we study two approaches to capture dislocation interactions in crystal plasticity that differ by the conservation of dislocations. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-convex rate dependent strain gradient crystal plasticity and deformation patterning

A rate dependent strain gradient crystal plasticity framework is presented where the displacement and the plastic slip fields are considered as primary variables. These coupled fields are determined on a global level by solving simultaneously the linear momentum balance and the slip evolution equation, which is derived in a thermodynamically consistent manner. The formulation is based on the 1D theory presented in Yalcinkaya et al. (2011), where the patterning of plastic slip is obtained in a system with non-convex energetic hardening through a phenomenological double-well plastic potential. In the current multi-dimensional multi-slip analysis the non-convexity enters the framework through a latent hardening potential presented in Ortiz and Repettto (1999) where the microstructure evolution is obtained explicitly via a lamination procedure. The current study aims the implicit evolution of deformation patterns due to the incorporated physically based non-convex potential.     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.