It is a tedious and excessive time consuming process to model the local area of crack tip part of structures in calculation of stress intensity factors by FEM. So, in this paper, the hybrid method of FEM and BEM approach was formulated to overcome this type of problems. The multi-domained BEM was adopted to simplify the modelling process of complex geometry and singularity characteristics of crack tip part and the ordinary FEM modelling was used in the rest part. The example calculations shows very good results compared with analytic solutions and other numerical method.