Source author record

Glaucio H. Paulino

Glaucio H. Paulino appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Discrete symmetries control mechanical response in parallelogram-based origami

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.

preprint2015arXiv

Unravelling Metamaterial Properties in Zigzag-base Folded Sheets

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.

preprint2013arXiv

A consistent operator splitting algorithm and a two-metric variant: Application to topology optimization

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.

preprint2013arXiv

Addressing Integration Error for Polygonal Finite Elements Through Polynomial Projections: A Patch Test Connection

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.

preprint2013arXiv

On the Virtual Element Method for Three-Dimensional Elasticity Problems on Arbitrary Polyhedral Meshes

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.

preprint2013arXiv

Topology Optimization Using Polytopes

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.

preprint2010arXiv

Dynamic Adaptive Mesh Refinement for Topology Optimization

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.