Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

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
Research Article

Resilient 3D hierarchical architected metamaterials

Lucas R. Meza, Alex J. Zelhofer, Nigel Clarke, Arturo J. Mateos, Dennis M. Kochmann, and Julia R. Greer
PNAS September 15, 2015 112 (37) 11502-11507; first published September 1, 2015; https://doi.org/10.1073/pnas.1509120112
Lucas R. Meza
aDivision of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alex J. Zelhofer
aDivision of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nigel Clarke
bNanotechnology Engineering, University of Waterloo, Waterloo, ON, Canada N2L 3G1
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Arturo J. Mateos
aDivision of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dennis M. Kochmann
aDivision of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julia R. Greer
aDivision of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: jrgreer@caltech.edu
  1. Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved August 11, 2015 (received for review May 8, 2015)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Fractal-like architectures exist in natural materials, like shells and bone, and have drawn considerable interest because of their mechanical robustness and damage tolerance. Developing hierarchically designed metamaterials remains a highly sought after task impaired mainly by limitations in fabrication techniques. We created 3D hierarchical nanolattices with individual beams comprised of multiple self-similar unit cells spanning length scales over four orders of magnitude in fractal-like geometries. We show, through a combination of experiments and computations, that introducing hierarchy into the architecture of 3D structural metamaterials enables the attainment of a unique combination of properties: ultralightweight, recoverability, and a near-linear scaling of stiffness and strength with density.

Abstract

Hierarchically designed structures with architectural features that span across multiple length scales are found in numerous hard biomaterials, like bone, wood, and glass sponge skeletons, as well as manmade structures, like the Eiffel Tower. It has been hypothesized that their mechanical robustness and damage tolerance stem from sophisticated ordering within the constituents, but the specific role of hierarchy remains to be fully described and understood. We apply the principles of hierarchical design to create structural metamaterials from three material systems: (i) polymer, (ii) hollow ceramic, and (iii) ceramic–polymer composites that are patterned into self-similar unit cells in a fractal-like geometry. In situ nanomechanical experiments revealed (i) a nearly theoretical scaling of structural strength and stiffness with relative density, which outperforms existing nonhierarchical nanolattices; (ii) recoverability, with hollow alumina samples recovering up to 98% of their original height after compression to ≥50% strain; (iii) suppression of brittle failure and structural instabilities in hollow ceramic hierarchical nanolattices; and (iv) a range of deformation mechanisms that can be tuned by changing the slenderness ratios of the beams. Additional levels of hierarchy beyond a second order did not increase the strength or stiffness, which suggests the existence of an optimal degree of hierarchy to amplify resilience. We developed a computational model that captures local stress distributions within the nanolattices under compression and explains some of the underlying deformation mechanisms as well as validates the measured effective stiffness to be interpreted as a metamaterial property.

  • hierarchical
  • nanolattices
  • structural metamaterial
  • recoverable
  • damage tolerance

Hierarchy is ubiquitous in the natural world; characterizing it, understanding its origins, and discovering its role in enhancing material properties are essential to designing new advanced materials (1⇓⇓–4). Natural structural materials, like Euplectella sponges, radiolarians, and bone, are exceptionally resilient against extreme mechanical environments and seem to draw their robustness from intricate mechanical networks that contain multiple levels of hierarchy (3⇓⇓⇓–7). Hierarchical engineered structures are used in modern architecture, with notable examples being the Eiffel tower and the Garabit viaduct (8); today, hierarchy is seen commonly in construction cranes and building scaffolding. Both natural and engineered structures use the concept of hierarchical design to minimize material use while optimizing structural integrity.

The hierarchical scale of a material is defined by its order, which represents the number of distinct structural length scales (2). Design principles and theories describing hierarchical structural materials exist (2, 9), and macroscopic second- and third-order 2D cellular solids, like honeycombs (10, 11) and corrugated core sandwich panels (12⇓–14), have been designed and tested experimentally. Theories that describe the design and optimization of 3D hierarchical trusses have been proposed (15⇓⇓–18); until recently, their fabrication presented a challenge. Modern fabrication breakthroughs have enabled the creation of 3D architected materials. These techniques include self-propagating photopolymer waveguides used to create Ni-based ultralight microlattices (19), microstereolithography capable of fabricating simultaneously strong and ultralight metal- and ceramic-based microlattices (20), and two-photon lithography direct laser writing, which enabled the creation of nanolattices with features sufficiently in the nanometer regime to capitalize on size-affected material properties in addition to structural effects (21⇓–23).

