Abstract
Mechanical metamaterials are usually designed to exhibit novel properties and functionalities that are rare or even unprecedented. What is common among most previous designs is the quasi-static nature of their mechanical behavior. Here, we introduce a previously unidentified class of strain rate-dependent mechanical metamaterials. The principal idea is to laterally attach two beams with very different levels of strain rate-dependencies to make them act as a single bi-beam. We use an analytical model and multiple computational models to explore the instability modes of such a bi-beam construct, demonstrating how different combinations of hyperelastic and viscoelastic properties of both beams, as well as purposefully introduced geometric imperfections, could be used to create robust and highly predictable strain rate-dependent behaviors of bi-beams. We then use the bi-beams to design and experimentally realize lattice structures with unique strain rate-dependent properties including switching between auxetic and conventional behaviors and negative viscoelasticity.
INTRODUCTION
The story of mechanical metamaterials started with researchers trying to achieve unusual properties through the rational design of architected materials. The quest that started with a relatively limited scope has, however, rapidly pullulated and now targets not only such well-established behaviors as auxeticity (1–3), negative thermal expansion (4, 5), snap-through instability (6, 7), and ultrahigh properties (8, 9) but also more novel functionalities including shape-morphing (10, 11), nonreciprocity in transferring motions (12), and even logic gate–like programmability (13, 14). In a way, the boundary between materials and machines is becoming more and more blurred, with some researchers even introducing concepts like “machine matter” (15, 16).
In this context, flexible metamaterials (2, 17–19) show the most promise. On the one hand, the large, reversible, and nonlinear deformations exhibited by these materials provide a fertile ground for cultivating unique functionalities. On the other hand, flexible mechanical metamaterials share multiple organic intersections with other current areas of research including soft robotics (9, 20, 21) and flexible electronics (22–24). They have, therefore, received the most attention during past few years.
Most previous studies have, however, focused on the quasi-static response of mechanical metamaterials. In the cases where the dynamic response of mechanical metamaterials is studied, the focus is usually on their acoustic and wave propagation properties [e.g., see (25–30)]. However, the viscoelastic and strain rate–dependent behaviors of mechanical metamaterials, particularly those exhibiting substantial flexibility, could present many novel opportunities. This is particularly important given the fact that most of the present designs of flexible mechanical metamaterials either are made from polymers that exhibit vastly strain rate–dependent behavior (31) or include other energy-dissipating mechanisms that often lead to strain rate dependency. For example, there is substantial friction and contact in origami-based flexible metamaterials (10, 24, 32) that causes energy dissipation and strain rate dependency. It is only natural that the strain rate–dependent behaviors of flexible mechanical metamaterials are taken into the account, or even better exploited to further expand the space of achievable properties and functionalities. However, very few [e.g., see (31, 33)] major attempts in that direction can be found in the literature.
Here, we used the concept of bi-beams to design, analyze, and fabricate metamaterials whose mechanical behavior either switches (e.g., from conventional to auxetic) depending on the applied strain rate or exhibits new types of viscoelastic behavior (e.g., negative viscoelasticity). Bi-beams are made by laterally attaching two beams made from two different materials of which one is both highly deformable and highly strain rate–dependent (i.e., visco-hyperelastic), while the other is highly deformable but largely strain rate–independent (i.e., hyperelastic) (Fig. 1A). Under quasi-static conditions, the stiffness of the hyperelastic beam is higher than that of the visco-hyperelastic beam. When subjected to high enough compressive forces, such a bi-beam construct will predictably buckle to either the left or the right side, depending on the rate of the applied force (strain) (Fig. 1, A and B). To study the behavior of a single bi-beam, let us assume that the hyperelastic layer is always positioned at the left side of a bi-beam that is clamped at its both ends (Fig. 1, A and B).
(A) A schematic drawing of a bi-beam made of two identical beams with limited viscoelasticity (the hyperelastic beam) and substantial viscoelasticity (the visco-hyperelastic beam). Because of the higher stiffness of the hyperelastic beam, the bi-beam buckles to the left when loaded slowly in compression. Under high strain rates, the apparent stiffness of the visco-hyperelastic beam may exceed that of the hyperelastic beam, causing the bi-beam to buckle to the right. (B) Visualization of the stress distributions corresponding to the three potential states of sufficiently compressed bi-beams (i.e., left buckling, no buckling, or right buckling). (C) Maps showing the computationally predicted buckling direction for a bi-beam with an aspect ratio of 8 and different values of the mechanical properties of the hyperelastic and visco-hyperelastic beams. The results are presented for different values of the applied nondimensional strain rate, and the color code indicates the regions where right buckling occurs. (D) Evolution of the stresses and strains of a visco-hyperelastic material as a function of the applied nondimensional strain rate for four different sets of material properties. Analytical relationships are used for drawing the stress-strain curves and compared with the corresponding hyperelastic materials.
RESULTS AND DISCUSSION
Analytical and computational models of perfect bi-beams
Neo-Hookean models were used for representing the nonlinear elastic behaviors of the hyperelastic (Neo-Hookean constant, C10H) and visco-hyperelastic (long-term Neo-Hookean constant,
The computationally predicted map of the buckling direction (for a beam with an aspect ratio of 8) is in agreement with the predictions of the analytical model and shows that when
Effects of bi-beam geometry
To study the effects of bi-beam geometry on its buckling behavior, we built computational models corresponding to bi-beams with smaller and larger aspect ratios (Fig. 2A). For aspect ratios larger than 8, the buckling behavior of the bi-beams was found to be less predictable. Moreover, the higher modes of buckling appeared in some cases (depending on the balance between the material parameters and the nondimensional strain rate). According to the Euler’s buckling theory, as the length of the beams, l, increases, the critical strain decreases quadratically. Because the lateral expansion of the beams is related to the axial strain through the incompressibility constraint, a quadratically smaller value of the critical strain also means a quadratically smaller value of the lateral deformation, w. The force exerted by an elastic foundation, Fel, is proportional to the lateral deformation of the beam it supports (Fel = kw). A much smaller value of the lateral deformation, therefore, decreases the force exerted by the elastic foundation. Consequently, the relative importance of other minor forces (e.g., the slight forces caused by asymmetric meshes or small asymmetries in the numerical solutions) will increase, making the prediction of the buckling direction more challenging. In the theory of beams on elastic foundations, the parameter
(A) Computationally predicted maps showing the buckling direction for bi-beams with four different aspect ratios. The color code indicates the regions where right buckling occurs as a consequence of applying the corresponding nondimensional strain rate. The gray regions are the regions with conflicting results, while the orange regions highlight the occurrence of some higher modes of buckling. (B) Sensitivity of the computational results to geometric imperfections (imp.). The geometric imperfections originating from the first buckling mode (calculated using a linear buckling analysis) are applied to the geometry of bi-beams either to the left side (negative values of the buckling mode) or to the right side (positive values of the buckling mode). The predictions of buckling to the right direction are presented for different magnitudes of the applied imperfection. (C) Modification of the bi-beams design to minimize their sensitivity to geometric imperfections.
Our computational models showed the diminished predictability of the buckling direction and the appearance of the higher modes of buckling for some higher values of the aspect ratio (i.e., 12 and 16, Fig. 2A). Regardless of whether global buckling or the higher modes of local buckling occur, we decided to discard aspect ratios above 8 for further analysis and the design of cellular structures. For a smaller value of the aspect ratio (i.e., 4), the sensitivity of the buckling direction to the applied nondimensional strain rate increased, meaning that the region within which right buckling occurs enlarges more prominently as a consequence of increasing the applied nondimensional strain rate (Fig. 2A). This is a positive development in terms of affording greater flexibility in the design of strain rate–dependent metamaterials. However, bi-beams with smaller aspect ratios show higher critical forces that may not be desirable, as we ultimately aim to use the bi-beams for designing cellular structures. If the forces required for the local buckling of the bi-beams are too high, the global buckling of the entire cellular solid starts to compete with the local buckling modes of individual bi-beams, rendering the buckling behavior of the cellular structure less predictable. An intermediate aspect ratio (e.g., 8 in the case of our simulations) balances the two abovementioned aspects to make the buckling behavior of the bi-beams as predictable as possible.
Imperfect bi-beams and purposefully introduced imperfections
Geometric imperfections and equivalent material imperfections are the other important parameters influencing the buckling direction of the bi-beams. In the buckling analyses presented so far, we assumed perfect beam geometries. To analyze the effects of imperfection size on the buckling direction predicted by our computational models, we performed a sensitivity analysis in which a gradually increasing percentage of the eigenvector corresponding to the first buckling mode was introduced as a geometric imperfection to the shape of the bi-beam. The imperfect geometries were then used to perform nonlinear buckling analyses. Only for very small values of the imperfection sizes (i.e., ≤0.001 mm or ≤0.02% of the width of the bi-beams), the buckling direction predicted by the perfect models remained fully valid (Fig. 2B). For an imperfection size of 0.01 mm (=0.2% of the width of the bi-beams), the buckling direction was already affected for some combinations of the material properties and nondimensional strain rates, while the entire switching effect was lost for an imperfection size as small as 0.1 mm (Fig. 2B). These results show the extremely high degree of sensitivity of the buckling direction to geometric (or equivalent material) imperfections. We used two strategies to mitigate the effects of imperfection dependency. As we aim to compare our computational results with experimental observations, we will use the dimensional version of the strain rate,
The first strategy was to change the design of our bi-beams. When studying the evolution of the stresses induced in both hyperelastic and visco-hyperelastic beams, we found that the nonuniform distribution of the stresses caused by geometric imperfections results in stress concentrations at the interface of both beams (Fig. 2C). For example, when the direction of the geometric imperfection is to the left (e.g., in a design with
The second strategy was the purposeful introduction of geometric imperfections into the shape of the bi-beams to predispose them such that they tend to buckle to the left (i.e., the side of the hyperelastic beam) under low strain rates. This predisposition ensures that the buckling direction of a bi-beam is robustly predictable under low strain rates, regardless of any random imperfections that may result from the applied fabrication process. However, the purposefully introduced imperfections change the shape of the regions within which switching to right buckling could take place (Fig. 3). There are two changes in the shape of the switching zones. First, the size of the switching zones decreases as the imperfection size increases from 0.1 to 0.5 mm (corresponding to 2 to 10% of the width of the bi-beam) (Fig. 3). A smaller size of the switching zone is caused by the fact that the elevated stiffness of the visco-hyperelastic beam should overcome not only the higher stiffness of the hyperelastic beam but also the forces created because of the presence of the purposefully introduced geometric imperfections. Second, the boundary between both zones of left and right buckling, which was initially diagonal, starts to rotate and becomes fully horizontal for the larger sizes of the geometric imperfection (e.g., a left imperfection size of 0.5 mm) (Fig. 3). This indicates the increased relative importance of the properties of the visco-hyperelastic beam in determining the switching behavior of the bi-beam as compared to those of the hyperelastic beam (Fig. 3). In addition to the effects of the imperfection size, it is important to note that higher values of τ1 minimize the size of the switching zone, because, at lower strain rates, the incidence of right buckling for high τ1 values is high.
The switching is from left (at lower strain rate) to right (at higher strain rates). These regions are determined as the intersections of the regions corresponding to right buckling at higher strain rate with the regions defining left buckling for lower strain rates.
A combined experimental and computational study of imperfect bi-beams
To experimentally evaluate our predictions, we fabricated four bi-beams from a hyperelastic silicone rubber and a 3D (three-dimensional)–printed visco-hyperelastic rubber and introduced geometric imperfections to predispose the bi-beams to left buckling (Fig. 4A). The dimensions of the purposefully introduced imperfections were measured to be <5% of the width of the bi-beams (visible in Fig. 4A). For low rates of the applied deformation (e.g., 100 mm/min, strain rate = 1.7 1/s), all four specimens reliably buckled to the left (Fig. 4A). For a higher rate of the applied deformation (e.g., 500 mm/min, strain rate = 8.3 1/s), all four specimens reliably buckled to the right (Fig. 4A). The experiments were repeated more than 10 times, and the buckling direction was consistent between all experiments and all four beams. The effects of the applied strain rate are demonstrated in movie S1. The measured stiffness values of the bi-beams clearly increased with the applied strain rate (Fig. 4A).
(A) Stress-strain curves of actual bi-beams compressed at different loading rates. The warm colors refer to left buckling, while the cold colors stand for right buckling. The bi-beams were predisposed to left buckling. Switching in the direction of buckling (from left to right) happens by increasing the loading speed (from 100 to 500 mm/min). (B) Comparison between the buckling direction and the stress-strain curves predicted by our computational models and the experimental observations of the same quantities for two different aspect ratios of the bi-beams (i.e., 6 and 8) and for bi-beams with a variation in the thickness of their hyperelastic beam. The lengths of the beams with aspect ratios of 8 and 6 (short bi-beam) were, respectively, 120 and 90 mm (beam thickness = 21 mm for all types of beams). The length of the thick bi-beams was 120 mm, while the width of the hyperelastic beam was 12.5 mm, and the width of the visco-hyperelastic beam was 7.5 mm. Two types of computational models were used with either single-term or seven-term Prony series representing the viscous behavior. Photo credit: Shahram Janbaz, Delft University of Technology.
To assess the validity of our computational predictions, we compared the buckling directions predicted by our computational models, where single-term Prony series were used for capturing the viscoelastic behavior of the bi-beams, with those predicted by more complex computational models (i.e., with seven-term Prony series) and experimental observations (Fig. 4B). For all considered design variations (i.e., different aspect ratios and a variation in the geometrical design of bi-beams), the strain rate–dependent buckling directions predicted by both types of computational models agreed with each other and with experimental observations (Fig. 4B). Moreover, the stress-strain curves predicted by the different types of computational models largely overlapped and were in reasonable agreement with the experimentally measured stress-strain curves (Fig. 4B). These observations confirm that the results presented in Figs. 1 to 3, which were predicted using a single-term Prony series, are predictive of the actual buckling direction of the bi-beams. Despite the general agreement between the computational and experimental results, there were some differences between the stress-strain curves predicted by our numerical models and those obtained experimentally. It is important to note that the estimated and measured stress-strain curves and their corresponding stiffness values highly agreed with each other at the beginning of the loading (i.e., for small strains) (Fig. 4B). For large strains, however, the experiments consistently showed a softer response than what was predicted by our computational models (Fig. 4B). This softer response could be attributed to imperfections, such as air bubbles and imperfect curing and some omissions in the computational models, such as the effects of the adhesive used for attaching the two beams on the mechanical response of the materials from which the beams are made.
Application of bi-beams for the design of strain rate–dependent metamaterials
The reliable and highly predictable switching of the buckling direction opens up many opportunities for the design of strain rate–dependent mechanical metamaterials whose structural elements are bi-beams. We designed and created two such types of cellular structures with the following properties: (i) strain rate–dependent switching between auxetic and conventional (i.e., nonauxetic) behaviors and (ii) negative viscoelasticity. We also explored the possibility of applying bi-beams for the design of soft mechanisms.
To achieve strain rate–dependent switching between a positive (i.e., conventional behavior) and a negative (i.e., auxetic behavior) value of the Poisson’s ratio, we connected two parallel but mirrored bi-beams (with purposefully introduced geometric imperfections) to each other through two highly stiff connectors to create the basic unit cell Fig. 5A, i). The unit cell was then repeated in different directions to create a lattice structure (Fig. 5A, iii). We created nonlinear finite-element models of the lattice structures with three different sets of material properties (
(A) (i) A near-orthogonal unit cell made from predisposed bi-beams was used to form a lattice structure whose Poisson’s ratio is highly dependent on the applied strain rate and could switch between an auxetic (i.e., negative Poisson’s ratio) and a conventional (i.e., a positive Poisson’s ratio) behavior. (ii) Three different sets of material properties were chosen to evaluate the evolution of the Poisson’s ratio with strain under different rates of the applied strain. (iii) Switching from an auxetic to a conventional behavior. Some unit cells from the middle part of the cellular lattice have been magnified for better visualization of the strain distributions. (iv) Experimental observation of strain rate–based switching from an auxetic to a conventional behavior. (B) A circular arrangement of bi-beams was used to transform axial compression to either clockwise or counterclockwise rotation, depending on the applied strain rate. Moreover, the stiffness of the rotational mechanism is highly dependent on the applied strain rate. Photo credit: Shahram Janbaz, Delft University of Technology.
A circular arrangement of the bi-beams was then used to transform axial (compressive) loads to either clockwise or counterclockwise rotation, depending on the applied strain rate (Fig. 5B and movie S3). Together with the direction of rotation, the stiffness of the mechanism changed as well (Fig. 5B), meaning that the characteristic curve of such a mechanism would be also dependent on the applied strain rate. Such possibilities in transforming axial deformation to rotation, in switching between clockwise and counterclockwise rotation depending on the applied strain rate, and in adjusting the characteristic curve of a soft mechanism could be instrumental in the development of dedicated powertrains for soft machines.
A third application of bi-beams is in creating mechanical metamaterials with negative viscoelastic behavior. We defined negative viscoelasticity as the property of a material that shows lower instantaneous stiffness under higher strain rates. Solid materials usually only show positive viscoelasticity, meaning that their stiffness increases with the strain rate. In some ways, negative viscoelasticity is similar to “thixotropy” (37) in non-Newtonian fluids where the viscosity and, thus, the stiffness of the viscoelastic fluid decrease upon the application of a high–strain rate load. Negative elasticity is, however, distinct from thixotropy, as it is defined for permanently solid materials and is not associated with large changes in the viscosity of the material. To achieve negative viscoelasticity, we changed the structural element of our cellular structure from one single bi-beam to two bi-beams that were positioned with a small distance, a, between them (Fig. 6A). The applied strain rate determines the buckling direction of each of the constituting bi-beams, creating three types of contact situations between them (Fig. 6A). The direction of the bi-beams and the material properties can be selected such that for very small strain rates (e.g., 10−4 1/s), both constituting bi-beams buckle toward each other, thereby giving rise to rapidly increasing contact forces that increase the stiffness of the cellular structure (Fig. 6A). For intermediate strain rates (e.g., 1 1/s), both bi-beams initially buckle toward each other followed by a deformation to the same direction after a second instability (Fig. 6A). For high strain rates (e.g., 103 1/s), the bi-beams buckle toward opposite directions and away from each other, avoiding any contacts and the associated increase in the stiffness (Fig. 6A). The stress-strain curves and instantaneous stiffness values calculated for all the three ranges of the applied strain rates showed that, for large enough strains (e.g., >0.15), the stiffness associated with a high strain rate is many times smaller than that of a low–strain rate case (Fig. 6A). When such double bi-beams were arranged in a cellular structure using highly stiff connecting elements, the same type of negative viscoelasticity was observed (Fig. 4B and movie S4). We also performed experiments with double bi-beams to explore their response to different rates of the applied strain (Fig. 6C). For a deformation rate of 500 mm/s, the bi-beams buckled away from each other, similar to what was predicted by the computational models (Fig. 6C). For a deformation rate of 100 mm/s, the expected stiffening behavior takes place where the bi-beams buckle toward each other (Fig. 6C). With an increasing level of the applied strain, a second instability mode is activated, where one of the bi-beams pushes the other bi-beam to straighten and buckle to the other side, ultimately resulting in the separation of the bi-beams (Fig. 6C). For high enough deformations (e.g., >7 mm), the instantaneous stiffness and the maximum force were much lower for a deformation rate of 500 mm/s as compared to 100 mm/s (Fig. 6C). Of course, the negative viscoelasticity described here is a post-buckling phenomenon. For infinitesimal strains, the stiffness of a double bi-beam is slightly higher in the case of higher strain rates, as a higher strain rate initially increases the stiffness of the visco-hyperelastic beam present in the construction of each of the bi-beams (Fig. 6C). Because such post-buckling behaviors are highly dependent on the contact between the bi-beams, geometric parameters such as the distance between both bi-beams, a, could be used to adjust the observed properties and their strain rate dependency to befit the specific application at hand.
(A) At lower strain rates, the parallel bi-beams buckle toward each other, resulting in contact forces that increase the stiffness of the cellular structure, once the threshold of bi-beam instability has been reached. At higher strain rates, the bi-beams buckle away from each other. The stress-strain curves and the evolution of the elastic modulus with strain are compared for different loading rates. (B) The arrangement of the unit cells in a lattice structure and the response of the lattice structure to different rates of the applied strain. (C) Experiments demonstrating the same behavior as predicted by our computational models. Initially (i.e., before the onset of instability), the unit cells exhibit a stiffer behavior under higher strain rates. This is reversed after the onset of instability beyond which point the stiffness corresponding to a higher strain rate is lower than the one corresponding to a lower strain rate. Photo credit: Shahram Janbaz, Delft University of Technology.
Practical consideration: Imperfect loading conditions and scalability
In the last step of the study, we examined a few aspects that are expected to be of importance for the practical applications of bi-beams. First, the real-world applications of bi-beams are likely to be associated with imperfections in the way the loads are applied even when the design of the metamaterials is based on fully orthogonal loading conditions. Moreover, the design of strain rate–dependent mechanical metamaterials or the devices that incorporate them may necessitate loading conditions that are not fully axial. We, therefore, examined the effects of two different types of deviations from the previously discussed loading conditions on the strain rate–dependent buckling behavior of bi-beams. In the first experiment, we applied 20-mm offsets to the left and to the right (Fig. 7A and movie S5). Compressive loads with different rates (corresponding to the same approximate strain rates as applied to the bi-beams tested in the previous experiments) were then applied to the bi-beams. The bi-beams continued to exhibit the desired strain rate–dependent behavior despite the applied offset (Fig. 7A and movie S5). In the second experiment, we applied ≈45° rotation to the bi-beams before applying the compressive loads (Fig. 7A and movie S5). Once more, the bi-beams continued to exhibit the desired strain rate–dependent behavior despite the applied rotation (Fig. 7A and movie S5). These experiments show the robustness of the strain rate–dependent buckling behavior of the bi-beams that remain unchanged despite such large deviations from the axial loading conditions and underscore the potential of the bi-beam concept for practical applications.
(A) The bi-beams (thick bi-beam, aspect ratio of 6, similar to those used in Fig. 4B) exhibited the same strain rate–dependent buckling behavior despite imperfect loading conditions, including 20-mm offsets to the left and right and ≈45° rotation in the bi-beams. (B) Two miniaturized bi-beams (left: length = 22 mm, aspect ratio = 3.7, and thickness = 10 mm; right: length = 32 mm, aspect ratio = 5.3, and thickness = 10 mm) were manufactured and tested under low and high loading speeds that were scaled with the length of the bi-beams so as to maintain the same approximate strain rates as was used in the previous experiments (i.e., Figs. 4 to 6). Scaling down the length of the bi-beams by up to ≈5.3 times did not change the strain rate–dependent buckling behavior of the bi-beams. Photo credit: Shahram Janbaz, Delft University of Technology.
The other important aspect is the size of the bi-beams and how that may affect their buckling behavior. To examine the effects of specimen size, we manufactured two smaller versions of bi-beams (i.e., 5.5 and 3.8 times smaller) with different aspect ratios (i.e., 3.7 and 5.3) (Fig. 7B and movie S6). Both smaller versions of the bi-beams exhibited the expected strain rate dependency in their buckling behavior (Fig. 7B and movie S6). This confirms the scalability of the presented approach at least in the range of the sizes considered here. However, the concept of bi-beams and their strain rate–dependent behavior could be of interest for a wide range of practical applications. While structural application in which the length scales of the bi-beams is larger and in the range of those demonstrated here are relatively straightforward to realize, many other applications particularly in microelectromechanical systems, (steerable) medical instruments, and microrobots will require manufacturing of (a large number of) bi-beams at much smaller length scales. While the mechanical principles that are used in the design of bi-beams and the demonstrated strain rate–dependent behavior are highly scalable, the realization of such bi-beam constructs at smaller scales could pose some challenges. In terms of manufacturing, the microscale fabrication techniques required for such applications are already available. In principle, (multi-material) additive manufacturing techniques such as material jetting (9) and direct laser writing (38) can be used to manufacture similar constructs at the microscale and submicrometer scale, respectively. The main challenge, however, is in finding photopolymerizable materials that can be processed using the currently available techniques and exhibit the differential viscoelastic properties required for creating bi-beams. Moreover, the adhesion between both types of polymers (i.e., visco-hyperelastic and hyperelastic) should be very strong so as to support the transfer of forces between them. These are very specific requirements that, according to our preliminary screening of the currently available materials, are difficult to simultaneously satisfy. Alternatively, the differential strain-rate response of the beams of a bi-beam construct could be created in other ways of which inertia forces are a prime example. The use of inertia forces is much easier to realize in practice, and such types of bi-beams would be highly scalable. However, there is a limit (i.e., the difference between the range of the density of processable materials and voids) that determines the maximum amount of the differential behavior achievable using inertia forces. The application of additive manufacturing techniques would also mean that it is, at least in principle, relatively easy to fabricate 3D version of bi-beam arrays. This is likely to constitute a worthy avenue for further research, as the mechanical behavior of such bi-beam arrays is expected to be quite rich.
CONCLUSIONS
In summary, we demonstrated how strain rate–dependent mechanical metamaterials can be realized using bi-beams, which consist of two laterally attached beams of which one is hyperelastic and the other is visco-hyperelastic. Our computational models show the importance of geometrical design including the shape of the bi-beams (e.g., cut-out patterns), aspect ratio, and purposefully introduced imperfections in the reliability of the strain rate–dependent control of their buckling direction. We also showed how cellular structures and soft mechanisms exhibiting novel strain rate–dependent behaviors could be realized using (double) bi-beams as their structural elements. The strain rate–dependent switching between auxetic and conventional behaviors, negative viscoelasticity, and strain rate–dependent transformation of axial deformation to (counter) clockwise rotation are expected to find potential applications in various areas, including soft robotics, the development of wearable medical devices, and flexible electronics.
MATERIALS AND METHODS
Nonlinear buckling analysis was performed using an implicit nonlinear solver (Abaqus 2016. Simulia, USA) to study the effects of material parameters on the buckling behavior of bi-beams. Quadrilateral elements (type CPE4H) were used for the discretization of the bi-beam geometries (length, 40 mm; width, 5 mm; thickness, 25 mm). On the basis of the results of a mesh convergence study, at least 12 elements were seeded through the width of the bi-beams. The aspect ratio of the beams was between 4 and 16 (varied in steps of 4). The bi-beams were clamped at both their ends, and a reference point was created to apply the loads to the top side of the bi-beams. A series of simulations were performed for applied strain rates ranging between
Given that
The incompressibility condition and uniaxial loading conditions (i.e., σ22 = σ33 = 0) could then be used to determine the arbitrary pressures pV and pH. The stress in the direction 1 can be determined by combining Eqs. 5 and 6 as
Assuming a single-term Prony series, the instantaneous elastic modulus in the direction 1 is given by
For the small axial deformations preceding buckling, we have λ ≅ 1. The elastic moduli of the visco-hyperelastic and hyperelastic beams could then be written as
To study the buckling behavior of the bi-beams, each of the beams (i.e., visco-hyperelastic and hyperelastic beams) can be modeled as a beam on an elastic foundation with each beam acting as the elastic foundation of the other beam. The buckling load for a clamped-clamped beam on elastic foundation, Pcr, is higher than the Euler’s buckling load for an equivalent unsupported beam, Pe, and can be obtained from the solution to the following implicit equation (34)
The condition for buckling to the right side (i.e., the side of the hyperelastic beam) can then be simplified from εcr,H> εcr,V to
Because the beams are geometrically identical and given that
Combining Eqs. 9 and 15 gives the condition for buckling to the right as
Using Eq. 2, Eq. 16 could be rewritten as
For extremely high strain rates (i.e., t/τ1 → 0), the condition for buckling to the right can be simplified as
One of the limitations of the beams on elastic foundation theory as applied here is that it neglects the global buckling of the beams. When both beams buckle simultaneously, the modulus of the elastic foundation, k, is different from the elastic modulus of the other beam. One way to estimate the elastic modulus of the elastic foundation would be to impose the requirement that the critical load predicted by the beams on the elastic foundation theory equals that of the Euler’s theory for a bi-beam made from two identical beams (i.e., with identical geometries and material properties). Such an approach would suggest that the modulus of the elastic foundation is a function of not only the elastic modulus of the beams but also their geometry, particularly their slenderness. This additional geometry-dependent coefficient will, however, not change the main requirement for buckling to the right (i.e., EV > EH).
In our computational analyses, we varied the ratio of the Neo-Hookean coefficient of the visco-hyperelastic layer (
To evaluate the effects of geometric imperfections on the buckling direction of bi-beams, we used an eigenvalue buckling analysis to predict the first mode of buckling for a single-material beam. Having added positive and negative values of the buckling mode as a geometric imperfection to the geometry of our bi-beams, we performed nonlinear buckling analyses. In the context of these analyses, the imperfection size refers to the maximum size of the offset from the centerline of the bi-beam geometry presented (dimension, mm).
To manufacture the bi-beams used in the experiments, visco-hyperelastic beams were separately fabricated using an advanced multi-jet printer (Stratasys Objet350 Connex3, USA) by using a ultraviolet-curable polymer (Agilus, Stratasys, USA). The surface of the visco-hyperelastic beams was then fully covered with a thin layer of a rubber silicone adhesive (Sil-Poxy, Smooth-On) to ensure good adhesion between Agilus layer and silicone rubber. We used a silicone rubber (Dublisil 30, Dreve Dentamid GmbH, Germany) to mold the hyperelastic part of the bi-beams. To fabricate the other half of the bi-beams, we printed PLA (poly lactic acid) molds using a fused deposition modeling (FDM) 3D printer (Ultimaker 3, the Netherlands). We then positioned the visco-hyperelastic beam in the PLA molds and used the same rubber silicone adhesive to ensure that they remained attached to the walls of the molds during the curing process of the silicone rubber used for creating the hyperelastic beam. Subsequently, we filled the empty parts of the molds with silicone rubber. The bi-beams were kept in the molds for 8 hours at room temperature after which they were removed from the molds. We characterized the Neo-Hookean parameter, C10H, of the hyperelastic rubber based on the compression tests of rubber disks (D = 28.5 mm and h = 12.5 mm) under quasi-static conditions, while the Neo-Hookean parameter,
To ensure well-defined clamping conditions during compression tests, the bi-beams were placed in a compression test setup, which was 3D-printed from PLA using the same FDM printer. The speed of the crosshead of the testing machine was allowed to stabilize before compressing the bi-beams by introducing a gap (10 mm) between the moving crosshead and the test setup. Before placing the bi-beams in the test setup, their end parts were cut in such a way that a small geometric imperfection toward the hyperelastic side was materialized once they were clamped in the test setup. The lattice structures and parallel bi-beams were fabricated by printing the additional parts from PLA using the same FDM 3D printer (Figs. 5 and 6). The bi-beams and the cellular designs were compressed using the same Lloyd test bench equipped with a 100-N load cell. A digital camera (Sony A7R with a Sony FE 90 mm f/2.8 macro OSS lens) was used to record the deformation patterns of the bi-beams and the cellular designs during the compression tests.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/25/eaba0616/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).