In this paper, the solution domain is divided into muti-patches on which B-Spline basis functions are used for approximation. The different B-Spline patches are connected by a transition region which is described by B-Spline basis functions. The basis functions in different patches are modified in the transition region to ensure the basic polynomial reconstruction condition and the compatibility of displacements and their gradients. The results of numerical examples show that the present method is accurate and stable.