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.     

Phase-field slip-line theory of plasticity

A variational approach to determine the deformation of an ideally plastic substance is proposed by solving a sequence of energy minimization problems under proper conditions to account for the irreversible character of plasticity. The flow is driven by the local transformation of elastic strain energy into plastic work on slip surfaces, once that a certain energetic barrier for slip activation has been overcome. The distinction of the elastic strain energy into spherical and deviatoric parts is used to incorporate in the model the idea of von Mises plasticity and isochoric plastic strain. This is a “phase field model” because the matching condition at the slip interfaces is substituted by the evolution of an auxiliary phase field that, similar to a damage field, is unitary on the elastic phase and null on the yielded phase. The slip lines diffuse in bands, whose width depends upon a material length-scale parameter.Numerical experiments on representative problems in plane strain give solutions with noteworthy similarities with the results from classical slip-line field theory, but the proposed model is much richer because, accounting for elastic deformations, it can describe the formation of slip bands at the local level, which can nucleate, propagate, widen and diffuse by varying the boundary conditions. In particular, the solution for a long pipe under internal pressure is very different from the one obtainable from the classical macroscopic theory of plasticity. For this case, the location of the plastic bands may be an insight to explain the premature failures that are sometimes encountered during the manufacturing process. This practical example enhances the importance of this new theory based on the mathematical sciences.     

Mode I and mixed mode crack-tip fields in strain gradient plasticity

Strain gradients develop near the crack-tip of Mode I or mixed mode cracks. A finite strain version of the phenomenological strain gradient plasticity theory of Fleck-Hutchinson (2001) is used here to quantify the effect of the material length scales on the crack-tip stress field for a sharp stationary crack under Mode I and mixed mode loading. It is found that for material length scales much smaller than the scale of the deformation gradients, the predictions converge to conventional elastic-plastic solutions. For length scales sufficiently large, the predictions converge to elastic solutions. Thus, the range of length scales over which a strain gradient plasticity model is necessary is identified. The role of each of the three material length scales, incorporated in the multiple length scale theory, in altering the near-tip stress field is systematically studied in order to quantify their effect.     

A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals

A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals is developed. The theory accounts for: an arbitrary number and arrangement of dislocation lines over a slip plane; the long-range elastic interactions between dislocation lines; the core structure of the dislocations resulting from a piecewise quadratic Peierls potential; the interaction between the dislocations and an applied resolved shear stress field; and the irreversible interactions with short-range obstacles and lattice friction, resulting in hardening, path dependency and hysteresis. A chief advantage of the present theory is that it is analytically tractable, in the sense that the complexity of the calculations may be reduced, with the aid of closed form analytical solutions, to the determination of the value of the phase field at point-obstacle sites. In particular, no numerical grid is required in calculations. The phase-field representation enables complex geometrical and topological transitions in the dislocation ensemble, including dislocation loop nucleation, bow-out, pinching, and the formation of Orowan loops. The theory also permits the consideration of obstacles of varying strengths and dislocation line-energy anisotropy. The theory predicts a range of behaviors which are in qualitative agreement with observation, including: hardening and dislocation multiplication in single slip under monotonic loading; the Bauschinger effect under reverse loading; the fading memory effect, whereby reverse yielding gradually eliminates the influence of previous loading; the evolution of the dislocation density under cycling loading, leading to characteristic ‘butterfly’ curves; and others.     

Well-posedness of a model of strain gradient plasticity for plastically irrotational materials

The initial boundary value problem corresponding to a model of strain gradient plasticity due to [Gurtin, M., Anand, L., 2005. A theory of strain gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations. J. Mech. Phys. Solids 53, 1624–1649] is formulated as a variational inequality, and analysed. The formulation is a primal one, in that the unknown variables are the displacement, plastic strain, and the hardening parameter. The focus of the analysis is on those properties of the problem that would ensure existence of a unique solution. It is shown that this is the case when hardening takes place. A similar property does not hold for the case of softening. The model is therefore extended by adding to it terms involving the divergence of plastic strain. For this extended model the desired property of coercivity holds, albeit only on the boundary of the set of admissible functions.     

A simple theory of strain gradient plasticity based on stress-induced anisotropy of defect diffusion

A new strain gradient plasticity theory is formulated to accommodate more than one material length parameter. The theory is an extension of the classical J₂ flow theory of metal plasticity to the micron scale. Distinctive features of the proposed theory as compared to other existing theories are the simplicities of mathematical formulation, numerical implementation and physical interpretation.     

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and micro-devices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either ‘microhard’ (impenetrable to dislocations) or ‘microfree’ (an infinite dislocation sink).     

Simulation of micro-indentation hardness of FCC single crystals by mechanism-based strain gradient crystal plasticity

