skip to main content

Title: Voronoi Spaghetti & VoroNoodles: Topologically Interlocked, Space-Filling, Corrugated & Congruent Tiles
In this work, we introduce an approach to model topologically interlocked corrugated bricks that can be assembled in a water-tight manner (space-filling) to design a variety of spatial structures. Our approach takes inspiration from recently developed methods that utilize Voronoi tessellation of spatial domains by using symmetrically arranged Voronoi sites. However, in contrast to these existing methods, we focus our attention on Voronoi sites modeled using helical trajectories, which can provide corrugation and better interlocking. For symmetries, we only use affine transformations based on the Bravais lattice to avoid self-intersections. This methodology naturally results in structures that are both space-filling (owing to Voronoi tessellation) as well as interlocking by corrugation (owing to helical trajectories). The resulting shapes of the bricks appear to be similar to a variety of pasta noodles, thereby inspiring the names, Voronoi Spaghetti and VoroNoodles.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
SA '22: SIGGRAPH Asia 2022 Technical Communications
Page Range / eLocation ID:
1 to 4
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. An approach for modeling topologically interlocked building blocks that can be assembled in a water‐tight manner (space filling) to design a variety of spatial structures is introduced. This approach takes inspiration from recent methods utilizing Voronoi tessellation of spatial domains using symmetrically arranged Voronoi sites. Attention is focused on building blocks that result from helical stacking of planar 2‐honeycombs (i.e., tessellations of the plane with a single prototile) generated through a combination of wallpaper symmetries and Voronoi tessellation. This unique combination gives rise to structures that are both space‐filling (due to Voronoi tessellation) and interlocking (due to helical trajectories). Algorithms are developed to generate two different varieties of helical building blocks, namely, corrugated and smooth. These varieties result naturally from the method of discretization and shape generation and lead to distinct interlocking behavior. In order to study these varieties, finite‐element analyses (FEA) are conducted on different tiles parametrized by 1) the polygonal unit cell determined by the wallpaper symmetry and 2) the parameters of the helical line generating the Voronoi tessellation. Analyses reveal that the new design of the geometry of the building blocks enables strong variation of the engagement force between the blocks.

    more » « less
  2. The recent development of three-dimensional graphic statics (3DGS) has greatly increased the ease of designing complex and efficient spatial funicular structural forms [1]. The reciprocal diagram based 3DGS approaches not only generate highly efficient funicular structures [2], but also result in planarity constraints due to the polyhedron nature of the reciprocal diagrams [3]. Our previous research has shown the feasibility of leveraging this planarity by using planar glass sheets to materialize the 10m-span, double-layer glass bridge [3]. This paper is framed as a proof of concept for the 10m bridge and explores the form-finding, detail configuration, fabrication constraints, and assembly logic by designing and constructing a small-scale bridge prototype with a span of 2.5m. The prototype is designed in a modular approach, where each polyhedral cell of the form is materialized using a hollow glass unit (HGU) (Figure 1a), which can be prefabricated and preassembled, and therefore, greatly simplifies the assembly of the whole bridge. The compression-only form of the prototype is generated using the PolyFrame beta [4] plug-in for Rhinoceros [5]. The form-finding is carried out with a comprehensive consideration of a variety of parameters, including fabrication constraints, assembly ease, construction cost, and practicality. To start the form-finding process, a group of closed convex force polyhedrons is aggregated, controlling the topology of the form diagram and the orientations of the form elements. By manipulating the face tilting angles of the force diagram, the supported edges at the end of the bridge are all made horizontal, reducing the difficulty of the support design. Then, vertex locations and edge lengths of the form diagram are constrained, determining the final dimensions of both the bridge and the cells. After getting the geometry of the bridge, the detail developments are streamlined. Each of the 13 HGUs consists of two flat deck plates and a series of side plates (Figure 1b). To interlock the adjacent cells and prevent possible sliding, a male-female connection mechanism is introduced to the conjoint side plates of the HGUs (Figure 1b). Additionally, to eliminate the direct contact of the glass parts and prevent the stress concentration, two softer transparent materials are involved for connecting purposes. Within each HGU, silicon-based binding agent is used to hold the glass parts together; between the neighboring HGUs, plastic sheets are placed as interface materials (Figure 1b). Figure 1. a) The 2.5m-span small-scale prototype dome, b) Exploded view showing deck plates, side plates, male-female connection, and interface material For the fabrication of the glass parts, 5-axis Waterjet cutting techniques are applied. While the glass sheets for the deck plates can be purchased from the market, the irregular side plates with male-female connections need to be made from kiln-cast glass. In terms of the Waterjet cutting constraints, there is a max cutting angle of 60 degrees from vertical. With respect to this, all the glass parts are examined during the design process to ensure they all satisfy the cutting angle requirements. Aiming to achieve a fast and precise assembly, several assistant techniques are developed. On the local HGU level, assembly connectors are designed and 3D-printed to help locate the glass parts. On the global prototype level, the assembly sequence of the HGUs are simulated to avoid interference. Besides, a labeling system is also established to organize the fabricated parts and guide the entire assembly process. The design and construction of this small-scale prototype provide important information for the future development of the full-scale bridge regarding the interlocking detail design, the fabrication constraints, and assembly logic. The actual structural performance of the prototype awaits further investigation through-loading experiments. 
    more » « less
  3. Summary Conditional density estimation seeks to model the distribution of a response variable conditional on covariates. We propose a Bayesian partition model using logistic Gaussian processes to perform conditional density estimation. The partition takes the form of a Voronoi tessellation and is learned from the data using a reversible jump Markov chain Monte Carlo algorithm. The methodology models data in which the density changes sharply throughout the covariate space, and can be used to determine where important changes in the density occur. The Markov chain Monte Carlo algorithm involves a Laplace approximation on the latent variables of the logistic Gaussian process model which marginalizes the parameters in each partition element, allowing an efficient search of the approximate posterior distribution of the tessellation. The method is consistent when the density is piecewise constant in the covariate space or when the density is Lipschitz continuous with respect to the covariates. In simulation and application to wind turbine data, the model successfully estimates the partition structure and conditional distribution. 
    more » « less
  4. Largely due to superior properties compared to traditional materials, the use of polymer matrix composites (PMC) has been expanding in several industries such as aerospace, transportation, defense, and marine. However, the anisotropy and nonhomogeneity of these structures contribute to the difficulty in evaluating structural integrity; damage sites can occur at multiple locations and length scales and are hard to track over time. This can lead to unpredictable and expensive failure of a safety-critical structure, thus creating a need for non-destructive evaluation (NDE) techniques which can detect and quantify small-scale damage sites and track their progression. Our research group has improved upon classical microwave techniques to address these needs; utilizing a custom device to move a sample within a resonant cavity and create a spatial map of relative permittivity. We capitalize on the inevitable presence of moisture within the polymer network to detect damage. The differing migration inclinations of absorbed water molecules in a pristine versus a damaged composite alters the respective concentrations of the two chemical states of moisture. The greater concentration of free water molecules residing in the damage sites exhibit highly different relative permittivity when compared to the higher ratio of polymer-bound water molecules in the undamaged areas. Currently, the technique has shown the ability to detect impact damage across a range of damage levels and gravimetric moisture contents but is not able to specifically quantify damage extent with regards to impact energy level. The applicability of machine learning (ML) to composite materials is substantial, with uses in areas like manufacturing and design, prediction of structural properties, and damage detection. Using traditional NDE techniques in conjunction with supervised or unsupervised ML has been shown to improve the accuracy, reliability, or efficiency of the existing methods. In this work, we explore the use of a combined unsupervised/supervised ML approach to determine a damage boundary and quantification of single-impact specimens. Dry composite specimens were damaged via drop tower to induce one central impact site of 0, 2, or 3 Joules. After moisture exposure, Entrepreneur Dr, Raleigh, North Carolina 27695, U.S.A. 553 each specimen underwent dielectric mapping, and spatial permittivity maps were created at a variety of gravimetric moisture contents. An unsupervised K-means clustering algorithm was applied to the dielectric data to segment the levels of damage and define a damage boundary. Subsequently, supervised learning was used to quantify damage using features including but not limited to thickness, moisture content, permittivity values of each cluster, and average distance between points in each cluster. A regression model was trained on several samples with impact energy as the predicted variable. Evaluation was then performed based on prediction accuracy for samples in which the impact energies are not known to the model.

    more » « less
  5. SUMMARY The uneven distribution of earthquakes and stations in seismic tomography leads to slower convergence of nonlinear inversions and spatial bias in inversion results. Including dense regional arrays, such as USArray or Hi-Net, in global tomography causes severe convergence and spatial bias problems, against which conventional pre-conditioning schemes are ineffective. To save computational cost and reduce model bias, we propose a new strategy based on a geographical weighting of sources and receivers. Unlike approaches based on ray density or the Voronoi tessellation, this method scales to large full-waveform inversion problems and avoids instabilities at the edges of dense receiver or source clusters. We validate our strategy using a 2-D global waveform inversion test and show that the new weighting scheme leads to a nearly twofold reduction in model error and much faster convergence relative to a conventionally pre-conditioned inversion. We implement this geographical weighting strategy for global adjoint tomography. 
    more » « less