Dissipation During Crack Growth in a Viscoelastic Material from a Cohesive Model for a Finite Specimen

In the present paper, we extend results recently given by Ciavarella, Zhang & McMeeking (2022, Journal of the Mechanics and Physics of Solids, 169, 105096) to show some actual calculations of the viscoelastic dissipation in a crack propagation at constant speed in a finite size specimen. It is usually believed that the cohesive models introduced by Knauss and Schapery and the dissipation-based theories introduced by de Gennes and Persson-Brener give very similar results for steady state crack propagation in viscoelastic materials, where usually only the asymptotic singular field is used for the stress. We show however that dissipation and the energy balance never reach a steady state, and we are therefore unable to use the de Gennes and Persson-Brener theories which suggested that the increase of effective fracture energy would go up to the ratio of instantaneous to relaxed modulus, at very fast rates. We show viscoelastic dissipation is in general a transient quantity, which can vary by orders of magnitude while the stress intensity factor is kept constant. Also, at intermediate rates dissipation can be orders of magnitude high er than work of fracture multiplied what is believed to be the “viscoelastic enhancement factor” at very large rates . Finally, the total work to break a specimen apart is found, in a realistic example, to be larger than this “limit”, and for quite a large range of realistic conditions at intermediate crack growth rates. This shows that the cohesive model of crack propagation in linear viscoelastic materials permits a more general understanding than models which assume steady state of energy fluxes and simplified energy balance of just the asymptotic singular stress field.


