Introduction

Many animals have evolved hard exoskeletons to resist predation or competitive attacks (e.g., crustaceans1, insects2, mollusks3, turtles4,5, seahorses6,7, and bony fish8,9,10). These ‘natural armors’ combine micro- and nano-scale materials design (e.g., materials selection, crystallography, composite architecture, porosity, surface chemistry) with macroscale geometrical design rules11,12,13,14,15 to provide additional functionalities16 such as enhanced mechanical properties17,18, transparency19, or flexibility20,21. Fish possess highly flexible armored exoskeletons which have attracted particular interest in the field of bio-inspired design22,23,24. The armored fish Polypterus senegalus (bichir) possesses an exoskeleton of imbricated, articulating, mineralized scales that provides penetration resistance from predatory and territorial attacks while granting the fish serpentine mobility25,26,27,28. While the complex geometry of the scales, their shape variation along the body length, and their protective capabilities have been studied20,27,29,30,31,32, there has yet been no comprehensive analysis of the geometric design rules, such as scale geometry, local articulation, and global assembly.

Here, we integrate advances in microtomographic imaging, parametric modeling, and multi-material 3D printing to achieve three goals: (1) translate the hierarchical geometric rules of assembly in the P. senegalus exoskeleton to a synthetic, flexible composite prototype, (2) experimentally and computationally assess the prototype’s mechanical behavior, and (3) elucidate the structure–function relationships that determine how shape of individual components can be used as a materials design parameter to tune the anisotropic behavior of composite materials. Parametric 3D modeling was used to design bio-inspired prototypes by generating abstracted models of scales from x-ray microtomography data and integrating the multi-material components into the flexible composite assembly. Multi-material 3D printing was used to fabricate the bio-inspired flexible armor prototypes. Anisotropic flexibility of the prototypes was examined under passive loading (self-weight). A mechanical tester was used to quantify the bending stiffness of the prototypes in active loading (bending), and finite element simulations were developed to correlate internal stress concentrations with mechanical response.

The results show how the complex scale shape contributes to local interscale mobility mechanisms that determine the bending response of the global prototype and generate anisotropic mechanical behavior. The most flexible orientations, including one in which the scales are ‘mechanically invisible’ without adding stiffness to the armor assembly, correspond to physiologically relevant bending modes in the fish. With one prototype design scheme, a wide array of mechanical behavior was generated with stiffness ranging over several orders of magnitude, thus showing how morphometry can tune the flexibility of protective, composite architectures without varying the constituent materials or their volume fractions. We anticipate that the introduction of geometric variation into synthetic prototypes will generate flexible, protective systems that are well adjusted to complex shapes, kinematics, and functional differentiations33.

Results

Translation of the geometric rules of assembly to a bio-inspired prototype

The P. senegalus exoskeleton covers the entire surface of the fish from head to tail. The scales are arranged in helical columns that wind around the body, with an angle (β) of 60° between the helical, paraserial axis and the horizontal plane in a straight body posture, shown in the schematic in Fig. 1a. Segmentation of the scales, their complex geometry, and their joint articulation mechanisms allow the fish to achieve large, bi-directional body curvatures in axial bending (Fig. 1b).

Fig. 1: The P. senegalus exoskeleton.
figure 1

a Schematic of the left side of a fish showing the helical arrangement of scales with angle β = 60° between the helical (paraserial) axis and the horizontal plane (sketch by the authors). b Anesthetized P. senegalus specimen (body length = 219 mm) exhibiting bi-directional body curvature in axial bending. Scale bar = 10 mm.

P. senegalus scales were scanned using x-ray microtomography (µCT) to generate reconstructed, digital 3D models. The shape and size of scales vary gradually across the fish exoskeleton, so in order to produce a tileable scale geometry, µCT data from an adjacent pair of scales were combined into a unitized scale shown in Fig. 2a. The individual scale geometry is complex with distinct features including the peg (P), socket (S), anterior process (AP), anterior shelf (AS), concave anterior margin (AM), and a thickened axial ridge (AR). Scales use two primary joints to interact with their neighbors and form an assembly in the exoskeleton (Fig. 2b–c): the articulated peg-and-socket joint in the paraserial direction within a column of scales, and the sliding overlap joint in the interserial direction between columns of scales situated at the angle β away from the paraserial axis.

