Size effect on structural strength

According to the classical theories of elastic or plastic structures made from a material with non-random strength ("f"t), the nominal strength ("σ"N) of a structure is independent of the structure size ("D") when geometrically similar structures are considered. Any deviation from this property is called the size effect. For example, conventional strength of materials predicts that a large beam and a tiny beam will fail at the same stress if they are made of the same material. In the real world, because of size effects, a larger beam will fail at a lower stress than a smaller beam.
The structural size effect concerns structures made of the same material, with the same microstructure. It must be distinguished from the size effect of material inhomogeneities, particularly the Hall-Petch effect, which describes how the material strength increases with decreasing grain size in polycrystalline metals.
The size effect can have two causes: 
Statistical Theory of Size Effect in Brittle Structures.
The statistical size effect occurs for a broad class of brittle structures that follow the weakest-link model. This model means that macro-fracture initiation from one material element, or more precisely one representative volume element (RVE), causes the whole structure to fail, like the failure of one link in a chain (Fig. 1a). Since the material strength is random, the strength of the weakest material element in the structure (Fig. 1a) is likely to decrease with increasing structure size formula_1 (as noted already by Mariotte in 1684).
Denoting the failure probabilities of structure as formula_2 and of one RVE under stress formula_3 as formula_4, and noting that the survival probability of a chain is the joint probability of survival of all its formula_5 links, one readily concludes that
The key is the left tail of the distribution of formula_4. It was not successfully identified until Weibull in 1939 recognized that the tail is a power law. Denoting the tail exponent as formula_7, one can then show that, if the structure is sufficiently larger than one RVE (i.e., if ), the failure probability of a structure as a function of formula_8 is
Eq. 2 is the cumulative Weibull distribution with scale parameter formula_9 and shape parameter formula_7; formula_11 = constant factor depending on the structure geometry, formula_12 = structure volume; formula_13 = relative (size-independent) coordinate vectors, formula_14 = dimensionless stress field (dependent on geometry), scaled so that the maximum stress be 1; formula_15 = number of spatial dimensions (formula_15 = 1, 2 or 3); formula_17 = material characteristic length representing the effective size of the RVE (typically about 3 inhomogeneity sizes).
The RVE is here defined as the smallest material volume whose failure suffices to make the whole
structure fail. From experience, the structure is sufficiently larger than one RVE if the equivalent number formula_18 of RVEs in the structure is larger than about formula_19 ; formula_20 = number of RVEs giving the same formula_2 if the stress field is homogeneous (always formula_22, and usually formula_23). For most normal-scale applications to metals and fine-grained ceramics, except for micrometer scale devices, the size is large enough for the Weibull theory to apply (but not for coarse-grained materials such as concrete).
From Eq. 2 one can show that the mean strength and the coefficient of variation of strength are obtained as follows:
(where formula_24 is the gamma function) The first equation shows that the size effect on the mean nominal strength is
a power function of size formula_1, regardless of structure geometry.
Weibull parameter formula_7 can be experimentally identified by two methods: 1) The values of formula_8 measured on many identical specimens are used to calculate the coefficient of variation of strength, and the value of formula_7 then follows by solving Eq. (4); or 2) the values of formula_29 are measured on geometrically similar specimens of several different sizes formula_1 and the slope of their linear regression in the plot of formula_31 versus formula_32 gives formula_33. Method 1 must give the same result for different sizes, and method 2 the same as method 1. If not, the size effect is partly or totally non-Weibullian. Omission of testing for different sizes has often led to incorrect conclusions. Another check is that the histogram of the strengths of many identical specimens must be a straight line when plotted in the Weibull scale. A deviation to the right at high strength range means that formula_18 is too small and the material quasibrittle.
Energetic Size Effect.
The fact that the Weibull size effect is a power law means that it is self-similar, i.e., no characteristic structure size formula_35 exists, and formula_17 and material inhomogeneities are negligible compared to formula_1. This is the case for fatigue-embrittled metals or fine-grained ceramics except on the micrometer scale. The existence of a finite formula_35 is a salient feature of the energetic size effect, discovered in 1984. This kind of size effect represents a transition between two power laws and is observed in brittle heterogenous materials, termed quasibrittle. These materials include concrete, fiber composites, rocks, coarse-grained and toughened ceramics, rigid foams, sea ice, dental ceramics, dentine, bone, biological shells, many bio- and bio-inspired materials, masonry, mortar, stiff cohesive soils, grouted soils, consolidated snow, wood, paper, carton, coal, cemented sands, etc. On the micro- or nano scale, all the brittle materials become quasibrittle, and thus must exhibit the energetic size effect.
A pronounced energetic size effect occurs in shear, torsional and punching failures of reinforced concrete, in pullout of anchors from concrete, in compression failure of slender reinforced concrete columns and prestressed concrete beams, in compression and tensile failures of fiber-polymer composites and sandwich structures, and in the failures of all the aforementioned quasibrittle materials. One may distinguish two basic types of this size effect.
Type 1: Structures that fail at crack initiation.
When the macro-crack initiates from one RVE whose size is not negligible compared to the structure size, the deterministic size effect dominates over the statistical size effect. What causes the size effect is a stress redistribution in the structure (Fig. 2c) due to damage in the initiating RVE, which is typically located at fracture surface.
A simple intuitive justification of this size effect may be given by considering the flexural failure of an unnotched simply supported beam under a concentrated load formula_39 at midspan (Fig. 2d). Due to material heterogeneity, what decides the maximum load formula_39 is not the elastically calculated stress formula_41 at the tensile face, where formula_42 = bending moment, formula_1 = beam depth, formula_44 and formula_45 = beam width. Rather, what decides is the stress value formula_46 roughly at distance formula_47 from the tensile face, which is at the middle of FPZ (2c). Noting that formula_46 = formula_49, where formula_50 = stress gradient = formula_51 and formula_52 = intrinsic tensile strength of the material, and considering the failure condition formula_46 = formula_54, one gets formula_55 = formula_56 where formula_57, which is a constant because for geometrically similar beams formula_58 = constant. This expression is valid only for small enough formula_59, and so (according to the first two terms of the binomial expansion) one may approximate it as
which is the law of Type 1 deterministic size effect (Fig. 2a). The purpose of the approximation made is: (a) to prevent formula_8 from becoming negative for very small formula_1, for which the foregoing argument does not apply; and (b) to satisfy the asymptotic condition that the deterministic size effect must vanish for formula_62. Here formula_63 = positive empirical constant; the values formula_64 = or 2 have been used for concrete, while formula_65 is optimum according to the existing test data from the literature (Fig. 2d).
A fundamental derivation of Eq. 5 for a general structural geometry has been given by
applying dimensional analysis and asymptotic matching to the limit case of energy release when the initial macro-crack length tends to zero. For general structures, the following effective size may be substituted in Eq. (5):
where formula_66 = strain gradient at the maximum strain point located at the surface, in the direction
normal to the surface.
Eq. 5 cannot apply for large sizes because it approaches for formula_67 a horizontal asymptote.
For large sizes, formula_8 must approach the Weibull statistical size effect, Eq. 3. This condition is satisfied by the generalized energetic-statistical size effect law:
where formula_69 are empirical constants (formula_70). The deterministic formula (5) is recovered as the limit case for formula_71. (Fig. 2d) shows a comparison of the last formula with the test results for many different concretes, plotted as dimensionless strength formula_72 versus dimensionless structure size formula_73.
The probabilistic theory of Type 1 size effect can be derived from fracture nano-mechanics. Kramer's
transition rate theory shows that, on the nano-scale, the far-left tail of the probability distribution of nano-scale strength formula_74 is a power law of the type formula_75. Analysis of the multiscale transition to the material macro-scale then shows that the RVE strength distribution is Gaussian but with a Weibull (or power-law) left tail whose exponent formula_7 is much larger than 2 and is grafted roughly at the probability of about 0.001.
For structures with formula_77, which are common for quasibrittle materials, the Weibull theory does not apply. But the underlying weakest-link model, expressed by Eq. (1) for formula_2, does, albeit with a finite formula_5, which is a crucial point. The finiteness of the weakest-link chain model causes major deviations from the Weibull distribution. As the structure size, measured by formula_18, increases, the grafting point of the Weibullian left part moves to the right until, at about formula_81, the entire distribution becomes Weibullian. The mean strength can be computed from this distribution and, as it turns out, its plot is identical with the plot of Eq. 5 seen in Fig. 2g. The point of deviation from the Weibull asymptote is determined by the location of the grafting point on the strength distribution of one RVE (Fig. 2g). Note that the finiteness of the chain in the weakest-link model captures the deterministic part of size effect.
This theory has also been extended to the size effect on the Evans and Paris' laws of crack growth in quasibrittle materials, and to the size effect on the static and fatigue lifetimes. It appeared that the size effect on the lifetime is much stronger than it is on the short-time strength (tail exponent formula_7 is an order-of-magnitude smaller).
Type 2: Structures in which a large crack or notch exists.
The strongest possible size effect occurs for specimens with similar deep notches (Fig. 4b), or for structures in which a large crack, similar for different sizes, forms stably before the maximum load is reached. Because the location of fracture initiation is predetermined to occur at the crack tip and thus cannot sample the random strengths of different RVEs, the statistical contribution to the mean size effect is negligible. Such behavior is typical of reinforced concrete, damaged fiber-reinforced polymers and some compressed unreinforced structures.
The energetic size effect may be intuitively explained by considering the panel in Fig. 1c,d,
initially under a uniform stress equal to formula_8 . Introduction of a crack of length formula_84, with a damage zone
of width formula_85 at the tip, relieves the stress, and thus also the strain energy, from the shaded undamaged triangles of slope formula_86 on the flanks of the crack. Then, if formula_86 and formula_88 are approximately the same for different sizes, the energy released from the shaded triangles is proportional to formula_89, while the energy dissipated by the fracture process is proportional to formula_90; here formula_91 = fracture energy of the material, formula_92 = energy density before fracture, and formula_93 = Young's elastic modulus. The discrepancy between formula_1 and formula_95 shows that a balance of energy release and dissipation rate can exist for every size formula_1 only if formula_8 decreases with increasing formula_1. If the energy dissipated within the damage zone of width formula_85 is added, one obtains the Bažant (1984) size effect law (Type 2):
(Fig. 4c,d) where formula_100 = constants, where formula_54 = tensile strength of material, and formula_102 accounts for the structure geometry.
For more complex geometries such an intuitive derivation is not possible. However, dimensional
analysis coupled with asymptotic matching showed that Eq. 8 is applicable in general, and that the dependence of its parameters on the structure geometry has approximately the following form:
where formula_103 half of the FPZ length, formula_104 = relative initial crack length (which is constant for geometrically similar scaling); formula_105 = dimensionless energy release function of linear elastic fracture mechanics (LEFM), which brings about the effect of structure geometry; formula_106, and formula_107 = stress intensity factor. Fitting Eq. 8 to formula_8 data from tests of geometrically similar notched specimens of very different sizes is a good way to identify the formula_91 and formula_110 of the material.
Size Effect in Cohesive Crack, Crack Band and Nonlocal Models.
Numerical simulations of failure by finite element codes can capture the energetic (or deterministic) size effect only if the material law relating the stress to deformation possesses a characteristic length. This was not the case for the classical finite element codes with a material characterized solely by stress-strain relations.
One simple enough computational method is the cohesive (or fictitious) crack model, in which it is assumed that the stress formula_111 transmitted across a partially opened crack is a decreasing function of the crack opening formula_112, i.e., formula_113. The area under this function is formula_91, and
is the material characteristic length giving rise to the deterministic size effect. An even simpler method is the crack-band model, in which the cohesive crack is replaced in simulations by a crack band of width formula_85 equal to one finite element size and a stress-strain relation that is softening in the cross-band direction as formula_116 where formula_117 = average strain in that direction.
When formula_85 needs to be adjusted, the softening stress strain relation is adjusted so as to maintain the correct energy dissipation formula_91. A more versatile method is the nonlocal damage model in which the stress at a continuum point is a function not of the strain at that point but of the average of the strain field within a certain neighborhood of size formula_85 centered at that point. Still another method is the gradient damage model in which the stress depends not only on the strain at that point but also on the gradient of strain. All these computational methods can ensure objectivity and proper convergence with respect to the refinement of the finite element mesh.
Fractal Aspects of Size Effect.
The fractal properties of material, including the fractal aspect of crack surface roughness and the lacunar fractal aspect of pore structure, may have a role in the size effect in concrete, and may affect the fracture energy of material. However, the fractal properties have yet not been experimentally documented for a broad enough scale and the problem has not yet been studied in depth comparable to the statistical and energetic size effects. The main obstacle to the practical consideration of a fractal influence on the size effect is that, if calibrated for one structure geometry, it is not clear how infer the size effect for another geometry. The pros and cons were discussed, e.g., by Carpinteri et al. (1994, 2001) and Bažant and Yavari (2005).
Practical Importance.
Taking the size effect into account is essential for safe prediction of strength of large concrete bridges, nuclear containments, roof shells, tall buildings, tunnel linings, large load-bearing parts of aircraft, spacecraft and ships made of fiber-polymer composites, wind turbines, large geotechnical excavations, earth and rock slopes, floating sea ice carrying loads, oil platforms under ice forces, etc. Their design depends on the material properties measured on much smaller laboratory specimens. These properties must be extrapolated to sizes greater by one or two orders of magnitude. Even if an expensive full-scale failure test, for example a failure test of the rudder of a very large aircraft, can be carried out, it is financially prohibitive to repeat it thousand times to obtain the statistical distribution of load capacity. Such statistical information, underlying the safety factors, is obtainable only by proper extrapolation of laboratory tests.
The size effect is gaining in importance as larger and larger structures, of more and more slender forms, are being built. The safety factors, of course, give large safety margins—so large that even for the largest civil engineering structures the classical deterministic analysis based on the mean material properties normally yields failure loads smaller than the maximum design loads. For this reasons, the size effect on the strength in brittle failures of concrete structures and structural laminates has long been ignored. Then, however, the failure probability, which is required to be formula_121, and actually does have such values for normal-size structures, may become for very large structures as low as formula_122 per lifetime. Such high failure probability is intolerable as it adds significantly to the risks to which people are inevitably exposed. In fact, the historical experience shows that very large structures have been failing at a frequency several orders of magnitude higher than smaller ones. The reason it has not led to public outcry is that the large structures are few. But for the locals, who must use the structures daily, the risk is not acceptable.
Another application is the testing of the fracture energy and characteristic material length. For quasibrittle materials, measuring the size effect on the peak loads (and on the specimen softening after the peak load) is the simplest approach.
Knowing the size effect is also important in the reverse sense—for micrometer scale devices if they are designed partly or fully on the basis of material properties measured more conveniently on the scale of 0.01m to 0.1m.