Abstract
Scaling effects on the creep crack growth behaviour are investigated by analyzing the results of compact tension (CT) tests on different-sized notched steel specimens appearing in the literature. Creep crack growth rate data are correlated to the elastic stress-intensity factor in terms of a Paris-type law, da∕dt = C 0 K q , where C 0 turns out to be a crack-size dependent coefficient of proportionality. Considering specimens with the same loading configuration (CT) and the same thickness, the observed crack-size effect on the creep crack growth rate is discussed on the basis of self-similarity considerations, and geometrically interpreted in terms of fractal tortuosity of the crack profile. A size-independent formulation of the creep crack growth law correlating renormalized quantities is finally deduced and confirmed by the experimental results.
Keywords
Introduction
Design and integrity assessment of structural components operating at high temperatures typically have to account for the creep behaviour under steady-state conditions, where the strain rate
In fact, failure of structural components operating under creep conditions can occur by either creep rupture or creep crack growth. The former is most likely in components which are initially flaw-free or containing benign defects. However, the presence of inherent sharp defects, which can grow and ultimately cause brittle fracture, has necessitated characterizing the creep crack growth as well, often relegating the “flaw-free” condition to a mere idealization.
Despite the predictive character of fracture mechanics, whereby the results from laboratory specimens can be extrapolated to make predictions on full-sized structural components, and the number of studies conducted about creep, only a few of them were directly aimed at investigating size effects. In the current paper, the approach of fracture mechanics is supplemented by self-similarity and fractal geometry considerations to re-examine the observed crack-size effects on the creep crack growth behaviour for metallic materials.
To properly identify a fracture parameter that is able to account for creep in the crack growth behaviour, a schematic representation of the creep deformation process in a cracked body has been considered. Immediately after loading, r −1∕2-type (being r the distance from the crack tip) and HRR [19,30] singular stress distributions are respectively generated in linear elastic and power-law hardening materials. In either case, if extensive creep deformation has time to increase in size before the crack propagation, neither the elastic stress-intensity factor, K, nor the J-integral [13,29] —characterizing the HRR solution— can adequately characterize the crack-tip stress field. This is due to the nonlinear and time-dependent character of creep deformation, which is not admitted in linear elastic fracture mechanics (LEFM) and elastic-plastic fracture mechanics (EPFM) formulations [16,34,46].
By analogy with the J-integral, Landes and Begley [20] and Nikbin et al. [24] defined the C*-integral as a path-independent crack-tip parameter characterizing stress and strain rate fields when creep strains dominate:
Several studies [17,20,25,31,33] demonstrated that the creep crack growth rate (CCGR) can be described by an expression of the form:
More recently, extensive work in the field of time-dependent fracture mechanics (TDFM) has been devoted to the analysis of the so-called constraint effects on CCGR. In analogy with EPFM [26,27,35], results of recent research have extended TDFM beyond the limits of single-parameter theory, where a two-parameter (C*-Q) approach has been proposed to characterize crack-tip constraint at elevated temperatures [4,6,23,35]. The Q parameter is introduced for quantifying the constraint effects during creep, in terms of deviations of the actual crack-tip stress-field—generally determined by numerical FEM analysis—from the HRR field in Eq. (3). Numerical and experimental studies reported in the literature [39,40,42,44,45,49–53] regard the dependence of these effects on loading configuration [44,49,52], as well as crack depth [39,45,51] (in-plane constraint effects) and specimen thickness (out-of-plane constraint effects) [40,42]. The general conclusion is that higher constraint levels at the crack tip provide higher creep cracking rates.
Within the single-parameter assumption, several parameters have been hystorically proposed for characterizing the creep crack growth behaviour, as reported for instance in [24,43]. It was emphasized that for creep-ductile materials, defined as those in which the creep crack growth rate is negligible compared to the creep zone expansion rate, the creep zone completely engulfs the uncracked ligament before the crack actually propagates, causing extensive stress redistribution [16,34,46]. This circumstance occurs in highly susceptible materials to creep deformation, i.e. when n ≫ 1 [38], whereby the crack-tip singularity tends to disappear according to Eq. (3a). That makes a description in terms of net section stress or reference stress 𝜎
N
appropriate [18,22]:
A relevant attempt [56–59] to consider creep as a thermally activated process led to the following Arrhenius-type relationship:
When a little creep deformation accompanies crack growth, the initial stress distribution remains virtually unaltered by creep, and K (or J) continues to characterize the stress state around the crack tip. Such creep-brittle behaviour, defined as that in which the creep crack growth rate da∕dt is comparable to the creep zone expansion rate [36,37], occurs in case of high CCGR values when the specimen always remains in small-scale creep condition. For this circumstance, occurring at short times (when the stress-strain response is predominantly elastic) or for materials designed to resist large-scale creep deformations (i.e. with creep sensitivity n ≅ 1), the creep crack growth rate is satisfactorily described in terms of stress-intensity factor K by means of a Paris-type law [12,28]:
Hereinafter, the attention will be focused on a set of CCGR data from an experimental campaign on CT steel specimens, originally correlated to both C* and K [41]. Despite the fact that CCGR is currently interpreted in terms of C* in the framework of constraint effects for most practical situations, on that occasion some success was found in correlating CCGR data to the elastic stress-intensity factor, K. In this case, the observed shift of the experimental curves is compatible with a crack-size effect, which can be interpreted in terms of self-similarity and fractal tortuosity of the crack profile, along the line of recent studies on fatigue crack growth [9–12,28].
The creep crack growth behaviour of CT 1% CrMoV steel specimens of various sizes —standard and ultra-large specimens respectively of 50.8 mm and 254 mm in width, with thickness ranging from 6.35 to 63.5 mm— was investigated by Tabuchi et al. [41]. After introducing a fatigue pre-crack of 3 mm for standard specimens and of 15 mm for ultra-large specimens at room temperature, the authors conducted creep crack growth tests at 538 °C. The tested geometry is outlined in Fig. 1, where the markers indicate the specimens under investigation. In this paper, the analysis is restricted to specimens with the same thickness (12.7 mm), and differing in width, i.e. in the initial crack length (3 and 15 mm), in order to isolate the crack-size effects on the CCGR.

