Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"Transient dynamic robust topology optimization methodology for continuum structure under stochastic uncertainties","authors":"Zeng Meng ,&nbsp;Zixuan Tian ,&nbsp;Yongxin Gao ,&nbsp;Matthias G.R. Faes ,&nbsp;Quhao Li","doi":"10.1016/j.cma.2025.118019","DOIUrl":"10.1016/j.cma.2025.118019","url":null,"abstract":"<div><div>Time-variant uncertainties are omnipresent in engineering systems. These significantly impact the structural performance. The main challenge in this context is how to handle them in dynamic domain response topology optimization. To tackle this challenge, a new transient dynamic robust topology optimization (TDRTO) method is proposed to optimize the topology of continuous structures. This method comprehensively considers the uncertainties of material property, loading directions, and time-variant stochastic parameters of loading amplitudes. The time-variant performance function is transformed into a set of independent instantaneous performance functions, where the stochastic processes are discretized by using the optimal linear estimation method to simulate the correlations among various time nodes. The mean and standard deviation of the structural compliance are approximated through a Taylor expansion. Moreover, the Hilber-Hughes-Taylor <em>α</em> method is employed to address the structural dynamic problem. The design and stochastic sensitivities are derived by the “discretize-then-differentiate” and the adjoint methods, thereby improving the computational efficiency. Three illustrative cases are tested to validate the efficacy of TDRTO method, which shows its superiority over the traditional robust topology optimization method for dealing with time-variant stochastic uncertainties.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118019"},"PeriodicalIF":6.9,"publicationDate":"2025-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143876669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Simultaneous topology and fiber path optimization for variable stiffness Double-Double laminates with strength control","authors":"Dan Wang ,&nbsp;Yucheng Zhong ,&nbsp;David W. Rosen ,&nbsp;Sridhar Narayanaswamy","doi":"10.1016/j.cma.2025.118021","DOIUrl":"10.1016/j.cma.2025.118021","url":null,"abstract":"<div><div>Variable stiffness laminates offer the advantage of tailoring structural performance by adjusting in-plane stiffness through curved fiber paths. Additionally, material distribution at the structural level can further fine-tune performance by varying the topology. If both the structural topology and curved fiber paths are optimized together, super-efficient composite laminates can be achieved. In this paper, a simultaneous topology and fiber path optimization method based on a coarse background mesh is proposed for variable stiffness Double-Double (DD) laminates with strength control. Firstly, elemental pseudo densities and nodal fiber orientations are selected as design variables to control topology and curved fiber paths, respectively. Due to the reduced design redundancy in DD laminates, only two independent fiber orientations are necessary for each node of a coarse background mesh. This ensures global fiber path continuity through interpolation of elemental values. Additionally, it helps decouple fiber path generation from the computationally expensive finite element analysis, significantly reducing the number of design variables and related constraints. Secondly, an optimization model is developed to maximize the laminate strength while satisfying weight and compliance constraints. The nested p-norm of Tsai-Hill failure indices is used to aggregate stresses across different elements and layers, enhancing the overall calculation efficiency. Additionally, minimum angle difference constraints are incorporated between the two groups of DD angles to ensure the laminate’s ability to resist secondary loads, thereby improving structural integrity. Notably, the proposed framework is the first to address the simultaneous optimization of both topology and curved fiber paths with strength considerations for multi-layer composite laminates. Sensitivities for both topology and fiber orientation design variables are efficiently calculated by solving a single adjoint problem, significantly improving computational efficiency. Finally, representative numerical examples demonstrate the effectiveness of the proposed method, achieving significant stress concentration reductions (over 40 %) compared to results from topology optimization alone. The optimized designs exhibit more compact topologies, improved load transfer, and enhanced resistance to secondary loading, validating the robustness of the proposed method.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118021"},"PeriodicalIF":6.9,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143875008","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An unconditionally stable variable time step scheme for two-phase ferrofluid flows","authors":"Aytura Keram ,&nbsp;Pengzhan Huang ,&nbsp;Yinnian He","doi":"10.1016/j.cma.2025.118018","DOIUrl":"10.1016/j.cma.2025.118018","url":null,"abstract":"<div><div>In this paper, a decoupled, linearized, unconditionally stable, and fully discrete numerical scheme is presented for simulating two-phase ferrofluid flows. This scheme is constructed by introducing two scalar auxiliary variables. It is based on the backward Euler scheme with variable time step and mixed finite element discretization. Nonlinear terms are treated explicitly to simplify the computational process. Meanwhile, without any restriction on time step, we show the stability of the proposed scheme. Finally, numerical examples are presented to check the accuracy and efficiency of the proposed scheme.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118018"},"PeriodicalIF":6.9,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143875009","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Overlapping Schwarz preconditioners for randomized neural networks with domain decomposition","authors":"Yong Shang ,&nbsp;Alexander Heinlein ,&nbsp;Siddhartha Mishra ,&nbsp;Fei Wang","doi":"10.1016/j.cma.2025.118011","DOIUrl":"10.1016/j.cma.2025.118011","url":null,"abstract":"<div><div>Randomized neural networks (RaNNs), characterized by fixed hidden layers after random initialization, offer a computationally efficient alternative to fully parameterized neural networks trained using stochastic gradient descent-type algorithms. In this paper, we integrate RaNNs with overlapping Schwarz domain decomposition in two primary ways: firstly, to formulate the least-squares problem with localized basis functions, and secondly, to construct effective overlapping Schwarz preconditioners for solving the resulting linear systems. Specifically, neural networks are randomly initialized in each subdomain following a uniform distribution, and these localized solutions are combined through a partition of unity, providing a global approximation to the solution of the partial differential equation. Boundary conditions are imposed via a constraining operator, eliminating the necessity for penalty methods. Furthermore, we apply principal component analysis (PCA) within each subdomain to reduce the number of basis functions, thereby significantly improving the conditioning of the resulting linear system. By constructing additive Schwarz (AS) and restricted AS preconditioners, we efficiently solve the least-squares problems using iterative solvers such as the Conjugate Gradient (CG) and generalized minimal residual methods. Numerical experiments clearly demonstrate that the proposed methodology substantially reduces computational time, particularly for multi-scale and time-dependent PDE problems. Additionally, we present a three-dimensional numerical example illustrating the superior efficiency of employing the CG method combined with an AS preconditioner over direct methods like QR decomposition for solving the associated least-squares system.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118011"},"PeriodicalIF":6.9,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868378","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An optimal-transport finite-particle method for driven mass diffusion","authors":"A. Pandolfi ,&nbsp;I. Romero ,&nbsp;M. Ortiz","doi":"10.1016/j.cma.2025.118013","DOIUrl":"10.1016/j.cma.2025.118013","url":null,"abstract":"<div><div>We formulate a finite-particle method of mass transport that accounts for general mixed boundary conditions. The particle method couples a geometrically-exact treatment of advection; Wasserstein gradient-flow dynamics; and a Kullback–Leibler representation of the entropy. General boundary conditions are enforced by introducing an adsorption/depletion layer at the boundary wherein particles are added or removed as dictated by the boundary conditions. We demonstrate the range and scope of the method through a number of examples of application, including absorption of particles into a sphere and flow through pipes of square and circular cross section, with and without occlusions. In all cases, the solution is observed to converge weakly, or in the sense of local averages.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118013"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143864541","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Parallel Algebraic Multigrid Solvers for Composite Discontinuous Galerkin Discretization of the Cardiac EMI Model in Heterogeneous Media","authors":"Edoardo Centofanti ,&nbsp;Ngoc Mai Monica Huynh ,&nbsp;Luca F. Pavarino ,&nbsp;Simone Scacchi","doi":"10.1016/j.cma.2025.118001","DOIUrl":"10.1016/j.cma.2025.118001","url":null,"abstract":"<div><div>In this paper, we develop and numerically study algebraic multigrid (AMG) preconditioners for the cardiac EMI (Extracellular space, cell Membrane, and Intracellular space) model, a recent and biophysically detailed framework for cardiac electrophysiology. The EMI model addresses the limitations of traditional homogenized cardiac models and leverages contemporary computational power to enable high-resolution simulations at the cellular scale. Using a composite Discontinuous Galerkin (DG) discretization, we introduce an AMG-EMI solver for the three dimensional EMI model. Our investigation includes the AMG-EMI scalability performance, both weak and strong, and evaluates its numerical robustness under ischemic conditions, addressing the challenges of heterogeneous media. Numerical tests exploit state-of-the-art pre-exascale supercomputers with hybrid CPU–GPU architectures. The results indicate better scalability performance of the AMG-EMI solver on CPUs compared to GPUs. However, the best solution times achieved using GPUs are up to 40x faster than those obtained on CPUs.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118001"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868380","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Monolithically coupled framework for mass and momentum balance: An open system approach","authors":"Samir El Masri,&nbsp;Barış Cansız,&nbsp;Johannes Storm,&nbsp;Michael Kaliske","doi":"10.1016/j.cma.2025.118017","DOIUrl":"10.1016/j.cma.2025.118017","url":null,"abstract":"<div><div>The finite element method (FEM) and its associated field have mainly been developed for adiabatic and closed systems. Nonetheless, open systems, which allow for the exchange of energy and mass with the surroundings, have gained increasing interest in applications where mass change occurs. For solving open systems two approaches can be undertaken. The first is the local approach, which incorporates mass change as an internal variable at the material level, while the second is the global approach, which treats mass change as an additional degree of freedom (DOF), solving the mass and momentum balance equations simultaneously. Although the global approach has been already developed, it has not yet incorporated a kinematic split of the deformation gradient. This split is necessary for modeling large strain deformations volume change (e.g. soft tissues). Hence, this study proposes a monolithic coupled mass-mechanical framework with a multiplicative split of the deformation gradient. The deformation gradient is multiplicatively split into mass-changing and mechanical components, with the mass-changing part accommodating orthotropic deformation and constraints enforcing density preservation. The study presents the complete finite element method from the kinematic foundations through to the discretization process. A sensitivity analysis is conducted to study the effects of various factors on the deformation and mass change. Moreover, a numerical example demonstrating the framework’s application to a general mass change problem is also conducted. The results show that the proposed framework effectively models mass-changing phenomena, offering a tool for future research in the field of open systems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118017"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868379","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Coalescing and break-up of viscous drops with surface tension through the Smoothed Particle Hydrodynamics","authors":"M. Antuono ,&nbsp;S. Marrone ,&nbsp;A. Colagrossi","doi":"10.1016/j.cma.2025.118014","DOIUrl":"10.1016/j.cma.2025.118014","url":null,"abstract":"<div><div>A study on the phenomena of drop coalescence and break-up is tackled for viscous flows with surface tension. The proposed analysis focuses on single-phase problems and is carried on by using the Smoothed Particle Hydrodynamics (SPH). This scheme allows for a straightforward treatment of the free surface and for an accurate description of the energy balance during the most delicate stages of coalescence and break-up. This latter aspect, combined with theoretical relations on the surface-tension energy, allowed for an in-depth characterization of the above phenomena in terms of energy evolution. A number of numerical test cases are considered in order to support and illustrate the main findings of the proposed work.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118014"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143864540","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A fourth-order reaction diffusion-based level set method for isogeometric topology optimization","authors":"He Li,&nbsp;Jianhu Shen,&nbsp;Xuyu Zhang,&nbsp;Shiwei Zhou","doi":"10.1016/j.cma.2025.118028","DOIUrl":"10.1016/j.cma.2025.118028","url":null,"abstract":"<div><div>This study presents a fourth-order reaction-diffusion isogeometric optimization method to effectively control curvature variations in minimum mean compliance optimization problems. Using isogeometric analysis with <em>k</em>-refinement technique, the level set function—parameterized using Non-Uniform Rational B-Splines (NURBS) to represent complex geometries while maintaining computational stability accurately—is updated to achieve smoother geometries with higher-order continuity. The elasticity equation is also solved using isogeometric analysis, which preserves precise geometric representation and eliminates the approximation errors associated with finite element analysis. Numerical examples show that the proposed method generates sharper, corner-free complex structures in significantly less computational time than traditional second-order reaction-diffusion methods. For instance, the proposed method produces a 2D quarter annulus under a 40 % volume constraint in just 13 iterations. At the same time, it only needs 20 iterations to yield an elegant 3D serpentine structure in an arbitrarily shaped design domain. The method demonstrates high efficiency, superior accuracy, and enhanced continuity, indicating its potential for various engineering applications.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118028"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868509","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"High-throughput machine learning framework for predicting neurite deterioration using MetaFormer attention","authors":"Kuanren Qian ,&nbsp;Genesis Omana Suarez ,&nbsp;Toshihiko Nambara ,&nbsp;Takahisa Kanekiyo ,&nbsp;Yongjie Jessica Zhang","doi":"10.1016/j.cma.2025.118003","DOIUrl":"10.1016/j.cma.2025.118003","url":null,"abstract":"<div><div>Neurodevelopmental disorders (NDDs) cover a variety of conditions, including autism spectrum disorder, attention-deficit/hyperactivity disorder, and epilepsy, which impair the central and peripheral nervous systems. Their high comorbidity and complex etiologies present significant challenges for accurate diagnosis and effective treatments. Conventional clinical and experimental studies are time-intensive, burdening research progress considerably. This paper introduces a high-throughput machine learning (ML) framework for modeling neurite deteriorations associated with NDDs, integrating synthetic data generation, experimental images, and ML models. The synthetic data generator utilizes an isogeometric analysis (IGA)-based phase field model to capture diverse neurite deterioration patterns such as neurite retraction, atrophy, and fragmentation while mitigating the limitations of scarce experimental data. The ML model utilizes MetaFormer-based gated spatiotemporal attention architecture with deep temporal layers and provides fast predictions. The framework effectively captures long-range temporal dependencies and intricate morphological transformations with average errors of 1.9641% and 6.0339% for synthetic and experimental neurite deterioration, respectively. Seamlessly integrating simulations, experiments, and ML framework can guide researchers to make informed experimental decisions by predicting potential experimental outcomes, significantly reducing costs and saving valuable time. It can also advance our understanding of neurite deterioration and provide a scalable solution for exploring complex neurological mechanisms, contributing to the development of targeted treatments.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"442 ","pages":"Article 118003"},"PeriodicalIF":6.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868381","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}