Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Programmable gear-based mechanical metamaterials

## Abstract

Elastic properties of classical bulk materials can hardly be changed or adjusted in operando, while such tunable elasticity is highly desired for robots and smart machinery. Although possible in reconfigurable metamaterials, continuous tunability in existing designs is plagued by issues such as structural instability, weak robustness, plastic failure and slow response. Here we report a metamaterial design paradigm using gears with encoded stiffness gradients as the constituent elements and organizing gear clusters for versatile functionalities. The design enables continuously tunable elastic properties while preserving stability and robust manoeuvrability, even under a heavy load. Such gear-based metamaterials enable excellent properties such as continuous modulation of Young’s modulus by two orders of magnitude, shape morphing between ultrasoft and solid states, and fast response. This allows for metamaterial customization and brings fully programmable materials and adaptive robots within reach.

## Main

Materials featuring tunable elastic properties1,2 offer tremendous possibilities for smart machines, robots, aircraft and other systems3,4,5. For example, robotic systems with variable stiffness can adapt to missions like grabbing6 and jumping, or maintain optimal performance in a changeable environment7. However, elastic properties of conventional materials are barely tunable even if phase changes are induced. Mechanical metamaterials8,9,10,11 are artificial architected materials that exhibit properties beyond those of classical materials12,13,14,15,16. Most existing metamaterials integrate monofunctional load-bearing elementary structures (such as rods, beams or plates) in specified topologies with fixed or hinged nodes (Fig. 1). Reconfigurable metamaterials open possibilities for drastic changes in properties17,18,19. When stimulated by stress, heat or electromagnetic fields, reconfigurations in these metamaterials are induced by the formation of new contacts, buckling20,21,22 or rotating hinges23,24,25. Due to node constraints, this permits the reshaping among only a few stable states26,27 and often includes unstable states, which limits the tunability. Reducing the connectivity28 or relaxing the constraints (for example, with chiral structures9 or by connecting elements with flexural traps29) can enable more states to improve the shape-changing capability, but this inevitably deteriorates the robustness and structural stability that are essential for most applications. Moreover, reconfiguration, including the shape-memory effect, usually involves large deformation that either leads to irreversible plastic deformation or adversely competes with the commonly required high stiffness18. Although the chemical-responsive materials30,31 enable some in situ tunability, the regulation process of their elastic properties is usually very slow, just like for thermal-responsive materials19. Assembling rods mounted with gears into special lattices can improve the stability while preserving the rotatable nodes32, but engineering practical and robust metamaterials with continuously tunable elasticity, especially with fast in situ tunability in service, remains a major challenge.

## Design concept

Overcoming these challenges requires an unprecedented design paradigm. First, tunability may be realized by assembling elements with built-in stiffness gradients. Second, the coupling between elements must comply with large deformation. Achieving tunable yet strong solids requires ensuring tunability under large force and robust controllability while avoiding plastic deformation in tuning. We find that such a mutable-yet-strong coupling can be realized with gear clusters. Gears provide an ideal mechanism to smoothly transmit rotation and heavy compressive loads thanks to the reliable gear engagement (meshing). Stiffness gradients can be built into an individual gear body or realized with hierarchical gear assemblies. Gear clusters can be assembled into manifolds and can, as metacells, be periodically arranged to form metamaterials (Fig. 1c). The proposed design concept is very general since there exist numerous architectures for gear assembly. Exotic functionality and flexible tunability can emerge from the diversity of gear types, built-in variability and cluster organization. We create several metamaterial prototypes with different gear clusters to demonstrate this.

## Metamaterial based on Taiji gears

The first prototype is created using compactly coupled periodic gears and two lattice frames (front and back) to arrange the gears into a simple quadratic pattern (Fig. 2a). The plane gears contain hollow sections. The outer part forms two elastic arms whose radial thickness smoothly varies with the rotation angle θ (Fig. 2b and Supplementary Figs. 1 and 2). Subject to compressive loading, deformation in the arms is dominated by bending (Fig. 2c). The effective stiffnesses of both an arm karm and a pair of arms Kp depend on the angle θ (Supplementary Fig. 3). The homogenized Young’s modulus of the metamaterial along the y axis Ey = Kp/B + Ef has contributions from the gears and frame (Methods, ‘Equivalent method’ section). Here, Ef is the stiffness of the frame and B denotes the gear width. Ey is continuously tunable by rotating the gears and dominantly depends on Kp since min(Kp/B) Ef = 2.06 MPa here.

The tunability relies on the shape of the built-in hollow section. Among diverse choices, the shape inspired by the Chinese Taiji diagram (Fig. 2b), characterized by a spiral direction, can give smooth variation and polarity. The angle difference between the two local coordinates is β (Fig. 2g). The spin rotations are opposite in any two meshing gears. Also, the spiral directions of the Taiji patterns on the front and the back faces are reversed. Therefore, the meshing mode of a pair of gears has two polarities. When the spiral directions of patterns are opposite (Fig. 2b), the polarity is positive, labelled as P+(β). The meshing pair in Fig. 2g features negative polarity, P(β).

