-
Triclinic metamaterials by tristable origami with reprogrammable frustration
Authors:
Ke Liu,
Phanisri P. Pratapa,
Diego Misseroni,
Tomohiro Tachi,
Glaucio H. Paulino
Abstract:
Geometrical frustration induced anisotropy and inhomogeneity are explored to achieve unique properties of metamaterials that set them apart from conventional materials. According to Neumann's principle, to achieve anisotropic responses, the material unit cell should possess less symmetry. Based on such guidelines, we present a triclinic metamaterial system of minimal symmetry, which is originated…
▽ More
Geometrical frustration induced anisotropy and inhomogeneity are explored to achieve unique properties of metamaterials that set them apart from conventional materials. According to Neumann's principle, to achieve anisotropic responses, the material unit cell should possess less symmetry. Based on such guidelines, we present a triclinic metamaterial system of minimal symmetry, which is originated from a Trimorph origami pattern with a simple and insightful geometry: a basic unit cell with four tilted panels and four corresponding creases. The intrinsic geometry of the Trimorph origami, with its changing tilting angles, dictates a folding motion that varies the primitive vectors of the unit cell, couples the shear and normal strains of its extrinsic bulk, and leads to an unusual Poisson's effect. Such effect, associated to reversible auxeticity in the changing triclinic frame, is observed experimentally, and predicted theoretically by elegant math formulae. The nonlinearities of the folding motions allow the unit cell to display three robust stable states, connected through snapping instabilities. When the tristable unit cells are tessellated, phenomena that resembles linear and point defects emerge as a result of geometric frustration. The frustration is reprogrammable into distinct stable and inhomogeneous states by arbitrarily selecting the location of a single or multiple point defects. The Trimorph origami demonstrates the possibility of creating origami metamaterials with symmetries that were hitherto non-existent, leading to triclinic metamaterials with tunable anisotropy for potential applications such as wave propagation control and compliant microrobots.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Experimental realization of tunable Poisson's ratio in deployable origami metamaterials
Authors:
Diego Misseroni,
Phanisri P. Pratapa,
Ke Liu,
Glaucio H. Paulino
Abstract:
Origami metamaterials are known to display highly tunable Poisson's ratio values depending on their folded state. Most studies on the Poisson effects in deployable origami tessellations are restricted to theory and simulation. Experimental realization of the desired Poisson effects in origami metamaterials requires special attention to the boundary conditions to enable deployable nonlinear deforma…
▽ More
Origami metamaterials are known to display highly tunable Poisson's ratio values depending on their folded state. Most studies on the Poisson effects in deployable origami tessellations are restricted to theory and simulation. Experimental realization of the desired Poisson effects in origami metamaterials requires special attention to the boundary conditions to enable deployable nonlinear deformations that give rise to tunability. In this work, we present a novel experimental setup suitable to study the Poisson effects in 2D origami tessellations that undergo simultaneous deformations in both the applied and transverse directions. The setup comprises a gripping mechanism, which we call a Saint-Venant fixture, to eliminate Saint-Venant end effects during uniaxial testing experiment. Using this setup, we conduct Poisson's ratio measurements of the Morph origami pattern whose configuration space combines features of the Miura-ori and Eggbox parent patterns. We experimentally observe the Poisson's ratio sign switching capability of the Morph pattern, along with its ability to display either completely positive or completely negative values of Poisson's ratio by virtue of topological transformations. To demonstrate the versatility of the novel setup we also perform experiments on the standard Miura-ori and the standard Eggbox patterns. Our results demonstrate the agreement between the theory, the simulations, and the experiments on the Poisson's ratio measurement and its tunability in origami metamaterials. The proposed experimental technique can be adopted for investigating other tunable properties of origami metamaterials in static and in dynamic regimes, such as finite-strain Poisson's ratios, elastic thermal expansion, and wave propagation control.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
Discrete symmetries control mechanical response in parallelogram-based origami
Authors:
James McInerney,
Glaucio H. Paulino,
D. Zeb Rocklin
Abstract:
Geometric compatibility constraints dictate the mechanical response of soft systems that can be utilized for the design of mechanical metamaterials such as the negative Poisson ratio Miura-ori origami crease pattern. Here, we develop a formalism for linear compatibility that enables explicit investigation of the interplay between geometric symmetries and functionality in origami crease patterns. W…
▽ More
Geometric compatibility constraints dictate the mechanical response of soft systems that can be utilized for the design of mechanical metamaterials such as the negative Poisson ratio Miura-ori origami crease pattern. Here, we develop a formalism for linear compatibility that enables explicit investigation of the interplay between geometric symmetries and functionality in origami crease patterns. We apply this formalism to a particular class of periodic crease patterns with unit cells composed of four arbitrary parallelogram faces and establish that their mechanical response is characterized by an anticommuting symmetry. In particular, we show that the modes are eigenstates of this symmetry operator and that these modes are simultaneously diagonalizable with the symmetric strain operator and the antisymmetric curvature operator. This feature reveals that the anticommuting symmetry defines an equivalence class of crease pattern geometries which possess equal and opposite in-plane and out-of-plane Poisson's ratios.
△ Less
Submitted 10 February, 2022; v1 submitted 12 August, 2021;
originally announced August 2021.
-
Speeding up Computational Morphogenesis with Online Neural Synthetic Gradients
Authors:
Yuyu Zhang,
Heng Chi,
Binghong Chen,
Tsz Ling Elaine Tang,
Lucia Mirabella,
Le Song,
Glaucio H. Paulino
Abstract:
A wide range of modern science and engineering applications are formulated as optimization problems with a system of partial differential equations (PDEs) as constraints. These PDE-constrained optimization problems are typically solved in a standard discretize-then-optimize approach. In many industry applications that require high-resolution solutions, the discretized constraints can easily have m…
▽ More
A wide range of modern science and engineering applications are formulated as optimization problems with a system of partial differential equations (PDEs) as constraints. These PDE-constrained optimization problems are typically solved in a standard discretize-then-optimize approach. In many industry applications that require high-resolution solutions, the discretized constraints can easily have millions or even billions of variables, making it very slow for the standard iterative optimizer to solve the exact gradients. In this work, we propose a general framework to speed up PDE-constrained optimization using online neural synthetic gradients (ONSG) with a novel two-scale optimization scheme. We successfully apply our ONSG framework to computational morphogenesis, a representative and challenging class of PDE-constrained optimization problems. Extensive experiments have demonstrated that our method can significantly speed up computational morphogenesis (also known as topology optimization), and meanwhile maintain the quality of final solution compared to the standard optimizer. On a large-scale 3D optimal design problem with around 1,400,000 design variables, our method achieves up to 7.5x speedup while producing optimized designs with comparable objectives.
△ Less
Submitted 26 April, 2021; v1 submitted 25 April, 2021;
originally announced April 2021.
-
The role of nonlinearities in topological protection: using magnetically coupled fidget spinners
Authors:
Javier Vila,
Glaucio H. Paulino,
Massimo Ruzzene
Abstract:
We investigate and experimentally observe the existence of topologically protected interface modes in a one-dimensional mechanical lattice, and we report on the effect of nonlinearities on topological protection. The lattice consists of a one-dimensional array of spinners with nearest neighbor coupling resulting from magnetic interactions. The distance between the spinners is spatially modulated t…
▽ More
We investigate and experimentally observe the existence of topologically protected interface modes in a one-dimensional mechanical lattice, and we report on the effect of nonlinearities on topological protection. The lattice consists of a one-dimensional array of spinners with nearest neighbor coupling resulting from magnetic interactions. The distance between the spinners is spatially modulated to obtain a diatomic configuration, and to produce a non-trivial interface by breaking spatial inversion symmetry. For small amplitudes of motion, the interactions are approximately linear, and the system supports topologically protected interface modes at frequencies inside the bulk bandgaps of the lattice. Nonlinearities induced by increasing amplitude of motion cause the interface modes to shift outside the bandgaps and merge with the bulk bands. The resulting edge-to-bulk transition causes the extinction of the topologically protected interface mode and extends it to the entire length of the chain. Such transition is predicted by analytical calculations and verified by experimental observations. The paper thus investigates the existence of topologically protected interface modes obtained through broken spatial inversion symmetry, and documents their lack of robustness in the presence of nonlinearities.
△ Less
Submitted 29 October, 2018; v1 submitted 10 October, 2018;
originally announced October 2018.
-
Stochastic Sampling for Structural Topology Optimization with Many Load Cases: Density-Based and Ground Structure Approaches
Authors:
Xiaojia Zhang,
Eric de Sturler,
Glaucio H. Paulino
Abstract:
We propose an efficient probabilistic method to solve a deterministic problem -- we present a randomized optimization approach that drastically reduces the enormous computational cost of optimizing designs under many load cases for both continuum and truss topology optimization. Practical structural designs by topology optimization typically involve many load cases, possibly hundreds or more. The…
▽ More
We propose an efficient probabilistic method to solve a deterministic problem -- we present a randomized optimization approach that drastically reduces the enormous computational cost of optimizing designs under many load cases for both continuum and truss topology optimization. Practical structural designs by topology optimization typically involve many load cases, possibly hundreds or more. The optimal design minimizes a, possibly weighted, average of the compliance under each load case (or some other objective). This means that in each optimization step a large finite element problem must be solved for each load case, leading to an enormous computational effort. On the contrary, the proposed randomized optimization method with stochastic sampling requires the solution of only a few (e.g., 5 or 6) finite element problems (large linear systems) per optimization step. Based on simulated annealing, we introduce a damping scheme for the randomized approach. Through numerical examples in two and three dimensions, we demonstrate that the stochastic algorithm drastically reduces computational cost to obtain similar final topologies and results (e.g., compliance) compared with the standard algorithms. The results indicate that the damping scheme is effective and leads to rapid convergence of the proposed algorithm.
△ Less
Submitted 10 September, 2016;
originally announced September 2016.
-
Unravelling Metamaterial Properties in Zigzag-base Folded Sheets
Authors:
Maryam Eidini,
Glaucio H. Paulino
Abstract:
Creating complex spatial objects from a flat sheet of material using origami folding techniques has attracted attention in science and engineering. In the present work, we employ geometric properties of partially folded zigzag strips to better describe the kinematics of the known zigzag/herringbone-base folded sheet metamaterials such as the Miura-ori. Inspired by the kinematics of a one-degree of…
▽ More
Creating complex spatial objects from a flat sheet of material using origami folding techniques has attracted attention in science and engineering. In the present work, we employ geometric properties of partially folded zigzag strips to better describe the kinematics of the known zigzag/herringbone-base folded sheet metamaterials such as the Miura-ori. Inspired by the kinematics of a one-degree of freedom zigzag strip, we introduce a class of cellular folded sheet mechanical metamaterials comprising different scales of zigzag strips in which the class of the patterns combines origami folding techniques with kirigami. Employing both analytical and numerical models, we study the key mechanical properties of the folded materials. Particularly, we show that, depending on the geometry, these materials exhibit both negative and positive in-plane Poisson's ratio. By expanding the design space of the Miura-ori, our class of patterns is potentially appropriate for a wide range of applications, from mechanical metamaterials to deployable structures at both small and large scales.
△ Less
Submitted 3 June, 2015; v1 submitted 20 February, 2015;
originally announced February 2015.
-
Topology Optimization Using Polytopes
Authors:
Arun L. Gain,
Glaucio H. Paulino,
Leonardo Duarte,
Ivan F. M. Menezes
Abstract:
Meshing complex engineering domains is a challenging task. Arbitrary polyhedral meshes can provide the much needed flexibility in automated discretization of such domains. The geometric property of the polyhedral meshes such as the unstructured nature and the facial connectivity between elements makes them specially attractive for topology optimization applications. Numerical anomalies in designs…
▽ More
Meshing complex engineering domains is a challenging task. Arbitrary polyhedral meshes can provide the much needed flexibility in automated discretization of such domains. The geometric property of the polyhedral meshes such as the unstructured nature and the facial connectivity between elements makes them specially attractive for topology optimization applications. Numerical anomalies in designs such as the single node connections and checkerboard pattern, which are difficult to manufacture physically, are naturally alleviated with polyhedrons. Special interpolants such as Wachspress, mean value coordinates, maximum entropy shape functions are available to handle arbitrary shaped elements. But the finite elements approaches based on these shape functions face some challenges such as accurate and efficient computation of the shape functions and their derivatives for the numerical evaluation of the weak form integrals. In the current work, we solve the governing three-dimensional elasticity state equation using a Virtual Element Method (VEM) approach. The main characteristic difference between VEM and standard finite element methods (FEM) is that in VEM the canonical basis functions are not constructed explicitly. Rather the stiffness matrix is computed directly utilizing a projection map which extracts the linear component of the deformation. Such a construction guarantees the satisfaction of the patch test (used by engineers as an indicator of optimal convergence of numerical solutions under mesh refinement). Finally, the computations reduce to the evaluation of matrices which contain purely geometric surface facet quantities. The present work focuses on the first-order VEM in which the degrees of freedom associated with the vertices. Utilizing polyhedral elements for topology optimization, we show that the mesh bias in the member orientation is alleviated.
△ Less
Submitted 25 December, 2013;
originally announced December 2013.
-
On the Virtual Element Method for Three-Dimensional Elasticity Problems on Arbitrary Polyhedral Meshes
Authors:
Arun L. Gain,
Cameron Talischi,
Glaucio H. Paulino
Abstract:
We explore the recently-proposed Virtual Element Method (VEM) for numerical solution of boundary value problems on arbitrary polyhedral meshes. More specifically, we focus on the elasticity equations in three-dimensions and elaborate upon the key concepts underlying the first-order VEM. While the point of departure is a conforming Galerkin framework, the distinguishing feature of VEM is that it do…
▽ More
We explore the recently-proposed Virtual Element Method (VEM) for numerical solution of boundary value problems on arbitrary polyhedral meshes. More specifically, we focus on the elasticity equations in three-dimensions and elaborate upon the key concepts underlying the first-order VEM. While the point of departure is a conforming Galerkin framework, the distinguishing feature of VEM is that it does not require an explicit computation of the trial and test spaces, thereby circumventing a barrier to standard finite element discretizations on arbitrary grids. At the heart of the method is a particular kinematic decomposition of element deformation states which, in turn, leads to a corresponding decomposition of strain energy. By capturing the energy of linear deformations exactly, one can guarantee satisfaction of the engineering patch test and optimal convergence of numerical solutions. The decomposition itself is enabled by local projection maps that appropriately extract the rigid body motion and constant strain components of the deformation. As we show, computing these projection maps and subsequently the local stiffness matrices, in practice, reduces to the computation of purely geometric quantities. In addition to discussing aspects of implementation of the method, we present several numerical studies in order to verify convergence of the VEM and evaluate its performance for various types of meshes.
△ Less
Submitted 4 November, 2013;
originally announced November 2013.
-
A consistent operator splitting algorithm and a two-metric variant: Application to topology optimization
Authors:
Cameron Talischi,
Glaucio H. Paulino
Abstract:
In this work, we explore the use of operator splitting algorithms for solving regularized structural topology optimization problems. The context is the classical structural design problems (e.g., compliance minimization and compliant mechanism design), parameterized by means of density functions, whose ill-posendess is addressed by introducing a Tikhonov regularization term. The proposed forward-b…
▽ More
In this work, we explore the use of operator splitting algorithms for solving regularized structural topology optimization problems. The context is the classical structural design problems (e.g., compliance minimization and compliant mechanism design), parameterized by means of density functions, whose ill-posendess is addressed by introducing a Tikhonov regularization term. The proposed forward-backward splitting algorithm treats the constituent terms of the cost functional separately which allows suitable approximations of the structural objective. We will show that one such approximation, inspired by the optimality criteria algorithm and reciprocal expansions, improves the convergence characteristics and leads to an update scheme that resembles the well-known heuristic sensitivity filtering method. We also discuss a two-metric variant of the splitting algorithm that removes the computational overhead associated with bound constraints on the density field without compromising convergence and quality of optimal solutions. We present several numerical results and investigate the influence of various algorithmic parameters.
△ Less
Submitted 18 July, 2013;
originally announced July 2013.
-
Addressing Integration Error for Polygonal Finite Elements Through Polynomial Projections: A Patch Test Connection
Authors:
Cameron Talischi,
Glaucio H. Paulino
Abstract:
Polygonal finite elements generally do not pass the patch test as a result of quadrature error in the evaluation of weak form integrals. In this work, we examine the consequences of lack of polynomial consistency and show that it can lead to a deterioration of convergence of the finite element solutions. We propose a general remedy, inspired by techniques in the recent literature of mimetic finite…
▽ More
Polygonal finite elements generally do not pass the patch test as a result of quadrature error in the evaluation of weak form integrals. In this work, we examine the consequences of lack of polynomial consistency and show that it can lead to a deterioration of convergence of the finite element solutions. We propose a general remedy, inspired by techniques in the recent literature of mimetic finite differences, for restoring consistency and thereby ensuring the satisfaction of the patch test and recovering optimal rates of convergence. The proposed approach, based on polynomial projections of the basis functions, allows for the use of moderate number of integration points and brings the computational cost of polygonal finite elements closer to that of the commonly used linear triangles and bilinear quadrilaterals. Numerical studies of a two-dimensional scalar diffusion problem accompany the theoretical considerations.
△ Less
Submitted 16 July, 2013;
originally announced July 2013.
-
Dynamic Adaptive Mesh Refinement for Topology Optimization
Authors:
Shun Wang,
Eric de Sturler,
Glaucio H. Paulino
Abstract:
We present an improved method for topology optimization with both adaptive mesh refinement and derefinement. Since the total volume fraction in topology optimization is usually modest, after a few initial iterations the domain of computation is largely void. Hence, it is inefficient to have many small elements, in such regions, that contribute significantly to the overall computational cost but co…
▽ More
We present an improved method for topology optimization with both adaptive mesh refinement and derefinement. Since the total volume fraction in topology optimization is usually modest, after a few initial iterations the domain of computation is largely void. Hence, it is inefficient to have many small elements, in such regions, that contribute significantly to the overall computational cost but contribute little to the accuracy of computation and design. At the same time, we want high spatial resolution for accurate three-dimensional designs to avoid postprocessing or interpretation as much as possible. Dynamic adaptive mesh refinement (AMR) offers the possibility to balance these two requirements. We discuss requirements on AMR for topology optimization and the algorithmic features to implement them. The numerical design problems demonstrate (1) that our AMR strategy for topology optimization leads to designs that are equivalent to optimal designs on uniform meshes, (2) how AMR strategies that do not satisfy the postulated requirements may lead to suboptimal designs, and (3) that our AMR strategy significantly reduces the time to compute optimal designs.
△ Less
Submitted 25 September, 2010;
originally announced September 2010.
-
Integral Equations with Hypersingular Kernels -- Theory and Applications to Fracture Mechanics
Authors:
Youn-Sha Chan,
Albert C. Fannjiang,
Glaucio H. Paulino
Abstract:
Hypersingular integrals of the type $I_α(T_n,m,r) = \int_{-1}^{1} \hpsngAbs \frac{T_n(s)(1-s^2)^{m-{1/2}}}{(s-r)^α}ds |r|<1$ and $I_α(U_n,m,r) = \int_{-1}^{1} \hpsngAbs \frac{U_n(s)(1-s^2)^{m-{1/2}}}{(s-r)^α}ds |r|<1$ are investigated for general integers $α$ (positive) and $m$ (non-negative), where $T_n(s)$ and $U_n(s)$ are the Tchebyshev polynomials of the 1st and 2nd kinds, respectively. Exac…
▽ More
Hypersingular integrals of the type $I_α(T_n,m,r) = \int_{-1}^{1} \hpsngAbs \frac{T_n(s)(1-s^2)^{m-{1/2}}}{(s-r)^α}ds |r|<1$ and $I_α(U_n,m,r) = \int_{-1}^{1} \hpsngAbs \frac{U_n(s)(1-s^2)^{m-{1/2}}}{(s-r)^α}ds |r|<1$ are investigated for general integers $α$ (positive) and $m$ (non-negative), where $T_n(s)$ and $U_n(s)$ are the Tchebyshev polynomials of the 1st and 2nd kinds, respectively. Exact formulas are derived for the cases $α= 1, 2, 3, 4$ and $m = 0, 1, 2, 3$; most of them corresponding to new solutions derived in this paper. Moreover, a systematic approach for evaluating these integrals when $α> 4$ and $m>3$ is provided. The integrals are also evaluated as $|r|>1$ in order to calculate stress intensity factors (SIFs). Examples involving crack problems are given and discussed with emphasis on the linkage between mathematics and mechanics of fracture. The examples include classical linear elastic fracture mechanics (LEFM), functionally graded materials (FGM), and gradient elasticity theory. An appendix, with closed form solutions for a broad class of integrals, supplements the paper.
△ Less
Submitted 9 June, 1999;
originally announced June 1999.