Next Article in Journal
On Coolant Flow Rate-Cutting Speed Trade-Off for Sustainability in Cryogenic Milling of Ti–6Al–4V
Previous Article in Journal
Photobiomodulation with Red and Near-Infrared Light Improves Viability and Modulates Expression of Mesenchymal and Apoptotic-Related Markers in Human Gingival Fibroblasts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Interfacial Dynamic Modulus of Biological Composites

1
Department of Mechanical Engineering, Ben-Gurion University of the Negev, Beer Sheva 84105, Israel
2
Nuclear Research Center-Negev, P.O. Box 9001, Beer-Sheva 84190, Israel
*
Author to whom correspondence should be addressed.
Authors with equal contributions.
Materials 2021, 14(12), 3428; https://doi.org/10.3390/ma14123428
Submission received: 11 May 2021 / Revised: 9 June 2021 / Accepted: 16 June 2021 / Published: 21 June 2021
(This article belongs to the Section Biomaterials)

Abstract

:
Biological composites (biocomposites) possess ultra-thin, irregular-shaped, energy dissipating interfacial regions that grant them crucial mechanical capabilities. Identifying the dynamic (viscoelastic) modulus of these interfacial regions is considered to be the key toward understanding the underlying structure–function relationships in various load-bearing biological materials including mollusk shells, arthropod cuticles, and plant parts. However, due to the submicron dimensions and the confined locations of these interfacial regions within the biocomposite, assessing their mechanical characteristics directly with experiments is nearly impossible. Here, we employ composite-mechanics modeling, analytical formulations, and numerical simulations to establish a theoretical framework that links the interfacial dynamic modulus of a biocomposite to the extrinsic characteristics of a larger-scale biocomposite segment. Accordingly, we introduce a methodology that enables back-calculating (via simple linear scaling) of the interfacial dynamic modulus of biocomposites from their far-field dynamic mechanical analysis. We demonstrate its usage on zigzag-shaped interfaces that are abundant in biocomposites. Our theoretical framework and methodological approach are applicable to the vast range of biocomposites in natural materials; its essence can be directly employed or generally adapted into analogous composite systems, such as architected nanocomposites, biomedical composites, and bioinspired materials.

1. Introduction

Load-bearing biocomposites are typically structured as arrays of rigid and predominantly elastic reinforcing elements (e.g., biominerals or crystalline biopolymers), which are connected by a more compliant and energy-dissipating matrix material (e.g., proteins or hemicellulose) through submicron length, compositionally graded, and irregularly-shaped interfacial regions [1,2,3,4,5,6,7,8]. The effective dynamic (viscoelastic) modulus of these interfacial regions provides the biocomposites’ diverse mechanical functions, including adsorbing impacts, detaining cracks, and filtering mechanical signals [9,10,11,12,13,14,15,16,17]. Identifying the interfacial dynamic modulus of biocomposites is a long-standing objective of biomaterial science research [18], it is considered the keystone toward understanding the fundamental structure–function relationships in various biocomposite systems [19,20,21,22,23,24,25].
Nanomechanical testing methods, e.g., nanoindentation and nanoscale dynamic mechanical analysis (DMA), are the benchmark approaches to identifying the mechanical characteristics of the interfacial regions in biocomposites [18,26,27,28,29,30,31,32]. These methods apply local contact loadings to certain locations within the interfacial region, analyze their mechanical response upon static or harmonic forces, and determine the elastic stiffness and viscous damping characteristics of the underlying reinforcement or matrix materials within the interfacial region. These underlying material characteristics link to the mechanical response of the interfacial region as a whole via shear-lag mechanisms, which transfer axial loads between adjacent reinforcements through the tensile-shear loadings of their intermediate matrix material [33,34,35,36]. Recent studies on planar interfacial morphologies (e.g., staggered, triangular, and trapezoid) have derived analytical relationships between the overall dynamic modulus of the interfacial region and those of its underlying reinforcement and matrix materials [37,38,39,40,41]. However, these analytical relationships cannot account for non-planar, irregularly-shaped, or unmarked interfacial morphologies (as commonly present in natural materials), which must be characterized through direct interfacial experiments. Practically, such direct interfacial experiments are nearly impossible due to the small dimensions and the confined locations of the interfacial regions within the biocomposite complex. The interfacial mechanical characteristics must be analytically extracted from far-field experiments on a larger-scale biocomposite segment [42,43,44,45,46]. Even small variations in the interfacial characteristics, i.e., material properties or relative content within the biocomposite, may substantially affect the mechanical response of the biocomposite segment [47,48,49]. Establishing feasible methodologies to approach the interfacial dynamic modulus of biocomposites from far-field experiments is a pending challenge of both biological and synthetic nanomaterial science [50,51,52,53,54].
Here, we employ composite-mechanics modeling, theoretical approximations, and numerical simulations to identify simple analytical relationships between the dynamic modulus (i.e., modulus magnitude and loss coefficient) of an interfacial region within a biocomposite to its larger-scale, enclosing biocomposite segment. With these relationships, we propose an analytical methodology that allows for the back-calculation (linear scaling) of the interfacial dynamic modulus from far-field DMA results on the biocomposite segment. Finally, we demonstrate the usability and adequacy of our methodology via numerical experiments on a class of sutural interfaces that are abundant in natural materials.

2. Analytical Relationships for the Interface–Biocomposite Dynamic Moduli

