Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"Structural transient dynamic topology optimization based on autoencoder-enhanced generative adversarial network and elitist guidance evolutionary algorithm","authors":"Haojie Ma ,&nbsp;Xiao Kang ,&nbsp;Yixing Huang ,&nbsp;Shengyu Duan ,&nbsp;Ying Li ,&nbsp;Daining Fang","doi":"10.1016/j.cma.2025.118417","DOIUrl":"10.1016/j.cma.2025.118417","url":null,"abstract":"<div><div>Structural transient dynamic optimization faces significant challenges stemming from material nonlinearities and geometric nonlinearities induced by large deformations. These nonlinear phenomena severely complicate gradient-based sensitivity analysis, while conventional non-gradient optimization approaches face limitations including prohibitive computational demands, suboptimal solution quality, and compromised robustness. To overcome these challenges, we present an integrated computational framework synergistically combining an autoencoder-enhanced generative adversarial network with an elitist guidance evolutionary algorithm for nonlinear dynamic optimization. The developed multi-fidelity surrogate modeling architecture achieves dual enhancement in computational efficiency and solution diversity, while the elitism-preserving mechanism in elitist guidance evolutionary algorithm ensures superior convergence characteristics. Furthermore, we introduce a self-supervised criterion noise rate metric for quantitatively evaluating structural performance under transient loads. Results demonstrate that the proposed method improves structural clarity and diversity by 18.56 and 21.55 times compared to conventional methods. Case studies with both cantilever and fixed-end beams across dynamic loading regimes confirm the method’s generalizability. This framework is easily transferable to other engineering fields, offering new insights for solving transient nonlinear problems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118417"},"PeriodicalIF":7.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094113","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 polytopal discontinuous Galerkin method for the pseudo-stress formulation of the unsteady Stokes problem","authors":"Paola F. Antonietti,&nbsp;Michele Botti,&nbsp;Alessandra Cancrini,&nbsp;Ilario Mazzieri","doi":"10.1016/j.cma.2025.118404","DOIUrl":"10.1016/j.cma.2025.118404","url":null,"abstract":"<div><div>This work aims to construct and analyze a discontinuous Galerkin method on polytopal grids (PolydG) to solve the pseudo-stress formulation of the unsteady Stokes problem. The pseudo-stress variable is introduced due to the growing interest in non-Newtonian flows and coupled interface problems, where stress assumes a fundamental role. The space-time discretization of the problem is achieved by combining the PolydG approach with the implicit <span><math><mi>θ</mi></math></span>-method time integration scheme. For both the semi- and fully-discrete problems we present a detailed stability analysis. Moreover, we derive convergence estimates for the fully discrete space-time discretization. A set of verification tests is presented to verify the theoretical estimates and the application of the method to cases of engineering interest.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118404"},"PeriodicalIF":7.3,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094114","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":"GALDS: A graph-autoencoder-based latent dynamics surrogate model to predict neurite material transport","authors":"Tsung Yeh Hsieh ,&nbsp;Yongjie Jessica Zhang","doi":"10.1016/j.cma.2025.118409","DOIUrl":"10.1016/j.cma.2025.118409","url":null,"abstract":"<div><div>Neurons exhibit intricate geometries within their neurite networks, which play a crucial role in processes such as signaling and nutrient transport. Accurate simulation of material transport in the networks is essential for understanding these biological phenomena but poses significant computational challenges because of the complex tree-like structures involved. Traditional approaches are time-intensive and resource-demanding, yet the inherent properties of neuron trees, which consists primarily of pipes with steady-state parabolic velocity profiles and bifurcations, provide opportunities for computational optimization. To address these challenges, we propose a Graph-Autoencoder-based Latent Dynamics Surrogate (GALDS) model, which is specifically designed to streamline the simulation of material transport in neural trees. GALDS employs a graph autoencoder to encode latent representations of the network’s geometry, velocity fields, and concentration profiles. These latent space representations are then assembled into a global graph, which is subsequently used to predict system dynamics in the latent space via a trained graph latent space system dynamic model, inspired by the Neural Ordinary Differential Equations (Neural ODEs) concept. The integration of an autoencoder allows for the use of smaller graph neural network models with reduced training data requirements. Furthermore, the Neural ODE component effectively mitigates the issue of error accumulation commonly encountered in recurrent neural networks. The effectiveness of the GALDS model is demonstrated through results on eight unseen geometries and four abnormal transport examples, where our approach achieves mean relative error of <span><math><mrow><mn>3</mn><mspace></mspace><mo>%</mo></mrow></math></span> with maximum relative error <span><math><mrow><mo>&lt;</mo><mn>8</mn><mo>%</mo></mrow></math></span> and demonstrates a 10-fold speed improvement compared to previous surrogate model approaches.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118409"},"PeriodicalIF":7.3,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094116","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-fidelity modelling of floor-borne vibrations in axisymmetric MRI magnets using hp-finite element method","authors":"Yashwanth Sooriyakanthan ,&nbsp;Antonio J. Gil ,&nbsp;Paul D. Ledger ,&nbsp;Michael J. Mallett","doi":"10.1016/j.cma.2025.118385","DOIUrl":"10.1016/j.cma.2025.118385","url":null,"abstract":"<div><div>Magnetic Resonance Imaging (MRI) relies on the stability of highly uniform fields from superconducting main coils and spatially varying fields from AC-driven gradient coils. Both types of coils are thermally separated, as the main coils are cryogenically cooled within a cryostat whilst gradient coils operate at room temperature. Externally generated floor-borne vibrations (FBV) can induce relative motion between radiation shields and coils, generating eddy currents in the shields. These in turn produce parasitic magnetic fields that compromise field homogeneity and degrade image quality. This paper presents a high-fidelity computational framework for simulating the magneto-mechanical effects of FBV in axisymmetric MRI scanners to inform the manufacturing design workflow. The approach introduces three key advancements: <strong>first</strong>, a nonlinear, fully coupled magneto-mechanical formulation solved using <span><math><mrow><mi>h</mi><mi>p</mi></mrow></math></span>-Finite Element Methods (<span><math><mrow><mi>h</mi><mi>p</mi></mrow></math></span>-FEM) in the open-source NGSolve framework, with a focus on optimal interpolation order <em>p</em> and time step size; <strong>second</strong>, explicit mechanical modelling of both main and gradient coils, moving beyond idealised Biot-Savart type current sources; and <strong>third</strong>, the use of realistic axisymmetric geometries with structural connectivity between coils and radiation shields in order to inform preliminary designs in Industry. A comprehensive series of numerical results is presented in order to validate the method against some benchmarked scenarios and highlight its potential for guiding vibration mitigation and improving MRI image fidelity.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118385"},"PeriodicalIF":7.3,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094117","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":"Thermodynamically consistent coupled chemo-thermo-mechanical model of interfaces in overmolded thermoplastic parts","authors":"Junhe Cui ,&nbsp;Tiansheng Liu ,&nbsp;Michele Valsecchi ,&nbsp;Martin Giersberg ,&nbsp;Hakan Çelik ,&nbsp;Jaan-Willem Simon ,&nbsp;Sanat Kumar ,&nbsp;Jan Petersen ,&nbsp;Jacob Fish","doi":"10.1016/j.cma.2025.118359","DOIUrl":"10.1016/j.cma.2025.118359","url":null,"abstract":"<div><div>Achieving reliable bonding between dissimilar semicrystalline polymers in overmolded components remains a critical challenge in advanced manufacturing, with significant implications for structural integrity, process efficiency, and material design. This work introduces a transformational, thermodynamically consistent multiphysics framework that, for the first time, captures the full coupling between heat conduction, crystallization, deformation, and nanoscale polymer diffusion during the cooling stage of the overmolding process. The framework rigorously links manufacturing conditions to the mechanical performance of the final product by integrating process-induced residual stresses, interfacial crystallinity, and polymer interpenetration into a cohesive zone model whose fracture properties evolve dynamically.</div><div>Unlike existing approaches, which rely on phenomenological models or decoupled analyses, our formulation provides predictive capability grounded in continuum thermodynamics and validated by experimental observations. This enables not only the detection of manufacturing-induced interfacial defects but also virtual process optimization through simulation. The resulting model serves as a digital twin for overmolded thermoplastics, offering a powerful new tool for engineering high-performance composite parts in automotive, aerospace, and biomedical applications.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118359"},"PeriodicalIF":7.3,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094120","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-order time-marching schemes for incompressible flow in particle methods","authors":"Takuya Matsunaga","doi":"10.1016/j.cma.2025.118395","DOIUrl":"10.1016/j.cma.2025.118395","url":null,"abstract":"<div><div>This paper presents novel high-order time-marching schemes for simulating incompressible flow in particle methods. The proposed schemes are based on a newly developed formulation that describes the time evolution of computational variables along particle trajectories, resulting in a new form of the pressure Poisson equation. This formulation enables the direct application of existing forward-advancing time integration schemes, such as explicit Runge–Kutta methods, to achieve high-order temporal accuracy. Furthermore, the proposed schemes are generalized for arbitrary particle movement, enabling the efficient incorporation of particle shifting without requiring additional particle movement or variable corrections. By applying Runge–Kutta methods, this study presents four single- or multistage schemes referred to as RK1–RK4, corresponding to the number of stages. The validity of the proposed schemes is rigorously evaluated through numerical investigations involving four test cases and three types of particle movement (Lagrangian, Eulerian, and quasi-Lagrangian). The results reveal that the proposed RK2, RK3, and RK4 schemes achieve second-, third-, and fourth-order temporal convergence, respectively, and exhibit substantially higher accuracy than conventional first-order schemes, leading to improved volume and energy conservation. In addition, the proposed schemes demonstrate high computational efficiency, indicating their practical value for the numerical analysis of incompressible flow.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118395"},"PeriodicalIF":7.3,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094121","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":"Isogeometric analysis for non-Newtonian viscoplastic fluids: challenges for non-smooth solutions","authors":"Nicolò Antonelli ,&nbsp;Andrea Gorgi ,&nbsp;Rubén Zorrilla ,&nbsp;Riccardo Rossi","doi":"10.1016/j.cma.2025.118386","DOIUrl":"10.1016/j.cma.2025.118386","url":null,"abstract":"<div><div>This work explores the application of high-order Isogeometric Analysis (IGA) to the numerical simulation of non-Newtonian viscoplastic fluids, particularly in the presence of yield surfaces and non-smooth solutions. While IGA has demonstrated superior accuracy in smooth problems due to its high-continuity basis functions, its performance in cases with sharp transitions, such as viscoplastic flows with localized singularities, presents unique challenges. To address this, we develop a stabilized isogeometric framework for viscoplastic Stokes flow using the Variational Multiscale (VMS) method, ensuring numerical stability and preventing spurious pressure oscillations in equal-order discretizations. Additionally, we integrate an embedded boundary approach based on the Shifted Boundary Method (SBM) to efficiently handle complex geometries without the need for body-fitted meshes. The effectiveness of this high-order stabilized IGA framework is assessed through numerical benchmarks. The results confirm that high-order B-Spline bases achieve optimal convergence in smooth regions, while their performance near yield surfaces is affected by localized oscillations due to the inherent continuity of the basis functions. Furthermore, we demonstrate that the SBM-IGA formulation successfully enforces boundary conditions in embedded domains while preserving high-order accuracy. These findings provide valuable insights into the role of basis smoothness, stabilization techniques, and embedded formulations in non-Newtonian flow simulations, offering a foundation for future advancements in isogeometric methods for complex fluids.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118386"},"PeriodicalIF":7.3,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094124","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":"Rigorous error analysis of ETDRK4P-SAV scheme for the Allen-Cahn equation","authors":"Xiaoyan Li ,&nbsp;Xinlong Feng ,&nbsp;Lingzhi Qian","doi":"10.1016/j.cma.2025.118398","DOIUrl":"10.1016/j.cma.2025.118398","url":null,"abstract":"<div><div>This paper proposes a novel numerical method for the Allen-Cahn (AC) equation. By combining the dimension-splitting technique with the fourth-order exponential time-differencing Runge-Kutta(ETDRK)-scalar auxiliary variable (SAV) extrapolation method, we construct a fourth-order accurate ETDRK4P-SAV scheme with energy decay property. In terms of spatial discretization, a fourth-order central difference combined with dimension-splitting technique is employed; for temporal discretization, a fourth-order ETDRK-SAV extrapolation method based on the Padé approximation is utilized. From a theoretical perspective, we rigorously prove that the fully discrete scheme preserves the maximum principle, energy decay property and unique solvability, while also establishing the optimal error estimation theory. Numerical experimental results show that this scheme not only has good convergence but also maintains the discrete energy dissipation property, verifying its effectiveness and reliability in solving the AC equation.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118398"},"PeriodicalIF":7.3,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094119","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":"Anant-Net: Breaking the curse of dimensionality with scalable and interpretable neural surrogate for high-dimensional PDEs","authors":"Sidharth S. Menon,&nbsp;Ameya D. Jagtap","doi":"10.1016/j.cma.2025.118403","DOIUrl":"10.1016/j.cma.2025.118403","url":null,"abstract":"<div><div>High-dimensional partial differential equations (PDEs) arise in diverse scientific and engineering applications but remain computationally intractable due to the curse of dimensionality. Traditional numerical methods struggle with the exponential growth in computational complexity, particularly on hypercubic domains, where the number of required collocation points increases rapidly with dimensionality. Here, we introduce <em>Anant-Net</em>, an efficient neural surrogate that overcomes this challenge, enabling the solution of PDEs in high dimensions. Unlike hyperspheres, where the internal volume diminishes as dimensionality increases, hypercubes retain or expand their volume (for unit or larger length), making high-dimensional computations significantly more demanding. Anant-Net efficiently incorporates high-dimensional boundary conditions and minimizes the PDE residual at high-dimensional collocation points. To enhance interpretability, we integrate Kolmogorov-Arnold networks into the Anant-Net architecture. We benchmark Anant-Net’s performance on several linear and nonlinear high-dimensional equations, including the Poisson, Sine-Gordon, and Allen-Cahn equations, as well as transient heat equations, demonstrating high accuracy and robustness across randomly sampled test points from high-dimensional spaces. Importantly, Anant-Net achieves these results with remarkable efficiency, solving 300-dimensional problems on a single GPU within a few hours. We also compare Anant-Net’s results for accuracy and runtime with other state-of-the-art methods. Our findings establish Anant-Net as an accurate, interpretable, and scalable framework for efficiently solving high-dimensional PDEs. The Anant-Net code is available at <span><span>https://github.com/ParamIntelligence/Anant-Net</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118403"},"PeriodicalIF":7.3,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094118","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":"Full-scale topology optimization for dynamic responses of functionally graded porous infill designs using Nitsche-type multi-patch isogeometric analysis","authors":"Zhen Yang ,&nbsp;Liang Gao ,&nbsp;Haibin Tang ,&nbsp;Jie Gao","doi":"10.1016/j.cma.2025.118365","DOIUrl":"10.1016/j.cma.2025.118365","url":null,"abstract":"<div><div>Porous structures, with their outstanding mechanical properties, play a crucial role in engineering applications and are an important consideration in material distribution optimization for structural dynamic performance. Recently, Isogeometric Analysis (IGA) has gained significant interest due to precise geometric representation, high-order continuity, and flexible topology evolution capabilities. Hence, this study proposes a novel infill design approach through a periodic constraint strategy in multiple Non-Uniform Rational B-Splines (NURBS) patches for two dynamic topology optimization problems, namely eigenfrequency maximization and dynamic compliance minimization. By coupling multiple NURBS patches in a conforming mesh, the complexity of the structural design domain is effectively enhanced. The Nitsche-type dynamic formulation is introduced within the IGA framework, and the theoretical analysis of the stabilization condition is performed. Furthermore, the periodic constraint strategy is imposed onto NURBS patches within the specified parameter direction, which controls the sensitivity update values of the objective function across these patches to generate a gradient porous structure. The global topology is described by the Density Distribution Function (DDF) to achieve full-scale topology optimization. The Multi-frequency Quasi-Static Ritz Vector (MQSRV) method is used to reduce the computational cost associated with dynamic problems. The mathematical models for the dynamic compliance minimization and the eigenfrequency maximization are established, where the sensitivity analysis is derived in detail. Finally, the optimized results produced by the work are fully applicable to complex structural design domains and exhibit well-defined boundaries and smooth gradient distributions. Several numerical examples are presented to demonstrate the effectiveness of the proposed multi-patch isogeometric topology optimization infill design method.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"447 ","pages":"Article 118365"},"PeriodicalIF":7.3,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094122","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}