Fig. 2: Translation of the biological design rules to a synthetic flexible composite material.
figure 2

a 3D µCT reconstruction of an averaged scale highlighting geometrical features: peg (P), socket (S), anterior process (AP), anterior shelf (AS), concave anterior margin (AM), and axial ridge (AR). Scale bar = 1 mm. b Homogeneous assembly of the biological scale geometry, interior view (scale bar = 1 mm), and c exterior view, showing the paraserial and interserial axes oriented at angle β = 60° from each other. Scale bar = 1 mm. d A photograph of the biological exoskeleton in a deceased P. senegalus specimen. Scale bar = 5 mm. e Abstracted geometry (3D model) of a single scale unit magnified to 20 mm length. Scale bar = 1 mm. f Associative 3D model of the scale assembly that incorporates the essential parts of the exoskeletal assembly, including the scales, substrate, paraserial connections, and scale–substrate attachment, interior view (scale bar = 20 mm), g exterior view (scale bar = 20 mm), and h cross-sectional view. Scale bar = 20 mm. i A photograph of a multi-material 3D-printed prototype of the scale assembly.

Collagenous Sharpey’s fibers run between the peg and socket of adjacent scales to reinforce the joint, and a collagenous, multi-layered fibrous stratum compactum attaches the scales to the underlying soft tissue dermis at the axial ridge27. A portion of the protective exoskeleton shown is in Fig. 2d.

Prototypes were designed using computer-aided design (CAD) software to render the individual scale geometry into an abstracted 3D model (Fig. 2e) of 20 mm length with an overall rhomboid shape allowing for a tapered overlap area between scales, a tetrahedral peg, and a corresponding inverted concave socket. Parametric modeling was used to replicate the abstracted scale into a 2D array which preserved the joint articulation mechanisms with β = 60° between the paraserial and interserial axes (Fig. 2f–h). An underlying substrate was added, and the two compliant connective tissue elements were then integrated into the model as paraserial scale connections and attachments between the scale and the substrate (Fig. 2f).

The translated 3D designs were fabricated by multi-material inkjet 3D printing, using VeroWhite photopolymer (hard plastic with elastic modulus (E) = 2.0 GPa34,35) for the scales and TangoPlus photopolymer (rubber-like elastomer with E = 0.63 MPa34,35) for the soft components (Sharpey’s fibers, stratum compactum, and underlying tissue) to yield synthetic, flexible composite prototypes (Fig. 2i). We choose this 3D printing technique due to its significant advantages over other fabrication techniques21,22, including high resolution and dimensional accuracy, ease of producing arbitrary multi-material 3D structures, and the large stiffness ratio between different photopolymers (order of ~104). This stiffness ratio reasonably agrees with the expected physiological value (order of ~103) based on representative stiffnesses of mineralized (bone, ~10 GPa) and soft (collagen, ~10 MPa) biological tissues. The distance between the scales in the 3D printed prototypes is 1 mm, approximately thirty times larger compared to the resolution of the printer (30 microns), which ensures the soft connections between the scales are resolved.

Anisotropic flexibility is governed by paraserial connectivity

The innate flexibility of the prototype under passive loading (self-weight) was examined by draping a prototype comprising a square array of seventy two 20 mm scales over a half-cylinder mold with radius of curvature (Rm) without the application of an external load (Fig. 3). To ensure that our biomimetic model can address the full range of flexibility of the fish, we chose the ratio of scale size to the mold radius of curvature to equal the ratio of physiological scale size to the maximum physiological radius of curvature calculated from the reported curvature value27.

Fig. 3: Anisotropy of prototype flexibility measured by radius of curvature.
figure 3

a Prototype with an array of 72 scales is draped over a half-cylinder mold at angle (α) = 0° between the prototype’s paraserial axis and the mold’s line of zero curvature. Lines show the curvature of the mold (white) and prototype (red). Inverse of curvature represents the radius of curvature of the prototype (Rp) and the radius of curvature of the mold (Rm). The prototype exhibits different curvatures when rotated to b α = 30°, c α = 60°, d α = 90°. e Relative radius of curvature (Rp/Rm) vs. α for the prototype, a variation without the paraserial connections, and the substrate only. Error bars represent standard deviation with N = 3 samples per prototype design. f Top view of the prototype in the α = 0° orientation. g Top view of the prototype in the α = 90° orientation.