We report the fabrication, mechanical characterization, and computational analysis of hierarchical nanolattices made out of three different materials: (i) solid polymer IP-Dip (Nanoscribe GmbH), (ii) a core-shell composite with a polymer core and a 20-nm-thick Al2O3 coating, and (iii) hollow 20-nm-thick Al2O3. In situ nanomechanical deformation experiments were performed on a number of different geometries and showed that hierarchically designed nanolattices exhibit enhanced recoverability over simple periodic ones (22) and have a close to linear scaling of yield strength and stiffness with relative density. These findings were confirmed by computational predictions, and they are consistent with theoretical predictions for stretching-dominated solids (24). Simulations were performed to further elucidate the local stress distributions within the nanolattices, which confirm the effective experimental nanolattice response and help shed light on the distribution of loadbearing components that are responsible for the overall observed nanolattice performance.

Various sample geometries were tested to quantify the effect of geometry and relative density on mechanical behavior (Materials and Methods and SI Appendix, Tables S1–S3). Samples were designed using a recursive method, where a unit cell is patterned along the length of the beams of a self-similar unit cell, resulting in a fractal-like geometry (Fig. 1A and SI Appendix, Fig. S1). This design concept can be extended to create hierarchical metamaterials of any order with previously unobserved combinations of properties across multiple length scales—for example, high strength to weight ratios, tunable mass density, near-infinite bulk to shear modulus ratios (25, 26), and negative Poisson’s ratios (27, 28). The versatility of this fabrication method is shown in Fig. 1 and SI Appendix, Fig. S2, which show examples of designed and fabricated hierarchical nanolattice geometries.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Computer-aided design (CAD) and scanning electron microscopy (SEM) images of various hierarchical nanolattices show the versatility of the nanolattice fabrication technique. (A) CAD images illustrating the process of making a third-order hierarchical nanolattice. A zeroth-order repeating unit, an elliptical beam, is arranged into a first-order octahedron; it becomes the repeating unit for a second-order octahedron of octahedra, which is then arranged to create a third-order octahedron of octahedra of octahedra. (B, Upper, C, Upper, D, Upper, and E, Upper) CAD and (B, Lower, C, Lower, D, Lower, and E, Lower) SEM images of the various second-order samples. (Scale bars: 20 µm.) (F) SEM image of a second-order octahedron of octahedra lattice. (Scale bar: 50 µm.) (G) A zoomed-in image of the second-order octahedron of octahedra lattice showing the first-order repeating units that make up the structure. (Scale bar: 10 µm.) (H) SEM image of a third-order octahedron of octahedra of octahedra. (Scale bar: 25 µm.)

Results

Strength and Stiffness.

All fabricated samples were compressed uniaxially in an in situ nanoindenter (InSEM; Nanomechanics Inc.) to 50% strain before unloading. The load displacement data for each of the samples showed an initial linear region from which the effective loading stiffness was determined followed by an inelastic region with behavior that varied depending on the constituent material (Figs. 2 and 3). Failure strength and modulus were calculated as defined in Materials and Methods. Linear elastic finite element simulations were performed on representative nanolattices to predict stiffness and local stress distributions.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Compression experiments on second-order octahedron of octet half-cells with N = 15 and L = 8. (A) Image of the hollow 20-nm walled Al2O3 sample before compression. (B) Load displacement data that show compression to 50% strain. Inset corresponds to 50% strain. (C) Postdeformation image of the hollow sample. (D) Image of the composite polymer and 20-nm Al2O3 sample before compression. (E) Load displacement data that show compression to 65% strain. Inset corresponds to the sample after the occurrence of a strain burst. (F) Postdeformation image of the composite sample. (G) Image of the polymer sample before compression. (H) Load displacement data that show compression to 50% strain. Inset corresponds to 50% strain. (I) Postdeformation image of the polymer sample. (Scale bars: 20 µm.)

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Compression experiments on third-order octahedron of octahedra of octahedra half-cells with N = 5 and L = 8. (A) Image of the hollow 20-nm walled Al2O3 sample before compression. (B) Load-displacement data that show cyclic compression to 50% strain. Insets correspond to 50% strain at various loading cycles. (C) Postcompression image of the hollow sample. (D) Image of the composite polymer and 20-nm Al2O3 sample before compression. (E) Load-displacement data that show compression to 65% strain. Inset corresponds to the sample after the occurrence of a strain burst. (F) Postdeformation image of the composite sample. (G) Image of the polymer sample before compression. (H) Load displacement data that show compression to 50% strain. Inset corresponds to 50% strain. (I) Postdeformation image of the polymer sample. (Scale bars: 50 µm.)

