Computational Algorithms for a Damage-Coupled Cyclic Viscoplasticity Material Model
The primary objective of the investigation is to develop efficient and robust computational schemes for a damage-coupled cyclic thermoviscoplasticity model for solder material. Three constitutive integration algorithms, Euler, modified Euler and semi-implicit algorithm for the model are examined. The three algorithms for the model are coded in the commercial finite element (FE) code ABAQUS (version 6.21) via the user-defined material subroutine UMAT. Two single-step algorithms of the substep scheme are applied for the modified Euler algorithm to control the error in the integration of constitutive laws. A semi-empirical formulation is established for an adaptive time stepping algorithm that is based on the Euler algorithm. Single-element, miniature specimen and notched specimen simulations have been conducted to compare with the test results, which include monotonic tensile, creep and fatigue tests of 63Sn-37Pb solder. It is observed that the explicit algorithm consistently requires much less CPU time than others. The modified Euler algorithm has shown on the other hand to be not only efficient but also accurate. The semi-implicit algorithm yields accurate solution. It is worth noting that the method is also effective when an appropriate integration scheme is chosen.