Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"Addressing concave boundaries in two-dimensional pointwise contact detection under the common-normal concept","authors":"Lucas da Silva,&nbsp;Marina Vendl Craveiro,&nbsp;Alfredo Gay Neto","doi":"10.1016/j.cma.2025.117865","DOIUrl":"10.1016/j.cma.2025.117865","url":null,"abstract":"<div><div>Contact search, the step where pairs of interacting points are identified, is crucial in computer methods for contact mechanics. This work explores the properties of contact pairs in a specific approach known as master-master method, combined with a hybrid-barrier enforcement method. The scope is on two-dimensional non-conformal contact, modeled as pointwise. Line-to-line and other instances of flat contact, for which a distribution of pressure over a region of finite size better represents their physics, are avoided. The main goal is to overcome the non-uniqueness of solutions when considering concave geometries. The bodies are defined by parameterized plane curves composed of strictly convex segments that represent either convex or concave boundaries. In the master-master approach, contact pairs are characterized by the common normal concept. Within this framework, contact pairs are classified into four types: convex-convex, matchable convex-concave, non-matchable convex-concave, and concave-concave. The Hessian of the squared distance function is analyzed for each type to further characterize them. Characterization using the Hessian matrix reveals that convex-convex and matchable convex-concave pairs are local minimizers of the squared distance function, while the other two types are either saddle points or maximizers. This enables a demonstration of the uniqueness of solutions for convex bodies. In the convex-concave case, projecting the concave boundary onto the convex one results in a univariate restricted objective function that distinguishes matchable pairs as minimizers and non-matchable pairs as maximizers. This function is used to propose a robust search algorithm that includes subdividing the domain into intervals with at most one minimizer, enabling the practical use of iterative minimization techniques to find all desired contact solutions. An algorithm for contact search that accommodates concave geometries is especially valuable in multibody applications.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117865"},"PeriodicalIF":6.9,"publicationDate":"2025-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143528722","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":"Forward and inverse simulation of pseudo-two-dimensional model of lithium-ion batteries using neural networks","authors":"Myeong-Su Lee ,&nbsp;Jaemin Oh ,&nbsp;Dong-Chan Lee ,&nbsp;KangWook Lee ,&nbsp;Sooncheol Park ,&nbsp;Youngjoon Hong","doi":"10.1016/j.cma.2025.117856","DOIUrl":"10.1016/j.cma.2025.117856","url":null,"abstract":"<div><div>In this work, we address the challenges posed by the high nonlinearity of the Butler–Volmer (BV) equation in forward and inverse simulations of the pseudo-two-dimensional (P2D) model using the physics-informed neural network (PINN) framework. The BV equation presents significant challenges for PINNs, primarily due to the hyperbolic sine term, which renders the Hessian of the PINN loss function highly ill-conditioned. To address this issue, we introduce a bypassing term that improves numerical stability by substantially reducing the condition number of the Hessian matrix. Furthermore, the small magnitude of the ionic flux <span><math><mi>j</mi></math></span> often leads to a common failure mode where PINNs converge to incorrect solutions. We demonstrate that incorporating a secondary conservation law for the solid-phase potential <span><math><mi>ψ</mi></math></span> effectively prevents such convergence issues and ensures solution accuracy. The proposed methods prove effective for solving both forward and inverse problems involving the BV equation. Specifically, we achieve precise parameter estimation in inverse scenarios and reliable solution predictions for forward simulations.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117856"},"PeriodicalIF":6.9,"publicationDate":"2025-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143527281","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":"Octree-based scaled boundary finite element approach for polycrystal RVEs: A comparison with traditional FE and FFT methods","authors":"Shiva Kumar Gaddam ,&nbsp;Sundararajan Natarajan ,&nbsp;Anand K. Kanjarla","doi":"10.1016/j.cma.2025.117864","DOIUrl":"10.1016/j.cma.2025.117864","url":null,"abstract":"<div><div>Finite Element Method (FEM) is one of the most widely used numerical techniques for solving partial differential equations. Despite its popularity, FEM faces challenges such as automatic mesh generation, handling stress singularities, and adaptive meshing. The recently developed Scaled Boundary Finite Element Method (SBFEM) overcomes these challenges by utilizing polyhedral elements, such as octree elements. SBFEM, combined with octree meshes, offer significant advantages over FEM, including rapid mesh transition, automatic mesh generation, adaptive meshing, and enhanced computational efficiency. Octree-based SBFEM has been successfully implemented and tested in various applications, such as homogenization, elastoplasticity, and adaptive phase-field fracture. However, its application to polycrystal representative volume elements (RVEs) remains unexplored. In this work, we implemented octree-based SBFEM for polycrystal RVEs and evaluated its performance for elasticity. A detailed algorithm is provided to generate balanced periodic octree meshes for polycrystal RVEs. The homogenized response and local stress fields are compared with those obtained from FEM and fast Fourier transforms (FFT). The results demonstrate that SBFEM closely matches with FEM and FFT while offering the added advantage of computational efficiency over FEM.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117864"},"PeriodicalIF":6.9,"publicationDate":"2025-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519698","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":"Optimal solutions employing an algebraic Variational Multiscale approach part I: Steady Linear Problems","authors":"Suyash Shrestha ,&nbsp;Marc Gerritsma ,&nbsp;Gonzalo Rubio ,&nbsp;Steven Hulshoff ,&nbsp;Esteban Ferrer","doi":"10.1016/j.cma.2025.117832","DOIUrl":"10.1016/j.cma.2025.117832","url":null,"abstract":"<div><div>This work extends our previous study from S. Shrestha et al. (2024) by introducing a new abstract framework for Variational Multiscale (VMS) methods at the discrete level. We introduce the concept of what we define as the optimal projector and present a discretisation approach that yields a numerical solution closely approximating the optimal projection of the infinite-dimensional continuous solution. In this approach, the infinite-dimensional unresolved scales are approximated in a finite-dimensional subspace using the numerically computed Fine-Scale Greens’ function of the underlying symmetric problem. The proposed approach involves solving the VMS problem on two separate meshes: a coarse mesh for the full PDE and a fine mesh for the symmetric part of the continuous differential operator. We consider the 1D and 2D steady advection–diffusion problems in both direct and mixed formulations as the test cases in this paper. We first present an error analysis of the proposed approach and show that the projected solution is achieved as the approximate Greens’ function converges to the exact one. Subsequently, we demonstrate the working of this method where we show that it can exponentially converge to the chosen optimal projection. We note that the implementation of the present work employs the Mimetic Spectral Element Method (MSEM), although, it may be applied to other Finite/Spectral Element or Isogeometric frameworks. Furthermore, we propose that VMS should not be viewed as a stabilisation technique; instead, the base scheme should be inherently stable, with VMS enhancing the solution quality by supplementing the base scheme.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117832"},"PeriodicalIF":6.9,"publicationDate":"2025-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143527280","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 robust and efficient rate-independent crystal plasticity model based on successive one-dimensional solution steps","authors":"B. Nijhuis,&nbsp;E.S. Perdahcıoğlu,&nbsp;A.H. van den Boogaard","doi":"10.1016/j.cma.2025.117815","DOIUrl":"10.1016/j.cma.2025.117815","url":null,"abstract":"<div><div>An efficient stress update algorithm for rate-independent crystal plasticity is presented. A series of successive one-dimensional solution (SODS) steps traces the hypersurfaces describing the slip state for which the yield criteria of individual slip systems are fulfilled to identify the intersection of all hypersurfaces. This provides both the active set and all slip components without requiring iterative active set search procedures or inducing spurious slip on inactive systems. The basic SODS algorithm is accelerated by tracking the evolution of the active set. A fast Newton–Raphson procedure enables to obtain the solution for an unchanging active set directly, while line search and extrapolation procedures direct the SODS steps towards the solution faster. A regularised tangent modulus is proposed that eliminates stiffness jumps upon changes in active set to improve the convergence behaviour of outer (equilibrium) iterations conducted with the algorithm. The resulting stress update algorithm is highly stable and efficient, making it an attractive candidate for use in large-scale crystal plasticity FE simulations and homogenisation algorithms.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117815"},"PeriodicalIF":6.9,"publicationDate":"2025-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143519688","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":"Based on purely physical information in deep learning optimizes soliton system parameter identification problem","authors":"Zhiyang Zhang ,&nbsp;Muwei Liu ,&nbsp;Xiaowei Xing ,&nbsp;Shuzhuang Zhang ,&nbsp;Zhenya Yan ,&nbsp;Wenjun Liu","doi":"10.1016/j.cma.2025.117852","DOIUrl":"10.1016/j.cma.2025.117852","url":null,"abstract":"<div><div>Solitons find widespread applications across diverse disciplines. Accurate identification of the internal parameters within soliton systems allows us for precise comprehension and effective regulation of these systems. The introduction of deep learning has revolutionized the way to address the issue of parameter identification in soliton systems. However, the lack of suitable weight initialization schemes leads to the identification outcomes being prone to blurriness and errors. Consequently, we propose a novel initialization method: physical meta-learning(PML). The unique approach which relies solely on the physical information related to the system allows us to obtain the initialization weights without relying on any labeled data. In basic soliton systems experiments, PML reduces the identification error by 25% to 80%. Regarding the parameter identification task of dissipative soliton system in mode-locked lasers, the PML method significantly reduces the identification error by 98.1%. In addition to the application scenarios, we also examine the effectiveness of the PML method in different parameter identification methods. Overall, our research provides a method for optimizing the identification and simulation of complex soliton systems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117852"},"PeriodicalIF":6.9,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143508337","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":"Implicit stabilized non-ordinary state-based peridynamics for finite deformation and fracture analysis of nearly incompressible materials","authors":"Chengxuan Li,&nbsp;Hanbo Zhang,&nbsp;Cunliang Pan,&nbsp;Hongfei Ye,&nbsp;Hongwu Zhang,&nbsp;Yonggang Zheng","doi":"10.1016/j.cma.2025.117879","DOIUrl":"10.1016/j.cma.2025.117879","url":null,"abstract":"<div><div>This work presents an implicit stabilized non-ordinary state-based peridynamic method (ImNSPD) to simulate the finite deformation and crack propagation in nearly incompressible hyperelastic materials. Firstly, a mixed displacement-pressure (<span><math><mi>u</mi></math></span><strong>-</strong><span><math><mi>p</mi></math></span>) formulation of the ImNSPD is derived to mitigate the volumetric locking issues expected in a purely displacement formulation near the incompressibility limit. Additionally, the shape-associated micromodulus and the dilatation-associated micromodulus are introduced as penalty factors in the modified force vector state and the modified hydrostatic pressure scalar state, respectively, to enhance the numerical stabilization. Subsequently, an incremental-iterative solution scheme based on the Newton-Raphson algorithm and the Newmark-beta method is proposed to capture the nonlinear response of hyperelastic materials in the time domain. The stability and robustness of the ImNSPD are assessed by studying several representative numerical examples involving the finite deformation of nearly incompressible hyperelastic material models. Moreover, an equivalent strain function is introduced as a failure criterion for nearly incompressible hyperelastic materials, and the good performance of the ImNSPD in predicting crack propagation is demonstrated through the fracture analysis of the twisting column test.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117879"},"PeriodicalIF":6.9,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143508449","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 coupled clustering FFT2 multiscale method for modeling the nonlinear behavior and failure of composites","authors":"Menglei Li ,&nbsp;Marco Magri ,&nbsp;Bing Wang ,&nbsp;Bing Wang","doi":"10.1016/j.cma.2025.117854","DOIUrl":"10.1016/j.cma.2025.117854","url":null,"abstract":"<div><div>We propose a novel FFT² parallel multiscale computational method to predict the nonlinear behavior and failure of composite materials. Unlike traditional multiscale methods, the proposed approach reformulates the mechanical boundary value problem into Lippmann-Schwinger type integral equations at both the micro- and macro-scale, thereby leveraging the numerical efficiency of the fast Fourier transform (FFT) method at both scales. The application of generic (e.g. non-periodic) boundary conditions at the macro-scale is carried out by using the virtual boundary technique and buffer zones. In addition, the introduction of a clustering algorithm further improves the computational efficiency of the numerical method during the information transfer between scales. To ensure accurate damage prediction and mitigate spurious strain localization at both scales, suitable regularization techniques are employed. The proposed multiscale method is applied to investigate the transverse tension of unidirectional composite dog-bone specimens. After experimental verification, the method is applied to simulate 2D and 3D brittle fracture, elasto-plastic damage, and examples with non-uniform material orientation. The results demonstrate the robustness and adaptability of the clustering approach, which achieves up to 65.90-fold speedup and 81.62-fold reduction in memory usage compared to non-clustered multiscale methods, while maintaining a comparable level of accuracy.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117854"},"PeriodicalIF":6.9,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143508450","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":"Hybrid quantum genetic algorithm for structural damage identification","authors":"Lianming Xu,&nbsp;Xiaojun Wang,&nbsp;Zhenghuan Wang,&nbsp;Geyong Cao","doi":"10.1016/j.cma.2025.117866","DOIUrl":"10.1016/j.cma.2025.117866","url":null,"abstract":"<div><div>Quantum computing, as an emerging technology, has experienced rapid development and attracted widespread attention across various disciplines over the past four decades. In the field of damage identification, as the demand for damage detection efficiency in real-world engineering continues to rise, the application of this fast computing technology also deserves attention. As a combination of genetic algorithm and quantum computing, the hybrid quantum genetic algorithm (HQGA) utilizes quantum superposition and entanglement to implement evolutionary optimization on quantum computers. Due to its flexibility in performing a global adaptive search without relying on the specific nature of the problem or the form of the objective function, the HQGA demonstrates the potential to solve complex optimization problems. In this paper, a mapping relationship is established between quantum bits and the stiffness reduction factors of elements, leading to the development of the HQGA for solving the damage identification problem based on measured strain responses. The performance of the proposed method is tested by damage identification for three distinct structures. The results show that the proposed method can enhance detection efficiency greatly while maintaining its effectiveness and stability.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117866"},"PeriodicalIF":6.9,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143510961","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 nonlocal mixed-mode fatigue crack growth model based on peridynamic differential operator theory","authors":"Jianrui Liu ,&nbsp;Junxiang Wang ,&nbsp;Zhaobo Song ,&nbsp;Liang Wang","doi":"10.1016/j.cma.2025.117855","DOIUrl":"10.1016/j.cma.2025.117855","url":null,"abstract":"<div><div>This study presents a novel peridynamics (PD) fatigue model for the fatigue crack growth analysis under mixed-mode loading conditions. The foundational aspect of this work involves the application of Peridynamic Differential Operator (PDDO) theory, based on which the analytical relationships between the non-local bond deformations and local strain/stress tensors are first established with the consideration of bond rotation kinematics. Furthermore, the correlations between the bond stretch and Stress Intensity Factors (SIFs) within the crack tip field are rigorously derived, which facilitates the description of fatigue damage in alignment with the classical Linear Elastic Fracture Mechanics (LEFM) theory. The PD fatigue model is implemented through a coupled PDDO and finite element (FE) approach to achieve higher numerical efficiency. Finally, the model's validity is demonstrated through high-fidelity simulation of several benchmark mixed-mode fatigue examples. A notable advantage of the proposed PD fatigue model is its seamless integration of peridynamic theory with classical fracture mechanics, and the model parameters can be rigorously and accurately calibrated for mixed-mode fatigue problems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117855"},"PeriodicalIF":6.9,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143488381","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}