Second-order half-cell samples of varying material compositions were tested with densities spanning over two orders of magnitude from ρ=0.30  kg m−3 to ρ=33.2  kg m−3. Individual material systems had relative densities that spanned more than one order of magnitude. Strength and modulus in architected materials scale with relative density asE=BEsρ¯β[1] andσy=Cσysρ¯γ,[2]where Es and σys are the constituent material’s Young’s modulus and yield strength, respectively; B and C are geometry-dependent proportionality constants, and β and γ are scaling constants (24). We found the experimentally measured stiffness in the hollow second-order half-cells to scale nearly linearly, with relative density as E=0.015 Ehρ¯ 1.04 and strength as σy=0.026 σyhρ¯ 1.17. We found similar scaling relations for polymer and composite second-order half-cell samples. Table 1 provides a full list of the constituent material properties used and the observed scaling parameters; polymer and composite properties can be found in SI Appendix, SI Materials, and Al2O3 properties were taken from refs. 29⇓⇓⇓–33.

View this table:
  • View inline
  • View popup
Table 1.

Material properties and proportionality and scaling constants for second-order half-cells as obtained from experiments and simulations

The definition of stiffness in experiments on structural metamaterials has been ambiguous because of various competing deformation mechanisms at each level of hierarchy ranging from elastic deformation to hinge-like phenomena to elastic–plastic buckling. We performed simulations on the elastic response of all tested second-order hierarchical nanolattices to calculate the stiffness of polymer, composite, and hollow ceramic nanolattices using an efficient finite element substructuring technique (Fig. 4 and SI Appendix, SI Materials) with the constituent properties listed in Table 1. Computed scaling exponents are included in Table 1 and show good agreement with experiments, differing by only 6.1%, 3.3%, and 3.9% for polymer, composite and hollow trusses, respectively. The absolute computed stiffnesses were, on average, 10.7% lower for polymer, 30.2% higher for composite, and 68.5% higher for hollow samples compared with experimental data (Fig. 5), which hints that geometric and/or material imperfections contribute significantly to a reduction in the effective stiffness.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Model flowchart showing truss and refined model generation. (A) Representative lattice geometry section. (B) Creation of a truss model lattice. (C) Example compression of truss model half-cell nanolattices. Stress is normalized by the maximum compressive stress in the sample, and stresses |σ|≤15% of the maximum stress have been grayed out to help illustrate the beams with high stresses. (D) Refined model creation process containing geometrically unique supernodes (SN) and superbeams (SB). (E) Example refined model half-cell nanolattice colored by unique geometry beam or node.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Comprehensive data plot of all tested hierarchical nanolattices. (A) Effective Young’s modulus of the hierarchical structures plotted against their relative density. Data are plotted for experimental (slope values are in bold) and refined node simulations (slope values are italicized) results. (B) Experimentally derived effective yield strength of the hierarchical nanolattices plotted against their relative density.

We found that the strength and stiffness of the equivalently dense hollow third-order half-cell samples were approximately a factor of two lower than those of second-order half-cells; under the same metric, polymer and composite samples had equivalent strength and stiffness for second- and third-order samples. Experimental and computational results for all second- and third-order half-cells are summarized in Fig. 5. Experiments on full second-order nanolattices revealed that the strength and stiffness align with second-order half-cell experiments (SI Appendix, SI Materials).

Failure, Deformation, and Recoverability.

The characteristic failure and postyield deformation of each sample were observed to correlate with (i) material system and (ii) architecture. We found that the material system (i.e., polymer vs. composite vs. hollow alumina) most strongly influenced the global deformation behavior and that the architecture directly affected the localization of failure within the hierarchical beams and the global recoverability of the samples. Examples of all three material systems for samples with two different orders of hierarchy are provided in Figs. 2 and 3 (Movies S1–S6).

The postyield deformation behavior of hollow ceramic samples was ductile-like with continuous serrated flow (Figs. 2 A–C and 3 A–C). Failure occurred through a combination of localized elastic beam buckling in the first-order beams, shell buckling in individual hollow tubes, and microcracking at the nodes. Of all material systems, the hollow ceramic samples had the highest average recovery, with samples recovering up to 85–98% of their original height after compressions exceeding 50% strain. Failure in the hollow samples initiated and localized primarily in individual first-order beams oriented along the axis of second-order beams, or referred to here as axially oriented beams. Buckling of these axially oriented first-order beams creates a compliant region in the higher-order beams that accommodates most of the ensuing displacement. We observed that first-order beams that did not buckle at the onset of failure remained completely intact during compression. On unloading, most of the beams within the cross-section of the higher-order beams remained intact, allowing significant global recovery. This effect is present in all hierarchical nanolattices and does not seem to correlate with hierarchical order. The recoverability mechanism described above is present in all material systems to varying degrees and most pronounced in hollow samples.

The postyield deformation of polymer samples was ductile and continuous, with samples recovering to ∼75–90% of their original height after unloading (Figs. 2 G–I and 3 G–I). Failure in polymer samples was governed primarily by buckling in first- and second-order beams. Most of the composite samples failed catastrophically by fracture of the second-order beams (Figs. 2 D–F and 3 D–F), except for the most slender samples, which showed some recoverability (SI Appendix, Fig. S6E).