The size effect observed in the micro-indentation of FCC single crystal copper is modelled by the employment of mechanism-based strain gradient crystal plasticity (MSG-CP). The total slip resistance in each active system is assumed to be due to a mixed population of forest obstacles arising from both statistically stored and geometrically necessary dislocations. The MSG-CP constitutive model is implemented into the Abaqus/Standard FE platform by developing the User MATerial subroutine UMAT. The simulation of micro-indentation hardness on (0 0 1) and (1 1 1) single crystal copper, with a conical indenter having a sharp tip, a conical indenter with a spherical tip and a three-sided Berkovich indenter, is undertaken. The phenomena of pile-up and sink-in have been observed in the simulation and dealt with by appropriate use of the contact analysis function in Abaqus. These phenomena have been taken into account in the determination of the contact areas and hence the average indentation depth for anisotropic single crystals. The depth dependence of the micro-indentation hardness, the size effect, is calculated. The micro-hardness results from the simulation are compared with those of the published experimental ones in the literature and a good agreement is found.     

A constitutive model for the formation of martensite in austenitic steels under large strain plasticity

A constitutive model for diffusionless phase transitions in elastoplastic materials undergoing large deformations is developed. The model takes basic thermodynamic relations as its starting point and the phase transition is treated through an internal variable (the phase fractions) approach. The usual yield potential is used together with a transformation potential to describe the evolution of the new phase. A numerical implementation of the model is presented, along with the derivation of a consistent algorithmic tangent modulus. Simulations based on the presented model are shown to agree well with experimental findings. The proposed model provides a robust tool suitable for large-scale simulations of phase transformations in austenitic steels undergoing extensive deformations, as is demonstrated in simulations of the necking of a bar under tensile loading and also in simulations of a cup deep-drawing process.     

The role of partial mediated slip during quasi-static deformation of 3D nanocrystalline metals

We present dislocation simulations involving the collective behavior of partials and extended full dislocations in nanocrystalline materials. While atomistic simulations have shown the importance of including partial dislocations in high strain rate simulations, the behavior of partial dislocations in complex geometries with lower strain rates has not been explored. To account for the dissociation of dislocations into partials we include the full representation of the gamma surface for two materials: Ni and Al. During loading, dislocation loops are emitted from grain boundaries and expand into the grain interiors to carry the strain. In agreement with high strain rate simulations we find that Al has a higher density of extended full dislocations with smaller stacking fault widths than Ni. We also observe that configurations with smaller average grain size have a higher density of partial dislocations, but contrary to simplified analytical models we do not find a critical grain size below which there is only partial dislocation-mediated deformation. Our results show that the density of partial dislocations is stable in agreement with in situ X-ray experiments that show no increase of the stacking fault density in deformed nanocrystalline Ni (Budrovic et al., 2004). Furthermore, the ratio between partial and extended full dislocation contribution to strain varies with the amount of deformation. The contribution of extended full dislocations to strain grows beyond the contribution of partial dislocations as the deformation proceeds, suggesting that there is no well-defined transition from full dislocation- to partial dislocation-mediated plasticity based uniquely on the grain size.     

Parameter identification of advanced plastic strain rate potentials and impact on plastic anisotropy prediction

In the work presented in this paper, several strain rate potentials are examined in order to analyze their ability to model the initial stress and strain anisotropy of several orthotropic sheet materials. Classical quadratic and more advanced non-quadratic strain rate potentials are investigated in the case of FCC and BCC polycrystals. Different identifications procedures are proposed, which are taking into account the crystallographic texture and/or a set of mechanical test data in the determination of the material parameters.     

On stress-state dependent plasticity modeling: Significance of the hydrostatic stress, the third invariant of stress deviator and the non-associated flow rule

It has been shown that the plastic response of many materials, including some metallic alloys, depends on the stress state. In this paper, we describe a plasticity model for isotropic materials, which is a function of the hydrostatic stress as well as the second and third invariants of the stress deviator, and present its finite element implementation, including integration of the constitutive equations using the backward Euler method and formulation of the consistent tangent moduli. Special attention is paid for the adoption of the non-associated flow rule. As an application, this model is calibrated and verified for a 5083 aluminum alloy. Furthermore, the Gurson-Tvergaard-Needleman porous plasticity model, which is widely used to simulate the void growth process of ductile fracture, is extended to include the effects of hydrostatic stress and the third invariant of stress deviator on the matrix material.     

A physically based model of temperature and strain rate dependent yield in BCC metals: Implementation into crystal plasticity

In this work, we develop a crystal plasticity finite element model (CP-FEM) that constitutively captures the temperature and strain rate dependent flow stresses in pure BCC refractory metals. This model is based on the kink-pair theory developed by Seeger (1981) and is calibrated to available data from single crystal experiments to produce accurate and convenient constitutive laws that are implemented into a BCC crystal plasticity model. The model is then used to predict temperature and strain rate dependent yield stresses of single and polycrystal BCC refractory metals (molybdenum, tantalum, tungsten and niobium) and compared with existing experimental data. To connect to larger length scales, classical continuum-scale constitutive models are fit to the CP-FEM predictions of polycrystal yield stresses. The results produced by this model, based on kink-pair theory and with origins in dislocation mechanics, show excellent agreement with the Mechanical Threshold Stress (MTS) model for temperature and strain-rate dependent flow. This framework provides a method to bridge multiple length scales in modeling the deformation of BCC metals.     

Constitutive modeling of temperature and strain rate dependent elastoplastic hardening materials using a corotational rate associated with the plastic deformation

