Next Article in Journal
Effect of Aggregate Size on Strength Characteristics of High Strength Lightweight Concrete
Next Article in Special Issue
Hybridized Love Waves in a Guiding Layer Supporting an Array of Plates with Decorative Endings
Previous Article in Journal
Assessment of the Effect of Laser Welding on the Properties and Structure of TMCP Steel Butt Joints
Previous Article in Special Issue
Cloaking In-Plane Elastic Waves with Swiss Rolls
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Planar Mechanical Metamaterials with Embedded Permanent Magnets

by
Viacheslav Slesarenko
1,2
1
Freiburg Institute for Advanced Studies, 79104 Freiburg im Breisgau, Germany
2
Cluster of Excellence livMatS @ FIT—Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, 79110 Freiburg im Breisgau, Germany
Materials 2020, 13(6), 1313; https://doi.org/10.3390/ma13061313
Submission received: 28 January 2020 / Revised: 5 March 2020 / Accepted: 11 March 2020 / Published: 13 March 2020
(This article belongs to the Special Issue Advances in Mechanical Metamaterials)

Abstract

:
The design space of mechanical metamaterials can be drastically enriched by the employment of non-mechanical interactions between unit cells. Here, the mechanical behavior of planar metamaterials consisting of rotating squares is controlled through the periodic embedment of modified elementary cells with attractive and repulsive configurations of the magnets. The proposed design of mechanical metamaterials produced by three-dimensional printing enables the efficient and quick reprogramming of their mechanical properties through the insertion of the magnets into various locations within the metamaterial. Experimental and numerical studies reveal that under equibiaxial compression various mechanical characteristics, such as buckling strain and post-buckling stiffness, can be finely tuned through the rational placement of the magnets. Moreover, this strategy is shown to be efficient in introducing bistability into the metamaterial with an initially single equilibrium state.

Graphical Abstract

1. Introduction

Mechanical metamaterials possess unique properties and demonstrate unconventional behavior owing to their intricate internal architecture [1]. In many cases, the mechanical properties of the constituents do not play a significant role in defining the overall performance of the metamaterial. Extreme strength can be achieved even in metamaterials made from relatively weak base components due to the rational design of their interior [2]. It resembles the design principles of natural materials and composites that demonstrate extremely high performance despite relatively weak constituents [3,4,5,6]. At the same time, natural materials, as well as some classical mechanical metamaterials, were initially treated as structural; therefore, their performance was usually evaluated based on classical mechanical criteria such as strength, stiffness, or toughness [2,7,8,9,10]. However, similar design principles can be employed for the creation of functional mechanical metamaterials that are responsive [11], programmable [12], adaptive [13], and smart [14].
Various mechanical phenomena have been harnessed to design mechanical metamaterials with intricate functionalities. For instance, it was shown that mechanical buckling could be exploited to amplify the applied deformation and trigger involved reconfiguration of the internal architecture [13,15,16,17,18]. The capacity to program the deterministic deformation sequence through the embedment of the viscoelastic constituents in mechanical metamaterials has been demonstrated both theoretically and experimentally [19]. A similar idea of combining buckling and viscoelasticity inspired the novel design of programmable mechanical metamaterials with memory [20]. Rational selection of the internal architecture enables the development of materials that can convert simple compression into a twisting motion [21,22]. Another emerging design strategy explores the capacity of the internal architecture to transfer local deformation through the metamaterial. In this case of so-called action-at-a-distance metamaterials, the internal architecture may define the involved geometrical changes observed on the surface of the metamaterial in response to applied mechanical stimuli [23,24,25,26].
Nevertheless, the above-mentioned mechanical metamaterials, despite very involved behavior, have one fundamental weakness—their mechanical properties usually cannot be altered after fabrication. To overcome this restriction, the employment of stimuli-responsive constituents has been proposed to provide finer control over the performance of the metamaterials. Indeed, the synergetic combination of a base material capable of reacting to external stimuli by changing its properties or shape, and the rational architecture of the metamaterial to transform these changes into the global deformation/reconfiguration, provides a powerful means to design new adaptive metamaterials [27,28]. For example, the employment of shape memory polymers enables control over the material’s performance by changing the temperature [29]. Moreover, current advances in additive manufacturing techniques promise very quick and simple production and programming of such metamaterials [28]. Another type of stimuli-responsive metamaterials utilizes the ability of the base materials to swell, causing a reconfiguration of the internal architecture, thereby altering the effective properties [30]. At the same time, these two actuation principles (temperature and swelling) are characterized by very high inertia, and therefore it is almost impossible to achieve a high frequency of the actuation.
An alternative external field/stimuli that can be employed to trigger deformation within a metamaterial is the electromagnetic field [11]. Embedding dielectric or piezoelectric elements and magnetorheological elastomers seems to be a feasible strategy to enable control over the performance of these mechanical metamaterials. Application of an external magnetic or electric field can almost immediately trigger the deformation and reconfiguration of the internal architecture, which directly affects the properties and overall performance of the metamaterial [31,32]. However, the current generation of dielectric elastomers/magnetorheological elastomers demonstrates a rather weak coupling between the electromagnetic field and the mechanical deformation [33,34]. Even for the above-mentioned magnetoactive metamaterial [31], a magnetic field of 0.24 T is required to trigger buckling in samples with dimensions not exceeding a couple of centimeters. Such field strength is simply unachievable on a larger scale, making this approach unfeasible for practical applications outside the laboratory. However, there is no doubt that additional magnetic interaction provides powerful means to enhance the functionality of mechanical metamaterials.
Fortunately, there is another strategy to harness the magnetic interaction by employing permanent magnets embedded into the architecture of the metamaterial. This approach was shown to be useful for the origami folding structures, where the employment of permanent magnets enabled the locking of certain states [35]. Other metamaterials employ the interaction between permanent magnets for elastic energy harvesting [36,37]. From a mechanical standpoint, permanent magnets provide a convenient and straightforward way to adjust the stiffness of the lattice structure. Moreover, the insertion of permanent magnets into the structure of a metamaterial can enhance the impact resistance [38]. In practice, mechanical metamaterials consisting of re-entrant unit cells with embedded magnets are able to exhibit negative stiffness and auxetic behavior simultaneously [39]. Moreover, the incremental stiffness of the system can be adjusted by changing the orientation of the magnets [40].
Even in a relatively simple system consisting of magnetic particles connected by elastic springs, advanced phenomena (e.g., multistability) have been observed [41]. More sophisticated systems can be designed by integrating the magnets into the auxetic metamaterials [14,42]. The classical auxetic metamaterial with Poisson’s ratio of −1, which supports the embedment of permanent magnets, consists of rotating squares [43,44] connected by soft hinges. The mechanical behavior of similar metamaterials with rotating units was thoughtfully investigated in multiple publications (e.g., [45,46,47]) in which analytical models that take into account the deformation of the hinges and units were presented and discussed. The embedment of magnetic particles into metamaterials requires introducing additional components into the models to account for a magnetic interaction between elements. For example, the bead-spring model, which takes into account the mechanical force between beads, the Lennard-Jones interaction, and the magnetic dipole–dipole interaction, was shown to be capable of describing the behavior of systems with positive and negative Poisson’s ratios on the nanoscale [48,49]. On a larger scale, similar systems with permanent magnets embedded into the interior of auxetic metamaterials with rotating units were proposed and discussed in [14]. In particular, it was shown that auxeticity and stiffness of mechanical metamaterials with permanent magnets could be tuned by rearranging the magnets and adjusting the external magnetic field. The proposed concepts [14] provide a solid foundation for further theoretical and experimental studies of magneto-mechanical metamaterials.
The present work focuses on mechanical metamaterials with rotating squares in which the behavior is programmed through the modification of the energy landscape by periodic embedment of permanent magnets in various configurations. Here, attention is devoted to the macroscopic properties of the metamaterials (e.g., stiffness and stability), rather than to local structural changes (e.g., rotation angles). The designed planar metamaterials are manufactured using three-dimensional printing and subjected to biaxial compression to obtain continuous force–displacement curves. The selected design enables efficient reprogramming of the mechanical properties by inserting the magnets into specific locations. Numerical simulations support the experimental study and provide the necessary explanation of the observed behavior.