The flexibility of the prototype, as shown by the calculated radius of curvature (Rp), varied with the orientation of the scales over the mold (Fig. 3a–d). The radius of curvature of the prototype relative to the radius of curvature of the mold (Rp/Rm) is plotted against orientation angle (α) in Fig. 3e (N = 3 per sample). The prototype attains maximum flexibility at α = 0° where the paraserial peg-and-socket axis is aligned with the mold’s axis of zero curvature (Fig. 3f); the overlap joint allows the columns of scales to slide away from each other without resistance, and the prototype completely conforms to the mold with Rp = Rm. The radius of curvature of the prototype increased with α until reaching a maximum rigidity at α = 90°, where the peg-and-socket joint was aligned along the axis of mold curvature (Fig. 3d, g); the articulated peg-and-socket joint restricts the motion of scales within the column of scales, and thus resists flexure of the prototype so that it lies almost flat above the mold with Rp = 4.2Rm. As the prototype was rotated past α = 90°, the peg-and-socket joint returns to alignment with the mold’s line of zero curvature, and the global flexibility returns. A second prototype was designed without the paraserial connective elements between the peg and socket of neighboring scales. This prototype assembly exhibited uniformly maximal flexibility in all orientations, mechanically equivalent to the flexible substrate without scales (Fig. 3e). The results suggest that collagenous Sharpey’s fibers are a critical structural component in restricting the ranges of motion between scales at the peg-and-socket joint, thereby increasing the resistance to flexure along the paraserial axis and introducing mechanical anisotropy to the system.

Interscale mobility mechanisms and the bending stiffness

We tested bioinspired flexible composite prototypes in active loading (bending) to examine how scale shape contributes to local interscale mobility mechanisms and generates anisotropic mechanical behavior (Fig. 4a–c).

Fig. 4: Anisotropic bending stiffness of the flexible composite prototypes.
figure 4

a Rendered front view of the prototype designed for bending with rigid rods at the top and bottom and with the paraserial axis oriented at an angle (φ) = 60° away from the loading direction. b Schematic drawing of the sample holder. c Side view of the φ = 60° sample loaded in bending at various vertical displacements (d): d = 20 mm (left) and d = 80 mm (right). d Force–displacement curves for prototypes with φ = 0°, 30°, 60°, 90°, 120°, and 150°. e Normalized stiffness (K) of each phase of each prototype’s loading curve (experiments and finite element simulations). Error bars represent standard deviation with N = 3 samples per prototype orientation.

Reaction force (F) vs. vertical displacement (d) for the prototypes (N = 3 per orientation) is plotted in Fig. 4d. Each orientation (φ) has a characteristic loading response which can be divided into phases. Stiffness (K) for each phase, plotted in Fig. 4e, is calculated as the slope of the loading curve and normalized by the stiffness of a control sample with no scales, comprising TangoPlus elastomer only. The scale-less prototype was chosen as the control as it represents an integument without protective scales, which would be expected to have maximum flexibility (lowest stiffness). Stiffness and interscale mechanisms observed in each loading phase of each orientation are tabulated in Table 1. Finite element (FE) models of the prototypes were also created to computationally simulate bending (Fig. 5). The interscale mobility mechanisms that contribute to the mechanical response of the global assembly are depicted in Fig. 6.

Table 1 Displacement range, stiffness (K), and interscale mobility mechanisms for each loading phase of each prototype orientation (φ).
Fig. 5: Orientation-dependent bending behavior of the scaled prototypes.
figure 5

