NURBS-based isogeometric analysis (IGA) has currently been applied as a new numerical method in a considerable range of engineering problems. Due to non-interpolatory characteristic of NURBS basis functions, the properties of Kronecker Delta are not satisfied in IGA, and as a consequence, the imposition of essential boundary condition needs special treatment. The main contribution of this study...