  This paper studies the convergence difficulty of cohesive zone models in
static analysis. It is shown that an inappropriate starting point of iterations
in the Newton-Raphson method is responsible for the convergence difficulty. A
simple, innovative approach is then proposed to overcome the convergence issue.
The technique is robust, simple to implement in a finite element framework,
does not compromise the accuracy of analysis, and provides fast convergence.
The paper explains the implementation algorithm in detail and presents three
benchmark examples. It is concluded that the method is computationally
efficient, has a general application, and outperforms the existing methods.