We employ finite element analysis (FEA) to simulate the contact problem in gear-based metamaterials (Methods and Supplementary Figs. 4 and 5). Contact nonlinearity becomes apparent for high Kp (Supplementary Fig. 3). Young’s modulus is evaluated from the slope of the uniaxial stress–strain curves at relatively large strain ε. An all-metallic prototype, consisting of 5 × 5 copper gears and steel frames, is manufactured and assembled with P+(3°) and P(15°) metacells, respectively (Supplementary Fig. 2). The gear has 60 teeth, with the tooth thickness tto = 0.35π mm, diameter D = 42 mm and width B = 20 mm. Measured cyclic loading–unloading curves show some hysteresis (Fig. 2e). This is ascribed to the sliding friction between the meshed teeth (Supplementary Fig. 6). Figure 2f,g demonstrates that experimental results of Ey(θ) are in excellent agreement with FEA. The modulation period is 180° in both the P+(3°) and P(15°) cases. The smooth Ey(θ) curve indicates that the obtainable stable states are dense and that continuous tunability is achieved. Both polarity and β affect the tunable range and the correlation between the tunable properties and θ (Fig. 2d). For P+(3°) in Fig. 2f, the zigzag curve of Ey(θ) reaches the maximum value Emax = 7.67 GPa at θ = 78°, where the solid parts of the gears are in contact, and a minimal value Emin = 0.102 GPa at θ = 114°, where the meshing connects the forearms. This experimentally obtained modulation range of Emax/Emin = 75 demonstrates the spectacular reconfigurability. For P(15°), Ey(θ) is sombrero-shaped with a tunable range of 33 from Emin = 0.156 GPa to Emax = 5.13 GPa. Since Ex(θ) = Ey(θ + 90°), the anisotropy in orthogonal directions also changes with θ. For example, the P+(3°) metamaterial can be continuously modulated from Ex = Ey to the maximum ratio (Ey/Ex)max = 24.8. The latter gives a metamaterial with negligible lateral expansion upon compression (Supplementary Fig. 7). Compared to existing designs, the node constraints in gear-based metamaterials are relaxed, but the connection stability and reconfiguration robustness are maintained at any θ even under large compressive loads (Supplementary Video 1). The design is also robust to accommodate manufacturing inaccuracies when regarding the angle β as an indicator of the alignment error of the gears. Figure 2d shows that the programmability is preserved even for large β.

The all-metallic metamaterial introduced above is manufactured by assembling individual gears. For scale-up and miniaturization, it is desirable to avoid the assembly of individual parts. Next, we demonstrate that integrated gear-based metamaterials can be directly manufactured with three-dimensional (3D) printing, even on the microscale. The major challenge for such integrated manufacturing is to guarantee that the meshing teeth are not fused together but still reliably engaged. To tackle this problem, a small clearance is reserved between the surfaces of the meshing teeth in the assembled digital model to overcome manufacturing errors (Methods). Here we manufacture an integrated micro metamaterial consisting of 5 × 6 Taiji gears (Fig. 2h) by adopting the projection micro-stereolithography 3D printing technique. The diameter and tooth thickness of the Taiji gear are 3.6 mm and 235 µm, respectively; the thickest arm is 75 µm (Supplementary Fig. 8). The micro gears are arranged with P+(0°), and the reserved minimal clearance between the teeth is 32 µm. The sample is made of a photosensitive resin with a Young’s modulus of 3.5 GPa. As experimentally demonstrated in Fig. 2i, the equivalent modulus Ey(θ) of this micro specimen can be smoothly tuned by 35 times (from 8.3 MPa to 295 MPa). Using this integrated design strategy, gear-based metamaterials could be scaled up in size and number of gears with appropriate high-resolution large-scale 3D printing facilities. Modulation of such integrated metamaterials can be achieved with distributed drives or motors (Supplementary Video 2 and Supplementary Fig. 9a).

## Metamaterial based on planetary gears

