1. Introduction
Jacques Heyman’s 1966 Article “The Stone Skeleton” is frequently cited at the beginning of articles on the assessment of masonry arches, including bridge arches and the arches used in ancient architecture [
1]. The limit analysis theorems, described in Heyman’s paper for the analysis of masonry arches, are still widely applied to the analysis of heritage structures. This form of analysis is very useful for determining reasonable bounds on behavior or for estimating the capacity of a structure. However, the simplifying assumptions adopted by Heyman and some of the effects in masonry structures that he has ruled out by assumption require that these methods for the analysis of arches be modified, elaborated, and sometimes corrected when analyzing real masonry structures. It can also be argued that many contemporary researchers overlook important effects in masonry arch analysis by focusing strictly on the type of limit analysis that Heyman considered.
The Stone Skeleton is an application of the theorems of plastic collapse. Heyman noted that masonry arches, similarly to steel frames, develop hinges under loading, and that their limit state is the formation of one hinge more than the degree of statical indeterminacy. In this observation, Heyman followed Kooharian [
2], although he expanded greatly on Kooharian’s findings. Heyman’s method, spurred on by a growing interest in the analysis of ancient structures, has been used extensively in the analysis of masonry structures built prior to 1900. Computer programs incorporating the principles of plastic analysis have been written and distributed (LimitState:RING [
3] and Archie-M [
4]), and the method has been widely employed in the assessment of masonry arch bridges [
5]. Some of these programs, such as LimitState:RING, incorporate failure modes beyond the hinging mechanism described by Heyman.
3. The Lasting Contributions of Heyman’s Method
Before the application of plasticity methods to the analysis of unreinforced masonry arches, there was a reliance on the “middle 1/3 rule”. That is that the thrust line of an arch must remain within the middle 1/3 to ensure the safety of the arch. In fact, as Heyman discussed, an excursion of the thrust line outside of the middle 1/3 causes neither a local failure nor a global failure of the arch. It simply indicates a migration of the compression zone towards the edge of the arch. For a voussoir with a notionally unlimited compressive capacity, the compression zone can extend to the face of the arch.
After its publication, Heyman’s method became the principal legitimate means of analyzing an arch available to most engineers. Although there are many previously published elastic methods, both analytical and graphical, they all require tedious calculations, and all depend on locating the thrust line in the middle third of the arch. The “middle third rule” is an unnecessary restriction that limits the apparent capacity of an arch, and does not explain why so many examples of arches still stand with a thrust line outside of the middle third. Heyman’s methods found instant applicability in the analysis of gothic architecture, and allowed the profession to make sense out of many of the pre-rational techniques and conceptions of medieval architecture.
5. Limit Analysis Assumptions
Heyman’s analysis is based on three assumptions:
Masonry has no tensile capacity;
Masonry has unlimited compressive capacity;
Masonry units do not slide with respect to each other (effectively, unlimited shear capacity).
In consequence of these assumptions, two adjacent stones can only “hinge”, with compression on one face (intrados or extrados), opening with zero stress on the other.
The “no tensile capacity” claim is a conservative assumption, although not always a realistic one. Fanning and Boothby [
16] and Boothby and Fanning [
17] noted that, to ensure the fidelity of the results of a masonry arch analysis, it is necessary to choose a small value of tensile capacity in the model used to replicate the behavior of a masonry arch bridge. Unlimited compressive capacity, and unlimited shear capacity, effectively mean that the crushing/sliding of the stones is not a significant factor in arch failure. Heyman justifies the second assumption by noting that stresses in masonry structures are far less intense than the compressive capacity; however, in the formation of a hinge, the compressive stress in the masonry becomes very large as the area resisting compressive force shrinks. Heyman does observe that “local spalling at the pressure points may sometimes be seen, but these are not signs of overall danger” [
18]. However, the images in
Figure 2 of a colonette at Milan Cathedral and a flying buttress at Reims cathedral show instances where the compressive stress may be a sign of overall danger. The compressive stresses within a masonry pier, arch or buttress may precipitate a failure by causing a compression failure in the material via compressive stress at a hinge.
Figure 3 shows a vault rib at the Cathedral of Milan, where a critical failure has occurred in the material at the base of a masonry arch. The failure mode is complicated, but the hinging of the arch is accompanied by a shear failure. A plan of the Cathedral in
Figure 4 shows the location of the defects depicted in
Figure 2a,
Figure 3,
Figure 5 and
Figure 6.
One can also note, in any study of actual failures of bridge or building arches, that the sliding of voussoirs, bricks, or portions of the arch occurs too commonly to be ruled out by assumption.
Figure 5a,b show the sliding of the blocks in an arch, also at Milan Cathedral. The plastic theory of masonry can be adapted to account for such sliding, as described above. For relatively low span-rise ratios and relatively thick arches, it is often necessary to consider this failure mode. A flat arch, for instance, is locked against a hinging mode and can only fail by sliding. Also from the Cathedral of Milan,
Figure 6 shows the sliding of a voussoir in a diagonal arch rib [
19]. In this case, as in many other instances, this sliding was partly the result of the abutment yielding.
6. What Plasticity Approaches to Unreinforced Structural Masonry May Be Missing
The elements of unreinforced masonry structures that are not modeled in Heyman’s analysis can be grouped into material behavior and structural behavior. The material behavior elements that we will consider are tensile capacity, compressive capacity, and creep. The structural behaviors of interest are failure modes beyond the classical four-hinged arch failure: creep instability, shear failure, compressive failure, or others.
There are some other types of structural actions that are excluded from the simple application of the plasticity approach to masonry structures. The most notable are foundation movements (beyond a simple self-limiting settlement that is accommodated by hinges in the arch), creep, cracking tensile capacity, and the influence of other parts of the structure, such as spandrel walls, fill, and haunching. Heyman’s justification for ignoring compressive failure in the masonry is discussed above.
Issues of creep, although noted by Heyman, do not really find a place in his form of analysis [
20]. Creep of masonry has had a destabilizing influence on many structures from the Gothic era, although it is often difficult to separate creep from the yielding of the foundation under abutments or piers. Beverley Minster, for instance, has vaults that are noticeably sagging, in some cases by over 10 cm (
Figure 7). This appears to be the result of a combined action of the rotation of the buttresses due to foundation conditions and creep within the fabric. The figure, from a survey taken above the vaults in 2008, describes the sagging of the nave vaults. Although it is true that creep stresses have no influence on the plastic capacity of a single hinge, creep cannot be dismissed as having no influence on capacity where a change in geometry is the result of the creep of a structure. The weakening of buttresses due to the creep of the vaults is described by Ochsendorf [
9], although as a problem in limit analysis, without specific reference to the expected magnitude of the creep.
In the apse of the church of San Bassiano (
Figure 8), a pattern of radial cracks has developed in the vault. The eventual failure of this structure would most likely be a combined result of cracking in the masonry, the subsequent modification of the load paths through the masonry, and the exacerbation of the cracks.
Heyman [
18] and others have described how, in finding an alternative load path that satisfies the lower bound condition, a pattern of cracking can be considered harmless. The examples that he cites are the “Sabouret cracks” in a Gothic vault, cracks along the axis of the nave in line with the upper nave walls, and the meridional cracking in the dome of St. Peter’s examined by Poleni [
21]. In reality, cracks in masonry structures are rarely harmless, and may propagate or grow larger, even if there are alternate load paths.
Figure 8.
Apse dome at San Bassiano (
a) from inside and (
b) from outside; Cardani and Massetti [
22].
Figure 8.
Apse dome at San Bassiano (
a) from inside and (
b) from outside; Cardani and Massetti [
22].
This damage in a vault or dome may also interact unfavorably with creep in the masonry or the settlement of the supports of the structures under investigation. In view of the foregoing, Heyman’s discussion of the stability through time of a historic masonry structure after its construction may be extended to include explicit discussions of masonry creep and foundation settlement.
On the other hand, even a very small permissible tensile strength within the masonry may have a significant effect on the capacity of the structure. A previous study of bridges [
23] concluded that an increase in the collapse load by a factor of 2 may result from the selection of a tensile strength of 0.5 N/mm
2 or less. Fanning and Boothby concluded similarly that the inclusion of a tensile strength of 0.5 MPa or less is required for fidelity to the results from testing [
16]. The same is probably true for the accuracy of models used to study buildings. Studies performed via diagonal compression tests on square specimens extracted from historic masonry buildings have found a significant tensile stress capacity of the building fabric [
24].
Figure 9, taken from [
25], shows the difference in collapse load between a case study arch with and without spandrel walls. The method used by these authors, Discontinuity Layout Optimization, is a computer-based implementation of upper bound analysis.
7. Continuing Reliance on Heyman’s Methods
Over 50 years after the publication of The Stone Skeleton, there have been many more studies of medieval masonry architecture, and significantly more failure modes of masonry arches have been noted beyond the four-hinged arch mechanisms described by Heyman. Heyman’s method clearly has enduring utility in establishing the bounds or approximate values of the strength of an arch. However, as an assessment tool, this method in general needs to be supplemented with a more realistic treatment of the capacity of an arch. Heyman’s method, however, often seems to be used to the exclusion of other methods, sometimes out of a lack of understanding of the other possible parameters. A brief look at some of the present-day literature, as presented in the paragraph below, shows that limit analysis, often two-dimensional limit analysis, considering hinging only continues in many cases to be used as a primary tool for the analysis of masonry arches, for the assessment of masonry arch bridges, for the assessment of masonry heritage buildings, and for the review of repairs to masonry arches. There is, of course, nothing wrong with this, if the analyst recognizes fully the limitations of these methods, and is prepared to modify the assumptions used in this analysis and to bring other methods to bear on some of the factors overlooked in this analysis.
The analysis of arches with yielding abutments essentially began with Ochsendorf’s dissertation [
9]. This is an original work on the application of theorems of plastic collapse to arches on spreading supports and the buttressing force required to resist this spreading. Only hinging mechanisms are considered. His section on future works describes a further need for investigations of sliding failures of the buttress, but not the arch. Ochsendorf carried out further investigations of hinging mechanisms in structures with failing abutments in [
8]. The approach of looking at abutment failures combined with strictly hinging mechanisms of arch collapse has continued to be applied since that time. Morer et al. [
26] and Riveiro et al. [
27], in looking at a complex, often rebuilt, multispan stone arch bridge, used Lidar scanning to determine the geometry of the bridge, and then applied a computer implementation of hinging analysis as their means of determining the capacity of the bridge. Given the two-dimensional nature of this analysis and the lack of accounting for other effects such as actual vehicle configurations, material failure, cracking, sliding, and settlement, this analysis would have to be considered incomplete. Galassi et al. [
28], in looking at failure mechanisms on moving supports, considered only hinging mechanisms on unloaded voussoirs without spandrel walls. Galassi showed settling arches as developing hinging mechanisms exclusively [
29]. Zampieri, Amoroso, and Pellegrino [
10] and Zampieri, Cavalagli, Gusella, Pellegrino [
30] offer similar studies of support settlement in arches that only consider hinging mechanisms.
As the above summary indicates, the analysis of arches with yielding abutments focuses on the mechanisms of an arch with abutment displacement, either horizontal or vertical, subject to the standard assumptions of no compression failure, no tensile stress, and no sliding. Continuing through the present time, a variety of authors invoke these assumptions in the investigation of the effects of abutment movements, and refer to Ochsendorf’s or similar experiments showing the result of such abutment movement to be the formation of a three-hinged mechanism in the arch.
Any arch with a spandrel wall above is likely to fail in a very different manner. The additional weight of masonry above the arch, applied at the same time to a larger portion of the back of the arch, is very likely to cause the upper part of the arch to fall in a wedge-like failure mode as the abutments of the arch spread. The propensity for this type of failure cannot be described effectively in energy or limit state terms because the frictional force between voussoirs increases in proportion to the weight of the segment of the arch that is subject to sliding. Nevertheless, this type of failure of an arch is common, as shown in the images in
Figure 10.
In the above discussion, we have looked at examples of the analysis of structures with spreading, yielding, or sinking supports, where one would expect to find sliding blocks, out of plane, fill, and other effects. This critique can be extended to the seismic analysis of masonry bridges, where, for instance Breccolotti et al. [
31] proposed a method for the seismic analysis of masonry arch bridges that considers only in-plane hinging mechanisms, omitting sliding mechanisms, abutment failure, and transverse failure. Zampieri et al. [
32] considered the transverse loading of multi-span high pier masonry bridges, but the analysis focused on hinging mechanisms resulting from such transverse loading. DeLorenzis, DeJong, and Ochsendorf [
33] applied rigid block theory, with assumptions intact, to the analysis of impulsive dynamic response, including the exclusive hinging mechanism response of the system to base motion. In a seismic event, the failure of an arch bridge in the transverse direction is much more likely than failure in the direction along the axis of the bridge, and the separation at joints and the spreading of supports make classic hinging mechanisms less likely.
Because sliding failure is a common failure mode, especially for low-rise arches, the “no sliding” assumption is overly restrictive at best, and misleading at worst. Although arches with spreading supports may fail by hinging when the bare arch ring is considered (except in some of the cases investigated by Rios et al. [
34]), the universal presence of a wall above the arch ring often shifts the failure mode to an internal sliding in the arch under yielding abutments. A comprehensive explanation for this frequently occurring sliding mechanism is not yet available, and a continuing focus on limit analysis methods may not produce an explanation for this failure mode.