Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations","authors":"Daniel O’Shea,&nbsp;Xiaoran Zhang,&nbsp;Shayan Mohammadian,&nbsp;Chongmin Song","doi":"10.1016/j.cma.2025.117831","DOIUrl":"10.1016/j.cma.2025.117831","url":null,"abstract":"<div><div>An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the authors on elastodynamics by presenting a new algorithm that eliminates the need for factorization of the mass matrix providing benefit for the solution of nonlinear problems. The improved algorithm directly obtains the acceleration at the same order of accuracy of the displacement and velocity using vector operations (without additional equation solutions). The nonlinearity is handled by numerical integration within a time step to achieve the desired order of accuracy. The new algorithm fully retains the desirable features of the previous works: 1. The order of accuracy is not affected by the presence of external forces and physical damping; 2. The amount of numerical dissipation in the algorithm is controlled by a user-specified parameter <span><math><msub><mrow><mi>ρ</mi></mrow><mrow><mi>∞</mi></mrow></msub></math></span>, leading to schemes ranging from perfectly nondissipative <span><math><mi>A</mi></math></span>-stable to <span><math><mi>L</mi></math></span>-stable; 3. The effective stiffness matrix is a linear combination of the mass, damping, and stiffness matrices as in the trapezoidal rule, leading to high efficiency for large-scale problems. The proposed algorithm, with its elegance and computational advantages, is shown to replicate the numerical results demonstrated on linear problems in previous works. Additional numerical examples of linear and nonlinear vibration and wave propagation are presented herein. Notably, the proposed algorithms show the same convergence rates for nonlinear problems as linear problems, and very high accuracy. It was found that second-order time integration methods commonly used in commercial software produce significantly polluted acceleration responses for a common class of wave propagation problems. The high-order time integration schemes presented here perform noticeably better at suppressing spurious high-frequency oscillations and producing reliable and usable acceleration responses. The source code written in <span>MATLAB</span> is available for download at: <span><span>https://github.com/ChongminSong/HighOrderTimeIngt_PartialFraction.git</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117831"},"PeriodicalIF":6.9,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143422193","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":"Global-local adaptive meshing method for phase-field fracture modeling","authors":"FengYu Cheng,&nbsp;Hao Yu,&nbsp;Quan Wang,&nbsp;HanWei Huang,&nbsp;WenLong Xu,&nbsp;HengAn Wu","doi":"10.1016/j.cma.2025.117846","DOIUrl":"10.1016/j.cma.2025.117846","url":null,"abstract":"<div><div>This work develops a global-local adaptive meshing method for the phase-field model of brittle fracture, offering flexible adjustment of mesh density to produce seamless and high-quality adaptive meshes. The method first establishes a direct mapping from phase-field values and displacement errors to a normalized nodal density field, which is used to control the computational accuracy. On this basis, a sampling procedure is performed by detecting the maximum value to progressively place sampling nodes, ensuring that first-level nodes are placed globally while preserving crack location information. Subsequently, a hexagonal seeding algorithm is used to multiply nodes, where the spacing of generated seeds (i.e., higher-level nodes) is adaptively adjusted based on local nodal density requirements to regulate element sizes. A spatial assessment algorithm is utilized to compare the expected nodal spacing of the newly generated node with its distance to existing nodes, which serves as a termination criterion for the loop of the seeding algorithm and effectively prevents the occurrence of low-quality elements. After the seeding process of all nodes is completed, all generated nodes are connected by constrained Delaunay triangulation. This method has been discussed under classical brittle fracture cases with various control parameters (e.g., the mapping function, the expected maximum/minimum element size, and the distance factor) to validate its advantage of reducing degrees of freedom and improving solution efficiency.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117846"},"PeriodicalIF":6.9,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143422191","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":"Prediction of evolving plasticity in rails under steady state rolling contact based on Reduced-Order Modeling","authors":"Caroline Ansin,&nbsp;Fredrik Larsson,&nbsp;Ragnar Larsson","doi":"10.1016/j.cma.2025.117828","DOIUrl":"10.1016/j.cma.2025.117828","url":null,"abstract":"<div><div>Predicting railhead damage due to multiple wheel passes in railway operations can be computationally demanding, especially when accounting for the rail’s inelastic material response. In this paper, we use a steady-state assumption within a convective coordinate system and employ Reduced-Order Modeling (ROM) through the Proper Generalized Decomposition (PGD) method to increase the computational efficiency. Our approach solves a nonlinear reduced-order problem for the displacements in a 3D railhead with elastic–plastic material properties considering various contact scenarios. The ROM framework includes domain decomposition and a parametric loading framework using PGD. It accounts for the elastic–plastic rail material through three key features: (1) treatment of the moving load under the assumption of steady state, by using a convective coordinate system along the railhead to convert the problem into a stationary contact load problem, (2) implementation of PGD to solve the 3D displacement field efficiently, and (3) use of fixed-point iterations to treat the coupling for solving plastic strains and displacements. In this iterative process, plastic strains are solved from displacements, and displacements are solved based on a loading scenario and updated plastic strains. The accuracy and computational efficiency are assessed by comparing our strategy with 3D finite element simulations for moving contact loads. The results show convergence with only a few fixed-point iterations for each over rolling, which results in a solution that is 63 times faster. This efficiency is crucial for assessing the accumulated plastic deformation from multiple wheel passes.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117828"},"PeriodicalIF":6.9,"publicationDate":"2025-02-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143422224","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":"Stabilization-free Virtual Element Method for 2D second order elliptic equations","authors":"Stefano Berrone,&nbsp;Andrea Borio,&nbsp;Davide Fassino,&nbsp;Francesca Marcon","doi":"10.1016/j.cma.2025.117839","DOIUrl":"10.1016/j.cma.2025.117839","url":null,"abstract":"<div><div>In this work, we present and analyze a Stabilization-free Virtual Element high order scheme for 2D second order elliptic equation. This method is characterized by the definition of new polynomial projections that allow the definition of structure-preserving schemes. We provide a necessary and sufficient condition on the polynomial projection space that ensure the well-posedness of the scheme and we derive optimal a priori error estimates. Several numerical tests assess the stability of the method and the robustness in solving problems characterized by anisotropies.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117839"},"PeriodicalIF":6.9,"publicationDate":"2025-02-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143422225","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":"An inter-system iteration technique in the whole time domain implemented by the coupling of FEM and TD-BEM","authors":"Xiaofei Qin ,&nbsp;Hongjun Li ,&nbsp;Changnv Zeng","doi":"10.1016/j.cma.2025.117840","DOIUrl":"10.1016/j.cma.2025.117840","url":null,"abstract":"<div><div>The paper introduces an innovative coupling technique, termed the Inter-System Iteration Technique (ISIT) in the whole time domain, which is effectively implemented by the coupling of the Finite Element Method (FEM) and Time Domain Boundary Element Method (TD-BEM) for dynamic analysis of intricate engineering problems. This approach divides the entire system into distinct FEM and TD-BEM parts, computed separately across the whole time domain using an optimal solution strategy. Its notable strength lies in achieving exceptional computational efficiency and accuracy by minimizing iterative calculations for each individual time step and employing iterative computations between the FEM and TD-BEM sub-systems. Notably, this method offers versatility, adaptable to various numerical algorithms and facilitating straightforward computation with commercial structural analysis software. The implemented ISIT in the whole time domain can also be used to couple arbitrarily two numerical algorithms. The paper validates its efficacy through two numerical applications, affirming its high computational efficiency and potential for addressing complex engineering challenges.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117840"},"PeriodicalIF":6.9,"publicationDate":"2025-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143418495","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":"Machine learning-accelerated peridynamics model for mechanical and failure behaviors of materials","authors":"Jiasheng Huang ,&nbsp;J.X. Liew ,&nbsp;Binbin Yin ,&nbsp;K.M. Liew","doi":"10.1016/j.cma.2025.117826","DOIUrl":"10.1016/j.cma.2025.117826","url":null,"abstract":"<div><div>Computational mechanics is essential for understanding and predicting complex material behaviors, particularly in areas such as material fracture mechanics and structural engineering. However, the high computational costs associated with traditional methods, especially for large-scale simulations, present significant challenges. Peridynamics (PD) offers a compelling alternative to classical continuum mechanics by effectively modeling discontinuities such as cracks. Despite its strengths, PD is computationally intensive, limiting its broader application. To address these challenges, we introduce a machine learning-accelerated PD model that significantly reduces computational time while maintaining high accuracy. Our method integrates a machine learning-based surrogate model trained on displacement field data, which efficiently approximates the behaviors of material points, bypassing the iterative processes of conventional PD simulations. This approach is validated through a series of benchmark tests, ranging from one-dimensional bars to three-dimensional beams, demonstrating speedups of over six times compared to traditional methods. The integration of machine learning with PD not only enhances computational efficiency but also expands the practical applicability of PD to large-scale engineering problems, making it a viable tool for a wide range of scientific and industrial applications.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117826"},"PeriodicalIF":6.9,"publicationDate":"2025-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143418496","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 discrete fiber dispersion model with octahedral symmetry quadrature for mechanical analyses of skin corrective surgeries","authors":"Riccardo Alberini ,&nbsp;Michele Terzano ,&nbsp;Gerhard A. Holzapfel ,&nbsp;Andrea Spagnoli","doi":"10.1016/j.cma.2025.117809","DOIUrl":"10.1016/j.cma.2025.117809","url":null,"abstract":"<div><div>Advanced simulations of the mechanical behavior of soft tissues frequently rely on structure-based constitutive models, including smeared descriptions of collagen fibers. Among them, the so-called Discrete Fiber Dispersion (DFD) modeling approach is based on a discrete integration of the fiber-strain energy over all the fiber directions. In this paper, we review the theoretical framework of the DFD model, including a derivation of the stress and stiffness tensors required for the finite element implementation. Specifically, their expressions for incompressible plane stress problems are obtained. The use of a Lebedev quadrature, built exploiting the octahedral symmetry, is then proposed, illustrating the particular choice adopted for the orientation of the integration points. Next, the convergence of this quadrature scheme is assessed by means of three numerical benchmark tests, highlighting the advantages with respect to other angular integration methods available in the literature. Finally, using the implemented model, we analyze the mechanical properties of the Z-plasty, a technique commonly used in reconstructive skin surgery, considering multiple geometrical configurations, orientations of the fibers, and levels of skin prestress. The results are presented in the form of mechanical quantities relevant to surgical practice.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117809"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386409","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":"A hybrid isogeometric collocation method on implicitly trimmed domains","authors":"Jingjing Yang ,&nbsp;Pei Zhou ,&nbsp;Lin Lan ,&nbsp;Chun-Gang Zhu","doi":"10.1016/j.cma.2025.117812","DOIUrl":"10.1016/j.cma.2025.117812","url":null,"abstract":"<div><div>We propose a novel isogeometric collocation method (IGA-C) for trimmed domains, using weighted extended B-splines (WEB-splines). Our approach employs the implicit representations of the trimming boundaries to construct the weighted basis, which allows for the subsequent calculations based on a parametrization over a single tensor-product patch, despite the nontrivial shape of the domain. The stabilization of the weighted basis is accomplished by means of extension. We present the classification criterion and the calculation procedure for the extension coefficients of the inner B-splines. The utilization of WEB-splines in analysis enables a natural application of the Dirichlet boundary conditions. Additionally, we adopt a hybrid collocation–Galerkin approach to impose the Neumann boundary conditions on the trimming boundaries. Our proposed method combines the advantages of WEB-splines and IGA-C in terms of straightforward implementation, well-conditioned system matrices and high computational efficiency, which we illustrate by numerical tests on 2D and 3D trimmed geometries. The numerical results further demonstrate that our methodology guarantees the same convergence rates as IGA-C.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117812"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386412","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":"Multiscale topology optimization of functionally graded lattice structures based on physics-augmented neural network material models","authors":"Jonathan Stollberg ,&nbsp;Tarun Gangwar ,&nbsp;Oliver Weeger ,&nbsp;Dominik Schillinger","doi":"10.1016/j.cma.2025.117808","DOIUrl":"10.1016/j.cma.2025.117808","url":null,"abstract":"<div><div>We present a new framework for the simultaneous optimization of both the topology as well as the relative density grading of cellular structures and materials, also known as lattices. Due to manufacturing constraints, the optimization problem falls into the class of mixed-integer nonlinear programming problems. Since no algorithm is capable of solving these problems in polynomial time, we obtain a relaxed problem from a multiplicative split of the relative density and a penalization approach. The sensitivities of the objective function are derived such that any gradient-based solver might be applied for the iterative update of the design variables. In a next step, we introduce a material model that is parametric in the design variables of interest and suitable to describe the isotropic deformation behavior of quasi-stochastic lattices. For that, we derive and implement further physical constraints and enhance a physics-augmented neural network from the literature that was formulated initially for rhombic materials. Finally, to illustrate the applicability of the method, we incorporate the material model into our computational framework and exemplary optimize two-and three-dimensional benchmark structures as well as a complex aircraft component.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117808"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386408","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":"Improving hp-variational physics-informed neural networks for steady-state convection-dominated problems","authors":"Thivin Anandh ,&nbsp;Divij Ghose ,&nbsp;Himanshu Jain ,&nbsp;Pratham Sunkad ,&nbsp;Sashikumaar Ganesan ,&nbsp;Volker John","doi":"10.1016/j.cma.2025.117797","DOIUrl":"10.1016/j.cma.2025.117797","url":null,"abstract":"<div><div>This paper proposes and studies two extensions of applying hp-variational physics-informed neural networks, more precisely the FastVPINNs framework, to convection-dominated convection–diffusion–reaction problems. First, a term in the spirit of a SUPG stabilization is included in the loss functional and a network architecture is proposed that predicts spatially varying stabilization parameters. Having observed that the selection of the indicator function in hard-constrained Dirichlet boundary conditions has a big impact on the accuracy of the computed solutions, the second novelty is the proposal of a network architecture that learns good parameters for a class of indicator functions. Numerical studies show that both proposals lead to noticeably more accurate results than approaches that can be found in the literature.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117797"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386407","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}