Simulations reveal the local stress distribution and expose the load-carrying members of the hierarchical nanolattices. One prominent feature revealed by computations is that the first-order axially oriented beams, which comprise ∼8.1% of all of the beams in an individual sample, carry an average of 91% of the load, which corroborates the observation of failure initiating in these beams. Fig. 4C illuminates the regions of high local stress revealed by computations in second-order samples with varying degrees of slenderness. Samples with low slenderness have stress concentrations highly localized to the topmost first-order beams and are more likely to experience failure there. Samples with high slenderness have a more even stress distribution throughout their length and are more likely to have failure initiate away from the topmost region. This stress localization trend agrees well with experimental observations (SI Appendix, Fig. S6).

Fig. 3 A–C shows cyclic experiments on a third-order hollow ceramic half-cell, which revealed, after the initial loading cycle, that the stiffness dropped from 420 to 39  N/m and that the applied load at yield decreased from 0.77 to 0.089 mN. In the first loading cycle, the sample recovered to 96% of the original height on unloading; all subsequent cycles showed nearly complete 100% recovery to this initial deformed height. The load displacement data quickly reached a stable hysteretic cycling behavior, with minimal degradation after the second loading cycle (Fig. 3 A–C and Movie S6).

Discussion

Strength and Stiffness Scaling.

Compression experiments revealed the strength and stiffness of the second-order hierarchical half-cells to scale nearly linearly with relative density, which follows analytical and computational predictions for stretching-dominated cellular solids (24). These results show a factor of 1.5 improvement in the scaling relationship for strength and a factor of 1.6 improvement for modulus over nonhierarchical hollow Al2O3 nanolattices, and they are brought about by increasing the hierarchy from first to second order (22) (Fig. 5). Simulations reveal that the load in hierarchical nanolattices is carried primarily through axially oriented first-order beams in uniaxial compression, with little contribution to the deflection from bending. Analytical models for architected materials predict a linear scaling of strength and stiffness with density for structures with members that are primarily in uniaxial compression and tension (34). This linear scaling behavior has not been previously observed in hollow tube stretching-dominated nanolattices (22). We postulate that the near-linear strength and stiffness scaling observed in hierarchical nanolattices in this work arises from the combination of axially loaded first-order beams and the reduced effects of bending on global compliance. One key observation that arises from the nearly linear scaling of strength and stiffness is that the density of the samples can be varied indefinitely with only minor changes in the stiffness to weight and strength to weight ratios.

Simulations show that axially oriented beams, which comprise 8.1% of the macroscopic sample mass, carry 91% of the axial load (Fig. 4C). This low percentage of load carrying beams suggests that the remaining beams within the structure are underused. The underutilization of nonaxially oriented beams is likely a major factor in the observed reduction of the proportionality constants B and C (Eqs. 1 and 2) from what is predicted analytically for an ideal stretching-dominated solid, which have B=C≈0.3 (24) (Table 1). The negative impact of underused beams could be improved through better optimization of the hierarchical geometries.

Recoverability.

The underutilization of nonaxially oriented beams plays a significant role in the ability of hierarchical nanolattices to recover. Axially oriented beams undergo failure at the point of highest stress in a sample; in the absence of a catastrophic failure event, nonaxially oriented beams are able to remain intact. Under global compression, the undamaged nonaxially oriented first-order beams either (i) efficiently distribute strain through bending or (ii) undergo local elastic buckling to accommodate large global deformation without failure. A large number of elastically deformed first-order beams enables the hierarchical samples to globally recover. In the absence of residual strain in the buckled beams, such as is the case in hollow nanolattices, samples show excellent recovery behavior. The residual plastic strain in the buckled beams of polymer nanolattices impedes their ability to fully recover. In samples that undergo multiple loading cycles, most of the first-order beam buckling modes are activated in the first loading cycle; multiple loading cycles serve to reactivate the same buckling modes, which leads to near-perfect recoverability of the initially deformed samples. It may be possible to remove underused beams through better optimization of the hierarchical geometries, but such a reduction in the nonloadbearing beams may reduce the postyield recoverability by impeding the recovery mechanism.

Deformation and Failure.

