A Review on Tailoring Stiffness in Compliant Systems, via Removing Material: Cellular Materials and Topology Optimization

Cellular Materials and Topology Optimization use a structured distribution of material to achieve specific mechanical properties. The controlled distribution of material often leads to several advantages including the customization of the resulting mechanical properties; this can be achieved following these two approaches. In this work, a review of these two as approaches used with compliance purposes applied at flexure level is presented. The related literature is assessed with the aim of clarifying how they can be used in tailoring stiffness of flexure elements. Basic concepts needed to understand the fundamental process of each approach are presented. Further, tailoring stiffness is described as an evolutionary process used in compliance applications. Additionally, works that used these approaches to tailor stiffness of flexure elements are described and categorized. Finally, concluding remarks and recommendations to further extend the study of these two approaches in tailoring the stiffness of flexure elements are discussed.


Introduction
Compliance is the term used when referring to low levels of stiffness, in other words, when flexibility is the characteristic feature. Compliant systems can be categorized in compliant mechanisms (Figure 1a), compliant joints (Figure 1b), and flexure elements (Figure 1c). Compliant mechanisms transmit motion, force, and energy through the elastic deformation of flexible members known as flexures. Compliant joints are the key-components in compliant mechanisms; they act analogously to the kinematic pairs in rigid-body mechanisms. The basic form of a compliant joint is a single flexure, which is usually formed by removing material from a bulk piece. A compliant joint can be composed of more than one flexure and could be considered as a compliant mechanism on its own. Flexures stretch, bend, or twist during deformation, which leads to the relative motion between parts in a compliant system. Contact-aided compliant systems act as stiffness switches during deformation and are suitable only in applications where abrupt changes of stiffness are needed. Approaches to tailor the stiffness of compliant systems based on selectively removing material from flexure elements maintaining their overall dimensions appear as the possible solution to deal with inaccuracies and limitations of compliant systems. These approaches termed as Cellular Materials and Topology Optimization allow an accurate modification of the stiffness of flexure elements avoiding the use of external elements and actuators. Tailoring stiffness via Cellular Materials is based on using predefined periodic arrays combining material and void spaces to substitute the flexure elements of compliant systems [22]. The process to tailor the stiffness of a flexure element via Cellular Materials is graphically explained in Figure 2a. Topology Optimization has been widely used in the design and synthesis of compliant systems [53,54], and even to obtain specific stiffness values [55]. However, the resulting compliant systems often have complex architectures and distributed compliance, thus leading to complicated analysis, manufacture, and further implementation. Therefore, recent efforts have been focused on applying Topology Optimization processes only in the flexible parts (flexures) of compliant systems allowing an acceptable characterization of the stiffness while maintaining original system macro-shape and thus manufacturability. Then, tailoring stiffness via Topology Optimization consists of removing material from flexure elements based on a goal function [56]. The process to tailor the stiffness of a flexure element via Topology Optimization is graphically explained in Figure 2b. Both techniques are intimately related, as shown in Figure 2, as to achieve the optimal configuration of the cellular structure, it is necessary to perform an optimization process commonly called size or shape optimization, and the result of a Topology Optimization process after a postprocessing refinement is often a cellular structure. Optimization via cellular materials is often done intuitively and manually, while Topology optimization is automated via computational algorithms.
While both Cellular Materials and Topology Optimization can be reviewed independently for other applications, here we focused on tailoring stiffness via these two approaches acting at flexure level that has applications in compliant systems. The con-trolled distribution of material within a certain structure often leads to several advantages; one of them is the customization of the resulting mechanical properties of it. Following these two approaches, this can be achieved; therefore, reviewing the related literature to aid the understating and implementation of the approaches gains relevance. The main effort of this work is in clarifying how the stiffness can be tailored through these approaches and reviewing the works found in the literature that used them on flexure elements. Comments on the mechanical characteristics that the compliant systems gained with these methods, such as weight reduction, stiffness ratios, stress relief and motion range, are also addressed. As this work focuses on tailoring stiffness, applications of these methods in the customization of properties, such as thermal [57], electrical [58], and acoustic [59], are beyond of the scope of this review. This review is structured as follows. Tailoring stiffness via Cellular Materials is presented in Section 2. Tailoring stiffness via Topology Optimization is detailed in Section 3. Concluding remarks and discussion on research trends of these two are given in Section 4.