Geometry and size of CT specimens tested (adapted from [41]).
The authors discussed the experimental data both in the framework of TDFM and LEFM, correlating satisfactorily creep crack growth rate by Eqs (4) and (7), respectively. Characterization of da∕dt in terms of the C*-integral pointed out increasing creep crack growth rate with increasing specimen thickness, whereas neither specimen-width nor crack-size effects were observed. The specimen-thickness effect was already interpreted in the framework of the constraint effects, whereby the achievement of plane-strain conditions in very thick sections result in less creep ductility and higher cracking rates. Instead, the lack of crack-size (or specimen-width) effects on the creep crack growth rate was considered surprising, given that the pre-crack of ultra-large specimens was five times longer than standard ones, and worthy of further investigation.
Actually, the characterization of creep crack growth rate in terms of K by Eq. (7) did not demonstrate thickness effects for specimens with the same width, but, rather demonstrated a right- and downward shift of CCGR data with increasing specimen width, or pre-crack length. The power-law exponent q of the two fitting curves varies over the range 14.1–14.7 (the average exponent is

Shift of experimental CCGR curves correlating da∕dt with K (adapted from [41]) for different initial crack lengths, accounting for the crack-size effect.
From a physical point of view, correlating creep crack growth rate da∕dt by an expression in the form of Eqs (4), (5) or (7) provides little information on the multiple factors that affect creep crack growth. In the case of the above results, where CCGR data are satisfactorily correlated to K, a functional dependence upon a series of LEFM parameters can be properly stated:
Considering 𝜎
u
, K
C
, X, T
M
as independent variables, the dimensional analysis [5] yields the following expression in terms of dimensionless parameters:
On the other hand, power-law relationships, such as those presented to describe creep crack growth rate, give the evidence of self-similarity, wherein a phenomenon reproduces itself over different time and space scales. It has been observed that self-similar solutions describe the intermediate asymptotic behaviour of fatigue crack growth when “the influence of the initial conditions has disappeared but the influence of the instability has not yet intruded” [1–3,32].
Accordingly, characterization of creep crack growth in terms of power laws is signature of an intermediate asymptotic regime, which is achieved in the mid-range of growth rates, i.e., when the crack has started to propagate although the material is sufficiently far from failure. The so-called incomplete self-similarity prevails at this stage, corresponding to a power-law dependence of creep crack growth rate on certain dimensionless parameters.
Let us consider a and K. Since the above experimental evidences have shown the dependence of the creep crack growth rate on a and K, incomplete self-similarity is assumed in the corresponding dimensionless parameters a∕a
0 and K∕K
C
. Thus, Eq. (9) takes the following power-law asymptotic form:
Deviations from the mid-range power-law behaviour, predictable within intermediate asymptotics, have been experimentally verified [54] when the conditions of non-propagating crack, for K → K th , and Griffith-Irwin instability, for K → K C , are approached, yielding a Paris-type curve illustrated in Fig. 3.

Schematization of CCGR diagram: power-law stage corresponds to the intermediate asymptotic behaviour.
Since the publication of the celebrated article by Mandelbrot and co-workers [21], the fractal nature of metal fracture surfaces has been widely recognized. Such a fractal character has been revealed also for metallic substructures under creep conditions [14,48]. Repeated observations at various magnifications have shown that the structure replicates itself in a self-similar way in a range of intermediate scales, falling between the micro-scale —influenced by grains, inclusions and dislocations— and the macro-scale —characterized by the size of the specimen and the notch from which the crack propagates. Accordingly, a geometric interpretation of crack-size effects on the creep crack growth rate is obtained in the context of fractal geometry [7,8], whereby the crack profile can be modelled as an invasive fractal set a*, with projected length a and topological dimension 1 + d
G
(0 < d
G
< 1):
Exploiting Eq. (12), an application of the derivation rule for composite functions gives:
The experimental CCGR curves (with da∕dt in mm h−1 and K in MPa m1∕2) are renormalized by a transformation of the axes K and da∕dt of the form, K → K
∗ = Ka
−d
G
∕2 and da∕dt → da
∗∕dt = (1 + d
G
)a
d
G
da∕dt, where a is equal to 3 and 15 mm. The renormalized data collapse onto a single line (shown in Fig. 4), illustrating the fulfilment of the crack-size independent creep crack growth law given by Eq. (16). The dimensional increment d
G
of the crack profile is adjusted to achieve the best data collapse, yielding d
G
= 0.58 ± 0.09 and, contextually, the following scale-invariant law:
Experimental data about creep crack propagation for specimens with a pre-existing sharp defect are re-examined, revealing a clear crack-size effect on cracking rates when correlated to the elastic stress-intensity factor, K. These findings are formalized in the framework of dimensional analysis and intermediate asymptotics, leading to a generalized formulation of the creep crack growth law, da∕dt = C 0 K q , where the effects of the initial crack length as well as multiple parameters are taken into account. An application of the fractal geometry leads to the definition of a scale-invariant law, da ∗∕dt = C 0 ∗ K ∗q , able to uniquely characterize the creep crack growth behaviour of the material, regardless of the size scale. Similarly to fatigue crack propagation (where analogous scaling laws are obtained by the change of variable t → N, i.e., time → number of cycles), fractal roughness of the crack surfaces—through the dimensional increment d G —is predicted to be responsible for the negative scaling of C 0, whereby short cracks propagate at higher rates according to the experimental evidence.

Collapse of the renormalized CCGR data onto the fitting line represented by Eq. (18) in the fractal K* vs da*/dt diagram.
Assuming for K
th
the same scaling effect as for K
C
and K,

Predicted nominal CCGR curves (a) and renormalized curve in the fractal diagram (b), where v
th
(or
Extensive experimental work has to be conducted to confirm the present fractal approach. In particular, the investigation of a threshold limit, below which creep crack growth does not occur (see dashed branch of Paris-type curves in Fig. 5), is related to the ability to inspect smaller and smaller cracks. Further studies could also concern the application of fractal geometry concepts to interpret crack-size effects on creep crack growth rate in terms of the C*-integral.
Footnotes
Conflict of interest
None to report.
