An official website of the United States government
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
Energetically consistent inelasticity for optimization time integration
In this paper, we propose Energetically Consistent Inelasticity (ECI), a new formulation for modeling and discretizing finite strain elastoplasticity/viscoelasticity in a way that is compatible with optimization-based time integrators. We provide an in-depth analysis for allowing plasticity to be implicitly integrated through an augmented strain energy density function. We develop ECI on the associative von-Mises J2 plasticity, the non-associative Drucker-Prager plasticity, and the finite strain viscoelasticity. We demonstrate the resulting scheme on both the Finite Element Method (FEM) and the Material Point Method (MPM). Combined with a custom Newton-type optimization integration scheme, our method enables simulating stiff and large-deformation inelastic dynamics of metal, sand, snow, and foam with larger time steps, improved stability, higher efficiency, and better accuracy than existing approaches.
more »
« less
- PAR ID:
- 10358324
- Date Published:
- Journal Name:
- ACM Transactions on Graphics
- Volume:
- 41
- Issue:
- 4
- ISSN:
- 0730-0301
- Page Range / eLocation ID:
- 1 to 16
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan’s law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test is first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. The simulations show that dislocations exhibit some self-organized structures as the strain is increased.more » « less
-
null (Ed.)Structural health monitoring of complex structures is often limited by restricted accessibility to locations of interest within the structure and availability of operational loads. In this work, a novel output-only virtual sensing scheme is proposed. This scheme involves the implementation of the modal expansion in an augmented Kalman filter. Performance of the proposed scheme is compared with two existing methods. Method 1 relies on a finite element model updating, batch data processing, and modal expansion (MUME) procedure. Method 2 employs a recursive sequential estimation algorithm, which feeds a substructure model of the instrumented system into an Augmented Kalman Filter (AKF). The new scheme referred to as Method 3 (ME-AKF), implements strain estimates generated via Modal Expansion into an AKF as virtual measurements. To demonstrate the applicability of the aforementioned methods, a rollercoaster connection was instrumented with accelerometers, strain rosettes, and an optical sensor. A comparison of estimated dynamic strain response at unmeasured locations using three alternative schemes is presented. Although acceleration measurements are used indirectly for model updating, the response-only methods presented in this research use only measurements from strain rosettes for strain history predictions and require no prior knowledge of input forces. Predicted strains using all methods are shown to sufficiently predict the measured strain time histories from a control location and lie within a 95% confidence interval calculated based on modal expansion equations. In addition, the proposed ME-AKF method shows improvement in strain predictions at unmeasured locations without the necessity of batch data processing. The proposed scheme shows high potential for real-time dynamic estimation of the strain and stress state of complex structures at unmeasured locations.more » « less
-
Abstract In this article we formulate a stable computational nonlocal poromechanics model for dynamic analysis of saturated porous media. As a novelty, the stabilization formulation eliminates zero‐energy modes associated with the original multiphase correspondence constitutive models in the coupled nonlocal poromechanics model. The two‐phase stabilization scheme is formulated based on an energy method that incorporates inhomogeneous solid deformation and fluid flow. In this method, the nonlocal formulations of skeleton strain energy and fluid flow dissipation energy equate to their local formulations. The stable coupled nonlocal poromechanics model is solved for dynamic analysis by an implicit time integration scheme. As a new contribution, we validate the coupled stabilization formulation by comparing numerical results with analytical and finite element solutions for one‐dimensional and two‐dimensional dynamic problems in saturated porous media. Numerical examples of dynamic strain localization in saturated porous media are presented to demonstrate the efficacy of the stable coupled poromechanics framework for localized failure under dynamic loads.more » « less
-
Abstract This article introduces a computational design framework for obtaining three‐dimensional (3D) periodic elastoplastic architected materials with enhanced performance, subject to uniaxial or shear strain. A nonlinear finite element model accounting for plastic deformation is developed, where a Lagrange multiplier approach is utilized to impose periodicity constraints. The analysis assumes that the material obeys a von Mises plasticity model with linear isotropic hardening. The finite element model is combined with a corresponding path‐dependent adjoint sensitivity formulation, which is derived analytically. The optimization problem is parametrized using the solid isotropic material penalization method. Designs are optimized for either end compliance or toughness for a given prescribed displacement. Such a framework results in producing materials with enhanced performance through much better utilization of an elastoplastic material. Several 3D examples are used to demonstrate the effectiveness of the mathematical framework.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1145/3528223.3530072
