A seamless analysis of material behavior incorporating complex geometry and crack- tip modeling is one of greatly interesting topics in engineering and computational fracture mechanics fields. However, there are still large gaps between the industrial applications and fundamental academic studies due to a time consuming detailed modeling. In order to resolve this problem, a numerical method to calculate an energy release rate by virtual crack closure technique was proposed in this paper. Both free mesh method and finite element method have been utilized and, thereafter, robust local and global elements for various geometries and boundary conditions were generated. A validity of the proposed method has been demonstrated through a series of fracture mechanics analyses without tedious crack-tip meshing.