We considered a biocomposite segment—an isolated specimen from a larger-scale biocomposite complex—that includes adjacent elastic reinforcements connected by an energy-dissipating, viscoelastic matrix in various possible structural forms (Figure 1a). We identified the interfacial region (of length L i ) within the biocomposite segment (of length L c ) as the region that is different from the pristine reinforcements in terms of material properties or architectural characteristics. We characterized the elastic behavior of the reinforcements via Young’s modulus E f and the viscoelastic behavior of the whole biocomposite segment via the dynamic (complex) modulus E c * = E c e j δ c , where E c is the modulus magnitude of the biocomposite, and tan δ c is the loss coefficient of the biocomposite. We characterized the viscoelastic behavior of the interfacial region by an effective dynamic modulus E i * = E i e j δ i that compiles the properties of its underlying reinforcement and matrix materials via the interfacial shear-lag mechanisms ( E i and tan δ i are the effective modulus magnitude and the effective loss coefficient of the interfacial region, respectively).
To connect the interfacial dynamic modulus ( E i * ) to that of the whole biocomposite segment ( E c * ), we modeled the biocomposite segment as a reinforcement-interface-reinforcement sequence (Figure 1b) and employed Reuss’s model from classical composite mechanics (adapted via the correspondence principle) [49,55], which aptly describes various analogous biocomposite configurations, both at the macromolecular level (e.g., nanofibrils) [44], the nanocomposite level (e.g., nanofibril arrays enriched with biominerals) [3], and the microcomposite level (e.g., lamellar architectures) [56]. Consequentially, we expressed E i * as follows [49,55]:
E i * = [ 1 / ( L i / L c )   E c * 1 / ( L i / L c ) 1 E f ] 1
Practically, the biocomposite dynamic modulus ( E c * ) is characterized by small-scale DMA experiments, e.g., atomic force microscopy or dynamic nanoindentation, which yield E c and tan δ c from the amplitude ratio and the phase shift between the harmonic stress and strain signals, respectively. Thus, we wished to draw direct analytical connections between these experimental measures ( E c and tan δ c ) to the corresponding characteristics of the interfacial region ( E i and tan δ i ). By employing standard analytical steps (see Section S1 in Supporting Information) into Equation (1), we obtained the following equivalent equations:
E i = E c L i L c 1 1 2 ( 1 L i L c ) ( E c E f ) 1 1 + tan 2 δ c + ( E c E f ) 2 ( 1 L i L c ) 2
tan δ i = tan δ c 1 1 E c E f ( 1 L i L c ) 1 + tan 2 δ c
Next, we considered the typical mechanical characteristics of biocomposites, E c / E f 1 / 4 and tan δ c 1 / 2 (which are mostly far below these bounds; see [18] and the references therein), introduced analytical approximations into Equations (2) and (3), and obtained the following relationships (see section S1 in Supporting Information):
E i = k E E c ;   k E = L i / L c 1 ( E c / E f ) ( 1 L i / L c )
tan δ i = k δ tan δ c ;   k δ = 1 1 ( E c / E f ) ( 1 L i / L c )
Equations (4) and (5) show that the interface–biocomposite characteristics, E i E c and tan δ i tan δ c , link via linear scaling. The scaling factors, k E and k δ , depend only on the biocomposite-to-reinforcement modulus ratio ( E c / E f ) and the relative length of the interface within the biocomposite ( L i / L c ), and exclude the loss coefficient of the biocomposite ( tan δ c ). Specifically, when the biocomposite-to-reinforcement modulus ratio is sufficiently small ( E c / E f < 1 / 10 ), or when the relative length of the interfacial region is sufficiently large ( 1 / 2 < L i / L c ), the denominators of k E and k δ approaches unity, and Equations (4) and (5) further simplify into E i L i / L c E c and tan δ i tan δ c . Notably, as k E and k δ are nondimensional, the E i E c and tan δ i tan δ c relationships in Equations (4) and (5) are independent of the absolute length scale and the mechanical characteristics of the biocomposite. Consequentially, these E i E c and tan δ i tan δ c relationships are generally applicable for the broad dimensional range of biocomposites in natural materials—including macromolecular, nanoscale, and microscale biocomposites—and their diverse mechanical characteristics, i.e., from highly rigid to substantially compliant and from nearly elastic to prominently viscous.
We employed numerical dynamic mechanical analysis via finite element (FE) simulations (Abaqus 6.12, T2D2H elements) to verify the interface–biocomposite relationships in Equations (4) and (5). We analyzed a biocomposite segment (length L c ) that was comprised of a pair of elastic reinforcements (elastic modulus E f ) and an intermediate viscoelastic interface (length L i and dynamic modulus E i * = E i e j δ m ). In each simulation, we applied harmonic strain loadings on the lateral edges of the biocomposite segment, probed the resulting harmonic stresses on these edges, and extracted the modulus magnitude and the loss coefficient of the biocomposite ( E c and tan δ c ) from the amplitude ratio and the phase shift between the stress and strain signals. We analyzed the simulation models with a wide range of input parameters that enclosed the typical characteristics of natural materials (see Tables S1 and S2 in the Supporting Information), we plotted their resultant pairs, E i E c pairs (Figure 2), and tan δ i tan δ c pairs (Figure 3). For each simulation model, we calculated the biocomposite-to-reinforcement modulus ratio ( E c / E f ), the relative length of the interfacial region within the biocomposite ( L i / L c ), and the resultant interface–biocomposite scaling factors ( k E and k δ , respectively). Then, we plotted the theoretical relationships in Equations (4) and (5) with the E i E c , and tan δ i tan δ c pair sets from the numerical simulations (Figure 2 and Figure 3) which showed excellent correspondence for the entire range of input parameters analyzed.

