We are developing a GeoFEM module for simulating earthquake cycle of interplate great earthquakes at a subduction zone in a 3-D laterally heterogeneous viscoelastic medium. Here, toward this final goal, we present a 2-D simulation system of quasi-static earthquake cycle following the rate and state friction law which is derived from laboratory experiments, employing a FEM code of ABAQUS instead of GeoFEM. In the 1st ACES Workshop, we presented mainly a simple earthquake cycle due to a subducting slab with artificially changing of friction on the plate interface using rather coarse FEM meshes in a laterally heterogeneous viscoelastic medium under gravitation. In this study, we reconstruct a FEM model with fine meshes in the region close to the plate interface on which we apply the rate and state friction law. We present the concept of our FEM modeling of earthquake cycle of interplate great earthquakes due to a subducting slab and the details of our system with several numerical simulation results.