2. Materials and Methods

A classical design consisting of rotating squares connected by soft hinges/connectors was adopted in the current study [14,43,44,50]. Figure 1a shows an elementary cell of the metamaterial that consists of four rotating squares. A continuous equibiaxial compression leads to the rotation of stiff squares until the closed state is reached (Figure 1b). Since the goal of this research was to harness magnetic interaction as a tool to program the performance of the metamaterial, some elementary cells were modified through the embedment of permanent magnets at the centers of rotating squares. Three configurations of the elementary cell were considered in this study. Neutral cell did not contain any embedded magnets (Figure 1a), repulsive cell contained four magnets with co-directional out-of-plane magnetic moments (Figure 1c), and attractive cell contained four magnets with counter-directional out-of-plane magnetic moments (Figure 1d). Among various auxetic models [51], this specific design was selected mainly due to the fabrication method and the experimental limitations.
Numerical simulations of magneto-mechanical metamaterial were performed using FE software COMSOL 5.4a. The plane strain formulation was employed, and out-of-plane thickness of 5 mm was selected to match the dimensions of experimentally tested specimens. The mechanical behavior of the metamaterial with the periodic arrangement of modified elementary cells, which is characterized by the repeated unit cell shown in Figure 1e, was studied. Only four central squares (shown by blue color) were modified by permanent magnets. Selected 4 × 4 unit cell was subjected to equibiaxial compression through the application of the displacement ε on the boundaries (Figure 1e). The corresponding reaction force f was measured through the integration over the boundary of the unit cell. To take into account the periodicity, the equibiaxial compression was implemented using pointwise boundary conditions on the edges of the unit cells as follows:
{ u | l e f t u | r i g h t = ε a v | l e f t v | r i g h t = 0   u | b o t t o m u | t o p = 0   v | b o t t o m v | t o p = ε a   ,
where u   and   v are horizontal and vertical displacements, respectively, ε —applied compressive strain, and a —width and height of the unit cell.
Rotating squares were treated as rigid non-deformable bodies, while soft connecting hinges were considered as neo-Hookean solids with shear modulus μ = 1.1 MPa and Poisson’s ratio of ν = 0.46 . While FE analysis is capable to take into account the magnetic interaction between permanent magnets, the out-of-plane orientation of the magnetic moments requires 3D formulation and extra meshing of empty space around the metamaterial. An assumption of the dipole–dipole interaction between embedded magnets helps to keep 2D formulation, avoid unnecessary meshing and decrease computation time. Therefore, the magnetic interaction between embedded magnets within one elementary cell was treated as the dipole–dipole interaction [52]. For two identical magnets with out-of-plane orientations of the magnetic moment, the force exerted on the second magnet by the first can be calculated as
F 2 = s 1 s 2 3 μ 0 m 2 4 π | r 12 | 5 r 12 ,
where μ 0 = 4 π × 10 7 H/m is the permeability of the vacuum, r 12 is the vector connecting centers of the first and second magnets, m is the value of the magnetic moment, and s 1 ,   s 2 = ± 1 define the orientation of the magnetic moment. For convenience s i = + 1 if the magnetic moment is codirectional with z-axis and s i = 1 otherwise. According to this definition, two magnets with codirectional out-of-plane magnetic moments experience repulsive forces, while counter-directional orientation corresponds to the attractive configuration. The resulted force exerted on each magnet in modified elementary cell can be obtained as the sum of all pairwise interactions, and it is equal to
F j = 3 μ 0 m 2 4 π i s i s j r i j | r i j | 5   for   i j .
Therefore, the magnetic energy of the elementary cell is equal to
U = μ 0 m 2 4 π i , j s i s j | r i j | 3   for   i j ,
where the sum is taken for all pairwise combinations of magnets (for each pair only once). Due to the large distance between modified elementary cells in considered periodic metamaterial, the magnetic interaction between magnets from different cells was neglected. The magnetic forces (Equation (3)) and the magnetic energy (Equation (4)) were manually integrated into COMSOL.
For the experimental study, mechanical metamaterials were produced with the help of multi-material 3D printer Objet Connex 260 (Stratasys, Eden Prairie, MN, USA). Rotating squares were printed by rigid polymeric VeroWhite (Stratasys, Eden Prairie, MN, USA) material with Young’s modulus around 2 GPa, while soft hinges were produced using soft digital material. A slight modification of the outer boundaries of rotating squares was made in order to ensure the contact between neighbor squares after 45° rotation and to accommodate the deformation of soft hinges without mutual contact. The interior of rotating squares was left hollow to facilitate quick insertion/removal of the magnets after printing. The cubic 5 mm × 5 mm × 5 mm Nd-magnets of N42 grade were used to modify 3D printed metamaterial. The samples were subjected to equibiaxial compression in transparent fixtures preventing out-of-plane deformation, which was attached to the universal testing machine Shimadzu EZ-LX (Shimadzu Corporation, Kyoto, Japan). The tests were performed at a relatively low rate of 10 mm/min in order to mitigate the contribution of the viscoelasticity [53]. A CCD camera was used to capture the deformation sequence. The displacement D b and the force F b acting on the boundaries were calculated as D d = D m 2 and F d = F m 2 , where D m and F m are the displacement and the force measured by the testing machine.

3. Results and Discussion

3.1. Validation of the Dipole–Dipole Interaction

Auxiliary experimental measurements were performed to explore the accuracy of used dipole–dipole simplification for describing the interaction between cubic magnets with out-of-plane magnetic moments. Two magnets were fixed in 3D printed fixtures attached to the opposite jigs of the universal testing machine. The magnets mounted precisely above each other were slowly moved towards each other, while the displacement of the jigs and the exerting force were measured by the testing machine. At least five different magnets in two configurations (codirectional/repulsive and counter-directional/attractive) were tested. Figure 2 shows the obtained force-distance curves for two tests. As expected, an increase in the distance between magnets led to a drastic decrease in the interaction force both for repulsive and attractive configurations. Solid symbols in Figure 2 correspond to the force-distance relation calculated according to Equation (2) for selected values of the magnetic moment m . As one can see, the estimation obtained for m = 0.105   Am 2 provided the best fit of the force-distance curves. Note, that this value was lower than the value that could be estimated based on the dimensions and the grade of the magnets as m = B r V / μ o 0.127   Am 2 , where V is the volume of the magnet and B r is the remanent magnetization in T. This discrepancy is probably caused by an inaccuracy of the dipole–dipole model for considered configuration, however, a decrease in m value from estimated m = 0.127   Am 2 to m = 0.105   Am 2 enables a very good fit of experimental curves, therefore the dipole–dipole interaction model with m = 0.105   Am 2 was adopted in the following numerical section.

3.2. Numerical Results for Periodic Metamaterials

Figure 3 shows the force–strain curves obtained numerically for periodic metamaterials with various configurations of elementary cells (neutral, repulsive, and attractive). Additional dashed curves are shown for the metamaterial with weaker embedded magnets ( m = 0.075   Am 2 ) . The metamaterial in neutral configuration underwent buckling upon reaching critical force f c r n that was accompanied by rotation of rigid squares in alternating directions. As a result, an initial compressive state in soft hinges was replaced by bending during post-buckling deformation, causing the drastic drop in the slope of the force–strain curve, as shown in Figure 3. The rotation continued with an increase in applied strain ε until all squares rotated by 45° degrees, and the metamaterial completed its transition to the closed state (Figure 1b), when ε = ε c l . Mechanical metamaterials with embedded magnets demonstrated a similar bucking behavior and underwent softening after buckling at the early stage of the deformation. At the same time, the critical force required to trigger the buckling depended on the configuration and strength of the magnets. Indeed, embedment of the magnets in the repulsive configuration led to an increase in the critical force, while an opposite effect was observed for the system with attractive elementary cells. Moreover, it is clear that the deviation from f c r n increased with an increase in the value of the magnetic moment.
Figure 3 reveals that the slope of the force–strain curve in the post-buckling regime was also highly affected by the configuration of embedded magnets. Through the comparison with neutral configuration, it is easy to notice that the contribution of the magnetic interaction in repulsive configuration led to a stiffer response of the metamaterial. Surprisingly, attractive configuration resulted in the negative stiffness after buckling, and therefore a lower force needed to be applied in order to continue equibiaxial compression. Numerical simulations revealed that for repulsive configuration of magnets and large enough m , the decaying force could reach zero at some value of strain ε s t prior to the achievement of the strain ε c l . In this case, a rapid snap-through transition to the closed state occurred, and the metamaterial remained in the closed state even after the removal of external load.
A closer examination of the energy of the unit cell provided a clear explanation of the observed behavior. The total energy of the system E t o t a l can be evaluated as the sum of magnetic energy ( E m a g ) and elastic energy stored in soft deformable hinges ( E e l ) . Elastic energy E e l was calculated by integrating the strain energy density function over soft hinges. Figure 4a shows the dependencies of the total energy of the system on applied strain ε , while Figure 4b shows both components of the total energy for repulsive and attractive configurations with m = 0.105   Am 2 . For neutral configuration, the total energy increased monotonically due to an increase in E e l during deformation. The addition of repulsive configuration led to an increase in the total energy, since E m a g , calculated according to Equation (4), had a positive sign. During the deformation, E m a g increased monotonically, contributing to the growth of E t o t . Therefore, E t o t as a function of the applied deformation had only one minimum corresponding to the open state of the metamaterial. An opposite observation could be made for the attractive configuration. The magnetic energy E m a g had a negative sign, and the biaxial compression led to a decrease in E m a g due to the decrease in the distance between the magnets in elementary cell, as shown in Figure 4b. As a result, the total energy of the system E t o t a l = E m a g + E e l became concave, and for large enough m , had a maximum for the strain ε s t < ε c l , marked by the blue circle in Figure 4. Therefore, under these conditions, E t o t a l ( ε ) had two local minima, separated by a potential barrier, that corresponded to the open and closed states of the metamaterial. As soon as the applied strain ε reached the critical value ε s t during equibiaxial compression, the metamaterial snapped into the second equilibrium state. Therefore, the rational embedment of interacting permanent magnets can be exploited to modify the energy landscape of the metamaterials and to create systems with multiple stable equilibrium states.

3.3. Experimental Results

For the experimental study, the metamaterials consisting of 16 × 16 rotating squares with two various arrangements of the magnets were considered. Arrangement I closely resembled the design used in the numerical study with the positions of the magnets shown in Figure 5a by blue squares. Arrangement II had magnets only at the diagonal locations, as shown by white circles in Figure 5a. Figure 6a shows experimentally obtained force–displacement curves (two for each arrangement with repulsive and attractive configurations, plus neutral metamaterial). Note, that the metamaterials were compressed only up to 17 mm ( ε = 0.12 ) to prevent possible failure of the hinges due to high local strains and relatively weak 3D-printable material. Dashed lines represent the additional numerical simulations performed on the finite 16 × 16 metamaterials with arrangements I and II. As one can see, experimental curves very closely match the numerical predictions with respect to the buckling force and the slope of the curves in the post-buckling regime. Similar to the numerical results, the significant softening and even negative effective stiffness were observed experimentally. Figure 6a also revealed that the increase/decrease in the critical buckling force for repulsive/attractive configurations relative to the neutral configuration was observed irrespectively of the magnet arrangement. At the same time, the metamaterial with the arrangement II demonstrated a lower deviation of the buckling force in comparison with the arrangement I due to a lower number of the modified elementary cells. The contribution of the magnetic energy towards the total energy of the system is directly proportional to the density of modified elementary cells in the metamaterial. Therefore, the employment of stronger magnets and the embedment of a larger number of magnets are two feasible strategies in amplifying the effect of the magnetic interaction on the overall performance of mechanical metamaterial.
During the equibiaxial loading of the specimens with attractive magnets, the force–displacement curves demonstrated the negative slope in the post-buckling regime. For arrangement II, the force remained positive during biaxial compression, and the removal of external load led to the restoration of the open state (Video S1 in Supporting Materials). On the contrary, the experimental study confirmed that for arrangement I (Figure 6b), when the measured force dropped to almost zero, the snap-through transition to the closed state was observed (Figure 6c and Video S2 in Supporting Materials). The sample remained in the closed state even after the removal from the fixtures proving that this closed state is another stable equilibrium. It is clear that the value of ε s t depends on the arrangement of the magnets and their strength. Therefore, by adjusting the position or the density of the modified elementary cells in the metamaterials, it seems possible not only to create the second stable state, but also to adjust the critical strain ε s t at which the snap-through transition between states occurs.

4. Conclusions

Here, the strategy to program and control the mechanical behavior of planar mechanical metamaterials through the embedment of the interacting permanent magnets was studied. On the example of the metamaterial with rotating squares, the advantages of this strategy were shown both experimentally and numerically. The planar metamaterial produced by three-dimensional printing enables efficient and fast adjustment of the mechanical properties by rearranging the magnetic cells. By utilizing this design, the control over the buckling strain and the stiffness of the system in the post-buckling regime was demonstrated experimentally. It was also experimentally confirmed that the rational placement of the attractive magnets could be harnessed to introduce multistability into the metamaterial with an initially single equilibrium state. By generalizing the observations, it can be concluded that the underlying strategy enabling the enhanced control over the behavior of mechanical metamaterials is based on the thoughtful modification of the energy landscape by employing intricate interplay between the elastic and magnetic energies. More sophisticated auxetic models enabling larger changes in the internal structure (e.g., [54]) can be considered as promising candidates to host permanent magnets, but the feasibility of their fabrication should be taken into account. With the rapid development of additive manufacturing techniques, similar design principles can be potentially applied for three-dimensional auxetic systems to design advanced mechanical metamaterials with programmable behavior.

Supplementary Materials

The following are available online at https://www.mdpi.com/1996-1944/13/6/1313/s1. Video 1: Equibiaxial compression of mechanical metamaterial with attractive magnets in the arrangement II, Video 2: The snap-through phenomenon during equibiaxial compression of mechanical metamaterial with attractive magnets in the arrangement I.

Funding

The research leading to these results has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie-Sklodowska-Curie grant agreement No.754340. The research received support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC-2193/1—390951807. Gefördert durch die Deutsche Forschungsgemeinschaft (DFG) im Rahmen der Exzellenzstrategie des Bundes und der Länder—EXC-2193/1—390951807. VS acknowledges support by the state of Baden-Württemberg through bwHPC (NEMO, RV bw19K001). The article processing charge was funded by the German Research Foundation (DFG) and the University of Freiburg in the funding program Open Access Publishing.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Zadpoor, A.A. Mechanical meta-materials. Mater. Horiz. 2016, 3, 371–381. [Google Scholar] [CrossRef] [Green Version]
  2. Zheng, X.; Lee, H.; Weisgraber, T.H.; Shusteff, M.; DeOtte, J.; Duoss, E.B.; Kuntz, J.D.; Biener, M.M.; Ge, Q.; Jackson, J.A.; et al. Ultralight, ultrastiff mechanical metamaterials. Science 2014, 344, 1373–1377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yang, W.; Chen, I.H.; Gludovatz, B.; Zimmermann, E.A.; Ritchie, R.O.; Meyers, M.A. Natural flexible dermal armor. Adv. Mater. 2013, 25, 31–48. [Google Scholar] [CrossRef] [PubMed]
  4. Munch, E.; Launey, M.E.; Alsem, D.H.; Saiz, E.; Tomsia, A.P.; Ritchie, R.O. Tough, Bio-Inspired Hybrid Materials. Science 2008, 322, 1516–1520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Slesarenko, V.; Volokh, K.Y.; Aboudi, J.; Rudykh, S. Understanding the strength of bioinspired soft composites. Int. J. Mech. Sci. 2017, 131–132, 171–178. [Google Scholar] [CrossRef]
  6. Gu, G.X.; Takaffoli, M.; Buehler, M.J. Hierarchically Enhanced Impact Resistance of Bioinspired Composites. Adv. Mater. 2017, 29, 1700060. [Google Scholar] [CrossRef]
  7. Ortiz, C.; Boyce, M. Bioinspired structural materials. Science 2008, 319, 1053. [Google Scholar] [CrossRef]
  8. Chen, Z.; Wang, Z.; Zhou, S.; Shao, J.; Wu, X. Novel Negative Poisson’s Ratio Lattice Structures with Enhanced Stiffness and Energy Absorption Capacity. Materials 2018, 11, 1095. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, Y.; Zhang, X.; Kadic, M.; Liang, G. Stiffer, Stronger and Centrosymmetrical Class of Pentamodal Mechanical Metamaterials. Materials 2019, 12, 3470. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, J.-H.; Singer, J.P.; Thomas, E.L. Micro-/Nanostructured Mechanical Metamaterials. Adv. Mater. 2012, 24, 4782–4810. [Google Scholar] [CrossRef]
  11. Jackson, J.A.; Messner, M.C.; Dudukovic, N.A.; Smith, W.L.; Bekker, L.; Moran, B.; Golobic, A.M.; Pascall, A.J.; Duoss, E.B.; Loh, K.J.; et al. Field responsive mechanical metamaterials. Sci. Adv. 2018, 4, eaau6419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Florijn, B.; Coulais, C.; van Hecke, M. Programmable Mechanical Metamaterials. Phys. Rev. Lett. 2014, 113, 175503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shim, J.; Shan, S.; Košmrlj, A.; Kang, S.H.; Chen, E.R.; Weaver, J.C.; Bertoldi, K. Harnessing instabilities for design of soft reconfigurable auxetic/chiral materials. Soft Matter 2013, 9, 8198–8202. [Google Scholar] [CrossRef]
  14. Grima, J.N.; Caruana-Gauci, R.; Dudek, M.R.; Wojciechowski, K.W.; Gatt, R. Smart metamaterials with tunable auxetic and other properties. Smart Mater. Struct. 2013, 22, 084016. [Google Scholar] [CrossRef]
  15. Vangelatos, Z.; Gu, G.X.; Grigoropoulos, C.P. Architected metamaterials with tailored 3D buckling mechanisms at the microscale. Extreme Mech. Lett. 2019, 33, 100580. [Google Scholar] [CrossRef]
  16. Gao, C.; Slesarenko, V.; Boyce, M.C.; Rudykh, S.; Li, Y. Instability-Induced Pattern Transformation in Soft Metamaterial with Hexagonal Networks for Tunable Wave Propagation. Sci. Rep. 2018, 8, 11834. [Google Scholar] [CrossRef]
  17. Li, J.; Slesarenko, V.; Rudykh, S. Auxetic multiphase soft composite material design through instabilities with application for acoustic metamaterials. Soft Matter 2018, 14, 6171–6180. [Google Scholar] [CrossRef]
  18. Li, J.; Slesarenko, V.; Galich, P.I.; Rudykh, S. Instabilities and pattern formations in 3D-printed deformable fiber composites. Compos. Part B Eng. 2018, 148, 114–122. [Google Scholar] [CrossRef]
  19. Che, K.; Yuan, C.; Wu, J.; Jerry Qi, H.; Meaud, J. Three-Dimensional-Printed Multistable Mechanical Metamaterials With a Deterministic Deformation Sequence. J. Appl. Mech. 2016, 84, 11004. [Google Scholar] [CrossRef]
  20. Berwind, M.F.; Kamas, A.; Eberl, C. A Hierarchical Programmable Mechanical Metamaterial Unit Cell Showing Metastable Shape Memory. Adv. Eng. Mater. 2018, 20, 1800771. [Google Scholar] [CrossRef]
  21. Frenzel, T.; Kadic, M.; Wegener, M. Three-dimensional mechanical metamaterials with a twist. Science 2017, 358, 1072–1074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Fernandez-Corbaton, I.; Rockstuhl, C.; Ziemke, P.; Gumbsch, P.; Albiez, A.; Schwaiger, R.; Frenzel, T.; Kadic, M.; Wegener, M. New Twists of 3D Chiral Metamaterials. Adv. Mater. 2019, 31, 1807742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Mukhopadhyay, T.; Ma, J.; Feng, H.; Hou, D.; Gattas, J.M.; Chen, Y.; You, Z. Programmable stiffness and shape modulation in origami materials: Emergence of a distant actuation feature. Appl. Mater. Today 2020, 19, 100537. [Google Scholar] [CrossRef]
  24. Coulais, C.; Teomy, E.; De Reus, K.; Shokef, Y.; Van Hecke, M. Combinatorial design of textured mechanical metamaterials. Nature 2016, 535, 529–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hedayati, R.; Mirzaali, M.J.; Vergani, L.; Zadpoor, A.A. Action-at-a-distance metamaterials: Distributed local actuation through far-field global forces. APL Mater. 2018, 6, 036101. [Google Scholar] [CrossRef] [Green Version]
  26. Mirzaali, M.J.; Janbaz, S.; Strano, M.; Vergani, L.; Zadpoor, A.A. Shape-matching soft mechanical metamaterials. Sci. Rep. 2018, 8, 965. [Google Scholar] [CrossRef]
  27. Wang, Q.; Jackson, J.A.; Ge, Q.; Hopkins, J.B.; Spadaccini, C.M.; Fang, N.X. Lightweight Mechanical Metamaterials with Tunable Negative Thermal Expansion. Phys. Rev. Lett. 2016, 117, 175901. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, C.; Boorugu, M.; Dopp, A.; Ren, J.; Martin, R.; Han, D.; Choi, W.; Lee, H. 4D printing reconfigurable, deployable and mechanically tunable metamaterials. Mater. Horiz. 2019, 6, 1244–1250. [Google Scholar] [CrossRef]
  29. Yuan, C.; Mu, X.; Dunn, C.K.; Haidar, J.; Wang, T.; Jerry Qi, H. Thermomechanically Triggered Two-Stage Pattern Switching of 2D Lattices for Adaptive Structures. Adv. Funct. Mater. 2018, 28, 1705727. [Google Scholar] [CrossRef]
  30. Liu, J.; Gu, T.; Shan, S.; Kang, S.H.; Weaver, J.C.; Bertoldi, K. Harnessing Buckling to Design Architected Materials that Exhibit Effective Negative Swelling. Adv. Mater. 2016, 6619–6624. [Google Scholar] [CrossRef]
  31. Yu, K.; Fang, N.X.; Huang, G.; Wang, Q. Magnetoactive Acoustic Metamaterials. Adv. Mater. 2018, 30, 1706348. [Google Scholar] [CrossRef] [PubMed]
  32. Harne, R.L.; Deng, Z.; Dapino, M.J. Adaptive magnetoelastic metamaterials: A new class of magnetorheological elastomers. J. Intell. Mater. Syst. Struct. 2018, 29, 265–278. [Google Scholar] [CrossRef]
  33. Gu, G.-Y.; Zhu, J.; Zhu, L.-M.; Zhu, X. A survey on dielectric elastomer actuators for soft robots. Bioinspir. Biomim. 2017, 12, 011003. [Google Scholar] [CrossRef] [PubMed]
  34. Ahamed, R.; Choi, S.-B.; Ferdaus, M.M. A state of art on magneto-rheological materials and their potential applications. J. Intell. Mater. Syst. Struct. 2018, 29, 2051–2095. [Google Scholar] [CrossRef]
  35. Fang, H.; Chang, T.-S.; Wang, K.W. Magneto-origami structures: Engineering multi-stability and dynamics via magnetic-elastic coupling. Smart Mater. Struct. 2019, 29, 015026. [Google Scholar] [CrossRef]
  36. Tan, X.; Chen, S.; Wang, B.; Zhu, S.; Wu, L.; Sun, Y. Design, fabrication, and characterization of multistable mechanical metamaterials for trapping energy. Extreme Mech. Lett. 2019, 28, 8–21. [Google Scholar] [CrossRef]
  37. Tan, X.; Wang, B.; Yao, K.; Zhu, S.; Chen, S.; Xu, P.; Wang, L.; Sun, Y. Novel multi-stable mechanical metamaterials for trapping energy through shear deformation. Int. J. Mech. Sci. 2019, 164, 105168. [Google Scholar] [CrossRef]
  38. Dudek, K.K.; Wolak, W.; Gatt, R.; Grima, J.N. Impact resistance of composite magnetic metamaterials. Sci. Rep. 2019, 9, 3963. [Google Scholar] [CrossRef] [Green Version]
  39. Dudek, K.K.; Gatt, R.; Dudek, M.R.; Grima, J.N. Negative and positive stiffness in auxetic magneto-mechanical metamaterials. Proc. R. Soc. Math. Phys. Eng. Sci. 2018, 474, 20180003. [Google Scholar] [CrossRef]
  40. Dudek, K.K.; Gatt, R.; Grima, J.N. 3D composite metamaterial with magnetic inclusions exhibiting negative stiffness and auxetic behaviour. Mater. Des. 2020, 187, 108403. [Google Scholar] [CrossRef]
  41. Schaeffer, M.; Ruzzene, M. Wave propagation in multistable magneto-elastic lattices. Int. J. Solids Struct. 2015, 56–57, 78–95. [Google Scholar] [CrossRef]
  42. Dudek, K.K.; Wojciechowski, K.W.; Dudek, M.R.; Gatt, R.; Mizzi, L.; Grima, J.N. Potential of mechanical metamaterials to induce their own global rotational motion. Smart Mater. Struct. 2018, 27, 055007. [Google Scholar] [CrossRef]
  43. Grima, J.N.; Evans, K.E. Auxetic behavior from rotating squares. J. Mater. Sci. Lett. 2000, 19, 1563–1565. [Google Scholar] [CrossRef]
  44. Ishibashi, Y.; Iwata, M. A Microscopic Model of a Negative Poisson’s Ratio in Some Crystals. J. Phys. Soc. Jpn. 2000, 69, 2702–2703. [Google Scholar] [CrossRef]
  45. Grima, J.N.; Alderson, A.; Evans, K.E. Auxetic behaviour from rotating rigid units. Phys. Status Solidi B 2005, 242, 561–575. [Google Scholar] [CrossRef]
  46. Grima, J.N.; Zammit, V.; Gatt, R.; Alderson, A.; Evans, K.E. Auxetic behaviour from rotating semi-rigid units. Phys. Status Solidi B 2007, 244, 866–882. [Google Scholar] [CrossRef]
  47. Dudek, K.K.; Gatt, R.; Mizzi, L.; Dudek, M.R.; Attard, D.; Evans, K.E.; Grima, J.N. On the dynamics and control of mechanical properties of hierarchical rotating rigid unit auxetics. Sci. Rep. 2017, 7, 46529. [Google Scholar] [CrossRef] [Green Version]
  48. Dudek, M.R.; Wojciechowski, K.W. Magnetic films of negative Poisson’s ratio in rotating magnetic fields. J. Non-Cryst. Solids 2008, 354, 4304–4308. [Google Scholar] [CrossRef]
  49. Dudek, M.R.; Grabiec, B.; Wojciechowski, K.W. Molecular dynamics simulations of auxetic ferrogel. Rev. Adv. Mater. Sci. 2007, 14, 167–173. [Google Scholar]
  50. Gatt, R.; Mizzi, L.; Azzopardi, J.I.; Azzopardi, K.M.; Attard, D.; Casha, A.; Briffa, J.; Grima, J.N. Hierarchical Auxetic Mechanical Metamaterials. Sci. Rep. 2015, 5, 8395. [Google Scholar] [CrossRef]
  51. Saxena, K.K.; Das, R.; Calius, E.P. Three Decades of Auxetics Research − Materials with Negative Poisson’s Ratio: A Review. Adv. Eng. Mater. 2016, 18, 1847–1870. [Google Scholar] [CrossRef]
  52. Griffiths, D.J. Introduction to Electrodynamics; Pearson: Boston, MA, USA, 2013; p. 255. [Google Scholar]
  53. Slesarenko, V.; Rudykh, S. Towards mechanical characterization of soft digital materials for multimaterial 3D-printing. Int. J. Eng. Sci. 2018, 123, 62–72. [Google Scholar] [CrossRef] [Green Version]
  54. Lakes, R.; Wojciechowski, K.W. Negative compressibility, negative Poisson’s ratio, and stability. Phys. Status Solidi B 2008, 245, 545–551. [Google Scholar] [CrossRef]
Figure 1. Design of studied mechanical metamaterial. (a) elementary cell in the open state; (b) elementary cell in the closed state; (c) modified (repulsive) elementary cell with co-directional magnets; (d) modified (attractive) elementary cell with counter-directional magnets; and (e) unit cell with embedded modified cell (shown by blue color) under biaxial compression.
Figure 1. Design of studied mechanical metamaterial. (a) elementary cell in the open state; (b) elementary cell in the closed state; (c) modified (repulsive) elementary cell with co-directional magnets; (d) modified (attractive) elementary cell with counter-directional magnets; and (e) unit cell with embedded modified cell (shown by blue color) under biaxial compression.
Materials 13 01313 g001
Figure 2. Dependencies of the force between two identical magnets with out-of-plane magnetic moments. Solid lines correspond to experimental measurements, solid symbols represent calculations according to Equation (2) for several values of the magnetic moment.
Figure 2. Dependencies of the force between two identical magnets with out-of-plane magnetic moments. Solid lines correspond to experimental measurements, solid symbols represent calculations according to Equation (2) for several values of the magnetic moment.
Materials 13 01313 g002
Figure 3. Numerically obtained force–strain curves for the metamaterial with unit cell (Figure 1e) modified by the embedment of the magnets in repulsive and attractive configurations. Solid lines correspond to the results for the magnets with m = 0.105   Am 2 , dashed lines for the magnets with m = 0.075   Am 2 . Blue circles on the horizontal axis denote the strain ε s t corresponding to the snap-through phenomena. The numbers near the snapshots represent applied strain ε .
Figure 3. Numerically obtained force–strain curves for the metamaterial with unit cell (Figure 1e) modified by the embedment of the magnets in repulsive and attractive configurations. Solid lines correspond to the results for the magnets with m = 0.105   Am 2 , dashed lines for the magnets with m = 0.075   Am 2 . Blue circles on the horizontal axis denote the strain ε s t corresponding to the snap-through phenomena. The numbers near the snapshots represent applied strain ε .
Materials 13 01313 g003
Figure 4. Numerically obtained dependencies of the energy of the unit cell on the applied strain ε . (a) Total energy vs. ε for various configurations of magnets. Solid lines correspond to the magnets with m = 0.105   Am 2 , dashed lines to the magnets with m = 0.075   Am 2 . (b) Magnetic energy for repulsive configuration (red), attractive configuration (blue), and elastic energy for neutral configuration (dashed). Note that all energy curves are shifted vertically to have a (0, 0) origin.
Figure 4. Numerically obtained dependencies of the energy of the unit cell on the applied strain ε . (a) Total energy vs. ε for various configurations of magnets. Solid lines correspond to the magnets with m = 0.105   Am 2 , dashed lines to the magnets with m = 0.075   Am 2 . (b) Magnetic energy for repulsive configuration (red), attractive configuration (blue), and elastic energy for neutral configuration (dashed). Note that all energy curves are shifted vertically to have a (0, 0) origin.
Materials 13 01313 g004
Figure 5. Arrangements of the magnets in 16 × 16 metamaterial used for experiments and simulations. (a) Blue squares mark the positions of the magnets for the arrangement I, white circles mark the positions for the arrangement II. (b) Example of 3D printed metamaterial with embedded magnets (arrangement II). Four magnets in the corners are used for the mounting in the fixtures and do not affect the behavior of the metamaterial during tests.
Figure 5. Arrangements of the magnets in 16 × 16 metamaterial used for experiments and simulations. (a) Blue squares mark the positions of the magnets for the arrangement I, white circles mark the positions for the arrangement II. (b) Example of 3D printed metamaterial with embedded magnets (arrangement II). Four magnets in the corners are used for the mounting in the fixtures and do not affect the behavior of the metamaterial during tests.
Materials 13 01313 g005
Figure 6. (a) Force–displacement curves obtained experimentally (solid lines) and numerically (dashed lines) for finite 16 × 16 metamaterials with arrangements I and II. Solid lines show the means for three consecutive loadings of the metamaterial, while gray areas correspond to the standard deviation. (b,c) Two stable states of the metamaterial (arrangement I), observed experimentally. The snap-through occurs at the displacement marked by the blue circle (Figure 6a).
Figure 6. (a) Force–displacement curves obtained experimentally (solid lines) and numerically (dashed lines) for finite 16 × 16 metamaterials with arrangements I and II. Solid lines show the means for three consecutive loadings of the metamaterial, while gray areas correspond to the standard deviation. (b,c) Two stable states of the metamaterial (arrangement I), observed experimentally. The snap-through occurs at the displacement marked by the blue circle (Figure 6a).
Materials 13 01313 g006

Share and Cite

MDPI and ACS Style

Slesarenko, V. Planar Mechanical Metamaterials with Embedded Permanent Magnets. Materials 2020, 13, 1313. https://doi.org/10.3390/ma13061313

AMA Style

Slesarenko V. Planar Mechanical Metamaterials with Embedded Permanent Magnets. Materials. 2020; 13(6):1313. https://doi.org/10.3390/ma13061313

Chicago/Turabian Style

Slesarenko, Viacheslav. 2020. "Planar Mechanical Metamaterials with Embedded Permanent Magnets" Materials 13, no. 6: 1313. https://doi.org/10.3390/ma13061313

APA Style

Slesarenko, V. (2020). Planar Mechanical Metamaterials with Embedded Permanent Magnets. Materials, 13(6), 1313. https://doi.org/10.3390/ma13061313

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

Article Metrics

Back to TopTop