Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"The Lagrangian-Eulerian described Particle Flow Topology Optimization (PFTO) approach with isogeometric material point method","authors":"Daji Lin ,&nbsp;Liang Gao ,&nbsp;Jie Gao","doi":"10.1016/j.cma.2025.117892","DOIUrl":"10.1016/j.cma.2025.117892","url":null,"abstract":"<div><div>Recently, several topology optimization methods in the Lagrangian or Eulerian description has accepted a wide of discussions, which have been applied to several design problems. In the current work, the primary intention is to propose the Particle Flow Topology Optimization (PFTO) approach in a Lagrangian-Eulerian description, which can sufficiently unify their unique characteristics and superiorities for the optimization with a higher effectiveness. Firstly, an Isogeometric Material Points Method (I-MPM) is applied to develop the numerical analysis framework for particles to solve the static equilibrium equation, which can maintain the consistency of geometrical model and analysis model using the same NURBS (non-uniform rational B-splines) basis functions. Secondly, a Lagrangian-Eulerian particle flow model is proposed for representing structural topology in the optimization, which consists of several critical components, namely the mappings of physical information (P2C: Particles to Control points, C2G: Control points to Gauss quadrature points), the inverse mappings of sensitivity information (G2C and C2P). Thirdly, the mathematical formulation for the maximization of structural loading-capability is developed using the PFTO approach, in which physical positions and material consumptions of particles are design variables to drive the evolvement of structural topology. The sensitivity analysis of the objective and constraint functions with respect to design variables at particles, namely physical information is derived in detail. Finally, several numerical design examples in 2D and 3D are performed to demonstrate the validity, effectiveness, and superiority of the proposed PFTO approach, and the indispensability of particles movement for the optimization is also sufficiently studied.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117892"},"PeriodicalIF":6.9,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143687270","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":"EM-based fast uncertainty quantification for bayesian multi-setup operational modal analysis","authors":"Wei Zhu ,&nbsp;Binbin Li ,&nbsp;Zuo Zhu","doi":"10.1016/j.cma.2025.117942","DOIUrl":"10.1016/j.cma.2025.117942","url":null,"abstract":"<div><div>The current Bayesian FFT algorithm relies on direct differentiation to obtain the posterior covariance matrix (PCM), which is time-consuming, memory-intensive, and hard to code, especially for the multi-setup operational modal analysis (OMA). Aiming at accelerating the uncertainty quantification in multi-setup OMA, an expectation-maximization (EM)-based algorithm is proposed by reformulating the Hessian matrix of the negative log-likelihood function (NLLF) as a sum of simplified components corresponding to the complete-data NLLF. Matrix calculus is employed to derive these components in a compact manner, resulting in expressions similar to those in the single-setup case. This similarity allows for the reuse of existing Bayesian single-setup OMA codes, simplifying implementation. The singularity caused by mode shape norm constraints is addressed through null space projection, eliminating potential numerical errors from the conventional pseudoinverse operation. A sparse assembly strategy is further adopted, avoiding unnecessary calculations and storage of predominant zero elements in the Hessian matrix. The proposed method is then validated through a comprehensive parametric study and applied to a multi-setup OMA of a high-rise building. Results demonstrate that the proposed method efficiently calculates the PCM within seconds, even for cases with hundreds of parameters. This represents an efficiency improvement of at least one order of magnitude over the state-of-the-art method. Such performance paves the way for a real-time modal identification of large-scale structures, including those with closely-spaced modes.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117942"},"PeriodicalIF":6.9,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143687271","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":"Residual stress-constrained space–time topology optimization for multi-axis additive manufacturing","authors":"Kai Wu ,&nbsp;Fred van Keulen ,&nbsp;Jun Wu","doi":"10.1016/j.cma.2025.117913","DOIUrl":"10.1016/j.cma.2025.117913","url":null,"abstract":"<div><div>Residual stresses and distortions are major barriers to the broader adoption of wire arc additive manufacturing. These issues are coupled and arise due to large thermal gradients and phase transformations during the directed energy deposition process. Mitigating distortions may lead to substantial residual stresses, causing cracks in the fabricated components. In this paper, we propose a novel method to reduce both residual stresses and distortions by optimizing the fabrication sequence. This approach explores the use of non-planar layers, leveraging the increased manufacturing flexibility provided by robotic arms. Additionally, our method allows for the concurrent optimization of the structural layout and corresponding fabrication sequence. We employ the inherent strain method as a simplified process simulation model to predict residual stresses and distortions. Local residual stresses are aggregated using a <span><math><mi>p</mi></math></span>-norm function, which is integrated into distortion minimization as a constraint. Through numerical examples, we demonstrate that the optimized non-planar fabrication strategies can effectively reduce both residual stresses and distortions.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117913"},"PeriodicalIF":6.9,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143672840","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The Junction Riemann Problem in 1D shallow water channels including supercritical flow conditions","authors":"Juan Mairal,&nbsp;Javier Murillo,&nbsp;Pilar Garcia-Navarro","doi":"10.1016/j.cma.2025.117919","DOIUrl":"10.1016/j.cma.2025.117919","url":null,"abstract":"<div><div>This work presents an advancement in solving the Shallow Water Equations (SWE) in one-dimensional (1D) networks of channels using the Junction Riemann Problem (JRP). The necessity for robust solvers for junctions in networks is evident from the extensive literature and the variety of proposed methods. While multidimensional coupled approaches that model junctions as two-dimensional spaces have shown success, they lack the computational efficiency of pure 1D methods that treat junctions as singular points. In this context, existing JRP-based methods have primarily been limited to subcritical flow regimes. For the first time, this paper demonstrates that the Junction Riemann Problem can be effectively used as an internal boundary condition across all flow regimes representing a junction of channels. The proposed JRP solution is both simple and robust, accommodating various flow regimes and an arbitrary number of channels without requiring additional information. Furthermore, it is shown that the JRP can be solved efficiently at internal boundaries and integrated with a standard first-order Godunov scheme to yield accurate results. The validation of this method is confirmed through a series of test cases, highlighting its effectiveness in modeling free-surface flows.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117919"},"PeriodicalIF":6.9,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143672838","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 multiscale design method using interpretable machine learning for phononic materials with closely interacting scales","authors":"Mary V. Bastawrous ,&nbsp;Zhi Chen ,&nbsp;Alexander C. Ogren ,&nbsp;Chiara Daraio ,&nbsp;Cynthia Rudin ,&nbsp;L. Catherine Brinson","doi":"10.1016/j.cma.2025.117833","DOIUrl":"10.1016/j.cma.2025.117833","url":null,"abstract":"<div><div>Manipulating the dispersive characteristics of vibrational waves is beneficial for many applications, e.g., high-precision instruments. architected hierarchical phononic materials have sparked promise tunability of elastodynamic waves and vibrations over multiple frequency ranges. In this article, hierarchical unit-cells are obtained, where features at each length scale result in a band gap within a targeted frequency range. Our novel approach, the “hierarchical unit-cell template method,” is an interpretable machine-learning approach that uncovers global unit-cell shape/topology patterns corresponding to predefined band-gap objectives. A scale-separation effect is observed where the coarse-scale band-gap objective is mostly unaffected by the fine-scale features despite the closeness of their length scales, thus enabling an efficient hierarchical algorithm. Moreover, the hierarchical patterns revealed are not predefined or self-similar hierarchies as common in current hierarchical phononic materials. Thus, our approach offers a flexible and efficient method for the exploration of new regions in the hierarchical design space, extracting minimal effective patterns for inverse design in applications targeting multiple frequency ranges.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117833"},"PeriodicalIF":6.9,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143672841","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":"Self-propelling, soft, and slender structures in fluids: Cosserat rods immersed in the velocity–vorticity formulation of the incompressible Navier–Stokes equations","authors":"Arman Tekinalp ,&nbsp;Yashraj Bhosale ,&nbsp;Songyuan Cui ,&nbsp;Fan Kiat Chan ,&nbsp;Mattia Gazzola","doi":"10.1016/j.cma.2025.117910","DOIUrl":"10.1016/j.cma.2025.117910","url":null,"abstract":"<div><div>We present a hybrid Eulerian–Lagrangian method for the direct simulation of three-dimensional, heterogeneous, active, and self-propelling structures made of soft fibers and operating in incompressible viscous flows. Fiber-based organization of matter is pervasive in nature and engineering, from biological architectures made of cilia, hair, muscles or bones to polymers, composite materials or soft robots. In nature, many such structures are adapted to manipulate flows for feeding, swimming or energy harvesting, through mechanisms that are often not fully understood. While simulations can support the analysis (and subsequent translational engineering) of these systems, extreme fibers’ aspect-ratios, large elastic deformations, two-way coupling with three-dimensional flows, and self-propulsion all render the problem numerically challenging. To address this, we couple Cosserat rod theory, where fibers’ dynamics is accurately captured in one-dimensional fashion, with the velocity–vorticity formulation of the Navier–Stokes equations, through a virtual boundary technique. The favorable properties of the resultant hydroelastic solver are demonstrated against a battery of benchmarks, and further showcased in a range of multi-physics scenarios, involving magnetic actuation, viscous streaming, biomechanics, multi-body interaction, and untethered swimming.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117910"},"PeriodicalIF":6.9,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644321","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Harnessing dynamic turbulent dynamics in parrot optimization algorithm for complex high-dimensional engineering problems","authors":"Mahmoud Abdel-Salam ,&nbsp;Saleh Ali Alomari ,&nbsp;Jing Yang ,&nbsp;Sangkeum Lee ,&nbsp;Kashif Saleem ,&nbsp;Aseel Smerat ,&nbsp;Vaclav Snasel ,&nbsp;Laith Abualigah","doi":"10.1016/j.cma.2025.117908","DOIUrl":"10.1016/j.cma.2025.117908","url":null,"abstract":"<div><div>The Parrot Optimization Algorithm (PO) is a nature-inspired metaheuristic algorithm developed based on the social and adaptive behaviors of Pyrrhura molinae parrots. PO demonstrates robust optimization performance by balancing exploration and exploitation, mimicking foraging and cooperative activities. However, as the algorithm progresses through iterations, it faces critical challenges in maintaining search diversity and movement efficiency diminishes, leading to premature convergence and a reduced ability to find optimal solutions in complex search space. To address these limitations, this work introduces the Dynamic Turbulent-based Parrot Optimization Algorithm (DTPO), which represents a significant advancement over the original PO by incorporating three novel strategies: a novel Differential Mutation (DM), Dynamic Opposite Learning (DOL), and Turbulent Operator (TO). The DM Strategy enhances exploration by introducing controlled variations in the population, allowing DTPO to escape local optima. Also, the DOL Strategy dynamically generates opposite solutions to refresh stagnated populations, expanding the search space and maintaining adaptability. Finally, the TO strategy simulates chaotic movements inspired by turbulence, ensuring a thorough local search while preserving population diversity. Together, these strategies improve the algorithm's ability to explore, exploit, and converge efficiently. Furthermore, the DTPO's performance was rigorously evaluated on benchmark functions from CEC2017 and CEC2022, comparing it against 23 state-of-the-art algorithms. The results demonstrate DTPO's superior convergence speed, search efficiency, and optimization accuracy. Additionally, DTPO was tested on seven engineering design problems, achieving significant improvements over the original PO algorithm, with superior performance gains compared to other algorithms in real-world scenarios. Particularly, DTPO outperformed competing algorithms in 37 out of 41 benchmark functions, achieving an overall success rate of 90.24%. Moreover, DTPO obtained the best Friedman ranks across all comparisons, with values ranging from 3.03 to 1.18, demonstrating its superiority over classical, advanced, and recent algorithms. These results validate the proposed enhancements and highlight DTPO's robustness and effectiveness in solving complex optimization problems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117908"},"PeriodicalIF":6.9,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644323","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":"Spatiotemporal modeling based on manifold learning for collision dynamic prediction of thin-walled structures under oblique load","authors":"Jian Xie ,&nbsp;Junyuan Zhang ,&nbsp;Hao Zhou ,&nbsp;Zihang Li ,&nbsp;Zhongyu Li","doi":"10.1016/j.cma.2025.117926","DOIUrl":"10.1016/j.cma.2025.117926","url":null,"abstract":"<div><div>Numerical simulation of the collision dynamics in thin-walled structures under oblique load involves complex spatiotemporal processes, including material, geometric, and contact nonlinearities, which often require significant computational resources and time. Moreover, predicting high-dimensional spatiotemporal responses remains a challenge for most surrogate-based models. This paper proposes a deep learning framework based on manifold learning for spatiotemporal modeling of collision dynamics in thin-walled structures under oblique load. The framework leverages multiple deep learning models, including Variational Autoencoders (VAE), Radial Basis Function Interpolation (RBFI), and regression Residual Network (ResNet18), to capture the complex nonlinearities inherent in structural deformation, stress distribution, and crush force, enabling continuous prediction of multimodal spatiotemporal responses. Using a rectangular thin-walled tube under oblique load as an example, the models are validated with simulation data, yielding average prediction errors of 5.80 % for structural deformation, 6.01 % for Energy Absorption (EA), 10.66 % for Peak Crush Force (PCF), and 16.66 % for crush force. Compared to traditional finite element (FE) simulations, prediction time is reduced by 98.6 % for structural deformation and stress distribution, and 97.4 % for crush force. Additionally, the method demonstrates stability and broad applicability across different design parameters and structural configurations, including rectangular and double-cell tubes. This work underscores the potential of deep learning techniques to enhance computational efficiency and predictive accuracy in the crashworthiness design of thin-walled structures.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117926"},"PeriodicalIF":6.9,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644320","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":"Mutual-information-based dimensional learning: Objective algorithms for identification of relevant dimensionless quantities","authors":"Lei Zhang,&nbsp;Guowei He","doi":"10.1016/j.cma.2025.117922","DOIUrl":"10.1016/j.cma.2025.117922","url":null,"abstract":"<div><div>The classical dimensional analysis provides powerful insights into underlying physical mechanisms, but has limitations in determining the uniqueness and measuring the relative importance of dimensionless quantities. To address these limitations, we propose a data-driven approach, called mutual-information-based dimensional learning, to identify unique and relevant dimensionless quantities from available data. The proposed method employs a novel information-theoretic criterion to measure the relative importance of dimensionless quantities, whereas the existing methodologies rely on sensitivity/derivative-based measures. This entropy-based measure provides two significant advantages: (1) invariance (objectivity) with respect to reparametrizations of variables, and (2) robustness against outliers. Numerical results show that our method outperforms the current state-of-the-art method in these aspects, and enables identifying dominant dimensionless quantities. Examples include the study of the friction factor in benchmark pipe flows, the eddy viscosity coefficients in turbulent channel flows and the vapor depression dynamics in laser–metal interaction.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117922"},"PeriodicalIF":6.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644319","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 novel energy-fitted hexagonal quadrature scheme enables low-cost and high-fidelity peridynamic computations","authors":"Emely Schaller ,&nbsp;Ali Javili ,&nbsp;Paul Steinmann","doi":"10.1016/j.cma.2025.117918","DOIUrl":"10.1016/j.cma.2025.117918","url":null,"abstract":"<div><div>In this contribution, we propose a novel hexagonal quadrature scheme for one-neighbor interactions in continuum-kinematics-inspired peridynamics equivalent to bond-based peridynamics. The hexagonal quadrature scheme is fitted to correctly integrate the stored energy density within the nonlocal finite-sized neighborhood of a continuum point subject to affine expansion. Our proposed hexagonal quadrature scheme is grid-independent by relying on appropriate interpolation of pertinent quantities from collocation to quadrature points. In this contribution, we discuss linear and quadratic interpolations and compare our novel hexagonal quadrature scheme to common grid-dependent quadrature schemes. For this, we consider both, tetragonal and hexagonal discretizations of the domain. The accuracy of the presented quadrature schemes is first evaluated and compared by computing the stored energy density of various prescribed affine deformations within the nonlocal neighborhood. Furthermore, we perform three different boundary value problems, where we measure the effective Poisson’s ratio resulting from each quadrature scheme and evaluate the deformation of a unit square under extension and beam bending. Key findings of our studies are: The Poisson’s test is a good indicator for the convergence behavior of quadrature schemes with respect to the grid density. The accuracy of quadrature schemes depends, as expected, on their ability to appropriately capture the deformation within the nonlocal neighborhood. Our novel hexagonal quadrature scheme, rendering the correct effective Poisson’s ratio of <span><math><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></math></span> for small deformations, together with quadratic interpolation consequently yields the most accurate results for the studies presented in this contribution, thereby effectively reducing the computational cost.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"440 ","pages":"Article 117918"},"PeriodicalIF":6.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644324","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}