a Experimental bending of the prototype with φ = 0°. Stress plots (Mises, linear elastic) from the FE simulation of the φ = 0° model. (i) Front view at d = 3 mm, (ii) d = 9 mm, and (iii) d = 22 mm. b Experimental bending of the prototype with φ = 30°. Stress plots, front view at (i) d = 3 mm, (ii) d = 14 mm, and (iii) d = 33 mm. c Experimental bending of the prototype with φ = 60°. Stress plots, front view at (i) d = 5 mm, (ii) d = 8 mm, and (iii) d = 38 mm. d Experimental bending of the prototype with φ = 90°. Stress plots, front view at (i) d = 12 mm, (ii) d = 43 mm, and (iii) d = 68 mm. e Experimental bending of the prototype with φ = 120°. Stress plots, front view at (i) d = 3 mm, (ii) d = 12 mm, and (iii) d = 52 mm. f Experimental bending of the prototype with φ = 150°. Stress plots, front view at (i) d = 20 mm, (ii) d = 38 mm, and (iii) d = 55 mm.

Fig. 6: Orientation-dependent bending behavior of the scaled prototypes.
figure 6

a Schematics of the six interscale mobility mechanisms observed during bending. b Experimental bending of the prototype with (i) φ = 0°. Paraserial bending at vertical displacement (d) = 3 mm and (ii) back view of sample showing anterior process pushing into the substrate at d = 15 mm. (iii) φ = 30°. Paraserial rotation at d = 15 mm. (iv) φ = 60°. Paraserial and interserial rotation at d = 50 mm. (v) φ = 90°. Interserial sliding at d = 15 mm, side view and (vi) further interserial sliding after interserial contacts at d = 60 mm (side view). (vii) φ = 120°. Paraserial and interserial rotations and (viii) interscale contacts at d = 20 mm. (ix) φ = 150. Paraserial and interserial rotations plus interserial splay at d = 14 mm.

In the φ = 0° orientation, the paraserial axis of the scales is aligned with the loading direction. In the first phase, the scales undergo paraserial bending as the compliant TangoPlus connection within the peg-and-socket joint resists deformation (K = 47.8, Table 1). In the second phase, the scales interlock paraserially (Fig. 5a.i), generating stress concentrations around the peg and socket (Fig. 5a.ii). The paraserial interlock causes paraserial rotation as the anterior margin rotates toward the substrate, coupled with interserial rotation when the anterior process pushes into the substrate (into plane) and the back of the anterior margin and socket lifts up (out of plane) (6a, 6b.i-ii). At large deformations, the high internal stresses are distributed throughout the body of the scales (Fig. 5a.iii).

In the φ = 30° orientation, the scales initially undergo paraserial bending (K = 126) until paraserial interlock (Fig. 6b.iii). The interlock causes small paraserial rotation alongside paraserial bending as the anterior margin rotates into plane (Fig. 6b.iii). The orientation of the scales gives tolerance to the anterior process so that it is not pushed into substrate; there is no coupling with interserial rotation, and the stiffness drops (K = 21.7). As the model bends, stresses form at the sites of paraserial, interscale contact around the peg and socket of adjacent scales; the stresses grow and become distributed over the scale surface with higher degrees of bending (Fig. 5b.i–iii). At large deformations, paraserial failures occur in parallel throughout the sample (Fig. 5b.iii), after which the prototype offers no resistance to global bending (K < 0).

In the φ = 60° orientation, the scales first undergo interserial sliding as the anterior shelf slides under the scale in the adjacent column without generating any stresses on the scales (Fig. 5c.i). The oblique angle of orientation then causes the scales to touch interserially and resist deformation (K = 31.3) while generating stresses at the site of contact (Fig. 5c.ii). As the sample continues to bend, the interserial interlock causes small paraserial rotation (anterior margin moves out of plane) coupled with large interserial rotation (anterior process moves out of plane; back of anterior margin and socket move into plane) to accommodate interserial sliding (Fig. 6b.iv). Stresses grow at the sites of interscale contact and spread throughout the body of the scales (Fig. 5c.iii). The compliant paraserial connections resist the rotations but do not break, so the prototype is able to bear load as it bends (K = 10.5).