Hierarchical nanolattices failed by a combination of different mechanisms: elastic beam buckling, shell wall bending, brittle fracture, and viscoplastic yielding. Failure initiates in beams through the weakest mechanism; equating the energy necessary to activate each mechanism enables the prediction of the dominant deformation mode. Previous work by Meza et al. (22) showed enhanced recoverability in hollow Al2O3 nanolattices when shell buckling acts as the preferential mechanism over fracture. Calculations performed on the hollow hierarchical samples in this work show that global (Euler) beam buckling acts as a dominant mechanism over shell buckling and that shell buckling acts preferentially over fracture (SI Appendix, SI Materials). This formulation implies that incipient failure initiates through beam buckling, and subsequent deformation is accommodated through shell buckling, leading to ductile-like behavior. This prediction agrees well with experimental observations (Figs. 2 A–C and 3 A–C). The brittleness of the Al2O3 causes large local strains to be relieved by either elastic buckling or fracture, and because shell buckling acts as a dominant mechanism, no residual strains are generated in the beams after unloading, which enables the samples to globally recover their original shape after unloading. In polymer samples, Euler buckling of the first-order beams is the dominant failure mechanism; the substantial postbuckled deformations generate nonnegligible residual stress in the beams, which impedes global recovery (Figs. 2 G–I and 3 G–I). The viscoelasticity of the polymer likely gives rise to the observed ductile load displacement behavior. The presence of the polymer–ceramic interface in composite samples suppresses buckling failure of either material and leads to a brittle fracture of the tube walls that results in the observed catastrophic collapse (Figs. 2 D–F and 3 D–F).

The large number of available failure modes in the hierarchical samples serves to reduce the yield strength. The activation of both beam and shell buckling causes the nanolattices to fail at a lower load than the yield limit, which decreases the global strength of samples and the proportionality constant C in Eq. 2 when these mechanisms are dominant. Shell wall bending in hollow samples generates stress concentrations near the nodes and along the beams. The negative impact of hollow nodes manifests itself in a twofold reduction in strength and stiffness in the third-order hollow samples, likely brought on by the higher-volume fraction of nodes in these samples.

The transition to buckling failure in nanolattices normally correlates with a significant reduction in the scaling of strength with relative density caused by the nonlinear dependence of buckling on slenderness in low-density samples (35). In this work, the length and aspect ratio of the first-order beams remains unchanged for different relative density samples; the relative density is modulated by the number of unit cells in a hierarchical beam N, which does not negatively affect the strength scaling with relative density.

Role of Imperfections.

The marked overestimation of the absolute stiffness obtained by the refined model simulations compared with experiments can be explained by the presence of geometric imperfections in the fabricated samples. Defects, like misaligned nodes, prebending of the beams, and variations in the wall thickness negatively impact the mechanical performance of nanolattices (21, 22). One dominant imperfection that was observed in all tested samples is sinusoidal waviness of the first-order beams caused by the external vibrations during the two-photon writing process. Simulations of the compression of beams with varying degrees of waviness and material compositions showed that, for a wave amplitude of 50 nm and a wavelength of 1 µm, the calculated effective stiffness of polymer beams decreased by 5%, that of the ceramic–polymer composites decreased by 32%, and that of the hollow ceramic beams decreased by 70%, which serve to illustrate the increased sensitivity to defects in the composite and hollow beams. The simulations overpredicted the stiffnesses of composite and hollow ceramic hierarchical lattices by 30.2% and 68.5%, respectively, which suggests that waviness-induced defects significantly contribute to this reduction. Model inaccuracy in underpredicting solid polymer stiffness by 10.7% can likely be attributed to the uncertainty in the polymer modulus.

We created hierarchical architected metamaterials with exceptional strength, stiffness, and damage tolerance and have architectures comparable with those of hard biomaterials, which derive their resilience from hierarchically arranged constituents. We show that adding a single level of hierarchy in nanolattices enhances the mechanical properties over simple periodic lattices, enables more efficient load transfer, and suppresses global failure but that the addition of additional levels of hierarchy does little to improve the mechanical properties and in some cases, degrades them. This ability to engineer material structure on the most fundamental length scales opens up a new design space where material properties—mechanical, thermal, electrical, photonic, etc.—can be controlled and tuned independently by properly choosing (i) atomic-level microstructure, (ii) critical material dimensions, and (iii) architecture.

Materials and Methods

Fabrication.

Hierarchical nanolattices were fabricated from solid polymer, ceramic–polymer core-shell composites, and hollow ceramic tubes. Solid polymer lattices were written in negative photoresist (IP-Dip 780) using the Photonic Professional TPL-DLW System (Nanoscribe GmbH). Core-shell composites were created by depositing a conformal coating of 20 nm Al2O3 onto the polymer nanolattices using atomic layer deposition in a Cambridge Nanotech S200 ALD System with H2O and trimethylaluminum (TMA) precursors. Hollow structures were made by removing the edges of the coated nanolattices using focused ion beam milling with the FEI Nova 200 Nanolab and then etching out the internal polymer structure with oxygen plasma using the Zepto Plasma Etcher (Diener GmbH), resulting in the hollow ceramic nanolattice. More details on this process can be found in ref. 22 (SI Appendix, SI Materials).