3. Assessing the Interfacial Dynamic Modulus from a Far-Field Dynamic Mechanical Analysis

3.1. Methodological Approach

We have used the interface–biocomposite relationships in Equations (4) and (5) to propose the following methodological approach to back-calculate the interfacial dynamic modulus from far-field DMA measurements on the biocomposite segment itself.
  • Step 1: Isolate a testing segment (length L c ) from the biocomposite complex and use microscopy observations to identify its underlying interfacial region (length L i ).
  • Step 2: Apply DMA testing on a biocomposite segment and quantify its modulus magnitude and loss coefficient ( E c and tan δ c ).
  • Step 3: Use nanomechanical testing (or the literature data) to determine the elastic modulus of the reinforcements ( E f ) outside the interfacial region.
  • Step 4: Calculate the interface–biocomposite scaling factors ( k E and k δ ), and use them to back-calculate the modulus magnitude and the loss coefficient of the interfacial region from the corresponding biocomposite characteristics ( E i = k E E c , and tan δ i = k δ tan δ c ).
Practically, the estimations of L i and E f (Steps 1 and 3 above) typically incorporate certain deviation ranges, which arise from the resolution limitation of the microscopy analysis, and from the inherent variability of the nanomechanical analysis, respectively, such that L i L i ± Δ L i , and E f E f ± Δ E f . Accordingly, the interface–biocomposite scaling factors also include deviation ranges k E k E ± Δ k E , and k δ k δ ± Δ k δ . To identify these deviations, we calculated the first-order differentials of Equations (4) and (5), namely Δ k E = d k E / d L i   Δ L i + d k E / d E f Δ E f , and Δ k δ = d k δ / d L i Δ L i + d k δ / d E f Δ E f , and we expressed them analytically as follows:
Δ k E = k E [ ( 1 A ) Δ L i L i B Δ E f E f ]
Δ k δ = k δ [ A Δ L i L i B Δ E f E f ]
where
A = L i L c E c E f 1 C   , B = E c E f [ 1 L i L c ] 1 C   , C = 1 ( E c E f ) ( 1 L i L c )
Notably, when the interfacial region occupies large portions of the biocomposite segment ( L i / L c 1 ), the coefficients B A , and the k E and k δ deviations reduce into Δ k E k E ( 1 A ) Δ L i / L i and Δ k δ k δ A Δ L i / L i ; both k E and k δ are insensitive to E f deviations. Moreover, when the biocomposite-to-reinforcement modulus ratio is sufficiently small ( E c / E f < 1 / 10 ), the coefficients A , B 1 , and the k E and k δ deviations, further simplify into Δ k E k E Δ L i / L i and Δ k δ k δ ; the k E deviations are proportional to the L i deviations, and the k δ deviations are approximately negligible.
Next, we demonstrate our approach on biocomposites with zigzag-shaped sutural interfaces and illustrate the effect of their zigzag angles on the modulus magnitude and the loss coefficient of these interfacial regions.

3.2. Example: Sutural Interfaces

