An algorithm is proposed for the finite element implementation of a rate-dependent crystal plasticity model. A 'planar double slip' two-dimensional crystal model is used. The constitutive equations are derived in detail, and it is shown that they reduce to a set of three ordinary differential equations: two for the resolved shear stresses on the slip planes, and one for the angle that defines the slip systems' orientation. The proposed algorithm, which is suitable for implementation in an explicit finite element code such as PRONTO2D, is based on a weighted residual method for the approximate solution of these equations. It turns out that the performance of the algorithm does not restrict the selection of the global time increment. The efficiency and performance of the algorithm is demonstrated through a comparison with the so-called rate-tangent method. Also, numerical examples are given to demonstrate the model's capability to capture shear localization phenomena at high strain rates, using the PRONTO2D explicit finite element code. © 1992.
Mendeley saves you time finding and organizing research
Choose a citation style from the tabs below