The paper has established a boundary element method with isoparameter elements for 3-D composite laminate by introducing boundary element subrigion theory. The interlaminar stresses around the hole in a [0ooooooooooooooooo/90ooooooooooooo]s composite laminates containing a circular holes and loaded in uniaxial tension were analyzed. The results by this 3-D BEM are in good agreement with that by 3-D FEM.