In the φ = 90° orientation, the paraserial axis is perpendicular to the loading direction. In bending, the columns of scales move relative to each other via interserial sliding (Fig. 5d.i, 6b.v), and the only resistance to bending comes from the tensile stresses in the substrate beneath the scales. Here, the rigid scales are ‘mechanically invisible’ as they do not contribute any resistance to bending, and the sample’s stiffness matches that of the substrate (K = 0.99). At large degrees of deformation, the scales begin to touch interserially (Fig. 6b.v, vi). The scale contacts cause the columns of scales to rotate about the axial ridge, further straining the compliant material at the site of attachment to the substrate and introducing stresses into the center of the scales (Fig. 5d.ii–iii). As a result, the stiffness of the sample increases (K = 4.23).

When the φ = 120° prototype bends, the scales first move via interserial sliding. However, the anterior process is pushed into the substrate immediately, and these contacts induce paraserial rotation (anterior margin moves into plane) coupled with interserial rotation (anterior process and peg move into plane; back of anterior margin and socket move out of plane) to generate additional paraserial and interserial scale contacts (Fig. 6b.vii–viii) that in turn generate high stiffness (K = 124). Stress concentrations form on the scales in regions surrounding the peg and socket where paraserial contacts occur (Fig. 5e.i), at the sites of interserial contact between the columns of scales (Fig. 5e.ii), and throughout the substrate where the anterior process pushes in. These scale interlocks contribute to the high stiffness with further bending (K = 46.8). At very large deformations, the scales splay interserially (Fig. 5e.iii) to further accommodate sliding with reduced stiffness (K = 24.1).

The φ = 150° prototype is the stiffest of all orientations. At the onset of bending, small degrees of paraserial bending causes the anterior process to push into the substrate immediately while the back end of the anterior margin impinges on the scale in the adjacent column. These scale contacts do not allow for further paraserial bending, provide no tolerance for interserial sliding, and cause high stiffness (K = 450). With further bending, the scale contacts cause paraserial rotation (anterior margin moves into plane) coupled with interserial rotation (Fig. 6b.ix); however, the interlocked scales continue to resist bending with high stiffness (K = 80.9) and high stresses are sustained at the sites of contact between scales and spread throughout selective columns of scales (Fig. 5f.i–ii). At very large deformations, the scale interlocks resist bending such that the rods which are inserted into the sample holders bow out to accommodate the deformation with high stresses throughout the model (Fig. 5f.iii). Eventually, the compliant paraserial interconnections start to tear, and small degrees of paraserial bending are observed, while the scale interlocks allow the prototype to sustain high loads (K = 21.4).

For all orientations, the measured stiffnesses K are less than what would be expected from a rule-of-mixtures combination (i.e., homogeneous blend) of the rigid and soft photopolymers. The volume percent of stiff (VeroWhite, 2.0 GPa) and compliant (TangoPlus, 0.63 MPa) materials in the printed prototypes were 46% and 54%, respectively, after removal of support material, so the rule-of-mixtures expected stiffness of a homogenous blend is 0.92 GPa (predicted K = 1460), similar to a semirigid plastic such as high density polyethylene, which would not be expected to provide measurable drapability in the passive loading test. This demonstrates the importance of the segmentation of the integument into discrete, morphometrically complex scales in order to achieve anisotropic flexibility while providing protection.

Discussion

A hierarchy of shape- and materials-based design principles were translated from the biological exoskeleton of P. senegalus and integrated into the bioinspired flexible composite prototypes, including the complex shape of rigid scales, interscale joint articulation structure, assembly of scales into an armored surface, and soft connective components (substrate, scale-to-substrate attachment, and paraserial connections). The prototypes were able to replicate the biomechanical behavior of the biological exoskeleton, where the complex scale shape and joint articulations contribute to local, interscale mobility mechanisms that in turn determine the bending response of the global sample, and also generate global anisotropic (orientation-dependent) mechanical behavior in the scale assembly.

