skip to main content


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
NSF-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. 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
  2. Abstract

    This article presents an approach for the topology optimization of frame structures composed of nonlinear Timoshenko beam finite elements (FEs) under time‐varying excitation. Material nonlinearity is considered with a nonlinear Timoshenko beam FE model that accounts for distributed plasticity and axial–shear–moment interactions through appropriate hysteretic interpolation functions and a yield/capacity function, respectively. Hysteretic variables for curvature, shear, and axial deformations represent the nonlinearities and evolve according to first‐order nonlinear ordinary differential equations (ODEs). Owing to the first‐order representation, the governing dynamic equilibrium equations, and hysteretic evolution equations can thus be concisely presented as a combined system of first‐order nonlinear ODEs that can be solved using a general ODE solver. This avoids divergence due to an ill‐conditioned stiffness matrix that can commonly occur with Newmark–Newton solution schemes that rely upon linearization. The approach is illustrated for a volume minimization design problem, subject to dynamic excitation where an approximation for the maximum displacement at specified nodes is constrained to a given limit, that is, a drift ratio. The maximum displacement is approximated using the p‐norm, thus facilitating the derivation of the analytical sensitivities for gradient‐based optimization. The proposed approach is demonstrated through several numerical examples for the design of structural frames subjected to sinusoidal base excitation.

     
    more » « less
  3. 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
  4. Purpose The purpose of this paper is as follows: to significantly reduce the computation time (by a factor of 1,000 and more) compared to known numerical techniques for real-world problems with complex interfaces; and to simplify the solution by using trivial unfitted Cartesian meshes (no need in complicated mesh generators for complex geometry). Design/methodology/approach This study extends the recently developed optimal local truncation error method (OLTEM) for the Poisson equation with constant coefficients to a much more general case of discontinuous coefficients that can be applied to domains with different material properties (e.g. different inclusions, multi-material structural components, etc.). This study develops OLTEM using compact 9-point and 25-point stencils that are similar to those for linear and quadratic finite elements. In contrast to finite elements and other known numerical techniques for interface problems with conformed and unfitted meshes, OLTEM with 9-point and 25-point stencils and unfitted Cartesian meshes provides the 3-rd and 11-th order of accuracy for irregular interfaces, respectively; i.e. a huge increase in accuracy by eight orders for the new 'quadratic' elements compared to known techniques at similar computational costs. There are no unknowns on interfaces between different materials; the structure of the global discrete system is the same for homogeneous and heterogeneous materials (the difference in the values of the stencil coefficients). The calculation of the unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of OLTEM at a given stencil width. The numerical results with irregular interfaces show that at the same number of degrees of freedom, OLTEM with the 9-points stencils is even more accurate than the 4-th order finite elements; OLTEM with the 25-points stencils is much more accurate than the 7-th order finite elements with much wider stencils and conformed meshes. Findings The significant increase in accuracy for OLTEM by one order for 'linear' elements and by 8 orders for 'quadratic' elements compared to that for known techniques. This will lead to a huge reduction in the computation time for the problems with complex irregular interfaces. The use of trivial unfitted Cartesian meshes significantly simplifies the solution and reduces the time for the data preparation (no need in complicated mesh generators for complex geometry). Originality/value It has been never seen in the literature such a huge increase in accuracy for the proposed technique compared to existing methods. Due to a high accuracy, the proposed technique will allow the direct solution of multiscale problems without the scale separation. 
    more » « less
  5. Summary

    We present an original method for multimaterial topology optimization with elastic and thermal response considerations. The material distribution is represented parametrically using a formulation in which finite element–style shape functions are used to determine the local material properties within each finite element. We optimize a multifunctional structure that is designed for a combination of structural stiffness and thermal insulation. We conduct parallel uncoupled finite element analyses to simulate the elastic and thermal response of the structure by solving the two‐dimensional Poisson problem. We explore multiple optimization problem formulations, including structural design for minimum compliance subject to local temperature constraints so that the optimized design serves as both a support structure and a thermal insulator. We also derive and implement an original multimaterial aggregation function that allows the designer to simultaneously enforce separate maximum temperature thresholds based upon the melting point of the various design materials. The nonlinear programming problem is solved using gradient‐based optimization with adjoint sensitivity analysis. We present results for a series of two‐dimensional example problems. The results demonstrate that the proposed algorithm consistently converges to feasible multimaterial designs with the desired elastic and thermal performance.

     
    more » « less