Delamination is one the most important defect that reduces strength of part. Many researchers have been studied delamination in drilling of composite materials and they tried to formulate this phenomenon by developing of analytical, numerical and experimental models. In this models up to now, only effect of modeΙ and modeΙΙΙ of crack propagation is considered and effect of modeΙΙ is neglected. The goal of this research is determination of modeΙΙ effect on propagation of interlaminar cracks during drilling of multilayered composites by means of finite element analysis (FEA).Thus the numerical analysis of delamination of unidirectional Carbon/epoxy composite during drilling is performed by modification of previous numerical models. The numerical method which is used for determination of strain energy release rates in modeΙ and modeΙΙ is Virtual Crack Closure Technique (VCCT). This analysis is performed for crack propagation under chisel edge when the drill have not come out from the workpiece and for crack propagation under cutting edges when the drill have come out from the workpiece. By determination of strain energy release rates in modeΙ and modeΙΙ and comparing with critical values G_Ιc andG_ΙΙc, the critical thrust force that causes to delamination is determined and contribution of each crack propagation mode in delamination is discussed. At the end of this research, it was found that the effect of modeΙ is more than modeΙΙ in all of cases such that more than 95percent of crack propagation parameter in power law criteria is due to strain energy release rates in modeΙ.