Mertcan Cihan, Robin Aichele, Dustin Roman Jantos, Philipp Junker
{"title":"Application of adaptive virtual element method to thermodynamic topology optimization","authors":"Mertcan Cihan, Robin Aichele, Dustin Roman Jantos, Philipp Junker","doi":"10.1002/nme.7575","DOIUrl":"https://doi.org/10.1002/nme.7575","url":null,"abstract":"In this work, a low‐order <jats:italic>virtual element method</jats:italic> (VEM) with adaptive meshing is applied for thermodynamic topology optimization with linear elastic material for the two‐dimensional case. VEM has various significant advantages compared to other numerical discretization techniques, for example the <jats:italic>finite element method</jats:italic> (FEM). One advantage is the flexibility to use arbitrary shaped elements including the possibility to add nodes on the fly during simulation, which opens a new variety of application fields. The latter mentioned feature is used in this publication to propose an adaptive mesh‐refinement strategy during the thermodynamic optimization procedure and investigate its feasibility. The thermodynamic topology optimization (TTO) includes a efficient gradient‐enhanced approach to regularization of the otherwise ill‐posed density based topology optimization approach and was already applied to various material models in combination with finite elements. However, the special numerical treatment leading to efficiency of the regularization approach within the TTO has to be modified to be applicable to VEM, which is presented in the publication. We show the performance of this framework by investigating several numerical results on benchmark problems.","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141934842","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A robust cut‐cell finite element method for Poisson's equation in three dimensions","authors":"Donghao Li, Panayiotis Papadopoulos","doi":"10.1002/nme.7577","DOIUrl":"https://doi.org/10.1002/nme.7577","url":null,"abstract":"SummaryThis article documents a cut‐cell finite element method for solving Poisson's equation in smooth three‐dimensional domains using a uniform, Cartesian axis‐aligned grid. Neumann boundary conditions are imposed weakly by way of a Delaunay triangulation, while Dirichlet boundary conditions are imposed strongly using a projection method. A set of numerical simulations demonstrates the proposed method is robust and preserves the asymptotic rate of convergence expected of corresponding body‐fitted methods.","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141934844","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Large‐scale elasto‐plastic topology optimization","authors":"Gunnar Granlund, Mathias Wallin","doi":"10.1002/nme.7583","DOIUrl":"https://doi.org/10.1002/nme.7583","url":null,"abstract":"This work presents large‐scale elasto‐plastic topology optimization for design of structures with maximized energy absorption and tailored mechanical response. The implementation uses parallel computations to address multi million element three‐dimensional problems. Design updates are generated using the gradient‐based method of moving asymptotes and the material is modeled using small strain, nonlinear isotropic hardening wherein the coaxiality between the plastic strain rate and stress is exploited. This formulation renders an efficient state solve and we demonstrate that the adjoint sensitivity scheme resembles that of the state update. Furthermore, the KKT condition is enforced directly into the path dependent adjoint sensitivity analysis which eliminates the need of monitoring the elasto‐plastic switches when calculating the gradients and provides a straight forward framework for elasto‐plastic topology optimization. Numerical examples show that structures discretized using several millions degrees of freedom and loaded in multiple load steps can be designed within a reasonable time frame.","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141969815","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Celine Lauff, Matti Schneider, John Montesano, Thomas Böhlke
{"title":"Generating microstructures of long fiber reinforced composites by the fused sequential addition and migration method","authors":"Celine Lauff, Matti Schneider, John Montesano, Thomas Böhlke","doi":"10.1002/nme.7573","DOIUrl":"10.1002/nme.7573","url":null,"abstract":"<p>We introduce the fused sequential addition and migration (fSAM) algorithm for generating microstructures of fiber composites with long, flexible, nonoverlapping fibers and industrial volume fractions. The proposed algorithm is based on modeling the fibers as polygonal chains and enforcing, on the one hand, the nonoverlapping constraints by an optimization framework. The connectivity constraints, on the other hand, are treated via constrained mechanical systems of d'Alembert type. In case of straight, that is, nonflexible, fibers, the proposed algorithm reduces to the SAM (<i>Comput. Mech.</i>, 59, 247–263, 2017) algorithm, a well-established method for generating short fiber-reinforced composites. We provide a detailed discussion of the equations governing the motion of a flexible fiber and discuss the efficient numerical treatment. We elaborate on the integration into an existing SAM code and explain the selection of the numerical parameters. To capture the fiber length distributions of long fiber reinforced composites, we sample the fiber lengths from the Gamma distribution and introduce a strategy to incorporate extremely long fibers. We study the microstructure generation capabilities of the proposed algorithm. The computational examples demonstrate the superiority of the novel microstructure-generation technology over the state of the art, realizing large fiber aspect ratios (up to 2800) and high fiber volume fractions (up to <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>32</mn>\u0000 <mo>%</mo>\u0000 </mrow>\u0000 <annotation>$$ 32% $$</annotation>\u0000 </semantics></math> for an aspect ratio of 150) for experimentally measured fiber orientation tensors.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7573","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141934911","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Accurate data‐driven surrogates of dynamical systems for forward propagation of uncertainty","authors":"Saibal De, Reese E. Jones, Hemanth Kolla","doi":"10.1002/nme.7576","DOIUrl":"https://doi.org/10.1002/nme.7576","url":null,"abstract":"Stochastic collocation (SC) is a well‐known non‐intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full‐field uncertainty propagation that characterizes the distributions of the high‐dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter‐to‐solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. This work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data‐driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC‐over‐dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full‐field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141934912","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An improved natural stabilized nodal integration for locking-related materials in meshfree methods","authors":"Huy Anh Nguyen, Satoyuki Tanaka, Tinh Quoc Bui","doi":"10.1002/nme.7584","DOIUrl":"10.1002/nme.7584","url":null,"abstract":"<p>An improved naturally stabilized nodal integration (NSNI) is presented for resolving displacement locking concerned with highly orthotropic and nearly incompressible materials in the linear setting. It is recognized that the original NSNI is susceptible to the locking when dealing with these types of materials. The proposed method utilizes spectral decomposition to split the elasticity matrix into stiff and nonstiff parts. The terms associated with the stiff modes in the bilinear form are sampled by nodal integration (NI) without stabilization, whereas the other terms are integrated with NSNI. This approach leads to a unified implementation to handle locking in both types of materials. The performance and convergence of the proposed formulation are verified through several two- and three-dimensional numerical examples, illustrating the advantages of the presented method over its standard counterpart.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7584","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141881117","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Angel A. Ciarbonetti, Sergio Idelsohn, Gisela L. Mazzieri, Ruben D. Spies
{"title":"Identification of an N-valued heterogeneous conductivity profile in an inverse heat conduction problem","authors":"Angel A. Ciarbonetti, Sergio Idelsohn, Gisela L. Mazzieri, Ruben D. Spies","doi":"10.1002/nme.7578","DOIUrl":"10.1002/nme.7578","url":null,"abstract":"<div>\u0000 \u0000 <p>In this article we deal with the problem of determining a non-homogeneous <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 </mrow>\u0000 <annotation>$$ N $$</annotation>\u0000 </semantics></math>-valued heat conductivity profile in a steady-state heat conduction boundary-value problem with mixed Dirichlet-Neumann boundary conditions over a bounded domain in <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mrow>\u0000 <mi>ℝ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>n</mi>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {mathbb{R}}^n $$</annotation>\u0000 </semantics></math>, from the knowledge of the temperature field over the whole domain. In a previous work we developed a method based on a variational approach of the PDE leading to an optimality equation which is then projected into a finite dimensional space. Discretization of the optimality equation then yields a linear although severely ill-posed equation which is then regularized via appropriate ad-hoc penalizers based upon a-priori information about the conductivities of all materials present. This process results in a generalized Tikhonov-Phillips functional whose global minimizer yields our approximate solution to the inverse problem. In our previous work we showed that this approach yields quite satisfactory results in the cases of two different conductivities. We considered here an appropriate extension of that approach for the <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 </mrow>\u0000 <annotation>$$ N $$</annotation>\u0000 </semantics></math> materials case and show a few numerical examples for the case <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 <mo>=</mo>\u0000 <mn>3</mn>\u0000 </mrow>\u0000 <annotation>$$ N=3 $$</annotation>\u0000 </semantics></math> in which the method is able to produce very good reconstructions of the exact solution.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141880941","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A novel phase-field monolithic scheme for brittle crack propagation based on the limited-memory BFGS method with adaptive mesh refinement","authors":"Tao Jin, Zhao Li, Kuiying Chen","doi":"10.1002/nme.7572","DOIUrl":"10.1002/nme.7572","url":null,"abstract":"<p>The phase-field formulation for fracture propagation is widely adopted due to its capability of naturally treating complex crack geometries. The challenges of the phase-field crack simulation include the non-convexity of the underlying energy functional and the expensive computational cost associated with the fine mesh required to resolve the phase-field length-scale around the crack region. We present a novel phase-field monolithic scheme based on the limited-memory Broyden–Fletcher–Goldfarb–Shanno (BFGS) method, or the L-BFGS method, to address the convergence difficulties usually encountered by a Newton-based approach because of the non-convex energy functional. Comparing with the conventional BFGS method, the L-BFGS monolithic scheme avoids to store the fully dense Hessian approximation matrix. This feature is critical in the context of finite element simulations. To alleviate the expensive computational cost, we integrate the proposed L-BFGS monolithic scheme with an adaptive mesh refinement (AMR) technique. We provide the algorithmic details about the proposed L-BFGS monolithic scheme, especially about how to handle the hanging-node constraints generated during the AMR process as extra linear constraints. Several two-dimensional (2D) and three-dimensional (3D) numerical examples are provided to demonstrate the capabilities of the proposed monolithic scheme, including the accuracy, the robustness, and the computational efficiency regarding the memory consumption and the wall-clock time. Particularly, we emphasize the importance of the appropriately chosen convergence criteria for brute crack propagation. The proposed L-BFGS phase-field monolithic scheme combined with the AMR technique offers an accurate, robust, and efficient approach to model brittle crack propagation in both 2D and 3D problems.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7572","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141866347","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Featured Cover","authors":"Bing-Bing Xu, Yi-Fan Wang, Peter Wriggers","doi":"10.1002/nme.7570","DOIUrl":"https://doi.org/10.1002/nme.7570","url":null,"abstract":"<p>The cover image is based on the Research Article <i>Stabilization-free virtual element method for 2D elastoplastic problems</i> by Bing-Bing Xu et al., https://doi.org/10.1002/nme.7490.\u0000 <figure>\u0000 <div><picture>\u0000 <source></source></picture><p></p>\u0000 </div>\u0000 </figure></p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-07-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7570","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141730009","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rezgar Shakeri, Leila Ghaffari, Jeremy L. Thompson, Jed Brown
{"title":"Stable numerics for finite-strain elasticity","authors":"Rezgar Shakeri, Leila Ghaffari, Jeremy L. Thompson, Jed Brown","doi":"10.1002/nme.7563","DOIUrl":"10.1002/nme.7563","url":null,"abstract":"<p>A backward stable numerical calculation of a function with condition number <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>κ</mi>\u0000 </mrow>\u0000 <annotation>$$ kappa $$</annotation>\u0000 </semantics></math> will have a relative accuracy of <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>κ</mi>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ϵ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mtext>machine</mtext>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ kappa {epsilon}_{mathrm{machine}} $$</annotation>\u0000 </semantics></math>. Standard formulations and software implementations of finite-strain elastic materials models make use of the deformation gradient <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>F</mi>\u0000 <mo>=</mo>\u0000 <mi>I</mi>\u0000 <mo>+</mo>\u0000 <mi>∂</mi>\u0000 <mi>u</mi>\u0000 <mo>/</mo>\u0000 <mi>∂</mi>\u0000 <mi>X</mi>\u0000 </mrow>\u0000 <annotation>$$ boldsymbol{F}=I+partial boldsymbol{u}/partial boldsymbol{X} $$</annotation>\u0000 </semantics></math> and Cauchy-Green tensors. These formulations are not numerically stable, leading to loss of several digits of accuracy when used in the small strain regime, and often precluding the use of single precision floating point arithmetic. We trace the source of this instability to specific points of numerical cancellation, interpretable as ill-conditioned steps. We show how to compute various strain measures in a stable way and how to transform common constitutive models to their stable representations, formulated in either initial or current configuration. The stable formulations all provide accuracy of order <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ϵ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mtext>machine</mtext>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ {epsilon}_{mathrm{machine}} $$</annotation>\u0000 </semantics></math>. In many cases, the stable formulations have elegant representations in terms of appropriate strain measures and offer geometric intuition that is lacking in their standard representation. We show that algorithmic differentiation can stably compute stresses so long as the strain energy is expressed stably, and give principles for stable computation that can be applied to inelastic materials.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":null,"pages":null},"PeriodicalIF":2.7,"publicationDate":"2024-07-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141738301","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}