skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Virtual element method (VEM)-based topology optimization: an integrated framework
We present a virtual element method (VEM)-based topology optimization framework using polyhedral elements, which allows for convenient handling of non-Cartesian design domains in three dimensions. We take full advantage of the VEM properties by creating a unified approach in which the VEM is employed in both the structural and the optimization phases. In the structural problem, the VEM is adopted to solve the three-dimensional elasticity equation. Compared to the finite element method, the VEM does not require numerical integration (when linear elements are used) and is less sensitive to degenerated elements (e.g., ones with skinny faces or small edges). In the optimization problem, we introduce a continuous approximation of material densities using the VEM basis functions. When compared to the standard element-wise constant approximation, the continuous approximation enriches the geometrical representation of structural topologies. Through two numerical examples with exact solutions, we verify the convergence and accuracy of both the VEM approximations of the displacement and material density fields. We also present several design examples involving non-Cartesian domains, demonstrating the main features of the proposed VEM-based topology optimization framework. The source code for a MATLAB implementation of the proposed work, named PolyTop3D, is available in the (electronic) Supplementary Material accompanying this publication.  more » « less
Award ID(s):
1663244
PAR ID:
10170741
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Structural and Multidisciplinary Optimization
ISSN:
1615-147X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Topology optimization has been proved to be an efficient tool for structural design. In recent years, the focus of structural topology optimization has been shifting from single material continuum structures to multimaterial and multiscale structures. This paper aims at devising a numerical scheme for designing bionic structures by combining a two-stage parametric level set topology optimization with the conformal mapping method. At the first stage, the macro-structural topology and the effective material properties are optimized simultaneously. At the second stage, another structural topology optimization is carried out to identify the exact layout of the metamaterial at the mesoscale. The achieved structure and metamaterial designs are further synthesized to form a multiscale structure using conformal mapping, which mimics the bionic structures with “orderly chaos” features. In this research, a multi-control-point conformal mapping (MCM) based on Ricci flow is proposed. Compared with conventional conformal mapping with only four control points, the proposed MCM scheme can provide more flexibility and adaptivity in handling complex geometries. To make the effective mechanical properties of the metamaterials invariant after conformal mapping, a variable-thickness structure method is proposed. Three 2D numerical examples using MCM schemes are presented, and their results and performances are compared. The achieved multimaterial multiscale structure models are characterized by the “orderly chaos” features of bionic structures while possessing the desired performance. 
    more » « less
  2. Hsia, KJ; Rogers, JA; Suo, Z; Zhao, X (Ed.)
    Topology optimization algorithms often employ a smooth density function to implicitly represent geometries in a discretized domain. While this implicit representation offers great flexibility to parametrize the optimized geometry, it also leads to a transition region. Previous approaches, such as the Solid Isotropic Material Penalty (SIMP) method, have been proposed to modify the objective function aiming to converge toward integer density values and eliminate this non-physical transition region. However, the iterative nature of topology optimization renders this process computationally demanding, emphasizing the importance of achieving fast convergence. Accelerating convergence without significantly compromising the final solution can be challenging. In this work, we introduce a machine learning approach that leverages the message-passing Graph Neural Network (GNN) to eliminate the non-physical transition zone for the topology optimization problems. By representing the optimized structures as weighted graphs, we introduce a generalized filtering algorithm based on the topology of the spatial discretization. As such, the resultant algorithm can be applied to two- and three-dimensional space for both Cartesian (structured grid) and non-Cartesian discretizations (e.g. polygon finite element). The numerical experiments indicate that applying this filter throughout the optimization process may avoid excessive iterations and enable a more efficient optimization procedure. 
    more » « less
  3. Summary This study focuses on the topology optimization framework for the design of multimaterial dissipative systems at finite strains. The overall goal is to combine a soft viscoelastic material with a stiff hyperelastic material for realizing optimal structural designs with tailored damping and stiffness characteristics. To this end, several challenges associated with incorporating finite‐deformation viscoelastic‐hyperelastic materials in a multimaterial design framework are addressed. This includes consideration of a thermodynamically consistent finite‐strain viscoelasticity model for simulating energy dissipation together with F‐bar finite elements for handling material incompressibility. Moreover, an effective multimaterial interpolation scheme is proposed, which preserves the physics of material mixtures in the context of density‐based topology optimization. A numerically accurate analytical design sensitivity calculation is also presented using a path‐dependent adjoint method. Furthermore, both prescribed‐load and prescribed‐displacement boundary conditions are considered in the optimization formulations, together with various strategies for controlling stiffness. As demonstrated by the numerical examples, the use of the stiffer hyperelastic material phase in a design not only improves stiffness but also increases energy dissipation capacity. Moreover, with the finite‐deformation theory, the effect of the loading magnitude on the optimized designs can be observed. 
    more » « less
  4. The level set method has been widely applied in topology optimization of mechanical structures, primarily for linear materials, but its application to nonlinear hyperelastic materials, particularly for compliant mechanisms, remains largely unexplored. This paper addresses this gap by developing a comprehensive level set-based topology optimization framework specifically for designing compliant mechanisms using neo-Hookean hyperelastic materials. A key advantage of hyperelastic materials is their ability to undergo large, reversible deformations, making them well-suited for soft robotics and biomedical applications. However, existing nonlinear topology optimization studies using the level set method mainly focus on stiffness optimization and often rely on linear results as preliminary approximations. Our framework rigorously derives the shape sensitivity analysis using the adjoint method, including crucial higher-order displacement gradient terms often neglected in simplified approaches. By retaining these terms, we achieve more accurate boundary evolution during optimization, leading to improved convergence behavior and more effective structural designs. The proposed approach is first validated with a mean compliance problem as a benchmark, demonstrating its ability to generate optimized structural configurations while addressing the nonlinear behavior of hyperelastic materials. Subsequently, we extend the method to design a displacement inverter compliant mechanism that fully exploits the advantages of hyperelastic materials in achieving controlled large deformations. The resulting designs feature smooth boundaries and clear structural features that effectively leverage the material's nonlinear properties. This work provides a robust foundation for designing advanced compliant mechanisms with large deformation capabilities, extending the reach of topology optimization into new application domains where traditional linear approaches are insufficient. The developed methodology is expected to provide a timely solution to computational design for soft robotics, flexible mechanisms, and other emerging technologies that benefit from hyperelastic material properties. 
    more » « less
  5. Topology optimization problems are typically non-convex, and as such, multiple local minima exist. Depending on the initial design, the type of optimization algorithm and the optimization parameters, gradient-based optimizers converge to one of those minima. Unfortunately, these minima can be highly suboptimal, particularly when the structural response is very non-linear or when multiple constraints are present. This issue is more pronounced in the topology optimization of geometric primitives, because the design representation is more compact and restricted than in free-form topology optimization. In this paper, we investigate the use of tunneling in topology optimization to move from a poor local minimum to a better one. The tunneling method used in this work is a gradient-based deterministic method that finds a better minimum than the previous one in a sequential manner. We demonstrate this approach via numerical examples and show that the coupling of the tunneling method with topology optimization leads to better designs. 
    more » « less