Explicit finite difference code was used to calculate the stability factors of shallow tunnels without internal support in limit state. The proposed method was formulated within the nonassociative plasticity. For the shallow tunnels in soft clay, without considering ...