The fish engages in both convex and concave bending during its normal undulatory motion27. We choose to look at concave bending, since scale-to-scale contacts generate greater resistance to bending. The two lowest stiffness orientations, φ = 90° and φ = 60°, correspond to the two commonly observed bending modes in the fish: axial bending and torsion, respectively27,28. In these orientations, the interserial sliding mechanism allows the assembly to bend under a small applied load without generating interscale contacts or introducing stress concentrations within the scales in the initial phase of deformation. In the 90° orientation, all strains are sustained in the substrate up to 36% deformation (40 vertical mm / 112 mm total sample height), and sample stiffness matches that of the control sample without any scales. Since the rigid scales are “mechanically invisible” and do not contribute any resistance to bending, we show that it is possible to use shape as a materials design parameter to create composite materials that provide added protection from the stiffer material (e.g., scales) while maintaining the low bending stiffness of the compliant material (e.g., substrate). In the 60° orientation, all strains are sustained in the substrate up to 4.5% deformation, after which the orientation of scales relative to the loading direction allows the interscale contacts to utilize low stiffness interscale mobility mechanisms (i.e., coupled paraserial and interserial rotation), to continue to enable interserial sliding of scales.

The high stiffness orientations correspond to bending modes in which the fish does not engage; for instance, 150° represents dorso-ventral bending about a horizontal plane through the middle of the fish. Thus, the complex geometry and orientation of scales enables flexibility of the integument in directions that it uses for axial bending and torsion, and restricts flexibility in the directions it does not need. Features such as the oblong anterior process and concave anterior margin that make contact with the substrate and interlock with neighboring scales to provide full-body coverage of scales over the dermis as the scales move apart from each other, e.g., in ventral scales that splay apart from each other since they are oriented away from the plane of axial bending or torsion29, and in scales that slide apart during convex bending.

Understanding complex materials-morphometric design rules in natural exoskeletons and translating them to synthetic designs holds tremendous application for bio-inspired flexible armor16,36. By creating flexible composite prototypes inspired by the scale armor of P. senegalus, we have generated one design scheme that exhibits a wide array of mechanical behavior with bending stiffness ranging over several orders of magnitude (K = 1–450), thus showing how morphometry can tune the flexibility of protective, composite architectures without varying the constituent materials or their volume fraction. Our concurrent and future work seeks to integrate morphometric heterogeneity29,33, ability to conform to arbitrary curved surfaces33, and intrascale material heterogeneity37 into the prototype design for a truly hierarchical design that replicates all aspects of the biological armor. Biomimetic armors utilizing a segmented design hold enormous potential for a wide variety of applications by allowing damage localization, flexibility, reduced cost of fabrication, and selective replacement of damaged units.

Methods

Fish scale extraction and 3D reconstruction

Scales were surgically excised from a live P. senegalus specimen (219 mm body length). The specimen was anesthetized with 0.03 wt.% solution of tricaine methanesulfonate (MS-222, Sigma Aldrich) in water neutralized with potassium hydroxide. Four scales were removed from the 49th column on the left flank of the specimen with sterile, surgical-grade scalpels. The fish was transferred to anesthetic-free water to recover, and then returned to a quarantine aquarium and treated with tetracycline antibiotics (250 mg per 10 gallons of water per day) for two weeks until the scales began to regenerate. All work with the live specimen was performed in accordance with MIT’s Committee on Animal Care and IACUC regulations. The excised scales were scanned using x-ray microtomography (µCT) and reconstructed as tessellated surface files (STL) following our previously published procedure29. Two neighboring scales were segmented into halves and reassembled into one “unitized” scale geometry (Fig. 1a) that can be assembled in tiled arrangements. The unitized scale geometry was tiled over the surface to study the geometrical principles of scale articulation (Fig. 1b). We then chose the main geometrical features to mimic in the 3D modeling of the biomimetic scale design: peg (P), socket (S), anterior process (AP), anterior shelf (AS), concave anterior margin (AM), and a thickened axial ridge (AR).

Computational 3D modeling of prototypes

Geometric morphometric analysis was used to define the geometry of the scale and its features from the 3D scale object following our previously published procedure29. Parametric CAD software (SOLIDWORKS®, Dassault Systèmes SolidWorks Corp., France) was used to design an abstracted 3D model of the scale geometry of 20 mm length with an overall rhomboid shape allowing for a tapered overlap area between scales, a tetrahedral peg, and a corresponding inverted concave socket. Associative modeling was used to replicate the individual scale geometry into square arrays with 1 mm spacing in the paraserial and interserial directions. The soft tissue substrate was modeled as a separate layer. Connective elements were modeled between the peg and socket of adjacent scales and between the scales and the substrate as additional layers. The prototypes used for measuring radius of curvature were designed as square arrays of 72 scales of 20 mm atop the substrate. The prototypes for the bending experiment were designed as 112 × 124 mm arrays of 20 mm scales with rigid rods of 10 mm diameter at the top and bottom of the sample. The prototype designs were exported as separate STL files for the rigid components (scale, rods) and soft components (substrate, connective elements).