Obviously, this first metamaterial is tunable only under compressive loading. The tensile load is carried by the frame, and the tensile modulus is Et = Ef. One may also aim at strong metamaterials whose compressive and tensile moduli are both tunable while preserving structural integrity. This can be achieved by organizing a planetary gear system as a metacell (Fig. 3a). In this example, the metacell contains six gears: an inner-toothed ring gear (Supplementary Fig. 10), a central sun gear and two pairs of planetary gears A1–A2 and B1–B2. Gear centres A1–O–A2 (and B1–O–B2) are colinear. Using this gear cluster, we create a hierarchical and strong metamaterial whose tunability emerges from the relative rotation of the gears inside the metacell. The thickness of the ring tr is uniform. Neighbouring rings are rigidly connected in a quadratic lattice, which ensures structural integrity. Planetary gears revolve along the ring when rotating the sun gear by θsun. Their position is given by the revolution angle θpr = θsunrsun/(Rin + rsun), where Rin and rsun are the radii of the ring and sun gears, respectively (Supplementary Table 1 for parameter values). The teeth prevent relative slippage between the two gears even under tension. The metamaterial elastic properties are given by the effective stiffness of the annulus ring supported by the planetary gears that act as fulcrums (Fig. 3b). The position θpr of the planetary gears determines the elastic properties. We adopt the orthogonal relation A1A2 B1B2, which gives a large tunable range and symmetric behaviour with a modulation period of 90°. The tunable range can be further modified using the angle A1OB1. For the assembled metamaterial, all sun gears are connected to transmission gears by shafts (Fig. 3a), and those transmission gears are compactly coupled. Thereby robust reconfiguration of all metacell patterns can be achieved by rotating transmission gears.

Fx and Fy denote the compressive loads in the x and y directions, respectively (Fig. 3a). Under uniaxial compression (Fy > 0, Fx = 0), only the pair of planetary gears with an angle smaller than that of the loading axial (min(YOA1, YOB1) < 45°) supports the load (Fig. 3b). Stress in the other pair is zero. Conversely, under uniaxial tension, only the other pair is load-bearing. The two pairs exchange roles at θpr = 45°, and the material is orthogonally isotropic. This metamaterial presents a more remarkable compressive nonlinearity because four pairs of meshing teeth in a metacell bear loads. Both the compressive and tensile moduli Ec and Et reach maxima at θpr = 0, but $$E_{{{\mathrm{c}}}}^{{{{\mathrm{max}}}}} \gg E_{{{\mathrm{t}}}}^{{{{\mathrm{max}}}}}$$, and thus the static compressive-tension symmetry is broken (Fig. 3g). Moreover, at θpr = 45°, no stress is transmitted to the planetary gears (Fig. 3b); both moduli reach minima there, and $$E_{{{\mathrm{c}}}}^{{{{\mathrm{min}}}}} = E_{{{\mathrm{t}}}}^{{{{\mathrm{min}}}}}$$.

We fabricate three kinds of specimens using this strategy. An all-steel macro metamaterial is manufactured by assembling 3 × 3 metacells (Fig. 3c) with lattice constants ax = ay = 27 mm. The steel gears have a small tooth thickness tto = 0.15π mm, and Rin = 12 mm, rsun = 6 mm and tr = 1 mm. Integrated manufacturing of this prototype is more challenging than that of the Taiji pattern because there are two layers and every metacell possesses eight pairs of meshing teeth. The integrated prototype can also be directly manufactured by 3D printing, at both macro and micro scales (Methods, Supplementary Figs. 11 and 12 and Supplementary Video 3 for more details). We print a 6 × 6 macro specimen (Fig. 3d) using a polymer with a Young’s modulus of 2.5 GPa, and print a 3 × 4 micro specimen (Fig. 3e) using a resin with a Young’s modulus of 3.5 GPa. The size and the number of metacells are limited by the capability of the 3D printer rather than the design strategy. The micro polymer sample (Fig. 3e) has Rin = 2.4 mm, rsun = 0.6 mm, tr = 0.3 mm and ax = ay = 5.4 mm, with a tooth width and height of 135 µm and 225 µm, respectively.

The experimental results are consistent with the FEA simulations for all specimens (Fig. 3g–i and Supplementary Fig. 11). In this strong hierarchical metamaterial, we can smoothly tune the compressive modulus Ec of the macro metallic specimen by 46 times (5.2–0.11 GPa), the macro polymer specimen by 55 times (69–1.25 MPa; Supplementary Fig. 12b) and the micro specimen by 25 times (100–4 MPa). Meanwhile, their tensile modulus Et can be tuned by 5 times (0.52–0.11 GPa), 5.6 times (7–1.25 MPa) and 5 times (20–4 MPa), respectively. In Fig. 3h, some differences between experiment and FEA near θpr = 45° arise from the boundary conditions (Supplementary Fig. 13). The in situ tunability combined with the reasonably large moduli in tension and compression as well as large shear rigidity makes this metamaterial design particularly robust and strong, yet tunable. Furthermore, the metamaterials can be synchronously controlled with distributed motors at both the macro and micro scales (Supplementary Video 4).

## Mechanisms for stability