Tailoring Stiffness via Cellular Materials
Cellular Materials are essentially patterns that, in contrast to their homogeneous counterparts, are heterogeneous materials that have a unit-cell and repetition which leads to a structured distribution of material and voids [60]. Cellular Materials can be random or periodic, the latter are defined by a unit-cell that is repeated (combining material and voids). Therefore, a cellular solid is one made up of an interconnected network of solid struts or plates which form the edges and faces of cells. Typical cellular structures are shown in Figure 3  Commonly, open-cell structures are those where the void space is interconnected with a network of slender struts (open-cell structures can also be composed of surface shells, such as gyroids [61,62]); these types of Cellular Materials are among the most popular choices in design and are commonly called lattice structures, or simply lattices [63], an example of a 3D lattice is shown in Figure 3b. On the other hand, closed-cell structures are those where the network consists of interconnected surfaces that separate voids by cell walls or membranes, an example of a 3D closed-cell Cellular Material is shown in Figure 3c.
Cellular Materials are designed to offer advantages that cannot be afforded by homogeneous structures, such as the ability to locally tune properties, add multi-functionality, and improve performance while maintaining structural integrity and decreasing weight [64]. Due to their great stiffness-weight ratio, lattice structures have been used extensively in lightweight applications to reduce mass while maintaining or even increasing stiffness levels [65][66][67][68]. Most of these works used an optimization process (often called size and shape optimization) to select the most efficient configuration of unit-cell, dimensions of struts elements of the lattice, and the best way to infill bodies with this lattice. All of these processes have the aim of achieving substantial mass reduction and high structural stability while ensuring manufacturability.
When the characteristic lattice dimensions are significantly smaller than the overall dimensions of a piece of solid, one can represent the properties of such structural solids by means of apparent properties considering them as homogeneous having equivalent continuum properties [69]. In general, an apparent property, P , of a Cellular Material is often be expressed as P = CPρ n , where P is the property of the parent material (e.g., Young's modulus), and n and C are constants related to the deformation mechanism of the lattice elements and topology. The exponent n, which depends on the micro-architecture response, defines the scaling relation of the apparent property with the relative density [70]. The relative density is equal to the volume fraction,ρ = V s /V c , where V s is the volume of the solid material and V c is the volume of the lattice material considering the external dimensions.
According to their deformation mechanism, cellular Materials are classified into two main groups: bending-dominated and stretch-dominated structures [71]. When loaded, the cellular elements of the structured material can respond by bending or stretchingshortening; thus, they are said to be bending-or stretch-dominated (Figure 4a,b), respectively. Therefore, customization of the stiffness via Cellular Materials is possible as their apparent mechanical properties show dependence on the following: (i) the mechanical properties of the parent material, (ii) the topology and cell geometry, and (iii) the relative density. Properly planned selection or modification on any of these three can lead to a meta-material with tailored properties. In this regard, the stiffness of a compliant system can be modified through the implementation of specific cellular arrays in the flexible parts of the system.  Tailoring stiffness has become more relevant with the emergence of modern manufacturing processes such as additive manufacturing techniques, capable of fabricating complex topologies and shapes [72][73][74][75]. Some studies [76][77][78] showed how the stiffness in beams can be customized by using different types of cellular structures while obtaining minimum weight designs. Other works [79][80][81] focused on obtaining optimal designs of lattice structures to achieve the desired stiffness along with exceptional load-bearing efficiency and minimum weight. Later, Maskery et al. [82] investigated the effect of cell topology, orientation, and volume fraction to provide a specified stiffness in the lattice structures. Further, Weisgraber et al. [83] conducted research works where they proved that by modifying cell parameters, it is possible to control the stress-strain curve of a Cellular Material to assess specific load requirements and performance. A similar study was developed by Wang et al. [84], where the authors optimized lattice structures to provide specific stiffness for given multi-loading scenarios in the low volume fraction limit. Most recently, Alvarez-Trejo et al. [85] synthesized a new type of Bézier-based metamaterial, characteristics of which enable the generation of structures with customized mechanical properties by adjusting the geometry of the unit cells.
Some studies [86][87][88][89][90] demonstrated the possibility to tailor specific stiffness depending on the loading direction, allowing the design of customized anisotropic Cellular Materials. Niknam et al. [87] introduced the concept of multi-directional functionally graded Cellular Materials made by assembling porous unit cells of dissimilar densities and cell topologies to improve the performance of lightweight structural elements. In [88], the authors evidenced that tailoring the apparent anisotropy in Cellular Materials at the macroscopic level is possible by changing the ratio of bending-to-axial stiffness of the constituent beams, as, depending on the direction of loading, the nature of the mechanical response of cubic elementary trusses changes smoothly from bending-dominated to stretching-dominated.
This concept was also studied in [86], where the authors performed customization of effective anisotropic material properties in micro-lattice structures by controlling the unit-cell topology. Recently, Chen et al. [90] designed isotropic Cellular Materials created by superimposing complementary anisotropic lattices, enhancing the potential of high-stiff isotropic Cellular Materials. Variable stiffness is also achieved by including contacts during deformation of Cellular Materials. An example is found in Tanaka et al. [91], where the authors synthesized a square-cell Cellular Material with the bi-stiffness property in response to compressive loading. The specific case of stiffness depends on the orientation of the unit cells leading to a specific deformation of cells and further overlapping of cell walls. Kim et al. [89] synthesized a 3D variable stiffness and dual Poisson's ratio Cellular Material by including internal contacts between slit surfaces during both tensile and compressive loading achieving engineered stiffness.
Applications of tailored Cellular Materials have been focused mainly on stiffness maximization. Examples of Cellular Materials customized to offer maximum possible stiffness using different cellular structures such as 2D honeycombs [92,93], 3D lattices [94,95], and 3D closed-cell foams [96] are found in the literature. Tailored Cellular Materials have been also used to maximize transverse shear stiffness in honeycomb cores [97,98]. Studies of applications including the use of Cellular Materials to enhance buckling behavior [99], tune stress localization, minimize stress levels [100], and engineer failure location for specific load conditions [101] are limited.
Optimization via Cellular Materials can also be used for match specific stiffness. This aspect is particularly attractive for the manufacturing of bone implants, where the elastic modulus of the implant should match that of the bone in order to avoid stress shielding effects, bone atrophy, and implant loosening. Additionally, offering advantages such as improving stability by promoting tissue ingrowth and high energy absorption to withstand impacts. Three-dimensional topologies are the most frequently encountered in this application. A variety of works for this application can be found in the literature. Parthasarathy et al. [102] presented a thorough review where they expose the capacity of additive manufacturing techniques to fabricate biomedical devices with mechanical properties of surrounding bone tissue with replacements parts.
Most recently, tailoring stiffness through Cellular Materials to enhance flexibility or compliance behavior has been studied. Some of the first examples can be found in the works of Kim et al. [103,104], where the authors proposed Cellular Materials with customized anisotropic stiffness. The proposed materials consisted of flexible mesostructures designed from the deformation of re-entrant auxetic (negative Poisson's ratio) honeycomb and compliant amplifiers mechanisms. They could tailor in-plane properties of the said mesostructures, obtaining low stiffness and high strain in one direction and high stiffness and low strain in the other direction. A complete review of auxetic cellular materials is found in [105]. After, Duoss et al. [106] demonstrated the possibility to independently tailor compression and shear response in 3D printed Cellular Materials offering the extremal property of negative shear stiffness in designed structures. Then, Neff et al. [107] used modified diamond lattice structures to create Cellular Materials with reduced stiffness, advantageous for energy absorption, vibration isolation, and reduction of stress. Later, Runkel et al. [108] modified the topology of hexagonal chiral lattices and demonstrated that purposely choosing the parameters permits to tailor the behavior of the chiral lattice, enabling extreme variations in stiffness, including the possibility to attain negative and zero-stiffness regimes over large compressive strains which exemplifies the potential to tune deformation of compliant systems. Additionally, Ion et al. [109] introduced the metamaterial mechanisms concept, which considers cellular materials as machines instead of materials and allows the creation of macroscopic movement through the well-defined interaction of cells in a single block of cellular material. Similarly, Ou et al. [110] presented a group of auxetic-inspired material structures that can transform into various shapes upon compressive loading. Most recently, in [111] a 2D kagome honeycomb was optimized for delivering high amplitude actuation. In the same year, Dunn et al. [112] demonstrate that changing the cell size in lattice materials might involve maximizing stiffness in one mode while minimizing it in another.
Regarding this, some applications have taken advantage of the compliance behavior offered by Cellular Materials. One of the most extensively studied applications is in shape morphing as the material used should not only be in-plane flexible enough to morph with actuation, but stiff enough not to deform excessively under out-of-plane loads. To this end, 2D cellular structures have been suggested for general morphing skin applications [113][114][115]. In [113], the authors presented a novel zero-Poisson's ratio honeycomb structure that can achieve deformations along two orthogonal directions and avoid the increase of effective stiffness in the morphing direction by restraining the Poisson's effect in the non-morphing direction. The same concept was explored by Liu et al. [114,115] in the development of zero-Poisson's ratio cellular structures with lightweight, low in-plane moduli and high strain capability. Results demonstrated the potential of the proposed structures for in-plane morphing. Moreover, the main application is in aircraft morphing, where efforts have been mainly focused in airfoil-core morphing [116][117][118][119] and airfoil-skin morphing [120][121][122][123][124][125][126]. Bornengo et al. [116] proposed the use of an auxetic hexagonal structure as the airfoil-core of wings as an attempt to take advantage of the high deformability of the structure thus providing the flexibility required for morphing, while preserving the structural integrity due to its high in-plane shear resistance. Dong et al. [117] designed airfoil-cores with a re-entrant cellular structure due to their ability to undergo large in-plane displacements with limited deformation in span-wise direction. Zhang et al. [118] investigated the use of a cross-shaped honeycomb structure in airfoil-core morphing, resulting in great in-plane deformation capacity and adequate outof-plane support. Heo et al. [119] investigated the airfoil-core morphing capacity of three types of Cellular Materials through the in-plane flexible properties. Hexagonal chiral, regular honeycomb, and re-entrant honeycomb were chosen for this analysis. All structures showed high levels of in-plane flexibility. The re-entrant honeycomb showed the highest flexibility and causes a lower stress in cell walls than the other cellular cores, thus suggesting potential to be used in real morphing applications. On the other hand, Olympio et al. [120] studied the use of hexagonal honeycomb cellular structures in airfoil-skin morphing. Results showed that Cellular Materials could easily undergo global strains several times greater that the parent material while maintaining good out-of-plane stiffness; however, when honeycomb Cellular Materials are stretched along the principal axis they geometrically stiffen, thereby reducing the maximum global strain achievable. Olympio et al. [121] proposed two types of zero-Poisson's ratio cellular structures (hybrid and accordion) to be used in one-dimensional airfoilskin morphing applications. Results showed that the axial stiffness in the motion direction of the zero-Poisson's ratio cellular structures remains low when the skins are restrained in the nonmorphing direction, making these cellular structures ideal for use in airfoil-skin morphing. Other works [122][123][124] have studied similar zero-Poisson's ratio cellular structures to be used in airfoil-skin morphing, obtaining high levels of flexibility in the actuation direction while maintaining the desired shape and strength in the others. Most recently, in [125,126] a new Cellular Material with zero-Poisson's ratio called Fish cell (Figure 5b) was introduced for applications in airfoil-skin morphing. The proposed cellular structures were parametrized to obtain desired stiffness values under specific conditions. Finally, 3D arrays have been reserved for complete wing morphing [127,128]. 2D cellular structures used in airfoil-core and airfoil-skin morphing can be subdivided in three main categories according to the value of their Poisson's modulus (ν): positive-Poisson (ν > 0), zero-Poisson (ν ≈ 0), and auxetic (ν < 0). A complete overview of airfoil morphing by 2D cellular structures is depicted in Figure 5a, graphically showing the main categories of 2D patterns: positive-Poisson, zero-Poisson, and auxetic. In general, 2D cellular structures used in airfoil morphing applications showed lower levels of in-plane stiffness than their parent material, and higher levels of out-of-plane stiffness than fully dense elements of comparable mass, as well as mass reduction percentages above 90% and higher strain capabilities [120]. In airfoil cores, auxetic cellular structures are preferred due to their higher in-plane flexibility which causes lower stress in cell walls than other types of cellular cores [119]. On the other hand, airfoil-skin morphing commonly uses sandwiched skins with cellular structure cores and thin, high-strain-capable face sheets that provide a smooth aerodynamic surfaces. However, this sandwich configurations also involve dealing with the boundary effects and Poisson's contractions or expansions. Therefore, zero-Poisson cellular structures are preferred to deal with this restriction in unidirectional morphing applications such as airfoil-skin morphing. Tables 1 and 2 present the principal works of airfoil-core and airfoil-skin morphing applications using 2D Cellular Materials. Works in these tables are categorized by author and cellular structure used. Stiffness reductions are presented by in-plane Young's modulus ratios, E cx /E and E cy /E, for airfoil-core morphing (see Table 1), and by main-direction Young's modulus ratio, E cx /E, and shear-stiffness ratio, G c /E, for airfoil-skin morphing (see Table 2), where E c and G c are the Young and shear moduli of the Cellular Material, while E is the Young's modulus of the parent material. Additionally, the percentage of weight reduction is reported in Table 1. Tables 1 and 2 are presented to show the potential to perform high reductions of the stiffness levels in flexure elements through the substitution of fully dense solids with cellular patterns. Table 1. Stiffness reduction via Cellular Materials: Airfoil-core morphing. Principal works found in literature by author and cellular structure used (P: positive Poisson, Z: zero-Poisson, A: auxetic). Young's modulus ratios E c /E between the obtained by the application of Cellular Material, E c , and that of the parent material E for in-plane E cx /E and E cy /E stiffness, and percentage of weight reduction. "−" indicates that there is no sufficient information reported to obtain that parameter. The percentage of weight reduction was obtained based on the reported relative densityρ of the cellular structure. Other important application of the compliance behavior of Cellular Materials is in energy absorption. Cellular Materials are suitable for energy absorption applications as their particular structures permit high levels of elastic deflections without permanent deformations or high values of stress, allowing the material to spring back to its original configuration after the load is released. There are many works in the literature that prove the tailorability of Cellular Materials by manipulating their structural topology to enhance energy absorption capacity e.g., general energy absorption [129,130] or impact absorption [131][132][133][134]. Some studies generate negative stiffness materials [135,136], increasing the amount of energy absorbed, and others include contacts during the deformation of Cellular Materials to ensure restorability [137,138]. Many authors extensively studied the use of Cellular Materials in energy absorption applications, and a review of this application is out of the scope of this work. Therefore, interested readers are referred to specialized reviews in [139][140][141][142].

Author Cellular Young's Modulus Ratios Weight Reduction
In the process of tailoring stiffness via cellular materials approach, certain topologies are preferred depending on the stiffness objective (increase stiffness or reduce stiffness). When stiffness maximization is the objective, a stretch-dominated lattice and closed cell structures are preferred. On the other hand, if flexibility is the goal, a bending-dominated auxetic and zero-Poisson's ratio honeycombs are the most commonly used. The unit-cells of some of the most commonly used cellular structures for tailoring stiffness are shown in Figure 6. Increasing stiffness via cellular materials has been widely studied; therefore, unit cells for this objective, shown in Figure 6, are categorized in 2D-lattice structures, 3D-lattice structures, and Closed-cell structures. On the other hand, reducing stiffness via cellular materials has been limited to explore 2D-lattice structures; thus, unit cells for reduce stiffness, shown in Figure 6, are only in the category of 2D-lattice structures, nevertheless, they are separated according to their Poisson's ratio behavior in positive Poisson's ratio, auxetic, and zero-Poisson's ration. Unit cells for increase stiffness are categorized in: 2D lattice structures (square [92], triangular [93], and kagome [100]), 3D-lattice structures (octet truss [94], and body centered cubic [81]), and closed-cell structures (body centered cubic [96]). Unit cells for reduce stiffness are only in the category of 2D lattice structures, and they are separated according to their Poisson's ratio behavior in: positive Poisson's ratio (hexagonal [119]), auxetic (hexagonal chiral [116], and re-entrant [117]), and zero Poisson's ration (chevron [122], and accordion [121]). Some of the unit cells were generated using nTopology ®.
Additionally, to address properly the problem of tailoring stiffness via cellular materials, it is convenient to follow the next guideline. First, define the stiffness objective (increase or reduce stiffness). Then, identify the appropriate cellular material, for which it is necessary to identify the desired deformation mechanism (bending-dominated or stretch-dominated) and select an adequate topology (honeycomb, closed-cell structure, or open-cell structure). After, adjust the unit-cell parameters (geometry and relative density) to modify the stiffness of the cellular material. Afterward, compare the obtainable stiffness with the stiffness objective and evaluate if the cellular material is adequate. Next, if the objective stiffness is achieved, then the final step is the implementation. On the other hand, if the objective stiffness is not be achieved with the selected cellular material, then return to the topology identification step and repeat the process. The above-mentioned guideline is summarized in the diagram shown in Figure 7.

Cellular Materials in Compliant Systems as Flexure Elements
There is a great amount of evidence for the capability of Cellular Materials to achieve desired levels of compliance in certain directions while maintaining high levels of stiffness in others. However, only Merriam and Howell [22] have applied a cellular lattice structure directly in the flexures elements of compliant joints, as shown in Figure 8. In [22], Merriam and Howell introduced the concept of lattice flexure, which was not common before additive manufacturing techniques were widely available. They demonstrated that these types of flexures have a reduced bending stiffness when compared to traditional rectangular-section blade flexures of similar dimensions while maintaining comparable stiffness in off-axis directions. Making use of additive-manufactured flexures, they found that both type of lattice flexures-X-type and V-type-shown in Figure 8 have the potential to reduce the bending stiffness by 60-80% compared to a traditional blade flexure of similar size, and at the same time to exhibit a torsional/bending stiffness ratio as much as 1.7 times higher than an equal aspect-ratio blade flexure, and a transverse bending/motion-direction bending stiffness ratio up to 6.5 times higher than an equal aspect-ratio blade flexure.   [22] to be used in cross-axis flexural joints. (Adapted with permission from ref. [22]. 2021, Elsevier.) Aiming to modify or allow compliance, it is possible to tailor the stiffness of flexure elements via Cellular Materials, while obtaining other advantageous mechanical conditions such as mass reduction, stress relief, and increments in stiffness ratios. However, when specific stiffness values are desired, tailorability is possible by using optimization processes that modified the final configuration of predefined cellular structures (see Figure 2). Then, one may express the effective stiffness as a function of cell parameters (strut or wall thickness, length, width, angle, etc.). Although this approach to design is powerful, its disadvantage is that the true optimization of the resulting materials architecture is ultimately limited by the designers initial choice of cellular structure. Therefore, other types of optimization processes that overcome this limitation are reviewed in the following section. This approach, termed as Topology Optimization, results in the optimal topology by removing or adding material to a design domain.

Tailoring Stiffness via Topology Optimization
Topology Optimization is a computational design method, which varies the material distribution within a structure or part, considering a predefined performance objective. In contrast to size and shape optimization in the use of Cellular Materials, topology optimization does not require a predefined initial design, and it can generate optimal geometries when intuitive design approaches fail [143]. The general Topology Optimization problem can be described as follows: To find the material distribution that minimizes an objective function F, subjected to a volume constraint g 0 ≤ V 0 and possibly n other constraints g i ≤ 0, i = 1, . . . , n. The material distribution is described by the density variable ρ(x) that can take any value among 0 (void) and 1 (solid material) at any point in the discretized design domain Ω. This optimization problem can be written in a mathematical form as : ρ(x) = value among 0 and 1 ∀x ∈ Ω. (1) Considering the density approach, this Topology Optimization problem is typically solved by discretizing the domain Ω into m finite elements and density can take either the value 0 or 1. In this case, the last line of Equation (1) must be rearranged as ρ j = 0 or 1 for j = 1, . . . , m. On the other hand, the continuous Topology Optimization problem allows the finite elements to take values between 0 (void) and 1 (solid material). In this case, the last line of Equation (1) is reformulated as 0 ≤ ρ j ≤ 1 j = 1, . . . , m.
In general, Topology Optimization approaches are density-based and level set-based. In density-based methods [144], which include the popular Solid Isotropic Material with Penalization (SIMP) method [145], the geometry composed of finite elements is described via a material distribution of two or more phases, of which usually one represents "no material" (i.e., the void phase) leading to jagged boundary geometries. On the other hand, level setbased methods [143] define the interfaces between material phases implicitly by iso-contours of a level set function allowing a crisp description of the boundaries and the accurate enforcement of boundary conditions. For further details on Topology Optimization methods, solution strategies, algorithms, regularization methods, etc. refer to the works in [146][147][148][149].
Thus, by selecting a specific objective function F in Equation (1), it is possible to create structures with customized properties, e.g., stiffness. When the design domain has a simple geometry and is free of specific boundary conditions, the topology-optimized result is often a cellular pattern [56] (Figure 2b). In addition, structures resulting from Topology Optimization processes need a postprocessing refinement in order to smooth the jagged curves generated from the discretized domain and to achieve manufacturability [150]. As with the use of Cellular Materials, the main efforts of Topology Optimization processes have been focused in reducing mass while maximizing stiffness [151][152][153][154]. Then, the basic process shown in Figure 2b is hereafter termed as Topology Optimization Resulting in Cellular Patterns (TORCP).
There are other processes that combine Topology Optimization and Cellular Materials concepts in order to achieve specific stiffness values. For instance, Topology Optimization can be adopted to guide the material distribution and final configuration of unit cells in Cellular Materials. After the optimized unit-cell is obtained, a process to obtain the properties of the macromaterial in terms of the effective properties of its unit cell, termed as homogenization [155], is applied to approximate the apparent properties of the Cellular Material [146,156]. This widely used methodology is termed here as Topology Optimization Applied to Unit Cells (TOAUC), and it is graphically explained in Figure 9. The first implementation of Topology Optimization applied to unit-cell design can be found in the works of Sigmund [151,157,158]. Most recently, Osanov and Guest [159] reviewed the key requirements to apply topology optimization process to generate unit cells of periodic structured materials for a variety of customizable properties. Topology Optimization applied to the unit-cells of Cellular Materials has been used to maximize stiffness under general conditions [160][161][162], shear forces [163], and predetermined loads [164,165]. Resulting structures from 3D maximization stiffness problems are commonly closed-cell cellular structures [159,166]. This type of structure is preferred because their wall elements spans and connects various faces of the unit-cell, carrying loads axially in more than one direction, thus exhibiting stretch behavior in all directions. However, an important disadvantage of closed-cell structures with located enclosed pores is that they are quite difficult to fabricated. Furthermore, TOAUC-resulting cellular materials have been used to match specific stiffness values in bone implants [167][168][169][170].  The maximization of stiffness by combination of the Topology Optimization and Cellular Materials techniques has also been explored. Two main concurrent methodologies were identified: The first begins with the typical Topology Optimization process to generate a conceptual design in the design domain. After that, the generated design is filled with a cellular structure especially selected to improve the optimization of the objective property. This concurrent method is termed here as Multiscale Topology Optimization with Cellular Materials (MTOCM), and it is graphically explained in Figure 10. Examples of relevant works using MTOCM approach can be found in [171][172][173][174][175]. The second concurrent methodology consists in modifying the relative density of the cellular design domain based on a Topology Optimization process, obtaining the socalled functionally graded Cellular Materials. Modification of the unit-cell relative density can be made by altering cell parameters [176] or by using different unit-cell topologies [177]. This methodology is termed as Topology Optimization of Functionally-graded Cellular Materials (TOFCM), and it is graphically explained in Figure 11. Other, examples of relevant works using TOFCM approach can be found in [178][179][180][181][182].
Topology optimization process Unit cell with three different relative densities Cellular design domain Functionally graded cellular material Figure 11. Process to tailor stiffness using Topology Optimization applied to cellular design domains that leads to functionally graded Cellular Materials.
Topology Optimization for tailoring stiffness with compliance objectives has been used to create compliant mechanisms [53], including micro-grippers [183], bridge-mechanisms [184], positioning stages [185], and path generating mechanisms [186]. However, using the Topology Optimization process to generate complete compliant mechanisms may lead to complex designs with distributed compliance throughout its topology, thus leading to complicated analysis, manufacturability, and further implementation [187]. Therefore, a new approach of Topology Optimization appeared as an alternative to deal with these drawbacks, maintaining the versatility of the process. This new approach is focused on applying the Topology Optimization process only on the flexures of compliant systems, taking the flexure section as the design domain of the Topology Optimization process. This could enhance the analysis of compliant systems through a proper characterization of their force deflection functions while maintaining their original macro-shape and lumped architecture.
In this regard, TORCP and TOAUC approaches are the most commonly employed when compliance is the goal property while applying Topology Optimization. They have been used to create cellular structures that are optimal in morphing applications [188][189][190][191], as well as in energy absorption applications [192][193][194]. Obtaining comparable levels of stiffness reductions and topologies with similar behaviors as the reported in Section 2. Moreover, the TOAUC method has been applied to create 2D [195][196][197] and 3D [198] extreme elastic materials that possess interesting properties such as a negative Poisson's ratio. The same process was used by Deng et al. [199] to create stretchable lattices with architected unit cells based on Bezier curves.

Topology Optimization in Compliant Systems as Flexure Elements
With the potential of Topology Optimization to design compliance materials already being stated, the works using this approach in compliance systems at flexure level are reviewed. For example, Topology Optimization has been used to design optimal compliant kinematic limbs that are part of larger compliant mechanisms [200][201][202]. A schematic description of this process applied to the limbs of 2D flexure-based compliant parallel mechanisms is shown in Figure 12. Lum et al. [200] performed a Topology Optimization process to synthesize and optimize the topology, shape, and size of the limbs in a flexure-based compliant parallel mechanism (see Figure 12a), resulting in a significant improvement in the mechanism's workspace dimensions and increments in the mechanism's stiffness ratios, i.e., ratios of off-axis to actuating stiffness. Another work that dealt with the optimization of compliant kinematic limbs was that of Jin et al. [201], where those of 2D compliant parallel mechanisms were regarded as separate design domains for a Topology Optimization problems (see Figure 12(b-i,c-i), and the Jacobian matrix was introduced into the field of Topology Optimization. This resulted in compliant parallel mechanisms with maximized workspaces that have compliant kinematic limbs accurately characterized facilitating the mechanism analysis. Later, Jin et al. [202] enhanced their previous research by including an improved stiffness evaluation generating new optimal compliant limbs and thus compliant parallel mechanisms (see Figure 12(b-ii,c-ii) that have better levels of off-axis stiffness, and are similar to the generated in their previous work. Most recently, a different technique, not applied at the flexure level, has been used to design the limbs of compliant parallel platforms; this is based on keeping the mobile and fixed platforms as rigid bodies while applying a Topology Optimization process on a design domain that connects both platforms, examples using this technique designing 2D and 3D compliant parallel platforms can be found in [203][204][205][206], respectively.  [202] (Adapted from [202] 2018, ASME).
Topology Optimization of single flexures, i.e., those that are not part of a compliant mechanism, has been also addressed. The most commonly studied single flexures are the Single-Axis Rotational Flexure Joints (SARFJ). Studies of Topology Optimization used in the design of SARFJ can be categorized according to the geometry of the design domain used in the process. Initial attempts were performed using 2D rectangular design domains as shown in Figure 13. Examples using this approach are the works of Zhu et al. [207][208][209]. In [207], Zhu et al. performed the design and optimization of SARFJ, searching to maximize the compliance in the desired direction while minimizing the compliance in the other directions. The design problem was modeled and analyzed as a small-displacement problem using linear Finite Element Analysis (FEA) resulting in the final design shown in Figure 13a. Compared with the constant rectangular cross section flexure hinge of similar mass, the optimized design showed more compliance, capable of undergoing a larger rotation angle when subjected to the same external load, while achieving rotation with better precision. However, the obtained design reaches higher levels of stress. Then, in [208], the same authors extended their previous work by including "fixing the coordinates of the rotational center O" as an additional objective in the Topology Optimization problem. This additional objective was attempted by setting the topology to be symmetric with respect to xand y-axes in the design domain resulting in the topology shown in Figure 13b. However, the resulting design showed lower precision than their previous design, demonstrating that minimizing the displacement of the point O must be added as a constraint in the Topology Optimization problem instead of an objective. Most recently, Zhu et al. [209] improved their work in [207] by maintaining the functional requirements of maximizing the compliance in the desired direction, while minimizing the compliance in the other directions and setting "minimizing the displacement of the rotational center O" and symmetry of the topology in both axes, as constraints of the process. The resulting design was similar to that of the cartwheel flexure joint (see Figure 13c), often regarded as a flexure joint with ultra-high rotation accuracy. Additional designs were created by changing the allowable material usage and the geometry of the design domain. In general, the resulting topologies showed great rotational precision and high levels of compliance in the motion direction.  Figure 13. Schematic description of the Topology Optimization process to design the optimal configuration of SARFJ using a 2D rectangular design domain. Optimal designs were obtained in (and adapted from) (a) [207]; (b) [208]; (c) [209]; (d) [210] i-design based on linear FEA and ii-design based on nonlinear FEA; (e) [211]; (f) [212] i-design with less stress allowed for a mesh of 50 × 50 and ii-design with less stress allowed for a mesh of 100 × 100; (g) [213]; (h) [214] i-design for general distributed flexure hinges, ii-design for distributed flexure hinges with prescribed compliances, and iii-design for distributed flexure hinges with minimal parasitic motion; and (i) [56] i-optimal topology, ii-simplified design named "Snowflake", and iii-simplified design named "'Half-Snowflake".
The Topology Optimization of SARFJ, using a 2D rectangular design domain, was addressed also by Liu et al. [210][211][212][213][214]. In [210], the authors dealt with the Topology Optimization of SARFJ undergoing large displacement by using a geometrically nonlinear elastic analysis in the Topology Optimization process. The optimization problem was formulated to set the stiffness of the flexure joint at the functional direction as low as possible and the stiffness at non-functional direction as high as possible. Additionally, minimizing the displacement of the center of rotation point and symmetry about x-axis were added as constraints. Two optimal designs were created: one based on linear FEA (see Figure 13(d-i) and other using nonlinear FEA (see Figure 13(d-ii). The resulting designs were compared with a filleted V-shape flexure hinge. Results showed that nonlinear design can reach the highest rotation angle and has the lowest level of stress concentration. On the other hand, the linear design has better precision but the highest level of stress. This was further extended in [211], with the design of a three-segment flexure joint named quasi-V-shaped flexure hinge (QVFH) shown in Figure 13e. The authors derived the in-plane compliance equations for small displacements based on Castigliano's second theorem. Compared with the filleted V-shaped flexure hinge, the QVFH showed less compliance, more precision in keeping the center of rotation with the minimal offset, and larger maximum stress level. After that, Liu et al. [212] designed flexure joints by using stress-constrained Topology Optimization, aiming to reduce the general stress level and remove the sharp corners of the flexure hinges, thereby reducing the stress concentration and improving the range of motion. Various topologies were created for several levels of stress limit. The resulting final topology without stress constraint is similar to the one obtained in [211]. It can be observed that when the stress constraints are added, the length and thickness of the middle segment become longer and thicker, respectively, thus removing sharp corners distributing the stress and reducing stress concentrations. Several designs were created varying the level of stress constraint and changing the mesh refinement of the design domain, only two final topologies are shown in Figure 13f, corresponding to the designs with less stress allowed for a mesh of 50 × 50 (see Figure 13(f-i) and for a mesh of 100 × 100 (see Figure 13(f-ii). As it can be seen, optimal designs with stress constraints tend to seem blade flexures, therefore losing precision.
In this regard, Liu et al. [213] presented a systematic method for designing flexure hinges with distributed stress by using Topology Optimization with the objective of maximizing the bending displacement while minimizing the axial displacement and at the same time to minimize the maximum level of stress. Additionally, a symmetry constraint in both xand y-axes was added, resulting in the non-conventional notch flexure hinge, shown in Figure 13g, that exhibited distributed stress. Finally, this work was extended in [214], where authors designed flexure hinges with a symmetrical distributed configuration that accomplish the desired performances of prescribed compliance and minimal parasitic motion. Several designs were obtained for each desired performance; only three of them are shown in Figure 13h, accounting for general distributed flexure hinges (see Figure 13(h-i), distributed flexure hinges with prescribed compliances (see Figure 13(h-ii), and distributed flexure hinges with minimal parasitic motion (see Figure 13(h-iii). Resulting configurations are similar to the cartwheel flexure hinge which, as was mentioned before, it is well known for being a high precision flexure joint. Finally, Pinskier et al. [56] used Topology Optimization to design SARFJ to account for the goals of motion-range, compliance, and precision simultaneously. The resulting main topology is shown in Figure 13(i-i). Based on this topology, a second design stage was performed to produce two simplified hinge designs named "Snowflake" and "Half-Snowflake" shown in Figure 13(i-ii,iii), respectively. After comparing their performance, the new type of flexure hinge showed higher precision and comparable levels of rotational compliance than blade and cross-axis flexures. Additionally, the resulting structures have extra members that could lead to an off-axis stiffness increment, reducing parasitic displacements and making them able to transfer forces.
Other optimal designs of SARFJ were created using 2D design domains, based on the geometry of a right circular flexure hinge as shown in Figure 14. For instance, Liu et al. [215] used Topology Optimization to remove material from a right circular flexure hinge in order to increase compliance while retaining high precision, resulting in the multinotched flexure hinge topology shown in Figure 14(a-i). Several topologies were obtained using different dimensions in the design domain and varying the percentage of removed material; however, they resulted in nearly the same for all examples. Then, a final shape of the multi-notched flexure hinge (see Figure 14(a-ii) was proposed after a postprocessing refinement stage and was used for comparison in its performance with various flexure hinges (right circular [216], corner-filleted [217], elliptical [218], parabolic, and hyperbolic [219]). Results showed that the multi-notched flexure hinges can have larger rotational compliance, higher rotational precision, and lower stress levels than any other type of flexure hinge when certain dimensional conditions are satisfied. Additionally, Zhu et al. [209] obtained a similar topology, shown in Figure 14b, using Topology Optimization with a multi-criteria objective function and a symmetry constraint. The resulting design showed great rotational precision and high levels of compliance in the motion direction.  Figure 14. Schematic description of the Topology Optimization process to design the optimal configuration of SARFJ using a 2D design domain based on the geometry of a right circular flexure hinge. Optimal designs obtained: (a) i-optimal topology and ii-final shape (Adapted with permission from ref. [215]. 2021, Elsevier) and (b) optimal topology obtained in [209] (Adapted from ref. [209]. 2019, Springer).
Most recently, 3D design domains were also used in Topology Optimization to generate SARFJ as shown in Figure 15. Qiu et al. [220] used 3D continuum Topology Optimization to find the optimal topology of a SARFJ in a 3D design domain. The resulting topology is the single-axis multi-cavity flexure hinge (MCFH) shown in Figure 15(a-i). After that, a postprocessing stage was used to simplify the configuration of the hinge resulting in the final design shown in Figure 15(a-ii). Dimensionless empirical equations of the stiffness performance were obtained based on FEA for this simplified MCFH. Then, the new flexure hinge was compared with other types of hinges demonstrating that the MCFH have greater compliance. Then, Qiu et al. [221] presented a dual-objective Topology Optimization process for maximization of compliance and precision performance of flexure hinges based on 3D continuum Topology Optimization. The resulting topology, shown in Figure 15(b-i), was further postprocessed and parameterized generating the single-axis quasi-leaf porous flexure hinge (QLPFH) shown in Figure 15(b-ii). Dimensionless empirical equations of the stiffness and rotational center performance were derived with the aid of FEA, and they were used for comparison with a leaf flexure hinge with the same envelope space dimensions. The results showed that the rotational compliance of QLPFH is much larger than of the leaf flexure. In addition, the precision performance of the QLPFH was also improved significantly. i ii Figure 15. Schematic description of the Topology Optimization process to design the optimal configuration of SARFJ using a 3D cubic design domain. Optimal designs obtained in (a) i-optimal topology and ii-final shape (Adapted with permission from ref. [220]. 2021, Elsevier), (b) i-optimal topology and ii-final shape (Adapted from ref. [221]. 2020, Springer).
Topology Optimization has also been used to find the optimal material distribution in a less common type of flexure joint called "compliant revolute joint" (Figure 16). The compliant revolute joint consists of a thin cylinder of compliant material with a inner hollow where its external face is fixed and its internal face is subject to a rotational displacement produced by a tangential force applied at the same face. Li et al. [222][223][224] presented a series of works dealing with the Topology Optimization of this type of flexure joint. In [222], the authors designed new compliant revolute joints reconstructed after removing material with a Topology Optimization process. The obtained designs are shown in Figure 16(i-iv). FEA was used to compare the performance of the optimized designs with the non-optimized compliant revolute joint. Results showed that optimized designs can achieve more rotational displacement. However, their stress levels were also increased. Further experimental validation in the rotational angle and stress levels of the optimal designs was performed in [223], obtaining a maximum error between experimental and FEA results of 8.7%. This was expanded in [224], by modeling the compliant revolute joint as a kinetoelastic model considering the kinematic properties and structural elastic properties simultaneously. This approach allowed obtaining optimized compliant revolute joints with desirable compliance properties that are quite similar to those obtained in their previous works. (i-iv) Optimal designs obtained in [222] (Adapted from ref. [222]. 2017, Springer).
Furthermore, Topology Optimization has been used to achieve the optimal material distribution of single-axis translational flexure joints (SATFJ) based on the geometry of blade flexures (see Figure 17). Pinskier et al. [2,52] investigated the design of this type of flexure to be used in parallelogram linear flexure guides using Topology Optimization to maximize their ratio of in-plane to out-of-plane stiffness. The resulting topology shown in Figure 17(a-i) resembles a X-type lattice structure similar to that used in [22] (see Figure 8). After that, a parameterized flexure design was generated based on the resulting topology (see Figure 17(a-ii), and it showed stiffness ratios of 1.5 to 6 times greater than standard blade structures with mass reductions up to 75%.

Rigid part
Design domain ii Figure 17. Schematic description of the Topology Optimization process to design the optimal configuration of SATFJ using a 2D design domain based on the geometry of blade flexures. (a) Optimal design obtained in [52] i-optimal topology and ii-final shape (reprinted with permission from [52]. 2019, Elsevier).

Remarks on Tailoring Stiffness via Cellular Materials and Topology Optimization
Cellular Materials and Topology Optimization can be used effectively for tailoring stiffness. The Cellular Materials approach is based on the proper selection of the unit-cell that constructs the cellular material; note that this approach works on the micro-scale. On the other hand, the Topology Optimization approach generates an optimal topology for specific external load conditions; then, it is said that this approach acts on the macroscale. Therefore, both approaches can be combined to take advantage of their particular characteristics at the same time. A summary that address advantages and disadvantages of Cellular Materials and Topology Optimization approaches as well as their combinations, is presented in the diagram in Figure 18.

Cellular Materials
Topology Optimization

TOAUC
Topology optimization applied to unit cells implements the optimal topology under specific conditions as the unit cell of cellular materials MTOCM Multiscale topology optimization with cellular materials extends the customization range by filling an optimal geometry with a cellular pattern TOFCM Topology optimization of functionally-graded cellular materials generates the optimal topology inside a cellular material by modifying the relative density of the unit cells in specific zones In addition, the implementation of Cellular Materials and Topology Optimization approaches for tailoring stiffness in compliant systems entails the modification of other mechanical characteristics such as weight reduction, stiffness ratios, stress relief, and motion range. Therefore, a brief discussion on how these characteristics are modified as the stiffness is modified is presented here.
Weight reduction Stiffness reduction via Cellular Materials and Topology Optimization approaches implies removing material. Therefore, a reduction in stiffness brings with it a weight reduction. The key aspect relies on a controlled removing of material (or distribution of material), so that this reduction in overall weight can be related to the resulting mechanical properties. A random removal of material will in turn reduce stiffness and weight; however, the control of the resulting mechanics becomes a trail and error process (non-systematic).
Motion range The motion range of a compliant system is related to the elastic displacement. This scales with the elastic modulus of the structure, i.e., its stiffness. Therefore, any attempt to tailor the stiffness of the cellular material or optimize the topology will lead to modifications in the motion range that can be achieved. Another aspect that can come with the incorporation of any of the two approaches is the anisotropy, as the proposed or resulting geometries may be directionally dependent. Structures made of anisotropic materials can have displacements in directions apart from the those aligned to the loading axis.
Stress relief The use of Cellular Materials and Topology Optimization approaches for stiffness reduction often leads to geometries with sharp corners, thus incrementing the stress level. However, by controlling the topology of a cellular material, it is possible to tune the localization and magnitude of the stress [100] in a structure. Moreover, certain cellular structures such as the contact-aided patterns have been used to provide stress relief [44][45][46], while others such as the auxetic patterns allow a better stress distribution due to their high flexibility [119].
Stiffness ratios Stiffness ratios allow one to compare the stiffness levels in different directions. A flexure element commonly needs great flexibility in certain directions while maintaining high levels of stiffness in the others, i.e., its stiffness ratios must be large. The use of Cellular Materials and Topology Optimization approaches for stiffness reduction in certain direction generally produce stiffness reductions also in the others. However, the stiffness reduction in the desired direction is significantly higher than in the non-desired directions, thus incrementing the stiffness ratios [22]. Moreover, the stiffness ratios can be maximized in the Topology Optimization process [52].

Conclusions
Cellular Materials and Topology Optimization can be considered as approaches to tailor properties by selectively removing material, e.g., stiffness, in a structure eliminating the need of additional elements and maintaining the overall dimensions of the piece. These approaches appear as a possible solution to deal with inaccuracies and limitations of compliant systems. In this work, these two were reviewed regarding their used in tailoring the stiffness of flexure elements. Interest in both has grown in recent years due to the increasing availability of additive manufacturing processes. These allow the fabrication of complex geometries with a controlled distribution of matter, comparable to the ones resulting from these approaches. It was then relevant to overview both approaches, so that the users and readers can compare and contrast prior selecting one of them.
Regarding Cellular Materials, stretch-dominated structures are preferred when maximizing stiffness is the objective, while bending-dominated structures are the best option for compliance applications. Two-dimensional bending-dominated lattice arrangements have been more frequently used in compliance applications. The use of bending-dominated 3D closed-cell structures has been studied recently in energy absorption applications, demonstrating their capacity to be used for compliance purposes. Therefore, tailoring stiffness of flexures via Cellular Materials using 3D and closed-cell structures appears to be a promising option, due to the design advantages such as high off-axis stiffness and greater load bearing capacity.
Auxetic architectures are the most widely used in compliance applications because of their high deformations under relatively low stresses. However, the use of these type of architectures in flexure elements has not been explored. In addition, cellular structures with a zero-Poisson ratio represent another type of non-explored architecture. This type of metastructure is characterized for providing deformation in one direction while maintaining unchanged its dimensions in the others. This particular deformation-mechanism could be useful in the design of flexure elements and compliant mechanisms to avoid parasitic motions due to non-desire displacements.
Numerous literature has focused on flexure joints for rotational motion, while very few works on joints for other types of motion were encountered. Research has been focused almost entirely on SARFJ leaving unexplored other types of joints such as prismatic or cylindric. The wide variety of Cellular Materials may allow the creation of unexplored compliant joints that commonly need more than one flexure to be conformed.
On the other hand, the use of Topology Optimization allows finding the best option of material distribution, and with the use of computational power, this could be done relatively fast. However, it often leads to a complex architecture that requires a postprocessing stage to obtain the final design. Furthermore, the resulting design could be difficult to parametrize and handle analytically, compromising the prediction of the behavior under any other conditions different from the ones used in formulation of the optimization process. As opposed to this, the use of Cellular Materials allows an accurate structure-property relation. Still, the potential tailorability is limited to the Cellular Material topologies. Therefore, combinations of Cellular Materials and Topology Optimization, i.e., TOAUC, MTOCM, and TOFCM, must be explored in the design of flexure elements making used of the advantages offered by both approaches.
The use of Topology Optimization in the design of flexure joints has been studied more than that of Cellular Materials. Looking at the evolution of the Topology Optimization of flexure joints, it is possible to see that the tendency is to create designs with distributed compliance through the implementation of repeatable members in their geometrical sapce, avoiding the concentration of deformation in a single element. Therefore, Topology Optimization can be used as a first stage in the design of flexure elements, after that it can be followed by adjusting the material distribution by means of Cellular Material topology, therefore having a controlled structure-property relation and perhaps a simplification of the design.
Finally, further trends in research in this area and recommendations to extend the study of tailoring stiffness of flexure elements via Cellular Materials and Topology Optimization are provided based on the available literature reviewed here.

1.
Due to their advantages in tailoring the mechanical properties, exploring the use of closed-cell structures, 3D arrays, auxetic, and zero-Poisson architectures in tailoring the stiffness of flexure elements could lead to benefits compliant systems.

2.
As flexures base their functioning on the elastic response of their elements, other material properties, i.e., Poisson's ratio, could complicate their use to build up other kinematic pairs, e.g., cylindrical. However, a correct selection of cellular pattern could be designed so that it restricts the Poisson's ratio to 0. This opens up a new range of possibilities in compliant systems applications.

3.
Combining both Cellular Materials and Topology Optimization approaches (TOAUC, MTOCM, and TOFCM) to expand the possibilities in tailoring stiffness of flexure elements. This combination leads to obtain the best out of each. An example of the use of these two could be by implementing Topology Optimization at the conceptual design stage for the proposal of a flexure element. Then, employ a Cellular Material to adjust the stiffness by modifying the parameters that define the structure. 4.
Analytical parameterization of the Cellular Material employed can result in an accurate and precise tracking of the stiffness. This allows a controlled tailoring of stiffness in compliant systems.

5.
Another aspect to keep in mind is that lattice honeycombs might have in-plane anisotropy. Therefore, when used in flexures, the proper selection of orientation becomes crucial. Failure in considering this could lead to deformations in non-desired directions, i.e., parasitic motions. Funding: The authors would like to acknowledge the financial support of Tecnologico de Monterrey, in the production of this work.

Acknowledgments:
The authors would like to acknowledge support from CONACyT for PhD studies of first author (scholarship number 863078) and nTopology for providing access to research and academic licenses.

Conflicts of Interest:
The authors declare no conflicts of interest.