This study is concerned with an application of the boundary element method on the crack problem. The stable and efficient analysis method of two dimensional elastostatic stress intensity factor on the mode I deformations is established from the result o stress analysis for the center cracked plates. In order to precisely analyse, The subelements of quadratic element, singular elements on the crack tip and interface and division into regions are applied to elastic stress. The usefulness of the method has been tested with a center cracked plates, a double edge cracked plate and a single edge cracked plate, and the results have turned out to be fairly satisfactory.