Interestingly, the metamaterial in Fig. 2a (a discrete gear lattice with a very soft frame) remains stable under compressive stresses and shows large rigidity in shear. One of the contributing factors underpinning the observed stability stems from the non-uniform loading of the meshing teeth at different points, which leads to bending deformations that tightly grip the teeth together (Supplementary Fig. 5). The relatively large shear modulus of the metamaterial, G = Gg + Gf, is composed of the shear moduli generated by gears (Gg) and by frames (Gf = 1.04 MPa). Shear force induces both spin and the planetary rotation of gears. For a pair of gears, the relative planetary rotation leads to zero shear resistance, Gg = 0, giving a highly unstable state (Fig. 4a). However, in a group of four gears (shown in Fig. 4b), shear stress τ induces mutual locking of the planetary rotation by the opposite spin of the neighbouring gears, which is referred to as shear interlock. We calculate the shear stiffness of the metacell with periodic boundary conditions and the finite n × n gear lattice (Supplementary Figs. 1416). Owing to the shear interlock, the shear modulus is large but only marginally tunable (Supplementary Fig. 16), which is demonstrated by the measured generalized shear stiffness Kshear/B of the finite 3 × 3 architecture in Fig. 4c.

## Gear metamaterial for shape morphing

The programmability of gear-based metamaterials is not limited to elastic constants. Removing every second gear in every second row of the metamaterial in Fig. 2 can release the shear interlock (inset in Fig. 4f) to generate a state with Gg = 0. The effective shear modulus is then determined solely by the low stiffness frame, and the metamaterial can be considered as ultrasoft matter. The vanishing shear modulus enables complex deformation modes (Fig. 4d,e and Supplementary Video 5), conducive to shape morphing. To verify this, an ultrasoft prototype consisting of 4 × 4 metacells with rubber frames is manufactured and tested (Fig. 4f). The gears are made of aluminium alloy (Supplementary Fig. 17). Shear tests on the prototype give a tiny modulus of G = 21.52 kPa. Independent measurement of the frames gives Gf = 21.11 kPa, so that Gg = GGf is indeed negligibly small. Moreover, the modulus G remains tiny until the shear strain γ reaches 25%, at which point the semi-free gear interacts with two other gears, which builds a new meshing connection among the three gears. The new connection supports high shear stresses and leads to a sharp rise in G, switching the soft matter to a stiff solid. The resulting solid represents a geometrically interlocked state (Methods). Oscillations of the shear stress in Fig. 4c arise from the critical meshing state among the three gears before they interlock (Supplementary Fig. 18). The shear strain at which the geometrical interlock occurs can be adjusted by the size of the neighbouring gears, which in turn determines the limiting (strong) states of a shape-morphing structure. Previous metamaterial designs with vanishing shear modulus, like pentamode metamaterials33, show vanishing moduli only at small strains and in extremely fragile structures.

## Potential applications

Conventional machines generally rely on materials with constant stiffness and therefore show constant stiffness themselves. The designed stiffness is then a compromise among stability, safety, efficiency and performance, thus hindering the pursuit of the best performance and efficiency in variable environments. Programmable materials featuring tunable elastic properties, including active mechanical metamaterials, are much anticipated in intelligent machines and systems1,2. Here we offer a comparison of typical material designs from the literature2,19.

The response time, stability, force and energy required for property changes are all critical attributes for variable-stiffness structures4. We take the strain and response time required to accomplish a tunable period of material stiffness as metrics to position our gear-based metamaterials among the existing active materials (Fig. 5). Shape-morphing metamaterials18,34 enable tunability between or among two or a few stable states. The achieved tunability is non-continuous and requires a large deformation (ε ≈ 30%). Thermal-responsive composites35 made of shape-memory alloys or polymers may give a continuously tunable modulus. However, they require a long response time, and some suffer from nearly 100% strain36,37. Chemical-responsive materials30,31 containing hydrogels can offer continuous and in situ tunability, but also require hours of response time. Conventional magneto- or electro-responsive metamaterials38 based on elastomers or magnetorheological fluids can give fast, continuous, but narrow tunability, which usually requires a high active voltage (~5 kV) and complex facilities19. Our gear-based metamaterials in Figs. 2 and 3 can offer a fast response and the desired broad-range, continuous and in situ tunability of stiffness.

We propose several scenarios to showcase the broad application potential of the proposed gear-based metamaterials in Supplementary Figs. 1921 and Supplementary Table 2. For robots, a tunable-stiffness leg/actuator can offer high stiffness to stably support a heavy load while walking and a low stiffness for shock protection while jumping or running4,39 (Supplementary Figs. 19 and 20). A similar tunable-stiffness isolator is desired in the aero-engine pylon system to maintain the best performance and efficiency at different flight stages (Supplementary Fig. 21). Moreover, the fast-response gear-based metamaterial may give rise to a sensitive variable-stiffness skin, which has been attracting wide attention40. Furthermore, resonators with tunable stiffness are critical components in programmable metamaterials for wave manipulation41,42. Therefore, gear-based programmable metamaterials can aid in the realization of extensive intelligent machines. In contrast with conventional methods, the programmability enabled by a gear-based metamaterial does not require large deformation and heavy controlling systems, such as hydraulic/pneumatic or magnetic systems, and thus benefits the miniaturization and integration of machines and can even be used in harsh environments such as outer space.

## Conclusions

