In this paper, the coupled model with hydrogen transport and elasto-plasticity behavior is introduced. This model is implemented to the general-purpose FE code, ABAQUS, via the user-defined subroutine UMAT and UMATHT. In UMAT, the spatial gradients of hydrostatic stress and hydrogen induced deformation are calculated, and then are passed into UMATHT. Heat transfer equation within UMATHT is substituted by hydrogen transport equation including the effects of stress states and strain hardening. To validate this model, the finite element analyses coupled with hydrogen transport and mechanical loading are performed for the boundary layer specimens with low and high strength steel properties. The FE results are compared with the previous studies by Taha and Sofronis (2001).