Introduction
Crack propagation in viscoelastic materials is rate dependent, and the load to induce crack growth typically increases with the rate of crack propagation (Mueller and Knauss 1971).Gent and Lai (1994) study the rupture of elastomer sheets adhered together and show that if we vary by 9 orders of magnitude the peel rate, the apparent toughness increases by a factor of 100.Gent (1996) summarizes results for tearing, peeling and cutting of 3 elastomers and the force to spans 5 orders of magnitude as the rate of rupture is increased, which is perhaps even higher than the increase of the elastic modulus between slow and fast rates.Master curves are typically used using the temperature-rate superposition principle, which confirms viscoelastic effects are the cause of the rate dependence of rupture propagation that is observed, although we cannot be sure if these viscoelastic effects are in the intrinsic bond scission near the crack tip, or dissipation in the bulk.Experiments of detachment of polymers from a substrate (Gent and Petrich (1969), Gent and Schultz (1972), Ondarçuhu (1997)) exhibit rate phenomena similar to those observed for cohesive fracture (Mueller and Knauss (1971), Gent andLai, (1994), Gent, (1996)).Kostrov and Nikitin (1970), McCartney (1977McCartney ( , 1980McCartney ( , 1981)), Christensen (1979Christensen ( , 1980Christensen ( , 1981) ) and Christensen and McCartney (1983), Graham (1969), Barenblatt, Entov and Salganik (1970), Rice (1978), Goleniewski (1990), Hui, Zhu and Long (2022) and others all make it clear that when there is no rupture process zone in the model for crack growth in a linear, viscoelastic material with a rate independent crack tip toughness, the applied Mode I stress intensity factor,   , is required to be equal to a critical value,   , where   is consistent with The other model of fracture involves adding viscoelastic energy dissipation that occurs during crack growth, to explain the increase of the apparent toughness measured during crack growth (de Gennes 1996, Saulnier et al. 2004, Persson and Brener 2005, Persson 2017, 2021).In these models, the crack tip process zone is included via a limiting, or cut-off stress at the crack tip.This second class of models seem to generate similar conclusions that the cohesive models, at least for a semi-infinite crack with a process zone that is very small compared to all component dimensions (Ciavarella, Cricrì andMcMeeking, 2021, Hui et al. (2022)).However, the issue is more controversial in the case of finite component size effects (Ciavarella and Papangelo, (2021), Persson, (2017)).It would seem that there is a limit to the possible size of the region of dissipation and Persson 2017 predicts that, at an intermediate crack propagation speed, the toughness will go through a maximum, and at higher rates of crack propagation the toughness will either decline or become insensitive to further increases in crack speed.This outcome has not yet been firmly established in treatments of crack growth in finite sized viscoelastic components.Indeed, in a solution for the double cantilever beam (Wang et al. (2016), Ciavarella (2021)), when the beams are thin in the bending direction there seems to be no reduction of the fracture load compared to the trend for a thick beam, contrary to the prediction of the de Gennes (1996) and Persson (2017) dissipation-toughening models.Also, for a finite specimen we have to consider that there are ample regions of nominal stress which can dissipate energy, which are not taken into account in the de Gennes (1996) and Persson (2017) theories.It is not clear if this dissipation should also be included in addition to the intrinsic fracture energy of the material.Hence, while it is demonstrated by Hui et al. (2022) the substantial equivalence of the Knauss (1973), Schapery (1975aSchapery ( , 1975bSchapery ( , 1975c) ) and Persson and Brener (2005) models for the highly idealized semi-infinite crack under constant applied stress intensify factor: we will see that practical realization of this model is not obvious.
In the Knauss (1973), Schapery (1975aSchapery ( , 1975bSchapery ( , 1975c) ) and Persson and Brener (2005) models of crack growth in viscoelastic materials, Hui et al. (2022) identify what they term the fraction of energy dissipated per unit area of crack extension, and attribute to it the role of toughening of the material during crack growth.They give it the symbol Λ.From their Eq.(8b), it enters the plane strain relationship where  is the length of the cohesive zone at the crack tip,   is the creep compliance retardation time, i.e., the time constant for viscoelastic deformation in a creep test, and   is the Young's modulus of the material in the relaxed state.The term outside of the square brackets on the right-hand side of Eq. ( 2) is the elastic energy release rate to the crack tip when the component behaves as an elastic one with Young's modulus equal to its value in the relaxed state of viscoelastic behaviour.Indeed, it is often observed that, during fracture tests and adhesion detachment tests, the overall behaviour of viscoelastic components is consistent with the relaxed state of the material (Gent and Kinloch (1971), Mueller and Knauss (1971), Knauss (1973), Maugis and Barquins (1980), Knauss (2015)).Hui et al. (2022) assert that when the rupturing viscoelastic component overall is in the relaxed state, Λ (1 −  2 )  2   ⁄ is the energy per unit area of crack advance that is dissipated so that the total crack tip energy release rate, (1 −  2 )  2   ⁄ , associated with the applied loads, is the sum of the dissipated fraction, Λ (1 −  2 )  2   ⁄ , and the intrinsic toughness, Γ.However, it is not obvious to us that Λ (1 −  2 )  2   ⁄ is dissipated in contrast to simply not being available for energy flux to the crack tip.In addition, we are of the view that (1 −  2 )  2   ⁄ is not always meaningful as a measure of the crack tip energy release rate associated with the applied loads.
To explore and illustrate this situation, and to address questions that we have raised, we obtain below some exact results from linear elastic fracture mechanics converted to linear viscoelasticity under constant applied stress intensity by use of the correspondence principle (Graham, 1968).Part of these results have been already illustrated in Ciavarella et al (2022), but are repeated for clarity, adding some actual calculation and graphical results which help clarifying the meaning of these results.From these results, we compute the energy stored and/or dissipated while crack growth is occurring.We show that a high toughness in a viscoelastic material at a high rate of crack propagation is not due solely to viscoelastic dissipation, and, in some cases, not due at all to viscoelastic dissipation.Further, we find possibly orders of magnitude higher work rates than expected from the literature so far, occurring for realistic elastomers at intermediate speeds rather than at very fast ones.For a Maxwell material, we obtain an explicit result for the rate of energy dissipation in a strip in tension with an edge crack, and show that, per unit area of crack advance, the dissipation rate falls as the rate of crack propagation is increased.As such, these results run counter to the idea that increased dissipation leads to increased viscoelastic toughness for crack growth.Whereas Hui et al. (2022) have shown that the Persson and Brener (2005) model is effectively one in which a rupture process zone is interacting with the viscoelastic properties of the material, just as in the case of the Knauss (1973) and Schapery (1975aSchapery ( , 1975bSchapery ( , 1975c) models, our results demonstrate that the Knauss (1973) and Schapery (1975aSchapery ( , 1975bSchapery ( , 1975c) ) models, that utilize a Barenblatt (1962) cohesive zone for rupture at the crack tip, cannot be viewed as ones in which viscoelastic dissipation causes all of the increase in applied loads required to drive cracks to grow at high rates of extension.

The crack tip and the energy release rate
In the absence of a model for the process zone at a crack tip, its stress field in a linear, isotropic, viscoelastic material in infinitesimal strain theory is the standard one of linear elastic fracture mechanics: in Mode I, where   is the stress tensor, (, ) is a cylindrical polar coordinate system with origin at the current crack tip and  ̃ () is a dimensionless tensorial function (Anderson, 1995).The correspondence principle (Lee, 1955;Graham, 1968) can be used to convert isotropic, linear elastic strain and displacement fields at the crack tip to those for isotropic linear viscoelastic materials.When a body with a Mode I crack is subjected to traction boundary conditions, we cam use the correspondence principle to analyse the strains and displacements throughout the body knowing the linear elastic solution for the body if the stress field is independent of the elastic moduli, as is the case for most plane strain problems (Graham 1968).
For the Dugdale cohesive zone model (Dugdale 1960), with the length of the cohesive zone very small compared to component dimensions, the cohesive zone length, c, is given by where   is the uniform, constant, rate independent, cohesive stress in the zone, and   is the magnitude of the stress intensity factor in the absence of the cohesive zone.In small scale cohesion, the stresses in the region of the crack tip are finite as stress singularities are suppressed.Following Knauss (1973) and Schapery (1975a), the tip opening displacement for a plane strain Mode I crack propagating steadily at the steady speed V is where we have assumed that Poisson's ratio, , is constant during viscoelastic response.The parameter   , the effective Young's modulus for the energy release rate, is given by In the integral in Eq. ( 6) the creep compliance function () gives the strain response of the material to uniaxial stress as and the function f, the elastic cohesive zone stretch for a homogeneous system, normalized to unity at the crack tip, is (Rice, 1968) Since the crack tip energy release rate is crack propagation occurs with     equal to Γ, the intrinsic crack tip toughness, thereby determining   , the critical value of   , so that We confine ourselves to materials in which the crack tip toughness, Γ, has a unique value independent of the rate of crack growth.By combination of Eq. ( 4) & ( 5) at the critical value of the crack tip opening displacement, plus Eq. ( 10) for the toughness, we deduce that the Irwin relationship between the crack tip toughness and the stress intensity factor is This relationship provides a means by which the rate of crack propagation, , can be predicted from the applied stress intensity factor,   .

Work done by the applied loads and viscoelastic dissipation
The correspondence principle of Graham (1968) can be used for the component shown in Fig. 1, namely a strip with a single edge crack subject to Mode I, uniform tension.The applied stress at the ends of the component shown in Fig. 1 is uniform and is   and the average displacement of one end of the specimen relative to the other end is ∆, while the specimen width is b, the crack length is a and the state of the component is plane strain.The cohesive zone length is assumed be very small throughout the process of crack growth, i.e.,   ⁄ ≪ 1,  ( − ) ⁄ ≪ 1 and  ℎ ⁄ < 1.

Figure 1.
The single edge cracked strip subject to uniform tensile stress.Tada, Paris and Irwin (2000), we obtain the elastic displacement of one end of the specimen relative to the other as

From
where  is Young's modulus of the linear elastic material.A good approximation for  is Tada et al. (2000) state that Eq. ( 13) is accurate to within 1% for any   ⁄ as long as ℎ  ⁄ ≥ 1.Because the cohesive zone is very small, it does not perturb the result in Eq. ( 13) to any significant extent.Furthermore, other than the small discrepancy introduced by the cohesive zone, Eq. ( 12) is exact, although obtaining an exact expression for (  ⁄ ) instead of the approximation in Eq. ( 13) is not a simple task.
We note that the elastic potential energy per unit thickness is As a consequence, the elastic crack tip energy release rate is where  ′ () = ()  ⁄ .From Eq. ( 15) and the elastic equivalent of Eq. ( 11) we deduce the stress intensity factor to be During steady state crack propagation with the cohesive zone model,   must be held constant.To achieve this, we treat   as a constant during crack growth and set This restriction may not be practical in the laboratory, but for our exploration we enforce it in our thought experiments.Therefore, Eq. ( 12) becomes Initially, the crack length is   and the component is free of applied load.At time zero, raise the applied stress from zero to the value in Eq. ( 16) with  =   .At time  =   the crack begins to grow at speed V. Below, we will identify this viscoelastic delay in the initiation of crack growth (Knauss (1970), Schapery (1975c)).For  ≥   the crack length is given by To set ideas in place, we first consider crack growth in a linear elastic component.The total work done on the specimen per unit thickness during crack growth is Since the system is elastic, the total energy per unit thickness is where (∆  , ) is the strain energy per unit thickness.Thus consistent with Eq. ( 20) and where we have equated the crack tip energy release rate and the fracture toughness.From Eq. ( 21) we have also Therefore, and as expected, the work done on the specimen per unit thickness, Eq. ( 23), is the sum of the change in strain energy and the increase in fracture energy.
Finally, we consider an elastic component that is unloaded to begin with, and is loaded to the critical stress intensity factor where Γ is the rate independent, constant fracture toughness of the material.At this stage the crack is permitted to grow from its initial length,   , to a final length,   , while   =   and then the component is unloaded to zero load.In this cycle of loading, crack growth and unloading, integration of Eq. ( 23) leads to as the initial and final strain energy in the component is zero.Thus, we have the well-known result that, for an elastic component, the work per unit thickness required in the cycle of loading, crack growth and unloading is just the work required to grow the crack.The result in Eq. ( 25) is, of course, independent of the magnitude of Young's modulus.

Extremely fast loading followed by crack growth
We now consider an initially unloaded component that is rapidly loaded to a specific value of   .However, in this treatment of rapid loading we neglect mass inertia so that dynamic effects are ignored.We define After loading is completed,   is held fixed, and so the result in Eq. ( 26) is then constant.Therefore, after rapid initial loading and during the time 0 ≤  ≤   , the viscoelastic crack opening for the nonpropagating crack, is We note that (0) = 1   ⁄ , and we choose   such that   (0) ≤   .The crack will begin to propagate when   () =   (  ) =   , thereby defining   , the delay time for crack growth.The result can also be written Note that if the right-hand side of Eq. ( 28) is equal to 1   ⁄ , then   = 0 and crack growth commences immediately after loading.We avoid loading above the level at which the right-hand side of Eq. ( 32) is equal to 1   ⁄ as that will lead to unstable brittle crack growth.If the right-hand side of Eq. ( 28) exceeds 1   ⁄ then crack growth will never commence as the critical condition   =   will never be reached and   will always remain below   .Therefore, we limit the applied stress intensity factor to but subject to the observation that at the lower limit   = ∞.
Note that once crack growth commences it immediately occurs at the rate determined by Eq. ( 11) as we are neglecting inertia and at  =   the cohesive zone has reached the configuration it will have for steady state crack growth.We note in passing that Willis (1967) considered a crack in a viscoelastic material in antiplane shear with mass inertia accounted for.
For the next step we define During the time 0 ≤  ≤   the crack is not propagating so that ∆ ̃() = ∆ ̃(  ) and is therefore constant.The viscoelastic deformation is After the crack begins to grow, Eq. ( 19) prevails, and This gives where  ′ () = ()  ⁄ .Therefore, the increment of work per unit area of crack advance is Analogous to Eq. ( 23), the increment of work, , is the sum of the work increments stored and dissipated in the material plus the increase in fracture energy, Γ, due to crack growth.

Work done during crack growth after rapid loading
We return to Eq. ( 34) and observe that the left-hand side is equal to the fracture work plus viscoelastic deformation work.Therefore, we rewrite Eq. ( 34 where    ⁄ is the increment of viscoelastic work that occurs in the specimen divided by the increment of crack length.Although this term is not purely dissipative, it contains the dissipated work as well as any work stored elastically.Just after  =   the crack begins to grow at its steady state rate.Therefore, immediately upon initiation of crack growth, Eq. ( 36 We note that the 2 nd term on the right-hand side depends on   but has no explicit dependence on the crack propagation speed.Furthermore, it does not depend on viscoelastic properties other than   and thus is not a dissipative term. To explore the result in Eq. ( 37), we consider the standard viscoelastic material for which and We now take advantage of an approximation introduced by Schapery (1975b) for the integration of Eq. ( 6), namely the use of  = 2 3 ⁄ in the integrand, justified by the typical relaxation spectrum exhibited by most viscoelastic materials.As a consequence, we obtain so that inversion of ( 3 ⁄ ) will allow the evaluation of the crack propagation rate, V, in term of the size of the cohesive zone and material parameters.We assume that Schapery's (1975b) approximation is valid for the standard material and obtain The result in Eq. ( 47) is such that  = 0 at the lower limit of (33) and  rises monotonically with   2 until it is infinite at the upper limit of (33).When the result from Eq. ( 42) is substituted into Eq.( 37), the term on the right-hand side of Eq. ( 37) that accounts for dissipation leads to This result is zero at each limit in (33), indicating the absence of dissipation at each of these limits, consistent with our previous observation.As   (  )∆ ̃(  ) is proportional to   2 , the quotient containing it is independent of   .As   2 increases from the lower limit of ( 33), the right-hand side of Eq. ( 43) increases, goes through a maximum and falls back to zero at the upper limit of (33).Thus, as  is increased from zero, the right-hand side of Eq. ( 37) rises from zero, goes through a maximum and then falls back to zero as  → ∞.The maximum value of the right-hand side of Eq. ( 43) occurs at As the ratio of   to   is typically 2 or 3 orders of magnitude, the right-hand side of Eq. ( 44) is well approximated by  ( − 1) ⁄ ≈ 1.58.Thus, the maximum of the right-hand side of Eq. ( 43) occurs at a value of   approximately 25% greater than the lower limit of (33).As the upper limit of ( 33) is typically at least 10 times greater than the lower limit of (33), the results show that the right-hand side of Eq. ( 43) rises with increasing   over a quite narrow range of   at the lower end of the range in (33) and falls with increasing   over almost all of the range in (33).
We note also that the term in square brackets in Eq. ( 43) is equal to the parameter Λ introduced by Hui et al. (2022) and featured in our Eq.( 2).As Λ rises monotonically with   within the limits of ( 33), the behaviour of the right-hand side of Eq. ( 43), rising and falling with increasing   , indicates that Λ is not a good measure of the fraction of energy dissipated per unit area of crack growth.We also observe that the right-hand side of Eq. ( 43) depends on ℎ, the half-length of the component.This reflects the fact that the rate of work done in a long specimen is greater than the rate of work done in a short specimen, simply because of the different volume of material that is being stretched viscoelastically by the applied stress.This feature does not factor out when the results are specialized to a short crack in the component shown in Fig. 1 for which the approximate result   = 1.122  √ can be used.The behaviour also indicates that the dissipation rate embedded in the right-hand side of Eq. ( 43) accounts for more than just the propagation of the crack, but also reflects the structural behaviour of a component that would be stretching even if the crack were not growing.
We should also point out that the right-hand side of Eq. ( 37) provides information about the component's viscoelastic behaviour only at the time of initiation of crack growth.It does not provide insights into the viscoelastic response of the component after some time has passed and the crack has grown beyond its initial length.Such information can be extracted from Eq. (35).As an example, we integrate Eq. ( 35) for a standard viscoelastic material having     ⁄ = 1000 for the component in Fig. 1 with ℎ   ⁄ =    ⁄ = 100, i.e. an initially short crack.We use the approximation in Eq. ( 47) to evaluate the rate of crack growth, , and plot results in Fig. 2 for the work rate,   ⁄ , normalized by the fracture energy Γ, from Eq. ( 35) for values of   √  Γ (1 −  2 ) ⁄ ⁄ = √1.58,16, 28.5 (black, blue, red curves).To assume reasonable values for the quantities involved, we notice that the maximum dissipation at the initial instant of propagation is 1.39(1 −  2 )  2 ℎ (    ) ⁄ which can also be written as , where  0 is the size of the cohesive zone at very low rates of crack growth and according to Hui et al.(2022) in many elastomers that have been tested, it is a very small quantity, ranging from 0.1nm to few nm.To obtain results consistent for elastomers with orders of magnitude suggested by Hui et al.(2022), we further assume Γ=50J/m 2 ,   =1MPa,   =1kPa,   =  ,   = 0.1, which leads to  0 =20nm.Fig. 2 shows that the work rate is at times close to the initial time t0 orders of magnitude larger than     ⁄ = 1000 for all the loads chosen, which cover almost the entire range (the load for   √  Γ (1 −  2 ) ⁄ ⁄ = 28.5 is 0.9 of the load needed for very fast crack propagation).In particular, the maximum work rate can be several orders of magnitude (in our case, almost 6 orders) larger than the expected "steady state" increase of   /  which in our case is a factor 1000.The reason for this important and unexpected result is that the factor (bh)/(  0  0 ) is easily a very large number.We have therefore found an "ultratough" regime of crack propagation, despite this occurs only in a transient.= √1.58,16, 28.5 (black, blue, red curves): solid lines refer to the standard material, while dashed lines refer to the equivalent Maxwell material with viscosity appropriate to obtain the same creep function of the standard material for not too long times, namely  =     .The load is suddenly applied at  = 0 and crack growth commences at  =   in each case where   has a different value in each case.The work rate,   ⁄ , is normalized by the fracture energy Γ.We assumed Γ=50J/m 2 ,   =1MPa,   =1kPa,   =  ,   = 0.1, which leads to  0 =20nm.The curves for standard and Maxwell materials become very close for fast propagation, as expected.
Looking in details to Fig. 2, it can be seen that there is an incubation period during which the crack tip is stationary.Initiation of crack growth occurs in each case at time  =   where there is a steep rise in the work rate in Fig. 2. Thereafter, the crack grows at a constant rate in each case, given by Eq. ( 42), with the case plotted in red with   √  Γ (1 −  2 ) ⁄ ⁄ = 28.5 growing most rapidly.In this case, the crack eventually breaks the component into 2 pieces within the time period plotted, with the steep rise in the work rate at approximately    ⁄ = 9 denoting the final fracture of the component.Similarly, the specimen breaks into two parts for the intermediate load (blue curve) at approximately    ⁄ = 200 The breaking of the component into 2 pieces also occurs for   √  Γ (1 −  2 ) ⁄ ⁄ = √1.58(black line), but at a time later than the interval plotted in Fig. 2. Notice that Fig. 2 plots results for the standard material as solid lines and for the equivalent Maxwell material as dashed lines, which we shall discuss later.Considering both materials, it can be seen that the external work per unit area of crack growth is positive and highest immediately after initiation of crack growth, and is highest for the most lightly loaded case having   √  Γ (1 −  2 ) ⁄ ⁄ = √1.58.As this case is the one with the slowest rate of propagation of the crack, the red case with   √  Γ (1 −  2 ) ⁄ ⁄ = 28.5 is that having the fastest rate of crack propagation, and the curve in blue for   √  Γ (1 −  2 ) ⁄ ⁄ = 16 is associated with an intermediate rate of crack growth, the results in Fig. 2 show that a faster rate of crack propagation is not associated with a larger rate of external work per unit area of crack growth.This was already noted immediately after initiation of crack growth in conjunction with Eq. ( 43).
It can be further seen in Fig. 2 that the work done per unit area of crack growth declines rapidly after initiation of crack growth and eventually becomes negative for the standard material case (solid lines).During the stage where   ⁄ < 0 the work of fracture is being provided by release of strain energy already stored in the component, and the component does work on the applied loads, leading to the negative value of   ⁄ .This indicates that, in this stage of crack growth, the applied load is dropping sufficiently rapidly that the component is experiencing a reduction in its stretch despite crack growth and viscoelasticity causing the component's compliance to increase, i.e., the displacement in Eq. ( 30) is declining and the viscoelasticity associated with Eq. ( 32) is insufficient to lead to increases in the displacement in Eq. ( 32).The conclusion from these observations is that, during the stage of crack growth where   ⁄ < 0, the higher external work per unit area of crack growth anticipated to be associated with higher rates of crack growth by de Gennes (1996) and Persson and Brener (2005) does not exist because the external work per unit area of crack growth is negative.
It can be seen in Fig. 2 that the external work per unit area of crack growth for the standard material becomes positive once again when the crack length approaches the width of the specimen (blue and red curves).This indicates that, when the crack tip approaches the distal free surface of the component, the viscoelastic compliance of the component increases rapidly and it stretches despite the applied load dropping as the crack continues to lengthen.As a result, the applied loads are forced once more to do work on the component.
We note further that at no time during crack growth is the external work per unit area of crack growth in steady state in Fig. 2. The value of   ⁄ is always in transient mode once crack growth commences, with the transient changes being greatest when   ⁄ > 0. Thus, one cannot claim that there is a steady state value of a positive external work rate per unit area of crack growth associated with a constant rate of crack growth, as the proposals of de Gennes (1996) and Persson and Brener (2005) would imply.To the extent that there is any hint of a steady state value for   ⁄ it is when   ⁄ < 0, as   ⁄ changes relatively little in this stage, but only for the low loads regime, and not in general.However, we have trouble rationalizing a negative rate of external work per unit area of crack growth in a viscoelastic body with the concepts of de Gennes (1996) and Persson and Brener (2005).Furthermore, we used ℎ =  to compute the results in Fig. 2. A value of ℎ >  will alter the results in Fig. 2 but without changing any of the features associated with crack growth.
As the work rate is interesting to compare with the dissipation-based theories, the total work for rupturing the specimen in two parts may be also a practical quantity of interest.The total work is integrated numerically from the results of Eq. ( 37), and then divided by the value Γ(b-a0), which is the pure work needed to create the two surfaces.This results in Fig. 3 where the total dimensionless work is plotted as a function of   /  where   = √  Γ (1 −  2 ) ⁄ .We are deliberately not taking into account the initial strain energy in loading the specimen, nor the possible final strain energy dissipated at the time of fracture, which would not quantitatively change the conclusions anyway, and would be very small quantities especially considering fast loading and unloading, compared to the quantities plotted in Fig. 3.As it is clear from Fig. 3, with the reasonable set of constants we used for an elastomer, we obtain even with standard material that the total work seems to be higher than the pure work of fracture amplified by the factor     ⁄ = 1000, due to the large work rate at the beginning of propagation, and the large one at the final instants.So we have effectively an ultratough interface, except for the very fast propagation limit, where the total work returns ultimately to the expected one without dissipation.Notice that the curve for the standard material case at very low   /  should return also to the same very low limit, but this is not reproduced here due to numerical difficulties in integrating the equations at very large times.The result for the Maxwell material instead grows unbounded for very low speed, consistently to what we will show analytically in the next paragraph.
The striking result of this example reproduced in Fig. 3 is that for most of the range of applied loads, the total work to propagate the crack across the entire section of the specimen is larger than the pure work of fracture by a factor bigger than the expected factor     ⁄ = 1000. .Other parameters are as in Fig. 2. Solid line is the standard material, and dashed line refer to the equivalent Maxwell material with viscosity appropriate to obtain similar creep function of the standard material for not too long times, namely  =     .For most of the range of applied loads, the total work to propagate the crack across the entire section of the specimen is larger than the pure work of fracture by a factor bigger than the expected factor     ⁄ = 1000.

Crack growth in a Maxwell material
To illustrate our insights more clearly, we consider a Maxwell material having creep response given by where  is the viscosity.To make a quantitative example for the case we are considering in Fig. 2,3, we plot in Fig. 4, the example of Maxwell material which best represents the standard material considered so far, having   =1MPa,   =1kPa (solid line), which is obtained with a viscosity  =     (dashed line).
It is clear that the Maxwell material and the standard material's creep functions depart significantly only for  >   , which corresponds to the fact that our results in Fig. 2 also depart at about that time for the two materials.The calculations follow the same steps as done for the standard material, and are not repeated here for brevity.

Discussion
Results represented as being derived from an energy balance (de Gennes (1996), Persson and Brener (2005), Persson (2017Persson ( , 2021)) look at the steady state, and we have shown in the present paper that viscoelastic propagation, looked from a cohesive model, can be intrinsically a transient process, despite the crack growth rate may be constant.Secondly, these models tend to neglect that strain energy could be cumulating in the bulk of the solid (and not just from the asymptotic, singular, region of stresses), and therefore in general are simplified models, which in some cases have been shown to give results compatible with the cohesive models (Ciavarella et al., 2021, Hui et al., 2022).However, the present paper has given an example of a quite more general picture of the energy fluxes during viscoelastic crack propagation which has given rise to interesting and quite unexpected results, in particular that the work rate can be much higher than what expected from the literature and happens not at large speeds but at intermediate speeds.Similarly, that the total work in propagating a crack up to a specimen failure can be also orders of magnitude higher than what expected from the literature, and being minimal only at loads extremely close to the load for very fast crack growth.
This shows that the cohesive model is more general in addressing viscoelastic crack propagation problems, which can show transient dissipation, where dissipation can come also from areas of nominal stress, and the de Gennes-Persson-Brener picture cannot be applied in general.

Figure 2 .
Figure 2. Absolute value of the external work per unit area of crack growth,   ⁄ , as a function of elapsed time for the edge crack in Fig.1 for a standard material having     ⁄ = 1000 and retardation time   for values of   √  Γ (1 −  2 ) ⁄ ⁄ = √1.58,16, 28.5 (black, blue, red curves): solid lines refer to the standard material, while dashed lines refer to the equivalent Maxwell material with viscosity appropriate to obtain the same creep function of the standard material for not too long times, namely  =     .The load is suddenly applied at  = 0 and crack growth commences at  =   in each case where   has a different value in each case.The work rate,   ⁄ , is normalized by the fracture energy Γ.We assumed Γ=50J/m 2 ,   =1MPa,   =1kPa,   =  ,   = 0.1, which leads to  0 =20nm.The curves

Figure 3 .
Figure 3. Value of the total external work done during propagation from crack size a=a0 up to a=b (complete breaking), for the edge crack in Fig.1 as a function of   /  where   = √  Γ (1 −  2 )⁄ .Other parameters are as in Fig.2.Solid line is the standard material, and dashed line refer to the equivalent Maxwell material with viscosity appropriate to obtain similar creep function of the standard material for not too long times, namely  =     .For most of the range of applied loads, the total work to propagate the crack across the entire section of the specimen is larger than the pure work of fracture by a factor bigger than the expected factor     ⁄ = 1000.

Figure 4 .
Figure 4.The creep compliance for standard material considered in our examples where for the solid material   =1MPa,   =1kPa (solid line), and for the "equivalent" Maxwell material with viscosity  =     (dashed line).It is clear that the Maxwell material and the standard material's creep functions depart significantly only for  >   .