This paper is concerned with a nonstationary Markovian chain of cascading damage that constitutes an iterated version of a classical damage model. The main problem under study is to determine the exact distribution of the total outcome of this process when the cascade of damages finally stops. Two different applications are discussed, namely the final size for a wide class of SIR (susceptible → infective → removed) epidemic models and the total number of failures for a system of components in reliability. The starting point of our analysis is the recent work of Lefèvre (2007) on a first-crossing problem for the cumulated partial sums of independent parametric distributions, possibly nonstationary but stable by convolution. A key mathematical tool is provided by a nonstandard family of remarkable polynomials, called the generalised Abel–Gontcharoff polynomials. Somewhat surprisingly, the approach followed will allow us to relax some model assumptions usually made in epidemic theory and reliability. To close, approximation by a branching process is also investigated to a certain extent.

In this paper we consider the problem of first-crossing from above for a partial sums process m+St, t ≥ 1, with the diagonal line when the random variables Xt, t ≥ 1, are independent but satisfying nonstationary laws. Specifically, the distributions of all the Xts belong to a common parametric family of arithmetic distributions, and this family of laws is assumed to be stable by convolution. The key result is that the first-crossing time distribution and the associated ballot-type formula rely on an underlying polynomial structure, called the generalized Abel-Gontcharoff structure. In practice, this property advantageously provides simple and efficient recursions for the numerical evaluation of the probabilities of interest. Several applications are then presented, for constant and variable parameters.