Sutural interfaces are abundant in natural materials, such as dermal armors [57,58], bird beaks [59], and seed coats [60], and they serve as locally flexible and energy-dissipating regions within a much harder and fragile bulk material. These sutural interfaces commonly view periodic zigzag geometries (period length λ i , zigzag angle θ ), which connect adjacent elastic reinforcements (elastic modulus E f ) via a viscoelastic matrix layer (length L m , dynamic modulus E m * = E m e j δ m ), and yield an interfacial region of length L i = L m + λ i / 2 tan θ (Figure 4). The zigzag angle of these sutural interfaces plays a major role in the various mechanical functions of the sutural interfaces, e.g., its elastic stiffness, failure strength, dumping capability, and impact resistance [11,36,61,62]. Here, we use numerical simulations to demonstrate the effect of the zigzag angle of the sutural interface on its effective dynamic modulus [41]. We focus on sutural interfaces with substantial matrix portions ( L m / λ i = 1 / 4 ) and small-to-moderate zigzag angles ( θ = 0 75 ° ), which complement the framework of corresponding analytical models and extract the effective modulus magnitude and loss coefficient of these interfacial regions from a far-field DMA analysis of a biocomposite segment via the above methodological approach (Section 3.1). We used plane stress FE simulations (Abaqus 6.12, CPS4 elements) with a substantially viscoelastic matrix (loss coefficient tan δ m = 1 / 2 ), and much stiffer, elastic reinforcements ( E f / E m = 10 ). For simplicity, we considered the matrix and reinforcements to be nearly incompressible ( E f / G f = E m * / G m * = 2.5 ). For each zigzag angle analyzed ( θ ), we realized the geometry of the sutural interface ( L m / λ i = 1 / 4 , and L i / L m = 1 + 2 tan θ ), and set the remaining length of the reinforcements to ensure that the overall interface-to-composite length ratio was kept constant ( L i / L c = 0.1 ). For each case, we applied harmonic strain loadings on the lateral edges of the biocomposite (upper and lower edges were free), measured the resulting harmonic stresses on these edges, and extracted the modulus magnitude and the loss coefficient of the biocomposite ( E c ( θ ) / E m and tan δ c ( θ ) ) from the amplitude ratio and the phase shift between the stress and strain signals. Then, we calculated the scaling factors k E and k δ factors, and back-calculated the modulus magnitude and the loss coefficient of the sutural interface, E i ( θ ) / E m and tan δ i ( θ ) ), via Equations (4) and (5), respectively (Figure 5a,b, Table S3 in Supporting Information). To verify our results, we performed a complementary analysis on the isolated sutural interfaces (i.e., L i / L c = 1 ) that directly yielded E i ( θ ) / E m  and tan δ i ( θ ) (i.e., without back-calculations); evidently, our estimations from the far-field DMA analysis on the biocomposite segments are in good agreement with the direct analysis on the isolated sutural interfaces for the entire θ range analyzed. Our results show that the zigzag angle strongly affects the modulus magnitude and the loss coefficient of the sutural interface (Figure 5a,b), and are in line with the theoretical trends of recent analytical studies [41]. As the zigzag angle approaches zero, the interfacial region is almost flat and mostly comprises the matrix layer ( L i / L m ~ 1 ), and its modulus magnitude and loss coefficient approach that of the pristine matrix material E i ( θ 0 ) ~ E m and tan δ i ( θ 0 ) ~ tan δ m . Conversely, as the zigzag angle increases, the interfacial region occupies greater portions of the stiffer elastic reinforcements (in addition to the matrix layer), which results in a progressive increase in the modulus magnitude of the interfacial region E i ( θ ) and a complementary decrease in its loss coefficient tan δ i ( θ ) . Thus, natural sutural interfaces with similar matrix and reinforcement characteristics can achieve different biomechanical functions merely through morphological adaptations of their zigzag shapes. Mild zigzag shapes will generate interfacial regions with greater flexibility and higher energy-dissipation capabilities, whereas acute zigzag shapes will generate interfacial regions with greater robustness and higher energy-storing capabilities.

4. Conclusions

The interfacial dynamic modulus of biocomposites dominates the structure–function relationships in various organisms, including vertebrates (e.g., mammals, birds, and fishes), invertebrates (e.g., insects, arachnids, and mollusks), and plants. Nevertheless, due to the small dimensions and the confined locations of these interfacial regions, measuring their direct mechanical characterizations is nearly impossible. In this study, we established compact analytical formulae that link the modulus magnitude and the loss coefficient of the interfacial region to those of its enclosing, large-scale, biocomposite segment. We used these formulae to propose an analytical and experimental methodology that yielded the interfacial characteristics via back-calculations, from a feasible, far-field DMA analysis of the biocomposite itself and demonstrated it on zigzag-shaped sutural interfaces. From a broader perspective, our approach can also be used to analyze the interfacial characteristics of advanced engineering materials, such as bioinspired composites, nanocomposites, and electromechanical devices [63,64,65,66].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14123428/s1, Table S1: Summary of DMA simulations results, shown in Figure 2 in the main text; Table S2: Summary of DMA simulations results, shown in Figure 3 in the main text; Table S3: Summary of the DMA simulation results for the zigzag-shaped interfaces and the corresponding back-calculations of E i  and tan δ i via Equations (4) and (5).

Author Contributions

Y.S., A.Y.U. and B.B.-O. conducted the analytical formulations, Y.S., A.Y.U. and O.B. performed the numerical simulations. All authors analyzed and discussed the results and wrote the manuscript. Y.S. and A.Y.U. provided equal contributions. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Israel Science Foundation (ISF) grant No. 1429/16 and by the German research foundation (DFG) grant GZ: PO 1725/7-1; AOBJ: 639691. The authors acknowledge the generous support of the Pearlstone Center for Aeronautical Engineering Studies.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or supplementary material.

Acknowledgments