In this paper, a constitutive model with a temperature and strain rate dependent flow stress (Bergstrom hardening rule) and modified Armstrong-Frederick kinematic evolution equation for elastoplastic hardening materials is introduced. Based on the multiplicative decomposition of the deformation gradient,new kinematic relations for the elastic and plastic left stretch tensors as well as the plastic deformation-dependent spin tensor are proposed. Also, a closed-form solution has been obtained for the elastic and plastic left stretch tensors for the simple shear problem.To evaluate model validity, results are compared with known experimental data for SUS 304 stainless steel, which shows a good agreement with the results of the proposed theoretical model.Finally, the stress-deformation curve, as predicted by the model, is plotted for the simple shear problem at room and elevated temperatures using the same material properties for AA5754-O aluminium alloy.     

Fracture of bicrystal metal/ceramic interfaces: A study via the mechanism-based strain gradient crystal plasticity theory

Two continuum mechanical models of crystal plasticity theory namely, conventional crystal plasticity theory and mechanism-based crystal plasticity theory, are used to perform a comparative study of stresses that are reached at and ahead of the crack tip of a bicrystal niobium/alumina specimen. Finite element analyses are done for a stationary crack tip and growing cracks using a cohesive modelling approach. Using mechanism-based strain gradient crystal plasticity theory the stresses reached ahead of the crack tip are found to be two times larger than the stresses obtained from conventional crystal plasticity theory. Results also show that strain gradient effects strongly depend on the intrinsic material length to the size of plastic zone ratio (l/R₀). It is found that the larger the (l/R₀) ratio, the higher the stresses reached using mechanism-based strain gradient crystal plasticity theory. An insight into the role of cohesive strength and work of adhesion in macroscopic fracture is also presented which can be used by experimentalists to design better bimaterials by varying cohesive strength and work of adhesion.     

Dislocation dynamics: from microscopic models to macroscopic crystal plasticity

In this paper we study the connection between four models describing dislocation dynamics: a generalized 2D Frenkel-Kontorova model at the atomic level, the Peierls-Nabarro model, the discrete dislocation dynamics and a macroscopic model with dislocation densities. We show how each model can be deduced from the previous one at a smaller scale.      

Influence of material parameters and crystallography on the size effects describable by means of strain gradient plasticity

In the context of single-crystal strain gradient plasticity, we focus on the simple shear of a constrained strip in order to study the effects of the material parameters possibly involved in the modelling. The model consists of a deformation theory suggested and left undeveloped by Bardella [(2007). Some remarks on the strain gradient crystal plasticity modelling, with particular reference to the material length scales involved. Int. J. Plasticity 23, 296–322] in which, for each glide, three dissipative length scales are considered; they enter the model through the definition of an effective slip which brings into the isotropic hardening function the relevant plastic strain gradients, averaged by means of a p-norm. By means of the defect energy (i.e., a function of Nye's dislocation density tensor added to the free energy; see, e.g., Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5–32]), the model further involves an energetic material length scale. The application suggests that two dissipative length scales may be enough to qualitatively describe the size effect of metals at the microscale, and they are chosen in such a way that the higher-order state variables of the model be the dislocation densities. Moreover, we show that, depending on the crystallography, the size effect governed by the defect energy may be different from what expected (based on the findings of [Bardella, L., 2006. A deformation theory of strain gradient crystal plasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 54, 128–160] and [Gurtin et al. 2007. Gradient single-crystal plasticity with free energy dependent on dislocation densities. J. Mech. Phys. Solids 55, 1853–1878]), leading mostly to some strengthening. In order to investigate the model capability, we also exploit a Γ-convergence technique to find closed-form solutions in the “isotropic limit”. Finally, we analytically show that in the “perfect plasticity” case, should the dissipative length scales be set to zero, the presence of the sole energetic length scale may lead, as in standard plasticity, to non-uniqueness of solutions.     

Analysis of worm-like locomotion driven by the sine-squared strain wave in a linear viscous medium

In this paper, a worm-like locomotion in a linear resistive medium is studied to achieve controlled shape changes of the worm-like body by choosing a kind of driving with low energy expended and high-velocity locomotion in certain condition. To this end, we first develop the full dynamic model of the system under consideration to obtain the mean velocity related to friction coefficient, wave speed, linear density, body length and wave width. Correspondingly, a quasi-static model is also given from which the velocity can be expressed analytically. In the case of the shape change driven by the sine-squared strain wave (SSSW), it is seen that these two velocities will tend to uniformity with the friction coefficient or length of the body increasing or the wave speed decreasing when keeping the other parameters unchanged. Thus, the inertia term is ignorable for a large friction, a long body-length but a small wave-speed of the SSSW, which implies that the dynamical model can be reduced to the quasi-static one. The relative criterion is approximately given. As a result, the corresponding quasi-static model is employed to consider two typical drives, namely, the SSSW and the square strain wave (SSW). The result shows the shape change driven by the SSSW has an advantage in both the mean velocity and the average energy expended over that by the SSW when the necessary condition is satisfied. The analytical results are verified by numerical simulation.