Two sets of samples were created: (i) half-unit cells (half-cells) and (ii) full nanolattices. Three different types of half-cells were fabricated and tested: a second-order octahedron of octahedra (Fig. 1B), a second-order octahedron of octets (Fig. 1C), and a third-order octahedron of octahedral of octahedra (Fig. 1H). For the second-order structures, two different base unit cell sizes (L) were used (8 and 12 µm), with three numbers of unit cells per fractal beam (N): 10, 15, and 20. For the third-order samples, two different octahedron of octahedral of octahedra were fabricated and tested with different geometric configurations: a unit cell with L = 3 µm and N = 10 and one with L = 8 µm and N = 5. Three different second-order octahedron of octahedra full nanolattices were fabricated and tested: a unit cell with L = 8 and N = 10, one with L = 6 and N = 15, and one with L = 4 and N = 20. The inherent limitations of the two-photon lithography fabrication methodology used to create the nanolattices render it impractical to create materials with more than three orders of structural hierarchy. A full list of samples and structural parameters can be found in SI Appendix, Tables S1–S3.

Experimental Methods.

In situ nanomechanical experiments were performed in a Quanta SEM (FEI Co.) using an InSEM Nanomechanical Module (Nanomechanics Inc.). All samples were compressed to 50% strain at a strain rate of 10−3 s−1. The structural stiffness was estimated based on the loading slope of the load displacement curve in the linear regime, and failure strength was taken to be the peak applied load before failure. We calculated the Young’s modulus by normalizing the measured loading stiffness by the sample height divided by the footprint area. Yield strength was determined by dividing the measured peak load by the sample footprint area. Scaling relations were obtained using an exponential best fit of the stiffness data; the scaling is computed as the average of the fits of the two second-order half-cell geometries: the octahedron of octets and the octahedron of octahedra.

Simulations.

Modeling of hierarchical nanolattices ultimately requires simulating thousands or millions of truss members while accurately resolving local deformation and failure mechanisms at the smallest scales. For elastic loading and stiffness predictions, simulations were performed based on a two-step computational strategy (using Abaqus and an in-house variational-based solid mechanics code). First, the characteristic elastic response of individual beams and representative lattice junctions was determined by finite element calculations using linear elastic shell and solid elements and a detailed representation of the CAD-generated geometry of individual truss members and junctions as shown in Fig. 4. This procedure was applied to monolithic, hollow, and composite beams and junctions. To dramatically reduce the number of degrees of freedom, we assume that the cross-sections of both ends of a beam deform rigidly through kinematic constraints. Thus, the deformation of each beam is reduced to depend only on the displacements and rotations of its end points. Similarly, truss junctions are reduced to depend only on the orientations of the attached truss members. Second, the thus-obtained load-displacement relations and stiffness matrices are input into an efficient simulation of the complete hierarchical structure based on the reduced degrees of freedom of all beam members and junctions in the lattice. The response of 12 second-order lattices of varying architecture was modeled by the aforementioned procedure. Constituent Young moduli were taken from Table 1. Poisson’s ratios of the polymer and alumina were chosen as νp=0.49 and νAl2O3=0.3, respectively.

To show the importance of the exact deformation mechanisms through fully resolved geometric models, we alternatively simulated the effective lattice response using bar elements (for solid beams) or shell elements (for hollow beams) without the aforementioned multiscale strategy and with the same base material properties. Calculations of first-order beam utilization in second-order beams were completed using the bar element model. Percentages of the axially aligned force carried in first-order beams relative to the whole second-order beam were calculated by sampling the middle 40% length of the second-order lattice beam. Axially aligned beam fractions were taken by dividing the number of axially oriented beams by the number of beams in the sample.

CAD Design/Relative Density.

CAD Solidworks models were used to determine the relative density of structures. The model uses measured voxel dimensions from each structure as design parameters in conjunction with base unit cell size and fractal number of structures. Each hierarchical repeating unit is patterned to construct higher-order structures, and similarly, higher-order unit cells were patterned to compose full nanolattices. The relative densities were calculated by evaluating the ratio of the volume of the structure by that of the solid from which the cell walls are made. A full illustration can be found in SI Appendix, Fig. S4.

Acknowledgments

The authors thank the Kavli Nanoscience Institute at Caltech for the availability of critical cleanroom facilities. Part of this work was carried out in the Lewis Group facilities at Caltech. The authors acknowledge financial support from the Defense Advanced Research Projects Agency under Materials with Controlled Microstructural Architecture (MCMA) Program Contract W91CRB-10-0305 (managed by J. Goldwasser), Institute for Collaborative Biotechnologies Grant W911NF-09-0001 from the US Army Research Office, and National Science Foundation Grant CMMI‐1234364.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: jrgreer{at}caltech.edu.
  • Author contributions: L.R.M., N.C., D.M.K., and J.R.G. designed research; L.R.M., A.J.Z., N.C., and A.J.M. performed research; A.J.Z. and D.M.K. contributed new reagents/analytic tools; L.R.M. and A.J.Z. analyzed data; and L.R.M., A.J.Z., A.J.M., D.M.K., and J.R.G. 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.1509120112/-/DCSupplemental.