We show that gear-based mechanical metamaterials provide in situ tunability while preserving stability, strength and high load-bearing capacity. The programmability is robust and easily implementable. Gear clusters provide a vast design space that permits customizable performance of the metamaterials. Besides the demonstarted Young’s modulus, shape morphing and shock protection, the tunability can be extended to other elastic properties like shear modulus, Poisson’s ratio, strength, deformation modes and even damping coefficient (Supplementary Fig. 21). One can also envision 3D metamaterials by using bevel gears, assembling planar gears into hierarchical configurations as in Fig. 3 or synthesizing different types of gear (Supplementary Fig. 22). Integrated manufacturing bridges these tunable properties to produce robust multipurpose devices43. With the example of micro metamaterials, further miniaturization and an extension of gear-based metamaterials are possible with high-resolution and large-scale 3D printing.

In conclusion, this work proposes and demonstrates an unconventional design paradigm for programmable dynamic metamaterials via the mutable-yet-strong coupling and built-in variability of gears. We establish the general concept, conceive prototypes, conduct mechanical analyses, demonstrate the flexible tunability and integrated manufacturing at both the macro and micro scales and showcase the broad potential applications. The proposed design paradigm broadens the horizon for designing fully programmable materials, thus offering an impetus to their exploration for practical applications.

## Methods

### Integrated manufacturing

The printer used for the projection micro-stereolithography micro metamaterial fabrication is a BMF NanoArch S130, with a precision of about 5 µm. The material used in microscale 3D printing is a photosensitive resin with a Young’s modulus of about 3.5 GPa. The manufacturing process for the integrated micro metamaterial sample consisting of 5 × 6 Taiji gears follows three steps. First, the assembled gears are printed on a baseplate; those gears are adhered to the plate. Second, the sample is wrapped in a box to constrain the motion of the gears (Supplementary Fig. 9a). Last, everything including the box is removed from the plate. The box with a frame helps maintain the relative angle of the assembly in the removal process.

For the metamaterial based on planetary gears, the layer of the planetary gears is printed first (Supplementary Fig. 12a and Supplementary Video 3). Except for the preserved clearance, the connection shaft between the transmission gear and the sun gear is conical at both the macro and micro scales (Supplementary Fig. 10), which ensures that every printing part, especially the teeth of the transmission gears, is tightly attached on the formed structure. Otherwise, the teeth could move and then fuse together during the printing. At the macro scale, the integrated model is printed with two photosensitive resins using polymer injection with the printer Stratasys Objet260, with a precision of about 50 µm. The stiff model material is wrapped in the soft, soluble support material. The metamaterial acquires the targeted tunability after removing the support material. At the microscale, the material is immersed in the fluid resin during the printing. No support material/structure is required for this model owing to the conical shaft and high precision of projection micro-stereolithography. The sample shown in Fig. 3d is printed with a resin (polymer) with a Young’s modulus of 2.5 GPa.

In integrated manufacturing, the clearance reserved between the surfaces of meshing teeth in the assembled digital model depends on the precision of the printer, the structure and the materials. The minimal clearance Δ should be higher than the printer’s precision p (manufacturing errors) but much smaller than the tooth height (ht = 2.25m for standard gears), where the gear module m denotes the ratio between the gear diameter D and the number of teeth z, m = D/z (see Supplementary Text). Here, the minimal clearance between meshing teeth in all macro specimens printed with Objet260 is 86 µm. The minimal clearance for the micro metamaterial consisting of Taiji gears is set to 32 µm, and that for the micro planetary gear-based metamaterial is 21 µm. These clearances are sufficient to alleviate the manufacturing uncertainties to keep the meshing teeth separated but reliably engaged. Based on our 3D printers and tests, we suggest Δ > 1.5p and Δ ≤ ht/10 = 0.225m. This requires p < 0.15m, which helps us determine the required precision scale with a specified gear size.

### Actuation

As shown in Supplementary Video 2, we prepare a microscale sample consisting of 5 × 5 Taiji gears to show its actuation process. They are embedded into a box, and those gears connect to the frames through micro shafts. The sample is synchronously driven by four d.c. brushless motors (8 mm diameter) connected to the 1 × 1st, 1 × 4th, 4 × 1st and 4 × 4th gears. Here n × m denotes the position at the nth row and mth column in the array. As shown in Supplementary Video 4, the macro metamaterial in Fig. 3d is synchronously actuated by four-step motors whose diameter is 20 mm. These motors are synchronously controlled by an electronic controller. The revolving speed of the step motor depends on the impulse frequency generated by the controller. Similarly, the micro sample in Fig. 3e is put in a box and actuated by five micro step motors whose diameter is 5 mm. The controller is identical to the one used for the macro sample.

### FEA

