A homogenization theory of strain gradient single crystal plasticity and its finite element discretization

In this study, a homogenization theory based on the Gurtin strain gradient formulation and its finite element discretization are developed for investigating the size effects on macroscopic responses of periodic materials. To derive the homogenization equations consisting of the relation of macroscopic stress, the weak form of stress balance, and the weak form of microforce balance, the Y-periodicity is used as additional, as well as standard, boundary conditions at the boundary of a unit cell. Then, by applying a tangent modulus method, a set of finite element equations is obtained from the homogenization equations. The computational stability and efficiency of this finite element discretization are verified by analyzing a model composite. Furthermore, a model polycrystal is analyzed for investigating the grain size dependence of polycrystal plasticity. In this analysis, the micro-clamped, micro-free, and defect-free conditions are considered as the additional boundary conditions at grain boundaries, and their effects are discussed.     

On the homogenization of metal matrix composites using strain gradient plasticity

The homogenized response of metal matrix composites（MMC） is studied using strain gradient plasticity.The material model employed is a rate independent formulation of energetic strain gradient plasticity at the micro scale and conventional rate independent plasticity at the macro scale. Free energy inside the micro structure is included due to the elastic strains and plastic strain gradients. A unit cell containing a circular elastic fiber is analyzed under macroscopic simple shear in addition to transverse and longitudinal loading. The analyses are carried out under generalized plane strain condition. Micro-macro homogenization is performed observing the Hill-Mandel energy condition,and overall loading is considered such that the homogenized higher order terms vanish. The results highlight the intrinsic size-effects as well as the effect of fiber volume fraction on the overall response curves, plastic strain distributions and homogenized yield surfaces under different loading conditions. It is concluded that composites with smaller reinforcement size have larger initial yield surfaces and furthermore,they exhibit more kinematic hardening.     

A deformation theory of strain gradient crystal plasticity that accounts for geometrically necessary dislocations

We propose a deformation theory of strain gradient crystal plasticity that accounts for the density of geometrically necessary dislocations by including, as an independent kinematic variable, Nye's dislocation density tensor [1953. Acta Metallurgica 1, 153-162]. This is accomplished in the same fashion as proposed by Gurtin and co-workers (see, for instance, Gurtin and Needleman [2005. J. Mech. Phys. Solids 53, 1-31]) in the context of a flow theory of crystal plasticity, by introducing the so-called defect energy. Moreover, in order to better describe the strengthening accompanied by diminishing size, we propose that the classical part of the plastic potential may be dependent on both the plastic slip vector and its gradient; for single crystals, this also makes it easier to deal with the “higher-order” boundary conditions. We develop both the kinematic formulation and its static dual and apply the theory to the simple shear of a constrained strip (example already exploited in Shu et al. [2001. J. Mech. Phys. Solids 49, 1361-1395], Bittencourt et al. [2003. J. Mech. Phys. Solids 51, 281-310], Niordson and Hutchinson [2003. Euro J. Mech. Phys. Solids 22, 771-778], Evers et al. [2004. J. Mech. Phys. Solids 52, 2379-2401], and Anand et al. [2005. J. Mech. Phys. Solids 53, 1789-1826]) to investigate what sort of behaviour the new model predicts. The availability of the total potential energy functional and its static dual allows us to easily solve this simple boundary value problem by resorting to the Ritz method.