New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Simultaneous improvements of strength and toughness in topologically interlocked ceramics
Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved July 23, 2018 (received for review April 27, 2018)

Significance
Topologically interlocked materials (TIMs) use frictional sliding to generate large deformations and build toughness in otherwise all brittle components. TIMs can be up to 10 times more impact resistant than their monolithic form, but this improvement usually comes at the expense of static strength. Here we report a TIM design based on octahedral blocks which is not only much tougher (50×) than monolithic plates of the same material, but also stronger (1.2×). With no evidence of upper bounds for strength and toughness, TIMs have a tremendous potential as high-performance structural materials. Based on our experiments we propose a nondimensional “interlocking parameter” which could guide the exploration of new TIM designs and other new architectured systems.
Abstract
Topologically interlocked materials (TIMs) are an emerging class of architectured materials based on stiff building blocks of well-controlled geometries which can slide, rotate, or interlock collectively providing a wealth of tunable mechanisms, precise structural properties, and functionalities. TIMs are typically 10 times more impact resistant than their monolithic form, but this improvement usually comes at the expense of strength. Here we used 3D printing and replica casting to explore 15 designs of architectured ceramic panels based on platonic shapes and their truncated versions. We tested the panels in quasi-static and impact conditions with stereoimaging, image correlation, and 3D reconstruction to monitor the displacements and rotations of individual blocks. We report a design based on octahedral blocks which is not only tougher (50×) but also stronger (1.2×) than monolithic plates of the same material. This result suggests that there is no upper bound for strength and toughness in TIMs, unveiling their tremendous potential as structural and multifunctional materials. Based on our experiments, we propose a nondimensional “interlocking parameter” which could guide the exploration of future architectured systems.
Specific microstructures, heterogeneities, or hybrid compositions are now widely used in modern materials to generate high performance (1). These concepts are now pushed to the extreme with architectured materials, which contain highly controlled structures and morphological features at length scales intermediate between traditional microstructure and the whole component (2). These information-rich materials can be programmed with specific architectures, geometries, and interfaces to generate unusual and attractive combinations of properties and functionalities, for example negative Poisson’s ratio (3), unusually high stiffness (4, 5), or high toughness (6⇓⇓–9). Architectured materials include the now well-studied lattice materials which contain only a small fraction of solid (10). In contrast, the less-studied dense architectured materials are fully solid and are made of building blocks of well-defined size and shape, arranged in two or three dimensions (11, 12). The geometry and arrangement of the blocks can be designed to generate interlocking in “topologically interlocked materials” (TIMs) without the need for adhesive or mortar, a powerful concept pioneered by Estrin and coworkers (13⇓⇓–16). In particular, Dyskin et al. (14) explored how regular, convex polyhedral (platonic) solids can interlock and serve as building blocks for TIMs, with examples including cubes (15, 17) and tetrahedra (18, 19). Nonplatonic geometries with planar surfaces such as truncated tetrahedra (20⇓–22) or blocks with nonplanar surfaces (8, 23⇓⇓⇓–27) have also been proposed. Interestingly the interlocked blocks can still slide, rotate, or separate to some extent, providing a wealth of tunable deformation mechanisms and properties (12, 16, 28). For example, the frictional sliding of the blocks on one another can dissipate energy and confer TIMs with very high impact resistance compared with monolithic panels of the same materials, a strategy which can be used to overcome brittleness of glasses and ceramics. However, improvements in impact resistance and energy absorption are achieved at the expense of 40–80% losses in strength (11, 20, 25, 27, 29). Despite recent efforts in unifying designs (30⇓–32) and optimization (33, 34), there are still no comprehensive guidelines to select optimum architectures for given applications and requirements. In this work we systematically explored the design of architectured panels based on platonic shapes and their truncated versions to gain insights into the structure–mechanics–performance relations of these systems.
Design and Fabrication of Architectured Ceramic Panels
The base material we used in this study was calcium sulfate, a brittle mineral which can be shaped into complex geometries at room temperature. We choose a brittle material for the building blocks to highlight the capability of architecture to generate toughness from brittle components, and to avoid plastic deformations within individual blocks which can complicate the interpretation of the results (35). We used a replica casting technique to make a large number of calcium sulfate building blocks of high density (36). An array of building blocks was first 3D printed (Fig. 1 A and B) using a high-resolution digital light-processing printer (Micro HiRes model; EnvisionTEC). These blocks served as a positive shape to make a silicone mold (Fig. 1 C and D). A ceramic slurry (calcium sulfate powder + 19 wt % water) was then pressure cast into the mold (Fig. 1E), which produced ceramic building blocks with high density and homogeneity, uniform geometries, and smooth surfaces (see Methods for details). We reused the silicone mold several times to make many building blocks. The building blocks were then tape-transferred into a rigid aluminum frame to hold the blocks in place and to serve as an external ligament (Fig. 1F). The gap between the edges of the panel and the frame was filled with calcium sulfate paste which was cured before testing. This homogeneous region ensured a uniform transfer of force from the frame to the peripheral blocks. We focused on convex polyhedral blocks, assembled into panels made of a 7 × 7 array of blocks. The number of blocks (and therefore their size) was fixed for all geometries to remove size effects (we recently showed that larger blocks lead to higher overall properties, SI Appendix, Fig. S1). Following the geometrical construction of Dyskin et al. (14), we required the medial sections of the blocks to form a surface-filling tessellation made of regular polygons (Fig. 2). There are only three regular shapes that can tile the plane: triangles, squares, and hexagons. We did not consider triangles as medial surfaces, because they cannot produce 3D interlocking geometries. Square medial section (Fig. 2A) was used to generate four flat faces that formed a prism capped with flat square faces at z = ±h/2. The simplest geometry was a prism with dimensions l × l × h which was assembled into panels held by friction only and without interlocking. The width of the medial square was set to l = 5 mm and the height was set to h ≈ 4.7 mm. Interlocking geometries were introduced by tilting the opposed side faces of the prism inward by an angle θ, and by tilting the other pair of faces outward by the same angle θ. This procedure generated truncated tetrahedra (20) with geometries defined by three independent parameters: the medial section size l, the height of the block h, and the angle θ. In an architectured panel, each block was geometrically constrained and interlocked by its four neighbors. Fig. 2A shows the effects of increasing θ from θ = 0° to 2.5°, 5°, 10°, 20° up to θ ≈ 35.2°, which corresponds to the limiting case of a tetrahedron. Increasing θ promotes interlocking, but also decreases the contact area between the blocks. In addition, for θ > 0 the assembly leaves pyramidal empty space on the surface of the panel. To isolate and focus on the effects of geometry, we set l = 5 mm and adjusted h to maintain a constant areal density ρV/A = 1.09 gr/cm2 across all of the designs (ρ = 2.3 g/cm3 is the mass density for calcium sulfate, V is the volume of an individual block, and A is the area of its medial section). Using the same approach, we also generated blocks based on hexagonal medial sections (Fig. 2 B and C). These blocks were assembled into 7 × 7 panels with individual block size a = 3.25 mm, and the thickness h was adjusted to maintain the same areal density across all designs. The case θ = 0° produced an eight-sided prism with hexagonal base (Fig. 2B). Tilting of the side faces by alternating angles (+θ/−θ/+θ/−θ/+θ/−θ on the six side faces) produced truncated octahedra (Fig. 2B). The special case θ ≈ 19.5° corresponds to a regular octahedron, and further increase of θ produced truncated rhombohedra, up to the limit case of θ ≈ 28°corresponding to a full rhombohedron. Once assembled into panels, each block had six geometrically interlocked neighbors. Finally, a last case based on the hexagonal medial section is the regular dodecahedron, composed of 12 pentagonal faces. Because of the regular geometry of the dodecahedron, the interlocking angle is fixed at θ ≈ 10.81°, and to keep the areal density constant we set
Fabrication steps for the architectured panels: (A) three-dimensional printing of polymeric building blocks; (B) the 3D printed blocks are transferred to a substrate; (C) silicone is poured to make a negative replica that serves as the mold; (D) the 3D printed blocks are removed from the cured silicone mold; (E) calcium sulfate (CaSO4) is pressure cast into the silicone mold; (F) the building blocks are assembled and tape-transferred into an aluminum frame.
Overview of the 15 architectured panel designs explored in this study: For each group the geometry of the individual block is shown, together with a schematic of the assembly and pictures of 7 × 7 panels made of CaSO4 blocks: (A) truncated tetrahedra based on a square medial section; (B) truncated octahedra based on a hexagonal medial section; and (C) dodecahedron.
Mechanical Performance
The response of the panels to a transverse concentrated force was measured in quasi-static and impact loading conditions (Fig. 3A). The panel was simply supported and confined laterally by an aluminum frame. A localized force (or a 1-m/s impact) was applied on the top face of the center block. Fig. 3B shows the responses of an architectured panel made of octahedral blocks in quasi-static and impact conditions together with the behavior of a monolithic calcium sulfate panel of the same areal density. As expected, the monolithic panel produced a linear elastic response immediately followed by a brittle failure at small deflections (∼0.1 mm). Failure was catastrophic, with multiple cracks and fragmentation that extended to the edge of the panels (Fig. 3C). In contrast, the architectured panel produced a bell-shaped response typical of tough materials. Compared with the monolithic panel the initial stiffness was in general lower, but the maximum deflection and the energy to failure were an order of magnitude greater. In addition, failure was localized with only the center block missing at the end of the test (Fig. 3C). We also used in situ stereoimaging, image correlation, and 3D scene reconstruction to measure the 3D displacements and rotations of the individual blocks during the tests (see Methods for details). The 3D reconstructed shape obtained from image correlation at six snapshots (points A–F in Fig. 3B) at various times during impact are shown in Fig. 3D. Typical results show that individual blocks rotate and slide on one another near the impact site (the deformation mode in quasi-static was identical). To characterize relative block motion for the entire plate, we computed the total sliding area at the interface and the average block rotation at different time points during the experiment. Taken together, these results point to a well-defined and repeatable deformation mode that follows two distinct stages (Fig. 3B): In the prepeak regime (region I), the average rotation and sliding area between the blocks rapidly increases with panel deflection. As the blocks slide on one another their contact area decreases, which is compensated by an increased geometrical locking so that an increasing force is produced. The slope of the force-deflection curve decreases gradually as the sliding increases, with fluctuations attributed to stick-slip mechanisms which are typical of dry friction (37). At the onset of transition from stage I to stage II, several sharp drops can be observed in the force-deflection curve due to surface cracking and chipping of individual blocks. The force however keeps increasing after each drop as the blocks reinterlock. In this example the maximum force (∼180 N in quasi-static and ∼190 N in impact) was achieved in this plateau-like region. In the postpeak region (stage II), sliding only occurs at the interface between the center block and its neighbors, so that the total sliding area increases more slowly with deflection. The average rotation of the blocks also decreases as the deformation of the panel partially recovers (Fig. 3B). The panel finally fails by complete push-out of the center blocks, in this case at a displacement ∼3–4.5 mm. The failure is therefore localized, and the panel largely retains its shape and structural integrity, in contrast with the monolithic panel which failed catastrophically with multiple cracks (Fig. 3C).
Mechanical response of the architectured panels. (A) Schematic of the experimental setup. (B) Force-deflection curves for architectured panels made of octahedral blocks tested in quasi-static and impact conditions. The response of monolithic plate with the same areal density is also shown for comparison. The other two plots show the total sliding area and average rotation of the blocks over the course of the test obtained from stereoimaging and 3D reconstruction. (C) Postmortem samples: Monolithic plates fail catastrophically and by fragmentation, while in architectured panels failure is localized to the central block. (D) Three-dimensional reconstruction of the panels showing the average vertical displacement of the blocks at six points (A–F) during loading.
The shape of the force-deflection curves and these mechanisms were similar for all of the panels (SI Appendix, Figs. S2 and S3). However, the performance of the panels in terms of stiffness, strength, and energy absorption (area under the force-deflection curve) varied across geometries. In general, the stiffness of the architectured panel was about 50–75% lower than the monolithic panel (SI Appendix, Fig. S4). The strength was on the same order, but the energy absorption was far greater. All properties improved significantly with increasing interlocking angle, up to
Effects of architecture on the mechanical performance. Effect of interlocking angle on (A) stiffness, (B) maximum force, and (C) energy absorption; Maximum force-energy absorption chart and for all of the architectured panels explored in this study tested in (D) quasi-static conditions and (E) impacts. The width of the colored regions is representative of the variations in the experimental results.
The performance of the panels can be conveniently plotted in an Ashby chart that shows the energy absorption vs. strength in quasi-static (Fig. 4D) and impact (Fig. 4E) conditions. These charts highlight the superiority of the architected designs over the monolithic panels, and the positive effects of increasing the interlocking angle. Square-based and hexagon-based tiling showed similar performance, and in terms of combined strength and energy absorption, panels made of octahedral blocks showed the overall best performance in both quasi-static and impact conditions. The energy absorption for octahedral-block-based panels was ∼20× and ∼25× more than that of the monolithic panels in quasi-static and impact, respectively. When considering the actual energy dissipated by formation of cracks (see Methods for details of this analysis), the architectured panels were 35–50× more energy absorbent than monolithic ones. In addition to dissipating more energy, octahedral-block-based panels were, respectively, ∼13% and ∼18% stronger in quasi-static and impact. To explain this unexpected and remarkable increase of strength, it is useful to consider the failure mode for the different panel designs. The failure of the monolithic panels was governed by flexural cracks, a detrimental mode of failure which creates long cracks and ruins the load-carrying capability of the panel after the first impact (Fig. 3C). In contrast, flexural fracture was completely absent in the architectured panels. The finite size of each block reduces their flexural span, so that compared with the monolithic plates the flexural stress in individual blocks is much smaller for the same applied force [we recently observed this strengthening effect for stiff plates of finite size supported by a compliant substrate (43)]. An alternate but equivalent explanation is given by considering that most critical flexural stresses occur in tension in the monolithic plate. In the architectured plate these tension lines are disrupted by the interfaces between the blocks, which cannot carry tensile stresses. In addition to this effect, the distributed forces over three faces and an optimum locking angle increased the stability of individual blocks over long sliding distances.
The overall performance of the panels is governed by collective sliding mechanisms that involve geometry and frictional contact. Capturing these phenomena can be extremely costly computationally when multiple contacting elastic bodies are considered. Here we explored the possibility of capturing the geometrical effects using a single nondimensional number conveying the efficacy of the interlocking mechanism for any arbitrary design. Since the locking mechanism largely relies on elastic contact between the blocks, our model started by considering a single elastic block interlocked by rigid and stationary neighbors (SI Appendix, Fig. S5). We then considered a deformation mode where the center elastic block is pushed downward. The sliding between the blocks is considered frictionless but the geometrical interference between the sliding block and its rigid neighbors generates elastic deformations in the center block. In this simplified model, the elastic block represents the center block of the panel, which is subjected to a point force, and which undergoes large amounts of sliding but little rotation (at most 1.8° obtained from the 3D reconstruction results). We approximated the elastic energy U stored in the elastic block as a function of pushout distance (see Methods for details) and then we computed its maximum value Umax for each geometry. We finally normalized the result to capture the effects of architecture only and to produce a nondimensional “locking parameter”
A single nondimensional interlocking parameter, based on the strain energy stored in individual blocks during interlocking, can be used to predict (A) maximum force and (B) energy absorption. This approach may be effective in predicting the performance of other designs without the need for experiments or costly computational models. The width of the colored regions is representative of the variations in the experimental results.
Summary
In this study we have systematically explored topologically interlocked panels made of convex ceramic blocks. Under impact or transverse forces individual blocks slide and rotate, providing large deformations and toughness which leads to highly localized failure. In contrast, monolithic panels fracture in a brittle fashion and catastrophically by fragmentation. We have identified an architecture based on octahedral blocks which not only produces a 50-fold improvement in toughness compared with the monolithic form, but which is also ∼20% stronger. The interfaces between the blocks cannot sustain tensile stresses and therefore disrupt tensile stresses which are the usual cause of failure in flexural panels. The octahedral design provides added stability because contact forces are distributed over three contact surfaces and because among the geometries explored here, this geometry provides the most efficient balance between interlocking strength and surface damage from frictional contact stresses. There is an infinite number of other architectures to be explored, with some architectures which could lead to even higher performance (concave blocks, nonplanar faces). The exploration of this large design space is however difficult because experiments are lengthy and numerical models are computationally costly, if at all possible. To guide this exploration, we proposed a nondimensional interlocking parameter which we show correlates well with the performance of the panel. This parameter can be easily calculated for any geometry, which could greatly accelerate the exploration of new designs (a formal proof that the proposed interlocking parameter is a performance predictor for any arbitrary geometry remains to be established). The variability of our experimental results was similar for the monolithic and architectured panels. However, it is not clear which type of distribution the strength of the TIM panels follows, because we do not have enough experiments for each configuration. The strength of the architectured panels is partially governed by the onset of sliding between the blocks, which may follow a statistical law which is different from the weakest link (Weibull) statistics (44). More experiments and models are needed in this area, and questions related to statistics of failure for TIMs remain largely open.
An important design parameter for the architectured panels is the coefficient of friction at the interfaces, which may be finely tuned by adding roughness on the surfaces of the blocks (24). Interestingly, nature is well ahead of engineers in making use of architectured materials. Materials such as bone, teeth, or mollusc shells are also made of stiff building blocks of well-defined sizes and shapes, bonded together by deformable bioadhesives. The remarkable mechanical performance of these materials (12, 45) can suggest new types of 3D architectures. In addition, the building blocks in natural materials do not simply interact through contact and friction, but also through complex polymers with sacrificial bonds, dynamic cross-links, and viscous behaviors (31, 46) which could also serve as inspiration for interfaces in synthetic architectured materials. The segmentation of load-carrying structures into smaller elements joined by weaker interfaces is a counterintuitive approach to generate mechanical performance, but biological materials and recent studies on architectured materials show that it is indeed a powerful strategy to overcome brittleness while retaining strength. New combinations of properties in these architectured materials and systems can make them attractive for a variety of application including protective panels and armors, structural panels, or high-temperature structures.
Methods
Details of the derivation of the interlocking parameter (strain energy model) as well as experimental methods related to fabrication, testing, and data analysis can be found in SI Appendix.
Acknowledgments
This work was supported by Natural Sciences and Engineering Research Council of Canada and by the Fonds Quebecois de la Recherche sur la Nature et les Technologies. T.Z. was partially supported by a Summer Undergraduate Research in Engineering Award from McGill University.
Footnotes
- ↵1To whom correspondence should be addressed. Email: francois.barthelat{at}mcgill.ca.
Author contributions: M.M. and F.B. designed research; M.M. and T.Z. performed research; M.M., T.Z., and F.B. analyzed data; M.M. and F.B. prepared the figures; and M.M. and F.B. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1807272115/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- ↵
- Brechet Y,
- Embury JD
- ↵
- Lakes R
- ↵
- ↵
- Jaglinski T,
- Kochmann D,
- Stone D,
- Lakes RS
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Siegmund T,
- Barthelat F,
- Cipra R,
- Habtour E,
- Riddick J
- ↵
- Barthelat F
- ↵
- Dyskin AV,
- Estrin Y,
- Kanel-Belov AJ,
- Pasternak E
- ↵
- Dyskin AY,
- Estrin Y,
- Kanel-Belov AJ,
- Pasternak E
- ↵
- Estrin Y, et al.
- ↵
- Estrin Y,
- Dyskin AV,
- Pasternak E
- ↵
- Dugue M,
- Fivel M,
- Brechet Y,
- Dendievel R
- ↵
- Khandelwal S,
- Siegmund T,
- Cipra RJ,
- Bolton JS
- ↵
- Khandelwal S,
- Siegmund T,
- Cipra RJ,
- Bolton JS
- ↵
- Mirkhalaf M,
- Tanguay J,
- Barthelat F
- ↵
- Brocato M,
- Mondardini L
- ↵
- ↵
- Autruffe A, et al.
- ↵
- Djumas L,
- Simon GP,
- Estrin Y,
- Molotnikov A
- ↵
- Javan AR,
- Seifi H,
- Xu S,
- Ruan D,
- Xie Y
- ↵
- Mirkhalaf M,
- Barthelat F
- ↵
- Dyskin AV,
- Estrin Y,
- Pasternak E,
- Khor HC,
- Kanel-Belov AJ
- ↵
- Brugger C,
- Bréchet Y,
- Fivel M
- ↵
- Krause T, et al.
- ↵
- Naleway SE,
- Porter MM,
- McKittrick J,
- Meyers MA
- ↵
- Barthelat F,
- Yin Z,
- Buehler MJ
- ↵
- ↵
- Barthelat F
- ↵
- Begley MR, et al.
- ↵
- Malik I,
- Mirkhalaf M,
- Barthelat F
- ↵
- ↵
- ↵
- Johnson HH,
- Paris PC
- ↵
- ↵
- ↵
- ↵
- ↵
- Chintapalli RK,
- Mirkhalaf M,
- Dastjerdi AK,
- Barthelat F
- ↵
- ↵
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Engineering
- Biological Sciences
- Applied Biological Sciences