FEA simulations are carried out with the commercial software ANSYS. We compare the accuracy of different finite element models, including two-dimensional (2D), 3D, linear and nonlinear models. The plane stress state is considered in the 2D model. In the linear models, the meshing points of gears are bonded by fixing together the two surfaces in contact, resulting in a linear stress–strain relationship. In the nonlinear models, the size of the contact area on the tooth surface at the meshing points depends on the load, and there is a relative sliding between the contact surfaces. The sliding induces frictional damping if the coefficient of friction is non-zero. We also use a simplified model by removing all teeth, where the contact between two gears becomes that between two cylinders. In principle, the 3D nonlinear model should be the most realistic representation of the experimental set-up. Supplementary Fig. 4 demonstrates that the 2D nonlinear model is in excellent agreement with the 3D nonlinear model. The two linear models produce a large discrepancy with the nonlinear ones, although they still can capture the general variation trend. The simplified model approximately presents the standard results. To enhance the simulation efficiency, we use the 2D nonlinear models in most cases. The 3D model is adopted only when considering the frictional contact.

Our metamaterials embrace a periodic architecture. To evaluate the homogenized elastic and shear moduli, ideal periodic boundary conditions are applied on the unit cell in the FEA. Boundary conditions depend on the deformation mode of the unit cell. The homogenized strain vector is ε = (εx, εy, γ). These strains are realized by enforcing the displacement fields (u, v) in the plane stress state.

As explained in Supplementary Fig. 14, two types of boundary condition are considered when calculating the shear modulus in the shear interlock state. To show the shear state of a finite n × n gear lattice, we fix the lower row/column of the gears and apply a displacement field to the upper row/column of the gears. As a second method, periodic boundary conditions are applied on a metacell to calculate the shear modulus. These periodic boundary conditions present the shearing state ε = (0, 0, γ). In both cases, the strain energy density W = 2/2 is extracted to evaluate the shear modulus G. For the n × n finite structure without periodic boundary conditions, although the equation of the generalized shear stiffness G′ = Kshear/B is the same as the formula for shear modulus G = τ/γ, the value of G′ may not equal the real shear modulus G (Supplementary Fig. 15) due to the free edge effects in finite structures (Supplementary Figs. 13a and 14a).

For the metamaterial based on a planetary gear system, the load is applied on the four blocks of the ring. For a metacell in FEA, we specify the uniaxial deformation v = εyay and make εx free for solving Ey.

### Equivalent method

For the metamaterial based on Taiji gears, the deformation mode for meshing gears can be represented by the overall stiffness of a pair of meshing elastic arms Kp = 1/(1/karm1 + 1/karm2 + 1/ktooth) (Supplementary Fig. 3 for their definitions). The stiffnesses of the two arms karm1 and karm2 are independent of the compressive deformation. As shown in Supplementary Fig. 5, the meshing of a pair of teeth features a line of contact on their surfaces. With compression, a small contact area is generated near the line where sliding occurs during the process. Therefore, the contact stiffness of the teeth ktooth depends on the contact pressure on the involute teeth. A high pressure leads to significant contact nonlinearity and results in a dependence of Kp on the displacement/load. By contrast, deformation mainly occurs in the elastic arm rather than the teeth if karmktooth, and Kp is constant in this case. The homogenized Young’s modulus in the y direction of the metamaterial is Ey = Kp/B + Ef. The equivalent methods for shear modulus are explained with Supplementary Fig. 14.

For the metamaterial consisting of a periodic planetary gear cluster, the Young’s modulus depends on the deformation of the ring. The influences of contact nonlinearity between teeth on Ey are the same as described above.

### Geometrical interlock

In a meshing pair, the rotation directions of the driving and driven gears are opposite. In a group of gears, if every meshing is viewed as a connection line, n gears form a closed polygon as shown in Fig. 4f. If n is odd, spin rotation is incompatible, leading to the locking among the gears. This meshing state is referred to as geometrical interlock.

### Mechanical tests for Young’s modulus

When measuring the Young’s modulus Ey in the metamaterial based on Taiji gears, a compressive load Fy is applied and released from the top of the prototype in Fig. 2a. We control the strain ε for different θ to overcome clearance nonlinearity while avoiding plastic deformation. The rotation angle θ is manually controlled. Similar cyclic loading–unloading tests are performed for the measurement of the shear modulus.

The experimental setting for the test on the metamaterial based on planetary gear systems is shown in Supplementary Fig. 9. When measuring the compressive modulus, a compressive load is applied on the top and the bottom blocks on the rings; when measuring the tensile modulus, we fix the tails on the sample to a pair of clamps and apply tensile loads through the tails.

As shown in Figs. 2e and 3g, the cyclic loading–unloading process features high repeatability, thus testifying to the experimental accuracy. Moduli Ey and G are both calculated as the slope around the maximum ε. The initial cycle is excluded when fitting Ey and G. The choice of the strain interval for the slope calculation affects the final modulus value. The error bars and the average values are evaluated by choosing different intervals along the curve.