The authors acknowledge the generous support of the Pearlstone Center for Aeronautical Engineering Studies.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Meyers, M.A.; Chen, P.Y.; Lin, A.Y.M.; Seki, Y. Biological materials: Structure and mechanical properties. Prog. Mater. Sci. 2008, 53, 1–206. [Google Scholar] [CrossRef] [Green Version]
  2. Dunlop, J.W.; Weinkamer, R.; Fratzl, P. Artful interfaces within biological materials. Mater. Today 2011, 14, 70–78. [Google Scholar] [CrossRef]
  3. Bar-On, B.; Wagner, H.D. Structural motifs and elastic properties of hierarchical biological tissues—A review. J. Struct. Biol. 2013, 183, 149–164. [Google Scholar] [CrossRef]
  4. Barthelat, F.; Yin, Z.; Buehler, M.J. Structure and mechanics of interfaces in biological materials. Nat. Rev. Mater. 2016, 1, 16007. [Google Scholar] [CrossRef]
  5. Shelef, Y.; Bar-On, B. Surface protection in bio-shields via a functional soft skin layer: Lessons from the turtle shell. J. Mech. Behav. Biomed. Mater. 2017, 73, 68–75. [Google Scholar] [CrossRef]
  6. Serrano, C.V.; Leemreize, H.; Bar-On, B.; Barth, F.G.; Fratzl, P.; Zolotoyabko, E.; Politi, Y. Ordering of protein and water molecules at their interfaces with chitin nano-crystals. J. Struct. Biol. 2016, 193, 124–131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Bar-On, B. On the form and bio-mechanics of venom-injection elements. Acta Biomater. 2019, 85, 263–271. [Google Scholar] [CrossRef]
  8. Ampaw, E.; Owoseni, T.A.; Du, F.; Pinilla, N.; Obayemi, J.; Hu, J.; Nigay, P.-M.; Nzihou, A.; Uzonwanne, V.; Kana, Z.; et al. Compressive deformation and failure of trabecular structures in a turtle shell. Acta Biomater. 2019, 97, 535–543. [Google Scholar] [CrossRef]
  9. Ji, B.; Gao, H. Mechanical properties of nanostructure of biological materials. J. Mech. Phys. Solids 2004, 52, 1963–1990. [Google Scholar] [CrossRef]
  10. Smith, B.L.; Schäffer, T.E.; Viani, M.; Thompson, J.B.; Frederick, N.A.; Kindt, J.; Belcher, A.; Stucky, G.D.; Morse, D.E.; Hansma, P.K. Molecular mechanistic origin of the toughness of natural adhesives, fibres and composites. Nature 1999, 399, 761. [Google Scholar] [CrossRef]
  11. Lee, N.; Williams, L.N.; Mun, S.; Rhee, H.; Prabhu, R.; Bhattarai, K.R.; Horstemeyer, M.F. Stress wave mitigation at suture interfaces. Biomed. Phys. Eng. Express 2017, 3, 035025. [Google Scholar] [CrossRef]
  12. Erko, M.; Younes-Metzler, O.; Rack, A.; Zaslansky, P.; Young, S.L.; Milliron, G.; Chyasnavichyus, M.; Barth, F.G.; Fratzl, P.; Tsukruk, T.; et al. Micro-and nano-structural details of a spider’s filter for substrate vibrations: Relevance for low-frequency signal transmission. J. R. Soc. Interface 2015, 12, 20141111. [Google Scholar] [CrossRef]
  13. Launey, M.E.; Buehler, M.J.; Ritchie, R.O. On the mechanistic origins of toughness in bone. Annu. Rev. Mater. Res. 2010, 40, 25–53. [Google Scholar] [CrossRef] [Green Version]
  14. Amini, S.; Tadayon, M.; Idapalapati, S.; Miserez, A. The role of quasi-plasticity in the extreme contact damage tolerance of the stomatopod dactyl club. Nat. Mater. 2015, 14, 943. [Google Scholar] [CrossRef]
  15. Huang, W.; Shishehbor, M.; Guarín-Zapata, N.; Kirchhofer, N.D.; Li, J.; Cruz, L.; Wang, T.; Bhowmick, S.; Stauffer, D.; Manimunda, P.; et al. A natural impact-resistant bicontinuous composite nanoparticle coating. Nat. Mater. 2020, 19, 1236–1243. [Google Scholar] [CrossRef]
  16. Xu, M.; An, B. Dynamic crack propagation in the turtle carapace. Mech. Mater. 2020, 151, 103614. [Google Scholar] [CrossRef]
  17. Jearanaisilawong, P.; Jongpairojcosit, N.; Glunrawd, C. Dynamic behaviors and protection mechanisms of sulcata tortoise carapace. Comput. Methods Biomech. Biomed. Eng. 2021, 1–13. [Google Scholar] [CrossRef] [PubMed]
  18. Zlotnikov, I.; Zolotoyabko, E.; Fratzl, P. Nano-scale modulus mapping of biological composite materials: Theory and practice. Prog. Mater. Sci. 2017, 87, 292–320. [Google Scholar] [CrossRef]
  19. Fratzl, P.; Weinkamer, R. Nature’s hierarchical materials. Prog. Mater. Sci. 2007, 52, 1263–1334. [Google Scholar] [CrossRef] [Green Version]
  20. Meyers, M.A.; McKittrick, J.; Chen, P.Y. Structural biological materials: Critical mechanics-materials connections. Science 2013, 339, 773–779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Naleway, S.E.; Taylor, J.R.; Porter, M.M.; Meyers, M.A.; McKittrick, J. Structure and mechanical properties of selected protective systems in marine organisms. Mater. Sci. Eng. C 2016, 59, 1143–1167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Politi, Y.; Bar-On, B.; Fabritius, H.O. Mechanics of Arthropod Cuticle-Versatility by Structural and Compositional Variation. In Architectured Materials in Nature and Engineering; Springer: Cham, Switzerland, 2019; pp. 287–327. [Google Scholar]
  23. Shtein, I.; Koyfman, A.; Eshel, A.; Bar-On, B. Autotomy in plants: Organ sacrifice in Oxalis leaves. J. R. Soc. Interface 2019, 16, 20180737. [Google Scholar] [CrossRef] [Green Version]
  24. Bentov, S.; Palmer, B.A.; Bar-On, B.; Shelef, Y.; Aflalo, E.D.; Sagi, A. Reinforcement of bio-apatite by zinc substitution in the incisor tooth of a prawn. Acta Biomater. 2021, 120, 116–123. [Google Scholar] [CrossRef]
  25. Shtein, I.; Bar-On, B.; Popper, Z.A. Plant and algal structure: From cell walls to biomechanical function. Physiol. Plant 2018, 164, 56–66. [Google Scholar] [CrossRef] [PubMed]
  26. Ebenstein, D.M.; Pruitt, L.A. Nanoindentation of biological materials. Nano Today 2006, 1, 26–33. [Google Scholar] [CrossRef]
  27. Moshe-Drezner, H.; Shilo, D.; Dorogoy, A.; Zolotoyabko, E. Nanometer-Scale Mapping of Elastic Modules in Biogenic Composites: The Nacre of Mollusk Shells. Adv. Funct. Mater. 2010, 20, 2723–2728. [Google Scholar] [CrossRef]
  28. Zhou, X.; Miao, H.; Li, F. Nanoscale structural and functional mapping of nacre by scanning probe microscopy techniques. Nanoscale 2013, 5, 11885–11893. [Google Scholar] [CrossRef]
  29. Li, T.; Zeng, K. Nanoscale elasticity mappings of micro-constituents of abalone shell by band excitation-contact resonance force microscopy. Nanoscale 2014, 6, 2177–2185. [Google Scholar] [CrossRef]
  30. Zlotnikov, I.; Shilo, D.; Dauphin, Y.; Blumtritt, H.; Werner, P.; Zolotoyabko, E.; Fratzl, P. In situ elastic modulus measurements of ultrathin protein-rich organic layers in biosilica: Towards deeper understanding of superior resistance to fracture of biocomposites. RSC Adv. 2013, 3, 5798–5802. [Google Scholar] [CrossRef] [Green Version]
  31. Labonte, D.; Lenz, A.K.; Oyen, M.L. On the relationship between indentation hardness and modulus, and the damage resistance of biological materials. Acta Biomater. 2017, 57, 373–383. [Google Scholar] [CrossRef] [PubMed]
  32. Shelef, Y.; Bar-On, B. Interfacial indentations in biological composites. J. Mech. Behav. Biomed. Mater. 2020, 114, 104209. [Google Scholar] [CrossRef] [PubMed]
  33. Ji, B.; Gao, H. Mechanical principles of biological nanocomposites. Annu. Rev. Mater. Res. 2010, 40, 77–100. [Google Scholar] [CrossRef]
  34. Bar-On, B.; Wagner, H.D. Mechanical model for staggered bio-structure. J. Mech. Phys. Solids 2011, 59, 1685–1701. [Google Scholar] [CrossRef]
  35. Bar-On, B.; Wagner, H.D. New insights into the Young’s modulus of staggered biological composites. Mater. Sci. Eng. C 2013, 33, 603–607. [Google Scholar] [CrossRef] [PubMed]
  36. Li, Y.; Ortiz, C.; Boyce, M.C. A generalized mechanical model for suture interfaces of arbitrary geometry. J. Mech. Phys. Solids 2013, 61, 1144–1167. [Google Scholar] [CrossRef]
  37. Zhang, P.; Heyne, M.A.; To, A.C. Biomimetic staggered composites with highly enhanced energy dissipation: Modeling, 3D printing, and testing. J. Mech. Phys. Solids 2015, 83, 285–300. [Google Scholar] [CrossRef] [Green Version]
  38. Qwamizadeh, M.; Zhou, K.; Zhang, Y.W. Damping behavior investigation and optimization of the structural layout of load-bearing biological materials. Int. J. Mech. Sci. 2017, 120, 263–275. [Google Scholar] [CrossRef]
  39. Wu, J.; Yuan, H.; Li, L.; Fan, K.; Qian, S.; Li, B. Viscoelastic shear lag model to predict the micromechanical behavior of tendon under dynamic tensile loading. J. Theor. Biol. 2018, 437, 202–213. [Google Scholar] [CrossRef]
  40. Liu, J.; Zhu, W.; Yu, Z.; Wei, X. Dynamic shear-lag model for understanding the role of matrix in energy dissipation in fiber-reinforced composites. Acta Biomater. 2018, 74, 270–279. [Google Scholar] [CrossRef]
  41. Yu, Z.; Liu, J.; Wei, X. Achieving outstanding damping performance through bio-inspired sutural tessellations. J. Mech. Phys. Solids 2020, 142, 104010. [Google Scholar] [CrossRef]
  42. Bar-On, B.; Wagner, H.D. Stiffness of the extrafibrillar phase in staggered biological arrays. Phys. Rev. Lett. 2012, 109, 078102. [Google Scholar] [CrossRef] [PubMed]
  43. Prapavesis, A.; Tojaga, V.; Östlund, S.; van Vuure, A.W. Back calculated compressive properties of flax fibers utilizing the Impregnated Fiber Bundle Test (IFBT). Compos. Part A Appl. Sci. Manuf. 2020, 135, 105930. [Google Scholar] [CrossRef]
  44. Khodayari, A.; Hirn, U.; Van Vuure, A.W.; Seveno, D. Inverse rule of mixtures at the nanoscale: Prediction of elastic properties of cellulose nanofibrils. Compos. Part A Appl. Sci. Manuf. 2020, 138, 106046. [Google Scholar] [CrossRef]
  45. An, B.; Sun, W. A theory of biological composites undergoing plastic deformations. J. Mech. Behav. Biomed. Mater. 2019, 93, 204–212. [Google Scholar] [CrossRef] [PubMed]
  46. Xu, M.; An, B. An analysis of fracture in staggered mineralized collagen fibril arrays. Int. J. Solids Struct. 2020, 193–194, 535–549. [Google Scholar] [CrossRef]
  47. Yu, L.; Ma, Y.; Zhou, C.; Xu, H. Damping efficiency of the coating structure. Int. J. Solids Struct. 2005, 42, 3045–3058. [Google Scholar] [CrossRef]
  48. Gusev, A.A.; Lurie, S.A. Loss amplification effect in multiphase materials with viscoelastic interfaces. Macromolecules 2009, 42, 5372–5377. [Google Scholar] [CrossRef]
  49. Lurie, S.; Minhat, M.; Tuchkova, N.; Soliaev, J. On remarkable loss amplification mechanism in fiber reinforced laminated composite materials. Appl. Compos. Mater. 2014, 21, 179–196. [Google Scholar] [CrossRef]
  50. Bar-On, B.; Bayerlein, B.; Blumtritt, H.; Zlotnikov, I. Dynamic response of a single interface in a biocomposite structure. Phys. Rev. Lett. 2015, 115, 238001. [Google Scholar] [CrossRef] [PubMed]
  51. Bayerlein, B.; Bertinetti, L.; Bar-On, B.; Blumtritt, H.; Fratzl, P.; Zlotnikov, I. Inherent Role of Water in Damage Tolerance of the Prismatic Mineral–Organic Biocomposite in the Shell of Pinna Nobilis. Adv. Funct. Mater. 2016, 26, 3663–3669. [Google Scholar] [CrossRef]
  52. Funari, M.F.; Greco, F.; Lonetti, P. Dynamic debonding in layered structures: A coupled ALE-cohesive approach. Frat. Integrita Strutt. 2017, 11, 524–535. [Google Scholar] [CrossRef] [Green Version]
  53. Cohen, Y.; Ronen, D.; Ya’akobovitz, A. Mechanical behavior of vertically aligned carbon nanotubes under electrostatic tension. Sens. Actuators A Phys. 2019, 292, 105–111. [Google Scholar] [CrossRef]
  54. Cohen, Y.; Ya’akobovitz, A. The influence of thermal loads on the physical properties of carbon nanotubes forests. Microelectron. Eng. 2021, 247, 111575. [Google Scholar] [CrossRef]
  55. Gibson, R.F. Principles of Composite Material Mechanics; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  56. Seto, J.; Gupta, H.S.; Zaslansky, P.; Wagner, H.D.; Fratzl, P. Tough lessons from bone: Extreme mechanical anisotropy at the mesoscale. Adv. Funct. Mater. 2008, 18, 1905–1911. [Google Scholar] [CrossRef]
  57. Krauss, S.; Monsonego-Ornan, E.; Zelzer, E.; Fratzl, P.; Shahar, R. Mechanical function of a complex three-dimensional suture joining the bony elements in the shell of the red-eared slider turtle. Adv. Mater. 2009, 21, 407–412. [Google Scholar] [CrossRef]
  58. Yang, W.; Chen, I.H.; Gludovatz, B.; Zimmermann, E.A.; Ritchie, R.O.; Meyers, M.A. Natural flexible dermal armor. Adv. Mater. 2013, 25, 31–48. [Google Scholar] [CrossRef]
  59. Lee, N.; Horstemeyer, M.F.; Rhee, H.; Nabors, B.; Liao, J.; Williams, L.N. Hierarchical multiscale structure–property relationships of the red-bellied woodpecker (Melanerpes carolinus) beak. J. R. Soc. Interface 2014, 11, 20140274. [Google Scholar] [CrossRef] [Green Version]
  60. Gao, C.; Hasseldine, B.P.; Li, L.; Weaver, J.C.; Li, Y. Amplifying strength, toughness, and auxeticity via wavy sutural tessellation in plant seedcoats. Adv. Mater. 2018, 30, 1800579. [Google Scholar] [CrossRef] [PubMed]
  61. Achrai, B.; Bar-On, B.; Wagner, H.D. Biological armors under impact—Effect of keratin coating, and synthetic bio-inspired analogues. Bioinspirat. Biomim. 2015, 10, 016009. [Google Scholar] [CrossRef] [PubMed]
  62. Liu, L.; Jiang, Y.; Boyce, M.; Ortiz, C.; Baur, J.; Song, J.; Li, Y. The effects of morphological irregularity on the mechanical behavior of interdigitated biological sutures under tension. J. Biomech. 2017, 58, 71–78. [Google Scholar] [CrossRef]
  63. Yin, Z.; Hannard, F.; Barthelat, F. Impact-resistant nacre-like transparent materials. Science 2019, 364, 1260–1263. [Google Scholar] [CrossRef]
  64. Serra, M.; Arenal, R.; Tenne, R. An overview of the recent advances in inorganic nanotubes. Nanoscale 2019, 11, 8073–8090. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Reddy, S.K.; Ya’akobovitz, A. Electromechanical behavior of graphene foams. Appl. Phys. Lett. 2019, 115, 211902. [Google Scholar] [CrossRef]
  66. Ben-Shimon, Y.; Ya’akobovitz, A. Flexible and bio-compatible temperature sensors based on carbon nanotube composites. Measurement 2021, 172, 108889. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic examples of the interfacial regions in biocomposites: graded matrix properties, shaped matrix-reinforcement connection, and staggered matrix-reinforcement connection. (b) Mechanical modeling of a biocomposite segment (length L c and dynamic modulus E c * ), which includes an effective viscoelastic interface (length L i and dynamic modulus E i * ) between adjacent elastic reinforcements (elastic modulus E f ).