Multi-material 3D printing

The prototypes were fabricated as a flexible array of scales via multi-material 3D printing (OBJET Connex500™, Stratasys, USA). The STL files for the prototype were imported into OBJET Studio software and assigned to commercially available UV-cured photopolymer materials: the rigid components were printed with VeroWhite (hard plastic with elastic modulus (E) = 2.0 GPa34,35), and the soft components were printed with TangoPlus (rubber-like elastomer with E = 0.63 MPa34,35). The print jobs were submitted using the digital printing mode at 30 µm resolution. Print support material was removed with a water jet and manual brushing.

Curvature experiment

The radius of curvature of the prototype under self-weight was examined by draping the prototype over a curved, half-cylinder mold (radius Rm = 120 mm) without the application of an external load, with a camera situated along the mold’s axis of zero curvature. The prototype was rotated over the mold by an angle α = 0–180°, where α = 0° corresponds to the paraserial peg-and-socket axis in line with the mold’s axis of zero curvature, and α = 90° corresponds to the paraserial axis in line with the mold’s axis of curvature. In each orientation, normal projection rods were inserted into three scales on a single line parallel to the mold’s axis of curvature, and the radius of curvature of the prototype (Rp) was measured by drawing a circle amongst points of connection between the normal rod and the scales. The experiment was repeated with three samples (N = 3).

Bending experiment

Prototypes were designed for mechanical testing in bending with rigid rods at the top and bottom of the assembly and with scales aligned in orientation angles (φ) of 0°, 30°, 60°, 90°, 120°, and 150°, where φ is defined as the angle between the peg-and-socket axis and the loading direction. A control sample consisting of a solid TangoPlus sheet with no scales was printed as a control sample. The prototypes were experimentally tested in bending induced by axial compression on a mechanical tester (Zwick Z010, Zwick Roell, Germany) using load cells ranging from 20 to 2500 N. Sample holders were designed, 3D-printed, and affixed to the load cell with Permacel (Nitto) tape for pin–pin boundary conditions that allow rod rotation about the x-axis and constrain rotation about the y- and z-axes to prevent global twisting of the sample. The samples were induced to deform concavely (scales facing in) by setting an initial lateral deflection of 1 mm and zeroing the force before displacement-controlled compressive loading at a strain rate of 1 mm/s. The reaction force (F) vs. vertical displacement (d) for the prototypes was measured, and the experiment performed with three samples (N = 3) per orientation. Sample stiffness (K) was calculated as the slope of the loading curve and normalized by the stiffness of the control sample consisting of a 4.4 mm sheet of TangoPlus without scales (K = 7.22 N/m).

Finite element modeling

Three parts were designed for the model using FE software (ABAQUS, Dassault Systemes, France) and meshed with C3D4 (standard, linear stress) elements: a simplified scale geometry, a substrate material, and rigid rods. VeroWhite was modeled as an isotropic linear elastic material (E = 2.0 GPa, ν = 0.43, and ρ = 1.175 g/cm3) and assigned to the scales and rigid rods. TangoPlus was modeled as a Neo-Hookean hyperelastic material (C11 = 0.63 MPa and D1 = 10-6) and assigned to the substrate. The TangoPlus interconnections between the peg and socket of adjacent scales were modeled as springs with a stiffness of 0.63 MPa. Pin–pin boundary conditions were applied to the model to match the experimental conditions, and concave bending was simulated as a displacement-controlled compressive loading to induce lateral bending of the sample. Force–displacement (F–d) curves were generated as a measure of the reaction force vs. vertical displacement of the top rod at every increment. Values for stiffness (K) were calculated as the slope of the F–d data and normalized by the stiffness of a control model consisting of a 4.4 mm sheet of TangoPlus without scales (K = 10.7 N/m). Stresses (Mises, linear elastic, averaged) were captured through the whole model at every 10 increments.