Numerical study for cohesive zone model in delamination analysis based on higher-order B-spline functions