Figure 1. (a) Schematic examples of the interfacial regions in biocomposites: graded matrix properties, shaped matrix-reinforcement connection, and staggered matrix-reinforcement connection. (b) Mechanical modeling of a biocomposite segment (length L c and dynamic modulus E c * ), which includes an effective viscoelastic interface (length L i and dynamic modulus E i * ) between adjacent elastic reinforcements (elastic modulus E f ).
Materials 14 03428 g001
Figure 2. The correspondence between E i and E c for various biocomposite configurations. The dashed lines indicate the theoretical results via Equation (4), and the symbols indicate the corresponding finite-element results (Table S1 in Supporting Information). The colors indicate data sets with the same k E parameter: black, red, and green correspond to k E = 0.1 ,     0.3 , and 0.8 , respectively. The circle and square symbols represent different viscoelastic biocomposite configurations ( tan δ c 0 ), while the plus and asterisk symbols represent the corresponding biocomposite configurations, but with completely elastic properties ( tan δ c = 0 ).
Figure 2. The correspondence between E i and E c for various biocomposite configurations. The dashed lines indicate the theoretical results via Equation (4), and the symbols indicate the corresponding finite-element results (Table S1 in Supporting Information). The colors indicate data sets with the same k E parameter: black, red, and green correspond to k E = 0.1 ,     0.3 , and 0.8 , respectively. The circle and square symbols represent different viscoelastic biocomposite configurations ( tan δ c 0 ), while the plus and asterisk symbols represent the corresponding biocomposite configurations, but with completely elastic properties ( tan δ c = 0 ).
Materials 14 03428 g002
Figure 3. The correspondence between tan δ i and tan δ c for various biocomposite configurations. The dashed lines indicate the theoretical results via Equation (4), and the symbols indicate the corresponding finite-element results (Table S2 in Supporting Information). The colors indicate data sets with the same k δ parameter: black, red, and green correspond to k δ = 1.05 ,   2.5 , and 5 , respectively. The range of tan δ i spans between completely elastic ( tan δ i = 0 ) and predominantly viscoelastic ( tan δ i = 1 / 2 ) interfaces.
Figure 3. The correspondence between tan δ i and tan δ c for various biocomposite configurations. The dashed lines indicate the theoretical results via Equation (4), and the symbols indicate the corresponding finite-element results (Table S2 in Supporting Information). The colors indicate data sets with the same k δ parameter: black, red, and green correspond to k δ = 1.05 ,   2.5 , and 5 , respectively. The range of tan δ i spans between completely elastic ( tan δ i = 0 ) and predominantly viscoelastic ( tan δ i = 1 / 2 ) interfaces.
Materials 14 03428 g003
Figure 4. Schematic descriptions of zigzag-shaped sutural interfaces in biocomposites and their geometrical parameters. (a) A biocomposite segment with an underlying sutural interface. (b) An isolated interfacial region (one period) from the biocomposite. (c) Schematic examples of sutural interfaces with the same matrix length ( L m ), but different zigzag angles ( θ = 10 ° ,   45 ° , and ~ 70 ° ).
Figure 4. Schematic descriptions of zigzag-shaped sutural interfaces in biocomposites and their geometrical parameters. (a) A biocomposite segment with an underlying sutural interface. (b) An isolated interfacial region (one period) from the biocomposite. (c) Schematic examples of sutural interfaces with the same matrix length ( L m ), but different zigzag angles ( θ = 10 ° ,   45 ° , and ~ 70 ° ).
Materials 14 03428 g004
Figure 5. Simulation results for the interfacial dynamic modulus of biocomposites with zigzag-shaped sutural interfaces. (a) the modulus magnitude and (b) the loss coefficient of the interfacial region for different zigzag angles ( θ ), achieved by back-calculations from far-field DMA testing on the biocomposite ( L i / L c = 0.1 ), and by direct DMA testing on the interfacial region ( L i / L c = 1 ).
Figure 5. Simulation results for the interfacial dynamic modulus of biocomposites with zigzag-shaped sutural interfaces. (a) the modulus magnitude and (b) the loss coefficient of the interfacial region for different zigzag angles ( θ ), achieved by back-calculations from far-field DMA testing on the biocomposite ( L i / L c = 0.1 ), and by direct DMA testing on the interfacial region ( L i / L c = 1 ).
Materials 14 03428 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shelef, Y.; Uzan, A.Y.; Braunshtein, O.; Bar-On, B. Assessing the Interfacial Dynamic Modulus of Biological Composites. Materials 2021, 14, 3428. https://doi.org/10.3390/ma14123428

AMA Style

Shelef Y, Uzan AY, Braunshtein O, Bar-On B. Assessing the Interfacial Dynamic Modulus of Biological Composites. Materials. 2021; 14(12):3428. https://doi.org/10.3390/ma14123428

Chicago/Turabian Style

Shelef, Yaniv, Avihai Yosef Uzan, Ofer Braunshtein, and Benny Bar-On. 2021. "Assessing the Interfacial Dynamic Modulus of Biological Composites" Materials 14, no. 12: 3428. https://doi.org/10.3390/ma14123428

APA Style

Shelef, Y., Uzan, A. Y., Braunshtein, O., & Bar-On, B. (2021). Assessing the Interfacial Dynamic Modulus of Biological Composites. Materials, 14(12), 3428. https://doi.org/10.3390/ma14123428

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop