Research Article  Open Access
Effects of Morphology and Topology on the Effective Stiffness of Chiral Cellular Materials in the Transverse Plane
Abstract
The present study investigates the influence of topology and morphology on the effective stiffness of chiral cellular materials in the transverse plane by means of a homogenization method. For this purpose, finite element models of representative volume elements for regular hexagonal and hexagonalchiral configurations are used and simulations are conducted to quantify how cell topology—that is, chirality inside the cell—and cell wall slenderness affect the effective stiffness. Closed form solutions for regular hexagonal square and triangular RVEs provided in the literature are then taken as a basis for model validation. The results indicate that there are drastic differences between regular hexagonal and hexagonalchiral configurations, which can be explained in terms of deformation mechanism transformations between bending and stretching. The investigations also reveal the positive impact of cell wall slenderness on stiffness due to volumetric increase in the cell wall material resisting the deformation.
1. Introduction
In recent years, with growing interest in energy efficient, robust, and lightweight material systems, cellular materials have found their way into numerous engineering applications. They are nowadays used as highstiffness sandwich panels, energy absorbers, thermal insulators, and vibration dampers. Due to their low densities and high potential as being tailored materials in terms of effective stiffness and strength properties, there are ongoing research and development activities on these materials [1–3].
Dependence of these effective properties on cell wall geometrical, topological, and mechanical properties has been investigated in the framework of computational homogenization methods [4, 5]. The framework is used to solve the associated boundary value problem with proper boundary conditions on representative volume elements (RVEs) that capture the prominent features of the material [6]. In the case of cellular materials, the term RVE refers to an individual cell or a group of cells, as depicted in Figure 1, that are assumed to represent the effective material properties at the mesoscale.
As a contribution to the current state of the art on cellular material characterization in the framework of computational homogenization, the current study aims at (i) understanding the effect of cell topology on the effective stiffness properties of cellular materials in the transverse plane by using regular hexagonal and hexagonalchiral RVEs shown in Figure 1, (ii) understanding the effect of cell wall slenderness—that is, cell wall thicknesstoheight ratio—on these effective properties, and (iii) comparing the extracted stiffness properties with analytical results given in the literature.
2. Materials and Methods
Computational homogenization has been a useful tool for determining the effective properties of materials that can be described through their microscopic properties. The basic idea behind this methodology is to evaluate the effective properties by solving the boundary value problem with appropriate boundary conditions defined on a representative volume element (RVE), which contains all necessary information [7, 8]. Some of the key advantages include being able to introduce all of the necessary geometrical and constitutive information at appropriate scales of the material behavior and being able to use requisitely detailed yet computationally effective idealization—such as beam and plate models—at the scale of the material cells [9, 10].
In the literature, numerous analytical and computational approaches have been offered to the homogenization problems involving cellular materials (see, e.g., [11–13]). Due to advances in theory and hardware, computational homogenization methods have become more frequently used for material characterization. In the present study, a computational homogenization method is thus used to investigate the effective stiffness properties of cellular materials, which is explained in detail in the following section.
2.1. A FirstOrder StrainDriven Computational Homogenization Method
In the homogenization approach adopted here, the macroscopic strain for is assumed to be uniform over the material and imposed at the RVE boundary as illustrated in Figure 2. Therefore, a firstorder microscopic displacement field decomposition for the RVE can be given asfor which the first term on the righthand side represents the macroscopic displacement contribution and the second one represents the displacement fluctuation field due to heterogeneities within the RVE [4]. Here, represents the position vector with respect to any described origin.
Since the overall body is assumed to be composed of periodic arrays of repeating RVEs, continuity conditions for the displacements and tractions must be satisfied at the boundaries of each neighbouring cell. The first condition can be satisfied by taking the relative positions of the node sets, for example, nodes at the opposite edges on , for which displacement fluctuations have the same value due to the periodicity condition [14]. Hence, (1) can be expanded aswhich forms the basis of the displacement boundary condition for periodic arrays. The traction boundary condition is satisfied with antiperiodicity of traction field in case of existence of traction on the boundaries [6]. However, the current study focuses on the displacement boundary conditions; therefore, only (2) is studied for the periodic boundary conditions.
In order to establish a consistent relationship between different scales, the HillMandel principle can be used such that [15]for which the superscript “m” represents the microscopic or RVE scale. The symbol “:” denotes an inner product, and for secondorder tensors it is defined asBy using Gauss’s theorem, (3) can be rewritten at the RVE boundary aswhere is the microscopic traction vector at . By plugging (1) into (5),which yields toHere, the symbol denotes a dyadic product. The second integrand at the righthand side vanishes in case of periodic boundary conditions, as elaborated in (2). Hence, the macroscopic stress can be expressed as the volume average of the microscopic stress such thatwhere is the total volume of the RVE.
Thereafter, the effective stiffness can be determined by means of the relationship between the macroscopic strain and stress under the linear elasticity assumption. This relationship can be constructed in terms of macroscopic compliance so that [16]Then, the Voigt notation for can be constructed asfor which is taken to be symmetric matrix. Superscripts 1, 2, and 3 represent three different loading conditions fulfilling the minimum requirement to obtain the compliance matrix. Under the assumption of orthotropic material properties,for which and are the effective Young’s moduli, are the shear modulus, and and are Poisson’s ratios with reciprocal relation [17].
2.2. Implementation of Computational Homogenization
In this study, regular hexagonal and hexagonalchiral RVEs are investigated. Therefore, the aforementioned equation sets are applied on matching node sets that form the periodicity conditions. The mapping is illustrated on a regular hexagonal RVE geometry as seen in Figure 3(a), whereas the geometrical parameters of regular hexagonal and hexagonalchiral RVEs are depicted in Figure 3(b).
(a)
(b)
3. Results and Discussion
3.1. Design of Experiments
In the present study, cell walls are taken to be transversely isotropic under the assumption of uniform distribution of cell wall material in the plane of interest. Mechanical parameters of the cell walls are selected as Pa and , which refer to typical cell wall Young’s modulus and Poisson’s ratio value, respectively. In relation to (9) and (10), three different inplane loading conditions—that is, tensile load along the  and directions and the shear load as listed in Table 1 and illustrated in Figure 4—are used to compute the effective stiffness properties and understand the effect of cell wall slenderness , for which , is cell wall thickness, and refers to cell wall height (as depicted in Figure 3), on the computed properties. The reason for taking the cell wall slenderness as is to be able to compare the numerical results with the closed form solutions provided for double cell walls in the literature [5, 9]. The slenderness is taken to be in the range so that with increment of whereas hub diameter is taken as for the chiral configuration illustrated in Figure 3(b).

(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
3.2. Results
For understanding how cell wall slenderness (namely, cell wall thicknesstoheight ratio) and cell topology (namely, chirality introduced inside the cell) affect the effective stiffness properties, simulation experiments are carried out on the RVEs. For this purpose, Abaqus/CAE is used, in which the periodic boundary conditions can be implemented. Thereafter, the solution domain is discretized into sets of finite elements—here, 6noded quadratic triangular CPS6 provided in Abaqus/CAE. The results are compared with closed form solutions provided in Table 2.
3.2.1. Effect of Cell Wall Slenderness via Hexagonal RVE Simulations
Computational and analytical investigations on hexagonal RVEs are conducted so as to understand the effect of cell wall slenderness on the effective stiffness properties in the transverse plane. In case of regular hexagonal RVEs, cell walls are assumed to behave as mechanisms rather than rigid structures where the bending moments are the dominant deformation mechanisms under both uniaxial and shear loadings [18].
The comparison graphs presented in Figure 5 indicate that there is a positive effect of on the effective Young’s moduli and and shear modulus . However, with increase in , there is a decrease in Poisson’s ratio values and . In addition to the comparative studies, it is also deduced that Surveys I and II for regular hexagonal unit cells are in accordance with the computational homogenization results. This validates the current model and can hence be implemented on the proposed hexagonalchiral configuration, which is elaborated in the following section.
3.2.2. Effect of Topology
For understanding the topological effects on the effective inplane elastic parameters, closed form equations for regular hexagonal, square, and equilateral triangular unit cells listed in Table 2 are compared with the current computational investigations on hexagonalchiral RVEs. The results plotted for effective moduli , , and in Figure 6 show that stiffness is favored as the topology changes from hexagonal to triangular RVE. As stated in [3, 11], the main reason is the increasing dominancy of stretching over bending with this type of topological change. As also seen in Figure 6, the proposed hexagonalchiral configuration results in much higher effective Young’s moduli compared to regular hexagonal RVE. For instance, the ratio between the effective Young’s moduli of two configurations for is ~100, favoring the use of chirality. According to Figure 6, the hexagonalchiral configuration yields the highest shear moduli among all of the investigated topologies. It is a promising result especially for the materials subject to shear loads—for example, packaging products and cargo containers.
Comparison of effective Poisson’s ratios and in Figure 7 indicates that the regular hexagonal RVE has the highest linear elastic response—in terms of stretching or contraction—to the loading in the perpendicular direction. However, this response is lower for the triangular RVE and almost negligible for square one. The proposed hexagonalchiral configuration has and values that are much lower than the regular hexagonal RVE, which are comparable with those for the triangular and square RVEs.
4. Conclusions
In the present study, topological and morphological effects on the effective stiffness of cellular materials in the transverse plane are investigated by means of a firstorder straindriven computational homogenization method. For this purpose, representative volume elements (RVEs) of regular hexagonal, square, triangular, and hexagonalchiral configurations are studied to understand how cell topology (i.e., nominal versus chiral) and cell wall slenderness affect the effective stiffness properties. Closed form solutions for regular hexagonal, square, and triangular RVEs provided in the literature are then taken as the basis for model validation and comparison. The results show that there is a positive effect of cell wall slenderness on the effective elastic parameters due to increase in the cell wall material volume. In addition to this, it is also deduced that there are drastic differences between regular hexagonal and hexagonalchiral configurations, which can be explained in terms of deformation mechanism transformations between bending and stretching due to topological changes. Hence, the introduced chirality can be used as a tailored solution for materials and structures subject to high shear loads.
Disclosure
Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the sponsors.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The first author acknowledges the financial support of Tekniikan edistämissäätiö (TES) through Foundations’ Post Doc Pool. The second author acknowledges support by the National Science Foundation under Grant no. 1031181 (Program Managers Dr. Glaucio Paulino and Dr. Martin Dunn) and by the California Energy Commission’s Energy Innovations Small Grants Program under Grant no. 0903TE (Contract no. 50098014, Program Managers Stacie Ishida and David Chambers).
References
 A. Rafsanjani, Multiscale Poroelastic Model: Bridging the Gap from Cellular to Macroscopic Scale, ETH Zurich, Zurich, Switzerland, 2013.
 R. Harb, E. Taciroglu, and N. Ghoniem, “Partitioning of elastic energy in opencell foams under finite deformations,” Acta Materialia, vol. 61, no. 5, pp. 1454–1468, 2013. View at: Publisher Site  Google Scholar
 A. Karakoç, K. Santaoja, and J. Freund, “Simulation experiments on the effective inplane compliance of the honeycomb materials,” Composite Structures, vol. 96, pp. 312–320, 2013. View at: Publisher Site  Google Scholar
 M. G. D. Geers, V. G. Kouznetsova, and W. A. M. Brekelmans, “Multiscale computational homogenization: trends and challenges,” Journal of Computational and Applied Mathematics, vol. 234, no. 7, pp. 2175–2182, 2010. View at: Publisher Site  Google Scholar
 J. Freund, A. Karakoç, and J. Sjölund, “Computational homogenization of regular cellular material according to classical elasticity,” Mechanics of Materials, vol. 78, pp. 56–65, 2014. View at: Publisher Site  Google Scholar
 V.D. Nguyen, E. Béchet, C. Geuzaine, and L. Noels, “Imposing periodic boundary condition on arbitrary meshes by polynomial interpolation,” Computational Materials Science, vol. 55, pp. 390–406, 2012. View at: Publisher Site  Google Scholar
 S. NematNasser and M. Hori, Micromechanics: Overall Properties of Heterogeneous Materials, NorthHolland Publishing, Amsterdam, The Netherlands, 1999.
 J. Sjölund, A. Karakoç, and J. Freund, “Effect of cell geometry and material properties on wood rigidity,” International Journal of Solids and Structures, vol. 62, pp. 207–216, 2015. View at: Publisher Site  Google Scholar
 L. Gibson and M. Ashby, Cellular Solids: Structure and Properties, Cambridge University Press, Cambridge, UK, 1999.
 A.J. Wang, R. S. Kumar, and D. L. McDowell, “Mechanical behavior of extruded prismatic cellular metals,” Mechanics of Advanced Materials and Structures, vol. 12, no. 3, pp. 185–200, 2005. View at: Publisher Site  Google Scholar
 M. Alkhader and M. Vural, “Mechanical response of cellular solids: role of cellular topology and microstructural irregularity,” International Journal of Engineering Science, vol. 46, no. 10, pp. 1035–1051, 2008. View at: Publisher Site  Google Scholar
 J. Sjölund, A. Karakoç, and J. Freund, “Accuracy of regular wood cell structure model,” Mechanics of Materials, vol. 76, pp. 35–44, 2014. View at: Publisher Site  Google Scholar
 E. Kahle and J. Woodhouse, “The influence of cell geometry on the elasticity of softwood,” Journal of Materials Science, vol. 29, no. 5, pp. 1250–1259, 1994. View at: Publisher Site  Google Scholar
 F. Fritzen and T. Bohlke, “Influence of the type of boundary conditions on the numerical properties of unit cell problems,” Technische Mechanik, vol. 30, pp. 354–363, 2010. View at: Google Scholar
 K. Pham, V. G. Kouznetsova, and M. G. D. Geers, “Transient computational homogenization for heterogeneous materials under dynamic excitation,” Journal of the Mechanics and Physics of Solids, vol. 61, no. 11, pp. 2125–2146, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Karakoç and J. Freund, “A direct simulation method for the effective inplane stiffness of cellular materials,” International Journal of Applied Mechanics, vol. 5, no. 3, Article ID 1350034, 2013. View at: Publisher Site  Google Scholar
 A. Laghdir, Y. Fortin, C. M. De La Cruz, and R. E. Hernández, “Development of a technique to determine the 3D elasticity tensor of wood as applied to drying stress modeling,” Maderas: Ciencia y Tecnologia, vol. 10, no. 1, pp. 35–44, 2008. View at: Google Scholar
 V. S. Deshpande, M. F. Ashby, and N. A. Fleck, “Foam topology: bending versus stretching dominated architectures,” Acta Materialia, vol. 49, no. 6, pp. 1035–1040, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Alp Karakoç and Ertuǧrul Taciroǧlu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.