skip to main content


Title: Simulation and optimization of magnetoelastic thin shells
Magnetoelastic thin shells exhibit great potential in realizing versatile functionalities through a broad range of combination of material stiffness, remnant magnetization intensity, and external magnetic stimuli. In this paper, we propose a novel computational method for forward simulation and inverse design of magnetoelastic thin shells. Our system consists of two key components of forward simulation and backward optimization. On the simulation side, we have developed a new continuum mechanics model based on the Kirchhoff-Love thin-shell model to characterize the behaviors of a megnetolelastic thin shell under external magnetic stimuli. Based on this model, we proposed an implicit numerical simulator facilitated by the magnetic energy Hessian to treat the elastic and magnetic stresses within a unified framework, which is versatile to incorporation with other thin shell models. On the optimization side, we have devised a new differentiable simulation framework equipped with an efficient adjoint formula to accommodate various PDE-constraint, inverse design problems of magnetoelastic thin-shell structures, in both static and dynamic settings. It also encompasses applications of magnetoelastic soft robots, functional Origami, artworks, and meta-material designs. We demonstrate the efficacy of our framework by designing and simulating a broad array of magnetoelastic thin-shell objects that manifest complicated interactions between magnetic fields, materials, and control policies.  more » « less
Award ID(s):
1919647
NSF-PAR ID:
10359113
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
ACM Transactions on Graphics
Volume:
41
Issue:
4
ISSN:
0730-0301
Page Range / eLocation ID:
1 to 18
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Elastic gridshell is a class of net-like structure formed by an ensemble of elastically deforming rods coupled through joints, such that the structure can cover large areas with low self-weight and allow for a variety of aesthetic configurations. Gridshells, also known as X-shells or Cosserat Nets, are a planar grid of elastic rods in its undeformed configuration. The end points of the rods are constrained and positioned on a closed curve—the final boundary—to actuate the structure into a 3D shape. Here, we report a discrete differential geometry-based numerical framework to study the geometrically nonlinear deformation of gridshell structures, accounting for non-trivial bending-twisting coupling at the joints. The form-finding problem of obtaining the undeformed planar configuration given the target convex 3D topology is then investigated. For the forward (2D to 3D) physically based simulation, we decompose the gridshell structure into multiple one-dimensional elastic rods and simulate their deformation by the well-established discrete elastic rods (DER) algorithm. A simple penalty energy between rods and linkages is used to simulate the coupling between two rods at the joints. For the inverse problem associated with form-finding (3D to 2D), we introduce a contact-based algorithm between the elastic gridshell and a rigid 3D surface, where the rigid surface describes the target shape of the gridshell upon actuation. This technique removes the need of several forward simulations associated with conventional optimization algorithms and provides a direct solution to the inverse problem. Several examples—hemispherical cap, paraboloid, and hemi-ellipsoid—are used to show the effectiveness of the inverse design process. 
    more » « less
  2. We introduce a computational framework that incorporates multiple scattering for large-scale three-dimensional (3-D) particle localization using single-shot in-line holography. Traditional holographic techniques rely on single-scattering models that become inaccurate under high particle densities and large refractive index contrasts. Existing multiple scattering solvers become computationally prohibitive for large-scale problems, which comprise millions of voxels within the scattering volume. Our approach overcomes the computational bottleneck by slicewise computation of multiple scattering under an efficient recursive framework. In the forward model, each recursion estimates the next higher-order multiple scattered field among the object slices. In the inverse model, each order of scattering is recursively estimated by a nonlinear optimization procedure. This nonlinear inverse model is further supplemented by a sparsity promoting procedure that is particularly effective in localizing 3-D distributed particles. We show that our multiple-scattering model leads to significant improvement in the quality of 3-D localization compared to traditional methods based on single scattering approximation. Our experiments demonstrate robust inverse multiple scattering, allowing reconstruction of 100 million voxels from a single 1-megapixel hologram with a sparsity prior. The performance bound of our approach is quantified in simulation and validated experimentally. Our work promises utilization of multiple scattering for versatile large-scale applications. 
    more » « less
  3. Recent developments in the computational automated design of electromagnetic devices, otherwise known as inverse design, have significantly enhanced the design process for nanophotonic systems. Inverse design can both reduce design time considerably and lead to high-performance, nonintuitive structures that would otherwise have been impossible to develop manually. Despite the successes enjoyed by structure optimization techniques, most approaches leverage electromagnetic solvers that require significant computational resources and suffer from slow convergence and numerical dispersion. Recently, a fast simulation and boundary-based inverse design approach based on boundary integral equations was demonstrated for two-dimensional nanophotonic problems. In this work, we introduce a new full-wave three-dimensional simulation and boundary-based optimization framework for nanophotonic devices also based on boundary integral methods, which achieves high accuracy even at coarse mesh discretizations while only requiring modest computational resources. The approach has been further accelerated by leveraging GPU computing, a sparse block-diagonal preconditioning strategy, and a matrix-free implementation of the discrete adjoint method. As a demonstration, we optimize three different devices: a 1:2 1550 nm power splitter and two nonadiabatic mode-preserving waveguide tapers. To the best of our knowledge, the tapers, which span 40 wavelengths in the silicon material, are the largest silicon photonic waveguiding devices to have been optimized using full-wave 3D solution of Maxwell’s equations. 
    more » « less
  4. Abstract Buckling, a phenomenon historically considered undesirable, has recently been harnessed to enable innovative functionalities in materials and structures. While approaches to achieve specific buckling behaviors are widely studied, tuning these behaviors in fabricated structures without altering their geometry remains a major challenge. Here, we introduce an inverse design approach to tune buckling behavior in magnetically active structures through the variation of applied magnetic stimuli. Our proposed magneto-mechanical topology optimization formulation not only generates the geometry and magnetization distribution of these structures but also informs how the external magnetic fields should be applied to control their buckling behaviors. By utilizing the proposed strategy, we discover magnetically active structures showcasing a broad spectrum of tunable buckling mechanisms, including programmable peak forces and buckling displacements, as well as controllable mechano- and magneto-induced bistability. Furthermore, we experimentally demonstrate that multiple unit designs can be assembled into architectures, resulting in tunable multistability and programmable buckling sequences under distinct applied magnetic fields. By employing a hybrid fabrication method, we manufacture and experimentally validate the generated designs and architectures, confirming their ability to exhibit precisely programmed and tunable buckling behaviors. This research contributes to the advancement of multifunctional materials and structures that harness buckling phenomena, unlocking transformative potential for various applications, including robotics, energy harvesting, and deployable and reconfigurable devices. 
    more » « less
  5. Abstract

    Smart materials are versatile material systems which exhibit a measurable response to external stimuli. Recently, smart material systems have been developed which incorporate graphene in order to share on its various advantageous properties, such as mechanical strength, electrical conductivity, and thermal conductivity as well as to achieve unique stimuli‐dependent responses. Here, a graphene fiber‐based smart material that exhibits reversible electrical conductivity switching at a relatively low temperature (60 °C), is reported. Using molecular dynamics (MD) simulation and density functional theory‐based non‐equilibrium Green's function (DFT‐NEGF) approach, it is revealed that this thermo‐response behavior is due to the change in configuration of amphiphilic triblock dispersant molecules occurring in the graphene fiber during heating or cooling. These conformational changes alter the total number of graphene‐graphene contacts within the composite material system, and thus the electrical conductivity as well. Additionally, this graphene fiber fabrication approach uses a scalable, facile, water‐based method, that makes it easy to modify material composition ratios. In all, this work represents an important step forward to enable complete functional tuning of graphene‐based smart materials at the nanoscale while increasing commercialization viability.

     
    more » « less