Influence of fracture criteria on dynamic fracture propagation in a discrete chain

The extent to which time-dependent fracture criteria affect the dynamic behavior of fracture in a discrete structure is discussed in this work. The simplest case of a semi-infinite isotropic chain of oscillators has been studied. Two history-dependent criteria are compared to the classical one of threshold elongation for linear bonds. The results show that steady-state regimes can be reached in the low subsonic crack speed range where it is impossible according to the classical criterion. Repercussions in terms of load and crack opening versus velocity are explained in detail. A strong qualitative influence of history-dependent criteria is observed at low subsonic crack velocities, especially in relation to achievable steady-state propagation regimes.


Introduction
Studying fracture propagating in discrete structures results in a tool capable of analyzing a broad range of phenomena which would not emerge in the settings of continuum mechanics. This approach has found fruitful applications when dealing with crystals, cellular materials, cracks in fiber-reinforced matrices and investigations at the atomic level (e.g. Abraham et al. 1998;Atrash and Sherman 2011;Bolander and Sukumar 2005;Tsai et al. 2010;Glassmaker et al. 2007). Lattice structure models become even more crucial in the framework of dynamic propagation. In this respect, many kinds of instabilities can be predicted by intuitive considerations without the need of ad hoc hypotheses (see Marder and Gross 1995;Bernstein and Hess 2003;Kessler and Levine 2003;Bitzek et al. 2015). Finally, discrete models can be also treated as the discretization of the corresponding continuum problems where a choice of the fracture criterion may play an important role when dynamic fracture propagation is in question. For example it may lead to different predictions on the stability of possible steadystate regimes (e.g. Morini et al. 2013;Mishuris et al. 2012).
A lattice structure, in the dynamic scenario, is composed of concentrated masses interacting via links characterized by an interaction potential. In this paper the latter is parabolic, relating only the closest neighbors, while non-local interactions have been studied in Gor-bushin and Mishuris (2016a). The analyzed structure is mono-dimensional: a chain of oscillators, which is detached from a substrate that reflects the problem symmetry.
The focus of this work is to investigate the influence of the fracture criteria of the links on the dynamic fracture propagation in such medium. Exclusively, cracks which advance at constant speed, in a steadystate, are analyzed. Such regimes, indeed, have traditionally been of extreme interest in the field of dynamic fracture and have repeatedly been observed experimentally. A few classical studies can be highlighted for example in Ravi-Chandar and Knauss (1984), Fineberg et al. (1991) and Hauch et al. (1999) while the topic gained new attention more recently in discrete structures such as bridges (Brun et al. 2014) or xyloexplosives (King https://www.youtube. com/watch?v=jiWxU3jXOFc).
Before addressing the problem of the propagation, though, the behavior of a single link is discussed hereafter. A linear elastic spring can be quasi-statically elongated until failure and its final elongation value, u s , is assumed to be a known constant. The simplest and most common failure criterion neglects dynamic effects on the spring resistance. It identifies the displacement u as critical when min t : u(t) = u s , by which the time t = t f when the fracture occurs can be found. A fracture event, though, in many materials turns out to be not simply determined by an instantaneous threshold value for some energy measure like the maximum elongation established above. We deal in the present work with non-instantaneous fracture criteria which nevertheless do not change the material stiffness. The rate with which a body is deformed or an integral measure of the deformation energy provided to a bond before it breaks are examined: the incubation time (IT) and the Tuler-Butcher (TB) criteria.
According to the first of those formulations, the average stress, or equivalently the average linear elastic stretch, over a period of time preceding the breakage is considered the cause of fracture. Such a period is actually called the incubation time τ . The criterion, originally formulated in terms of stresses in Petrov and Utkin (1989), can be written for the elastic bond in object in terms of the elongation as min t : Notice that u s is still the threshold elongation of the spring when measured statically. The main idea of this approach is that a transient process does not occur instantly but includes a more complicated breakage process. Its realization demands a time period depending on the intensity and shape of loading and on the internal structure of the fractured media. The introduction of the incubation time τ as strength parameter in addition to u s allows one to predict the stress level at the instant of fracture for a variety of loading pulses with different intensities and shapes (Volkov et al. 2017). So the physical meaning of the fracture incubation time is a characteristic time determining the material's ability to resist to dynamic loading. The IT approach has shown to be reliable in different branches of mechanics and physics, such as the dynamic fracture of rocks and concretes, the dynamic yielding of metals, the acoustic ultrasonic cavitation of liquids, etc (Petrov 2004;Gruzdkov and Petrov 2008). Particularly, this criterion was successfully applied to problems of fracture in materials with pre-existing cracks. In this case the criterion can be reformulated in terms of stress intensity factor. The steady-state fracture propagation in a lattice structure is indeed crack growth and this analogy allows the application of the criterion to the problem considered here. On the other hand, it has been observed that cumulative damage can also be the cause of fracture. A way to quantify this is via the Tuler-Butcher criterion discussed in Tuler and Butcher (1968). Again, one linear elastic bond is statically tested until it breaks at u = u s . In a TB material, u s represents the elongation to be exceeded in dynamics as before the fracture occurs at t = t f . Here, H(u − u s ) is the Heaviside step function by which it is possible to write that only the work of the overstretch u − u s contributes to damage. Note that also this criterion was originally formulated in terms of stresses and that the exponent two was left general in the original formulation, but it turns out to be such in most experiments. In this way the physical meaning of the criterion is that a maximum work has to be done by an external overload on the spring before it collapses. Looking at Eq. (3), it turns out that, as for the IT criterion, TB materials can be regarded as one possible extension of ideal brittleness. The latter can be retrieved indeed by setting the cumulated energetic damage D to zero. The criterion has found fruitful applications in analyzing fracture in rocks, glass, aluminum and copper caused by spallation, impact loading and thermal shock. (see Boustie and Cottet 1991;Wei et al. 2009;Grady 2010). In order to illustrate the peculiarities of the criteria, one can imagine applying a ramp displacement of rate r at t = 0 to three links with the same static strength u s but different failure behavior. The ideally brittle one would break at t f = u s /r as soon as its elongation reaches u s . The IT spring would break, accordingly to Eq. (2), at t f = u s /r + τ/2 if r < 2u s /τ or at t f = √ 2u s τ/r alternatively, thus establishing a distinction between low and high deformation rates. The criterion shows a delayed failure causing an ultimate elongation bigger than the static one in this loading condition. In the case of a non-monotonic load, though, such delay might result in an elongation at failure which is smaller than the static one or during unloading (see Volkov et al. 2017). The TB criterion also predicts a delayed failure, but at t f = u s /r +u s 3 3D 2 /r 2 according to (3). The difference with the IT case is that, now, an oscillating load which is strong enough to break the spring in statics will also do it in dynamics. Notice in fact that a constantly increasing cumulated damage would sooner or later surpass D [see left-hand side in Eq. (3)].
The rest of the paper is devoted to the model of a fracture in a structured medium subjected to the aforementioned criteria and their effects on the stable regimes of propagation. The results are expressed in terms of trapped lattice energy, applied remote force and crack tip opening.

Background
Consider an infinite number of masses M linked to each other and to a rigid substrate via linear elastic springs of length 1 and stiffness k. A force F, applied at infinite distance, introduces energy into the system and finally breaks the links causing the crack to propagate to the right in Fig. 1. If x n (t) is the link that fails at the time t, the equilibrium for a generic oscillator at x = x i in terms of its displacement u i (t) holds as where the Heaviside step function H (i − n) allows for the combination of the equation for the detached (i < n) and intact part (i ≥ n) of the chain.
Only fracture with a constant speed v, i.e. steadystate fracture, is analyzed here. In such a way the problem is reduced to the long known settings of Slepyan and Troyankina (1984) and Slepyan (2012). The fracture can travel slower than sound in the broken structure, that is As a result of the steady-state assumption we search for a solution in terms of the unknown function for any i and t > 0. Further on we adopt a coordinate system which moves together with the crack tip in a way that the crack tip sits conveniently always at η = 0. By using such new moving frame, the coordinate η accounts for time and position simultaneously. Thus, the equation of motion Eq. (4) for u(η) can be written in the broken (η < 0) and intact (η ≥ 0) sides of the tip as: With the help of the mathematical tools of Fourier transform and Wiener-Hopf technique, such an equation has been repeatedly solved for this and more complex structures, for example in Slepyan et al. (1999), Marder and Gross (1995), Kresse and Truskinovsky (2004), Nieves et al. (2016) and Gorbushin and Mishuris (2016a), such to give the displacement profile u(η) which travels along the structure at a given steady-state crack speed v. If one intends to describe the trajectory of a single mass during time, one has to just apply Eq. (6).
The concept of energetic lattice trapping of a structured material was originally introduced in Thomson et al. (1971) and it can be quantified by the ratio G/G 0 not smaller than unity. Such energy may be introduced into the system in different ways. Analytical relations for the energy release rate G for every crack velocity have been retrieved in Slepyan and Troyankina (1984). The quantity is the link strain energy which is released locally at the crack tip at the moment of fracture, where u 0 = u(0). For the examined chain, G 0 /G has been plotted in Fig. 2a. In our work we assume that the energy derives from a constant force F far away and the consequent crack speed has been recently derived in Gorbushin and Mishuris (2017) as and plotted in Fig. 2b. Note that the same curve also implies that the limiting velocity v c can not be reached via a finite force besides requiring an infinite energy release rate (as from Fig. 2a).
The assumption that the crack propagates at a constant speed also requires some additional consideration. In particular it means, for a given oscillator i sitting at x i , that it is not allowed to break before all the links situated at x < x i do (links on the left-hand side of be clarified that a regime which involves nucleation of daughter cracks ahead of the mother crack tip (η > 0) is non-admissible. The detachment of the chain has to progress continuously. We shall discuss in the next sections how drastically the failure criteria change the admissible scenarios of stable detachment velocity.

Ideally brittle links
If the links are ideally brittle, the critical condition to be reached at the crack tip at the instant of fracture before further propagation is independent of the fracture propagation speed. With such a condition, that is when u 0 and u s are interchangeable, the diagrams of force and energy release rate in Fig. 2 are directly applicable. Furthermore, in this case, the condition of admissibility is easily checked i.e. that no points for η > 0 are lifted higher than the crack tip. Looking at Fig. 3, one can notice that for such materials, configurations occurring at low v are unphysical since there are points ahead of the crack tip where the failure criterion has been encountered already before the arrival of the fracture front itself and thus must be labelled as not admissible. Discussions on the matter have been dealt in Marder and Gross (1995), Mishuris (2016a, 2017). For example, the speed 0.2v c does not fulfill such requirement, then this must be discarded as non-admissible. On the contrary, the speeds 0.3v c and 0.47v c are admissible. Observing all the u(η) profiles, for the isotropic chain the minimum velocity of the crack corresponds to about 0.27v c and all larger subsonic velocities are admissible (see Fig. 4). It is perhaps worth pointing out that such a limit is smaller than the minimum energy release rate, which sits at about 0.38v c . This implies that a single G may correspond to two possible steady-states like 0.3v c and 0.47v c . Such speeds, anyway, correspond to two different loads (see Fig. 2b). The highest of the two speeds is achieved uniquely by means of a larger force.

Problem and methods
When dealing with non-instantaneous criteria for fracture, the crack opening before fracture depends on the crack speed. In general, must be determined accordingly to the new fracture parameters and does not simply equal u s as for the ideal brittle criterion Eq. (11). In a similar way, also the energy released locally at the crack tip G 0 (v) = ku 2 0 (v)/2 is a function of the crack speed. Further on we express F as a multiple of the force required to break the spring in a static test via the function which differs from the general Eq. (10) where the denominator incorporates the elongation at failure, independently on the particular fracture criterion adopted. The reason is conceptual and follows from the possibility of conducting experiments. For obtaining the same crack speed, the loading condition, indeed, must be accurately designed depending on how the failure happens (i.e. within the context of this paper, which fracture criterion better describes the constituent material). We stick to the ratio G/G 0 , instead of introducing a hypothetical G s = ku 2 s /2, because the energy release rate incorporates the type of the structure and its deformation properties (in the present case a linear elastic chain) and can hardly be measured. Moreover, in this way, as we will discuss further in the continuation, the dependency of G/G 0 remains untouched by the particular fracture criterion characterizing the links, while the force versus velocity relation depends on the particular criterion.

Incubation time criterion
The incubation time failure criterion Eq. (2) can be applied to the chain by a change of variable according to Eq. (7) which leads to The normalization by the displacement at the crack tip u 0 is convenient because, in the ideally brittle case, the crack opening u 0 before fracture was known in every case and given by the maximum elongation criterion, now it is unknown and dependent on velocity.
The shape u(η, v)/u 0 (v) of the deformation profile, though, is given once and for all as it does not depend on the particular value of the crack opening. The advantage is that, once one calculates the shape at a certain velocity from the solution of Slepyan and Troyankina (1984), this can be used for all the possible steady-state fracture criteria. If τ goes to zero, that is, the material is ideally brittle, Eq. (14) returns u 0 = u s coherently. By this respect, one can say that IT materials are an extension of ideally brittle ones by means of τ . Moreover, for steady-state propagation, the length vτ is constant in time and thus incubation time can be considered as a non-local criterion as well as a non-instantaneous one.

Tuler-Butcher criterion
In order to deal with the usual moving coordinate frame Eq. (7) and a steady-state regime of velocity v, the TB criterion Eq. (3) can be transformed into the equation Taking advantage of the invariance of u(η, v)/u 0 (v) with respect to η, the function can be obtained as solution of Eq. (15).

Results
In order to have dimensionless strength parameters, from this section on we express u s in units of the distance between the masses, whereas D and τ are expressed in units of the same distance divided by the sound velocity v c .

Incubation time criterion
The steady-state analogue of the incubation time criterion in Eq. (14), which also defines the function Ψ (v, τ ), solves the issue of calculating the crack opening. Given τ and u s , .
The force to apply for achieving a certain steadystate velocity is obtained by substituting Eq. (17) into Eq. (13): The behavior of the function Ψ (v, τ ) influences the way τ modifies the crack opening with respect to the ideally brittle one (τ → 0) and this is shown in Fig. 5a. The crack speed dependent opening u(0) = u(0, v) expressed through Ψ (v, τ) changes the prediction of the force to apply in order to cause a certain velocity A linear elastic bond which exhibits a non-zero incubation time will in general allow a bigger crack opening at the instant of fracture. Crack openings smaller than u s , though, are admissible at low velocities due to rapid oscillations and negative ∂u/∂t close to the tip (see Fig. 6). The influence of τ on the force is plotted in Fig. 5. Given the result of the static test on the spring u s , if the goal is achieving a certain velocity v, an IT type material predicts that the steady-state regime would be reached in general via a bigger force than one could expect if τ is neglected. The region where the relation between force and velocity is not bijective is stretched to the right and the difference in velocities for the same force decreases steadily while raising the Fig. 6 Crack opening for small velocities and high incubation times. With τ = 5, the crack tip opening u 0 for the admissible speed v = 0.14v c is predicted to be smaller than the static strength u s . For comparison, v = 0.2v c shows the most common situation of u 0 > u s . See also Fig. 5a incubation time. Note that τ does not play any role in the limiting case of zero velocity of propagation, where the force produces the same results as for the ideally brittle spring (see Eq. (14) for v → 0). In order to verify the analytical solution, we searched for the steady-state regimes via solving Eq. (8) by a finite difference scheme similarly as in Gorbushin and Mishuris (2017) where the same numerical procedure is extensively explained. A chain of 2000 masses was loaded with a distant vertical constant force. The instant of fracture for a link was identified according to the condition Eq. (2). After iterating, the next failures tended to occur at constant time interval and such an interval was used for calculating the stable crack speed for a given force. The analytical solution is perfectly matched and the numerical approach confirms that a steady-state propagation is not achievable at low velocities (results not shown).
This kind of computation is quite heavy because of the algorithm adopted to identify the time of fracture at every location x i . Before a steady fracture propagation is reached, indeed, during the transient regime, the history u i (t) in the last interval τ must be recorded and the integral Eq. (3) updated for every x i . This marks a principal difference with the instantaneous traditional criterion in which case the quick check u ≤ u s is sufficient. We have tried to simplify the procedure via making the criterion pseudo-instantaneous. Theoretically, Eq. (17) should be valid only for steady-state fracture.
Nevertheless, if one estimates the instantaneous crack velocityṽ(t) from the last two failures, the criterion can be artificially reduced to u ≤ u 0 (ṽ, τ ) instead of calculating the integral Eq. (17) at all. Such an attempt has been proven to be effective beyond expectations for the particular studied problem in achieving the same steady-states as in theory and the rigorous numerical simulation for same applied load.
Another crucial effect of τ > 0 on the crack propagation is that it monotonically enlarges the regions of achievable steady-states as illustrated in Fig. 4. For instance, the speed 0.2v c , that is non-admissible for an ideally brittle material with critical elongation parameter u s , can be reached with τ = 3 and bigger.

Tuler-Butcher criterion
By means of Eq. (16), one can retrieve the crack opening associated to all the combinations of crack speed and D.
As consequence, keeping the force proportional to the material property u s instead of the velocity dependent crack opening, Eq. (13) becomes .
The plots in Fig. 7 allow us to visualize how the dynamic strength parameter D affects the chain behavior. As observed for IT materials, the immediate impact of TB damage accumulation results in an augmented crack opening at equal crack speed as an ideally brittle material showing the same static strength u s , at least in the range of medium or high v/v c (see Fig. 7a). The force needed for obtaining a desired velocity is depicted in Fig. 7b. It is evident that the capability of the material to bear a certain work of the overstretch before failure, i.e. bigger D, makes the chain detachment increasingly slower for same F/ku s . A structure of TB bonds can be predicted to be dynamically tougher than its ideally brittle counterpart. For D → 0 and low fracture speed the complex structure does not respond like an ideally brittle one. Looking at Fig. 8 The crack speed dependent opening u 0 expressed through Λ(v, D) changes the prediction of the force to apply in order to cause a certain velocity decrease monotonically ahead of the crack tip, the criterion in the form Eq. (15) may return u 0 < u s . Such a feat differs from all the cases analyzed in this work: single bonds as well as the ideally brittle and IT crack tips always fail for u ≥ u s . However, in the cases where the latter inequality is violated in a TB link, we obtain nonadmissible steady-state propagation regimes. A point can be made, therefore, that a theoretical limit for the crack opening is for a TB material to satisfy u 0 ≥ u s . Like for incubation time materials, as it can be seen in Fig. 4, new zones of admissibility appear in the low velocity region for larger D. Nevertheless, there is a significant qualitative difference with the incubation time situation: such admissible intervals emerge small and scattered, but then, with increasing D, expand gradually and merge until every subsonic crack speed can be obtained for D close to unity.

Discussion and conclusions
The dynamic fracture propagation in discrete structures has been investigated in a considerable number of possible scenarios (see references above) but the influence of failure criteria different from a threshold stress has not gained the attention that it deserves despite the fact that non-instantaneous criteria have already been shown to be reliable in continuum mechanics [e.g. as recently discussed in Alves and Lobo (2017)]. As a first step to fill this gap, two time-dependent criteria have been analyzed in detail when applied to the dynamic fracture propagation of a chain of oscillators and they have been compared to the classical ideal brittle fracture. In both cases, enhanced admissibility have appeared at low crack speeds and mapped in Fig. 4. An increasing incubation time τ enlarges the admissibility continuously but never covers all the subsonic crack speeds. More than that, a TB material which requires a larger work of the overstretch for fracture and characterized by a bigger D also creates completely new zones of achievable steady-states and it is predicted that all the subsonic range is possible if D 0.13.
Speaking of the steady-state crack opening, the timedependent criteria cause a delay in fracture after reaching the static strength of the bonds. This means that in most cases one should expect u 0 > u s like it would happen when monotonically elongating a single spring. At low v, though, this is not the behavior caused by a constant force applied on a complex structure. Ample and rapid oscillations ahead of the crack tip cause the delayed fracture to happen at u 0 < u s . While such propagation regimes are admissible at high τ for incubation materials, the same is not true for TB ones (see Figs. 6,7 and 8). The mathematical form of the latter failure criterion indeed excludes such steady-states on the grounds that daughter cracks would jeopardize the steady-state assumption. In short, a theoretical limit has been found which states that for a TB chain a dynamic fracture can propagate at constant speed only if u 0 ≥ u s .
As intuition suggests, the two examined non-instantaneous criteria make the structure tougher than the corresponding ideal brittle one with the same static strength u s . Thus, for obtaining a certain velocity v one needs a bigger force if τ or D increase. In this way, the curves in Figs. 5b, 6b and 7b also result in an important practical application. With a few experiments on materials whose u s has been independently measured, the couples F-v allow for the material characterization in terms of the second fracture parameter τ or D at least if stable crack speeds are retrieved in the monotonic interval of the curves (medium and high v/v c ).
The velocity dependent energy release rate ratio G/G 0 is a solution which is irrespective of the particular fracture criterion adopted. It is also valid regardless of the way the energy is introduced into the system. In the present work we use a constant force as the external load. However, in case one prefers to implement different kind of loading, like for instance in Marder and Gross (1995) when dealing with lattices, or for example to facilitate a specific experimental procedure, the relation between the new load and the crack speed has to be evaluated, while the energy-crack speed diagram remains the same.
Two numerical integration schemes have been used to solve the set of governing equations in the IT scenario and compared with the analytical solution derived in the present paper. The first of these verified the criterion condition in its integral form Eq. (2) at every time step for all the unbroken bonds. Such an approach enables one to simulate also the transient regimes before a steady-state is reached. The steady-states were achieved only in the admissible regions of Fig. 4 and there the force-velocity relations agreed perfectly with the ones in Fig. 5. The second test was performed by adopting a pseudo-static failure criterion: a dynamic threshold elongation u 0 (ṽ, τ ) was established based on the instantaneous crack velocityṽ and Eq. (17). This simplified algorithm performed much faster than the first one, and still returned correct results. For this specific structure at least then, it seems that many of the conclusions drawn here can still be valid in the transient regimes.
Beyond the particular scope of this paper, various propagation regimes, in absence of crack arrest, can appear: steady-state, other regular ones (clustering or forerunning as discussed in Mishuris et al. (2008Mishuris et al. ( , 2009, Slepyan et al. (2010Slepyan et al. ( , 2015, Nieves et al. (2016Nieves et al. ( , 2017 or chaotic regimes. The realization of one or another heavily depends on the loading type, its intensity and on the structure itself. However, when the problem is faced from a mathematical point of view, assumptions of steady-state regimes have always been made in order to obtain simple solutions. With the analytical results in one hand, an a posteriori examination is required which identifies where the solution fulfills the assumptions and constraints: only that part of the solution is labeled as admissible. Generally speaking, though, it does not mean that all those regimes will emerge in practice as steady-state. In most of the cases it happens; nevertheless, as it has been shown in numerical simulations, other ordered regimes of propagation may arise such as clustering. In such circumstances, the steady-state velocity predicted theoretically reveals as the average speed at which the cluster moves (see Mishuris et al. (2009), Nieves et al. (2016 on the matter).
We have not treated the problem of branching in the present settings of history-dependent criteria. It has turned out already in Marder and Gross (1995), Mishuris (2016b, 2017) that such instabilities can become relevant at high crack speeds. In the considered geometry, loading condition and material parameters, the branching mostly happens along the crack surfaces but not on the crack line ahead. That is evident from the solution profile u(η) for v/v c = 0.2 in Fig. 3. If the horizontal springs show the same dynamic resistance as the vertical ones, their fracture can pre-cede the chain detachment from the substrate, making a steady-state propagation impossible and sensibly reducing the limiting speed with respect to v c . The admissibility check would imply that, for none of the consecutive oscillators, the difference u i+1 − u i ) does reach the condition imposed by the fracture criterion. More complex scenarios may occur with structures characterized by flexural stiffness, heterogeneities or localized feeding waves (Mishuris et al. 2009;Nieves et al. 2016Nieves et al. , 2017. Furthermore, a steady-state regime can be unattainable, resulting in unstable or alternating velocities, when the structure is not loaded far from the crack tip, but via accumulated energy in the form of residual stresses of the bonds (Ayzenberg-Stepanenko et al. 2014). Complications have also been the object of investigation in the framework of bridged cracks (Mishuris et al. 2008).
In conclusion, the fracture criteria considered here sensibly affect the dynamic propagation of cracks in discrete structures. The effects are particularly important both in terms of force vs velocity relations and in new regimes of admissibility at low crack speeds. Succinctly, by the present results, we are able to underline two main messages: • it is at low speed regimes that an experimental investigation should be carried out more carefully for understanding whether it is necessary to incorporate history-dependent fracture criteria in the dynamic fracture model; • the energy release rate ratio and shapes of the displacement profiles as functions of the velocity are invariants, in linear theory, and can promptly be used and adapted to the most suitable fracture criterion for the analyzed problem.
The possible outlook of this research is the application of the approach to (a) more complex lattice structures (inhomogeneous, triangle ones) as in Mishuris et al. (2007), Nieves et al. (2013); (b) highly ordered bidimensional lattices, for instance to crack propagation in graphene layers (Xiao et al. 2009;Tsai et al. 2010) or (c) to the unbinding of long protein chains whose analysis has been made feasible by the improvements in the field of atomic force microscopy and for which the bonds strength has already shown to be eminently dependent on the strain rate (Merkel et al. 1999;Hinterdorfer and Dufrêne 2006). The model considered here can be useful for modeling the peel test of flexible films (see e.g. Wei and Hutchinson 1998;Nase et al. 2016).