Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approxi-mately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper.In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the cor-ner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingularboundary integral equation numerically in a non-regularized form and in a local manner by using conforming C0 quadratic boundary ele-ments and standard Gaussian quadratures similar to those employed in the conventional displacement-BIE formulations. The approxi-mate formulation is very convenient to use because the corner information is comprised naturally in the representations of those ap-proximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results canbe achieved in comparison with those of the conventional BIE formulations.