### Mechanical tests for shear stiffness

For the metamaterial based on Taiji gears, a sample consisting of 3 × 3 gears and steel frames is manufactured for the measurement of the shear stiffness in the shear interlock state, as shown in Supplementary Fig. 15. A fixture apparatus is fabricated to obtain the shearing state. For the shape-morphing metamaterial, the sample is put in two right-angle grooves, and the load from the testing machine directly transfers to the sample.

## Data availability

The main data and models supporting the findings of this study are available within the paper and Supplementary Information. Further information is available from the corresponding authors upon reasonable request.

## References

1. McEvoy, M. A. & Correll, N. Materials that couple sensing, actuation, computation, and communication. Science 347, 1261689 (2015).

2. Levine, D. J., Turner, K. T. & Pikul, J. H. Materials with electroprogrammable stiffness. Adv. Mater. 33, 2007952 (2021).

3. Oztemel, E. & Gursev, S. Literature review of Industry 4.0 and related technologies. J. Intell. Manuf. 31, 127–182 (2020).

4. Wolf, S. et al. Variable stiffness actuators: review on design and components. IEEE ASME Trans. Mechatron. 21, 2418–2430 (2016).

5. Cully, A., Clune, J., Tarapore, D. & Mouret, J. Robots that can adapt like animals. Nature 521, 503–507 (2015).

6. Shintake, J., Cacucciolo, V., Floreano, D. & Shea, H. Soft robotic grippers. Adv. Mater. 30, 1707035 (2018).

7. Barbarino, S., Bilgen, O., Ajaj, R. M., Friswell, M. I. & Inman, D. J. A review of morphing aircraft. J. Intell. Mater. Syst. Struct. 22, 823–877 (2011).

8. Pham, M. S., Liu, C., Todd, I. & Lertthanasarn, J. Damage-tolerant architected materials inspired by crystal microstructure. Nature 565, 305–311 (2019).

9. Frenzel, T., Kadic, M. & Wegener, M. Three-dimensional mechanical metamaterials with a twist. Science 358, 1072–1074 (2017).

10. Kadic, M., Milton, G. W., Van Hecke, M. & Wegener, M. 3D metamaterials. Nat. Rev. Phys. 1, 198–210 (2019).

11. Fernandes, M. C., Aizenberg, J., Weaver, J. C. & Bertoldi, K. Mechanically robust lattices inspired by deep-sea glass sponges. Nat. Mater. 20, 237–241 (2021).

12. Zheng, X. et al. Ultralight, ultrastiff mechanical metamaterials. Science 344, 1373–1377 (2014).

13. Berger, J. B., Wadley, H. N. G. & McMeeking, R. M. Mechanical metamaterials at the theoretical limit of isotropic elastic stiffness. Nature 543, 533–537 (2017).

14. Fernandez-Corbaton, I. et al. New twists of 3D chiral metamaterials. Adv. Mater. 31, 1807742 (2019).

15. Xu, X. et al. Double-negative-index ceramic aerogels for thermal superinsulation. Science 363, 723–727 (2019).

16. Fang, X., Wen, J., Bonello, B., Yin, J. & Yu, D. Ultra-low and ultra-broad-band nonlinear acoustic metamaterials. Nat. Commun. 8, 1288 (2017).

17. Faber, J. A., Arrieta, A. F. & Studart, A. R. Bioinspired spring origami. Science 359, 1386–1391 (2018).

18. Chen, T., Pauly, M. & Reis, P. M. A reprogrammable mechanical metamaterial with stable memory. Nature 589, 386–390 (2021).

19. Qi, J. et al. Recent progress in active mechanical metamaterials and construction principles. Adv. Sci. 9, 2102662 (2022).

20. Javid, F. et al. Mechanics of instability-induced pattern transformations in elastomeric porous cylinders. J. Mech. Phys. Solids 96, 1–17 (2016).

21. Yang, Y., Terentjev, E. M., Wei, Y. & Ji, Y. Solvent-assisted programming of flat polymer sheets into reconfigurable and self-healing 3D structures. Nat. Commun. 9, 1906 (2018).

22. Fu, H. et al. Morphable 3D mesostructures and microelectronic devices by multistable buckling mechanics. Nat. Mater. 17, 268–276 (2018).

23. Coulais, C., Sabbadini, A., Vink, F. & van Hecke, M. Multi-step self-guided pathways for shape-changing metamaterials. Nature 561, 512–515 (2018).

24. Lipton, J. I. et al. Handedness in shearing auxetics creates rigid and compliant structures. Science 360, 632–635 (2018).

25. Silverberg, J. L. et al. Using origami design principles to fold reprogrammable mechanical metamaterials. Science 345, 647–650 (2014).

26. Frenzel, T., Findeisen, C., Kadic, M., Gumbsch, P. & Wegener, M. Tailored buckling microlattices as reusable light-weight shock absorbers. Adv. Mater. 28, 5865–5870 (2016).

