A penalty function based delamination model for postbuckling analysis of composite plates with delamination