View Abstract

References

  1. ↵
    1. Mandelbrot BB
    (1982) The Fractal Geometry of Nature (Freeman, San Francisco), 1st Ed
    .
  2. ↵
    1. Lakes R
    (1993) Materials with structural hierarchy. Nature 361(6412):511–515
    .
    OpenUrlCrossRef
  3. ↵
    1. Aizenberg J, et al.
    (2005) Skeleton of Euplectella sp.: Structural hierarchy from the nanoscale to the macroscale. Science 309(1985):275–278
    .
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Fratzl P,
    2. Weinkamer R
    (2007) Nature’s hierarchical materials. Prog Mater Sci 52(8):1263–1334
    .
    OpenUrlCrossRef
  5. ↵
    1. Robinson WJ,
    2. Goll RM
    (1978) Fine skeletal structure of the radiolarian Callimitra carolotae Haeckel. Micropaleontology 24(4):432–439
    .
    OpenUrlAbstract
  6. ↵
    1. Meyers MA,
    2. McKittrick J,
    3. Chen P-Y
    (2013) Structural biological materials: Critical mechanics-materials connections. Science 339(6121):773–779
    .
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Nalla RK,
    2. Kinney JH,
    3. Ritchie RO
    (2003) Mechanistic fracture criteria for the failure of human cortical bone. Nat Mater 2(3):164–168
    .
    OpenUrlCrossRefPubMed
  8. ↵
    1. Sundaram MM,
    2. Ananthasuresh GK
    (2009) Gustave Eiffel and his optimal structures. Resonance 14(9):849–865
    .
    OpenUrlCrossRef
  9. ↵
    1. Farr RS,
    2. Mao Y
    (2008) Fractal space frames and metamaterials for high mechanical efficiency. Europhys Lett 84(1):14001
    .
    OpenUrlCrossRef
  10. ↵
    1. Ajdari A,
    2. Jahromi BH,
    3. Papadopoulos J,
    4. Nayeb-Hashemi H,
    5. Vaziri A
    (2012) Hierarchical honeycombs with tailorable properties. Int J Solids Struct 49(11):1413–1419
    .
    OpenUrlCrossRef
  11. ↵
    1. Oftadeh R,
    2. Haghpanah B,
    3. Vella D,
    4. Boudaoud A,
    5. Vaziri A
    (2014) Optimal fractal-like hierarchical honeycombs. Phys Rev Lett 113(10):104301
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Wadley HNG
    (2006) Multifunctional periodic cellular metals. Philos Trans A Math Phys Eng Sci 364(1838):31–68
    .
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Kooistra GW,
    2. Deshpande V,
    3. Wadley HNG
    (2007) Hierarchical corrugated core sandwich panel concepts. J Appl Mech 74(2):259–268
    .
    OpenUrlCrossRef
  14. ↵
    1. Li G,
    2. Fang Y
    (2014) Failure mode analysis and performance optimization of the hierarchical corrugated truss structure. Adv Mech Eng 6:1–11
    .
    OpenUrl
  15. ↵
    1. Farr RS
    (2007) Fractal design for an efficient shell strut under gentle compressive loading. Phys Rev E Stat Nonlin Soft Matter Phys 76(5 Pt 2):056608
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Rayneau-Kirkhope D,
    2. Mao Y,
    3. Farr R,
    4. Segal J
    (2012) Hierarchical space frames for high mechanical efficiency: Fabrication and mechanical testing. Mech Res Commun 46:41–46
    .
    OpenUrlCrossRef
  17. ↵
    1. Rayneau-Kirkhope D,
    2. Mao Y,
    3. Farr R
    (2013) Optimization of fractal space frames under gentle compressive load. Phys Rev E Stat Nonlin Soft Matter Phys 87(6):063204
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Rayneau-Kirkhope D,
    2. Mao Y,
    3. Farr R
    (2012) Ultralight fractal structures from hollow tubes. Phys Rev Lett 109(20):204301
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. Schaedler TA, et al.
    (2011) Ultralight metallic microlattices. Science 334(6058):962–965
    .
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Zheng X, et al.
    (2014) Ultralight, ultrastiff mechanical metamaterials. Science 344(6190):1373–1377
    .
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Meza LR,
    2. Greer JR
    (2014) Mechanical characterization of hollow ceramic nanolattices. J Mater Sci 49(6):2496–2508
    .
    OpenUrlCrossRef
  22. ↵
    1. Meza LR,
    2. Das S,
    3. Greer JR
    (2014) Strong, lightweight, and recoverable three-dimensional ceramic nanolattices. Science 345(6202):1322–1326
    .
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Bauer J,
    2. Hengsbach S,
    3. Tesari I,
    4. Schwaiger R,
    5. Kraft O
    (2014) High-strength cellular ceramic composites with 3D microarchitecture. Proc Natl Acad Sci USA 111(7):2453–2458
    .
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Fleck NA,
    2. Deshpande VS,
    3. Ashby MF
    (2010) Micro-architectured materials: Past, present and future. Proc R Soc Lond A Math Phys Sci 466(2121):2495–2516
    .
    OpenUrlCrossRef
  25. ↵
    1. Kadic M,
    2. Bückmann T,
    3. Stenger N,
    4. Thiel M,
    5. Wegener M
    (2012) On the practicability of pentamode mechanical metamaterials. Appl Phys Lett 100(19):191901
    .
    OpenUrlCrossRef
  26. ↵
    1. Bückmann T,
    2. Thiel M,
    3. Kadic M,
    4. Schittny R,
    5. Wegener M
    (2014) An elasto-mechanical unfeelability cloak made of pentamode metamaterials. Nat Commun 5:4130
    .
    OpenUrlPubMed
  27. ↵
    1. Krödel S,
    2. Delpero T,
    3. Bergamini A,
    4. Ermanni P,
    5. Kochmann DM
    (2013) 3D auxetic microlattices with independently controllable acoustic band gaps and quasi-static elastic moduli. Adv Eng Mater 16(4):357–363
    .
    OpenUrl
  28. ↵
    1. Bückmann T, et al.
    (2012) Tailored 3D mechanical metamaterials made by dip-in direct-laser-writing optical lithography. Adv Mater 24(20):2710–2714
    .
    OpenUrlCrossRefPubMed
  29. ↵
    1. Berdova M, et al.
    (2014) Mechanical assessment of suspended ALD thin films by bulge and shaft-loading techniques. Acta Mater 66:370–377
    .
    OpenUrlCrossRef
  30. ↵
    1. Ilic B,
    2. Krylov S,
    3. Craighead HG
    (2010) Young’s modulus and density measurements of thin atomic layer deposited films using resonant nanomechanics. J Appl Phys 108(4):1–11
    .
    OpenUrlFREE Full Text
  31. ↵
    1. Jen S-H,
    2. Bertrand JA,
    3. George SM
    (2011) Critical tensile and compressive strains for cracking of Al2O3 films grown by atomic layer deposition. J Appl Phys 109(8):084305
    .
    OpenUrlCrossRef
  32. ↵
    1. Tripp MK, et al.
    (2006) The mechanical properties of atomic layer deposited alumina for use in micro- and nano-electromechanical systems. Sens Actuators A Phys 130-131:419–429
    .
    OpenUrl
  33. ↵
    1. Bauer J, et al.
    (2015) Push-to-pull tensile testing of ultra-strong nanoscale ceramic–polymer composites made by additive manufacturing. Extreme Mech Lett 3:105–112
    .
    OpenUrlCrossRef
  34. ↵
    1. Deshpande VS,
    2. Fleck NA,
    3. Ashby MF
    (2001) Effective properties of the octet-truss lattice material. J Mech Phys Solids 49(8):1747–1769
    .
    OpenUrlCrossRef
  35. ↵
    1. Valdevit L,
    2. Godfrey SW,
    3. Schaedler TA,
    4. Jacobsen AJ,
    5. Carter WB
    (2013) Compressive strength of hollow microlattices: Experimental characterization, modeling, and optimal design. J Mater Res 28(17):2461–2473
    .
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Resilient 3D hierarchical architected metamaterials
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Hierarchy in 3D metamaterials brings resilience
Lucas R. Meza, Alex J. Zelhofer, Nigel Clarke, Arturo J. Mateos, Dennis M. Kochmann, Julia R. Greer
Proceedings of the National Academy of Sciences Sep 2015, 112 (37) 11502-11507; DOI: 10.1073/pnas.1509120112

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Hierarchy in 3D metamaterials brings resilience
Lucas R. Meza, Alex J. Zelhofer, Nigel Clarke, Arturo J. Mateos, Dennis M. Kochmann, Julia R. Greer
Proceedings of the National Academy of Sciences Sep 2015, 112 (37) 11502-11507; DOI: 10.1073/pnas.1509120112
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 112 (37)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Engineering

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Large piece of gold
News Feature: Tracing gold's cosmic origins
Astronomers thought they’d finally figured out where gold and other heavy elements in the universe came from. In light of recent results, they’re not so sure.
Image credit: Science Source/Tom McHugh.
Dancers in red dresses
Journal Club: Friends appear to share patterns of brain activity
Researchers are still trying to understand what causes this strong correlation between neural and social networks.
Image credit: Shutterstock/Yeongsik Im.
White and blue bird
Hazards of ozone pollution to birds
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490