27. Jenett, B. et al. Discretely assembled mechanical metamaterials. Sci. Adv. 6, c9943 (2020).

28. Overvelde, J. T. B., Weaver, J. C., Hoberman, C. & Bertoldi, K. Rational design of reconfigurable prismatic architected materials. Nature 541, 347–352 (2017).

29. Shaw, L. A., Chizari, S., Dotson, M., Song, Y. & Hopkins, J. B. Compliant rolling-contact architected materials for shape reconfigurability. Nat. Commun. 9, 4512–4594 (2018).

30. Auletta, J. T. et al. Stimuli-responsive iron-cross-linked hydrogels that undergo redox-driven switching between hard and soft states. Macromolecules 48, 1736–1747 (2015).

31. Li, T. et al. ‘Freezing’, morphing, and folding of stretchy tough hydrogels. J. Mater. Chem. B 5, 5726–5732 (2017).

32. Meeussen, A. S., Paulose, J. & Vitelli, V. Geared topological metamaterials with tunable mechanical stability. Phys. Rev. X 6, 41029 (2016).

33. Bückmann, T., Thiel, M., Kadic, M., Schittny, R. & Wegener, M. An elasto-mechanical unfeelability cloak made of pentamode metamaterials. Nat. Commun. 5, 4130 (2014).

34. Fang, H., Chu, S. A., Xia, Y. & Wang, K. Programmable self-locking origami mechanical metamaterials. Adv. Mater. 30, 1706311 (2018).

35. Shan, W., Lu, T. & Majidi, C. Soft-matter composites with electrically tunable elastic rigidity. Smart Mater. Struct. 22, 85005 (2013).

36. Yang, C. et al. 4D printing reconfigurable, deployable and mechanically tunable metamaterials. Mater. Horiz. 6, 1125–1244 (2019).

37. Xin, X., Liu, L., Liu, Y. & Leng, J. 4D printing auxetic metamaterials with tunable, programmable, and reconfigurable mechanical properties. Adv. Funct. Mater. 30, 2004226 (2020).

38. Jackson, J. A. et al. Field responsive mechanical metamaterials. Sci. Adv. 4, u6419 (2018).

39. Vanderborght, B. et al. Variable impedance actuators: a review. Robot. Auton. Syst. 61, 1601–1614 (2013).

40. Lin, X. et al. Ultra-conformable ionic skin with multi-modal sensing, broad-spectrum antimicrobial and regenerative capabilities for smart and expedited wound care. Adv. Sci. 8, 2004627 (2021).

41. Jones, M. R., Seeman, N. C. & Mirkin, C. A. Programmable materials and the nature of the DNA bond. Science 347, 1260901 (2015).

42. Wang, W. et al. Active control of the transmission of Lamb waves through an elastic metamaterial. J. Appl. Phys. 128, 65107 (2020).

43. Begley, M. R., Gianola, D. S. & Ray, T. R. Bridging functional nanocomposites to robust macroscale devices. Science 364, eaav4299 (2019).

## Acknowledgements

This research was funded by the National Natural Science Foundation of China (projects no. 12002371 and no. 11991032), the Hong Kong Scholars Program, the Fraunhofer Cluster of Excellence ‘Programmable Materials’ and the Excellence Cluster EXC 2082 ‘3D Matter Made to Order’ (3DMM2O) in Germany.

## Funding

Open access funding provided by Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V.

## Author information

Authors

### Contributions

X.F. and P.G. designed the study. X.F. conceived the idea and performed the experiments. X.F. and P.G. carried out the numerical simulations. L.C., D.Y. and H.Z. analysed the data. All authors interpreted the results. X.F., L.C., J.W. and P.G. wrote the manuscript with input from all authors. P.G. supervised the study.

### Corresponding authors

Correspondence to Xin Fang, Jihong Wen or Peter Gumbsch.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

Nature Materials thanks Amir A. Zadpoor and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

### Supplementary Information

Supplementary Figs. 1–23 with explanatory text, Tables 1 and 2 and legends for Videos 1–5.

### Supplementary Video 1

Introduction of the metamaterial based on Taiji gears. This video shows the robust tunability and high stability of the proposed metamaterial under a large force.

### Supplementary Video 2

Actuation of the micro metamaterial consisting of Taiji gears.

### Supplementary Video 3

Integrated 3D printing process of the macro metamaterial consisting of a planetary gear system.

### Supplementary Video 4

Actuation of the macro metamaterial consisting of planetary gears.

### Supplementary Video 5

Shape-morphing metamaterial. This video shows the structure, manufacturing process and protected shape morphing of this metamaterial.

## Rights and permissions

Reprints and Permissions

Fang, X., Wen, J., Cheng, L. et al. Programmable gear-based mechanical metamaterials. Nat. Mater. (2022). https://doi.org/10.1038/s41563-022-01269-3

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41563-022-01269-3