首页 > 最新文献

Computer Methods in Applied Mechanics and Engineering最新文献

英文 中文
Real-time forecasting of chaotic dynamics from sparse data and autoencoders 基于稀疏数据和自编码器的混沌动态实时预测
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-13 DOI: 10.1016/j.cma.2025.118600
Elise Özalp , Andrea Nóvoa , Luca Magri
The real-time prediction of chaotic systems requires a nonlinear-reduced order model (ROM) to forecast the dynamics, and a stream of data from sensors to update the ROM. Data-driven ROMs are typically built with a two-step strategy: data compression onto a lower-dimensional latent space, and prediction of the temporal dynamics on it. Although these methods have proven effective, to achieve real-time prediction, there are two key challenges to overcome: (i) ROMs of chaotic systems can become numerically unstable; and (ii) sensors’ data are sparse, i.e., partial, and noisy. To overcome these challenges, we propose an augmented sequential data assimilation (DA) framework based on the Ensemble Kalman filter (EnKF) that updates the latent state of ROM by assimilating noisy and sparse measurements. We demonstrate the proposed DA-ROM framework using a ROM that consists of a convolutional autoencoder (CAE) to compresses the system’s state onto a lower-dimensional latent space, and an echo state network (ESN) formulated as a state-space model to forecast the temporal evolution on the latent space. The DA-CAE-ESN provides a numerically stable and real-time adaptive ROM. The DA-CAE-ESN is tested on spatio-temporally chaotic partial differential equations: the Kuramoto–Sivashinsky equation, and a two-dimensional Navier-Stokes equation (Kolmogorov flow). We show that the method provides accurate and stable forecasts across different levels of noise, sparsity, and sampling rates. As a by-product, the DA-CAE-ESN acts as a localization strategy that mitigates spurious correlations, which arise when applying the EnKF to high-dimensional systems. The DA-CAE-ESN provides a numerically stable method to perform real-time predictions, which opens opportunities for deploying data-driven latent models.
混沌系统的实时预测需要一个非线性降阶模型(ROM)来预测动力学,并需要来自传感器的数据流来更新ROM。数据驱动的ROM通常采用两步策略构建:将数据压缩到低维潜在空间,并预测其上的时间动力学。虽然这些方法已被证明是有效的,但要实现实时预测,有两个关键挑战需要克服:(i)混沌系统的rom可能变得数值不稳定;(ii)传感器的数据是稀疏的,即部分的和有噪声的。为了克服这些挑战,我们提出了一种基于集成卡尔曼滤波器(EnKF)的增强顺序数据同化(DA)框架,该框架通过同化噪声和稀疏测量来更新ROM的潜在状态。我们使用由卷积自编码器(CAE)组成的ROM来将系统状态压缩到低维潜在空间,并使用回声状态网络(ESN)作为状态空间模型来预测潜在空间上的时间演变,从而演示了所提出的DA-ROM框架。DA-CAE-ESN提供了一个数值稳定的实时自适应ROM。DA-CAE-ESN在时空混沌偏微分方程:Kuramoto-Sivashinsky方程和二维Navier-Stokes方程(Kolmogorov流)上进行了测试。我们表明,该方法在不同的噪声水平、稀疏度和采样率下提供了准确和稳定的预测。作为副产品,DA-CAE-ESN作为一种定位策略,可以减轻在将EnKF应用于高维系统时产生的虚假相关性。DA-CAE-ESN提供了一种数值稳定的方法来执行实时预测,这为部署数据驱动的潜在模型提供了机会。
{"title":"Real-time forecasting of chaotic dynamics from sparse data and autoencoders","authors":"Elise Özalp ,&nbsp;Andrea Nóvoa ,&nbsp;Luca Magri","doi":"10.1016/j.cma.2025.118600","DOIUrl":"10.1016/j.cma.2025.118600","url":null,"abstract":"<div><div>The real-time prediction of chaotic systems requires a nonlinear-reduced order model (ROM) to forecast the dynamics, and a stream of data from sensors to update the ROM. Data-driven ROMs are typically built with a two-step strategy: data compression onto a lower-dimensional latent space, and prediction of the temporal dynamics on it. Although these methods have proven effective, to achieve real-time prediction, there are two key challenges to overcome: (i) ROMs of chaotic systems can become numerically unstable; and (ii) sensors’ data are sparse, i.e., partial, and noisy. To overcome these challenges, we propose an augmented sequential data assimilation (DA) framework based on the Ensemble Kalman filter (EnKF) that updates the latent state of ROM by assimilating noisy and sparse measurements. We demonstrate the proposed DA-ROM framework using a ROM that consists of a convolutional autoencoder (CAE) to compresses the system’s state onto a lower-dimensional latent space, and an echo state network (ESN) formulated as a state-space model to forecast the temporal evolution on the latent space. The DA-CAE-ESN provides a numerically stable and real-time adaptive ROM. The DA-CAE-ESN is tested on spatio-temporally chaotic partial differential equations: the Kuramoto–Sivashinsky equation, and a two-dimensional Navier-Stokes equation (Kolmogorov flow). We show that the method provides accurate and stable forecasts across different levels of noise, sparsity, and sampling rates. As a by-product, the DA-CAE-ESN acts as a localization strategy that mitigates spurious correlations, which arise when applying the EnKF to high-dimensional systems. The DA-CAE-ESN provides a numerically stable method to perform real-time predictions, which opens opportunities for deploying data-driven latent models.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118600"},"PeriodicalIF":7.3,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145753448","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}
引用次数: 0
Large-scale topology optimisation of time-dependent thermal conduction using space-time finite elements and a parallel space-time multigrid preconditioner 基于时空有限元和平行时空多网格预调节器的大规模时变热传导拓扑优化
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-13 DOI: 10.1016/j.cma.2025.118605
Joe Alexandersen, Magnus Appel
This paper presents a novel space-time topology optimisation framework for time-dependent thermal conduction problems, aiming to significantly reduce the time-to-solution. By treating time as an additional spatial dimension, we discretise the governing equations using a stabilised continuous Galerkin space-time finite element method. The resulting large all-at-once system is solved using an iterative Krylov solver preconditioned with a parallel space-time multigrid method employing a semi-coarsening strategy. Implemented in a fully parallel computing framework, the method yields a parallel-in-time method that demonstrates excellent scalability on a distributed-memory supercomputer, solving problems up to 4.2 billion degrees of freedom. Comparative studies show up to 52 ×  speed-up over traditional time-stepping approaches, with only moderate increases in total computational cost in terms of core-hours. The framework is validated on benchmark problems with both time-constant and time-varying designs, and its flexibility is demonstrated through variations in material properties. These results establish the proposed space-time method as a promising approach for large-scale time-dependent topology optimisation in thermal applications.
本文提出了一种新的时变热传导问题的时空拓扑优化框架,旨在显著缩短求解时间。通过将时间作为一个额外的空间维度,我们使用稳定连续伽辽金时空有限元方法离散控制方程。采用半粗化策略的平行时空多重网格法作为前置条件,利用迭代Krylov求解器求解得到的大型一次性系统。该方法在完全并行计算框架中实现,产生了一种并行实时方法,在分布式内存超级计算机上展示了出色的可扩展性,可解决高达42亿自由度的问题。比较研究表明,与传统的时间步进方法相比,加速高达52 × ,而以核心小时计算的总计算成本仅略有增加。该框架在时间常数和时变设计的基准问题上得到了验证,并通过材料性能的变化证明了其灵活性。这些结果表明,所提出的时空方法是热应用中大规模时变拓扑优化的一种有前途的方法。
{"title":"Large-scale topology optimisation of time-dependent thermal conduction using space-time finite elements and a parallel space-time multigrid preconditioner","authors":"Joe Alexandersen,&nbsp;Magnus Appel","doi":"10.1016/j.cma.2025.118605","DOIUrl":"10.1016/j.cma.2025.118605","url":null,"abstract":"<div><div>This paper presents a novel space-time topology optimisation framework for time-dependent thermal conduction problems, aiming to significantly reduce the time-to-solution. By treating time as an additional spatial dimension, we discretise the governing equations using a stabilised continuous Galerkin space-time finite element method. The resulting large all-at-once system is solved using an iterative Krylov solver preconditioned with a parallel space-time multigrid method employing a semi-coarsening strategy. Implemented in a fully parallel computing framework, the method yields a parallel-in-time method that demonstrates excellent scalability on a distributed-memory supercomputer, solving problems up to 4.2 billion degrees of freedom. Comparative studies show up to 52 ×  speed-up over traditional time-stepping approaches, with only moderate increases in total computational cost in terms of core-hours. The framework is validated on benchmark problems with both time-constant and time-varying designs, and its flexibility is demonstrated through variations in material properties. These results establish the proposed space-time method as a promising approach for large-scale time-dependent topology optimisation in thermal applications.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118605"},"PeriodicalIF":7.3,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731701","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}
引用次数: 0
Learning a hyperelastic constitutive model from 3D experimental data 从三维实验数据中学习超弹性本构模型
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-12 DOI: 10.1016/j.cma.2025.118592
M. Bourdyot , M. Compans , R. Langlois , B. Smaniotto , E. Baranger , C. Jailin
A hyperelastic constitutive law is experimentally learned and validated from three-dimensional full-field kinematic data. 3D-printed thermoplastic polyurethane specimens were subjected to monotonic uniaxial tension within a laboratory micro-computed tomography system. Global Digital Volume Correlation yields volumetric displacement fields, reaching up to 35 % axial strain. A Physics-Augmented Neural Network (PANN) architecture, embedding polyconvexity constraints through an input-convex neural network, is trained in an unsupervised manner by enforcing mechanical equilibrium via the EUCLID loss, using only measured displacements and global reaction forces. This three-dimensional formulation eliminates the need for depth-related assumptions inherent in 2D measurements and provides improved representation of boundary conditions, which is critical for equilibrium-based training. For comparison, Neo-Hookean and Saint-Venant-Kirchhoff models are identified in parallel. The PANN achieves superior equilibrium satisfaction, outperforming both linear and Neo-Hookean benchmarks. In addition to the original specimen, a second experiment on a different sample was performed under comparable loading conditions. This configuration provides an independent validation dataset acquired under slightly varying experimental conditions, thereby testing the robustness and transferability of the learned constitutive model. This work presents the first experimental demonstration of unsupervised PANN model training and establishes a practical protocol for data-based hyperelastic characterization from volumetric measurements.
通过实验学习并验证了三维全场运动数据的超弹性本构规律。3d打印的热塑性聚氨酯样品在实验室微计算机断层扫描系统中受到单调单轴张力。全球数字体积相关产生体积位移场,达到35%的轴向应变。物理增强神经网络(PANN)架构通过输入-凸神经网络嵌入多凸约束,通过EUCLID损失强制机械平衡,仅使用测量的位移和全局反作用力,以无监督的方式进行训练。这种三维公式消除了对二维测量中固有的深度相关假设的需要,并提供了改进的边界条件表示,这对于基于平衡的训练至关重要。为了比较,新hookean和Saint-Venant-Kirchhoff模型是并行识别的。PANN实现了卓越的均衡满意度,优于线性和新胡克基准。除原始试样外,在相同的加载条件下对不同的试样进行了第二次实验。这种配置提供了在略有变化的实验条件下获得的独立验证数据集,从而测试了学习的本构模型的鲁棒性和可移植性。这项工作提出了无监督PANN模型训练的第一个实验演示,并建立了一个基于数据的基于体积测量的超弹性表征的实用协议。
{"title":"Learning a hyperelastic constitutive model from 3D experimental data","authors":"M. Bourdyot ,&nbsp;M. Compans ,&nbsp;R. Langlois ,&nbsp;B. Smaniotto ,&nbsp;E. Baranger ,&nbsp;C. Jailin","doi":"10.1016/j.cma.2025.118592","DOIUrl":"10.1016/j.cma.2025.118592","url":null,"abstract":"<div><div>A hyperelastic constitutive law is experimentally learned and validated from three-dimensional full-field kinematic data. 3D-printed thermoplastic polyurethane specimens were subjected to monotonic uniaxial tension within a laboratory micro-computed tomography system. Global Digital Volume Correlation yields volumetric displacement fields, reaching up to 35 % axial strain. A Physics-Augmented Neural Network (PANN) architecture, embedding polyconvexity constraints through an input-convex neural network, is trained in an unsupervised manner by enforcing mechanical equilibrium via the EUCLID loss, using only measured displacements and global reaction forces. This three-dimensional formulation eliminates the need for depth-related assumptions inherent in 2D measurements and provides improved representation of boundary conditions, which is critical for equilibrium-based training. For comparison, Neo-Hookean and Saint-Venant-Kirchhoff models are identified in parallel. The PANN achieves superior equilibrium satisfaction, outperforming both linear and Neo-Hookean benchmarks. In addition to the original specimen, a second experiment on a different sample was performed under comparable loading conditions. This configuration provides an independent validation dataset acquired under slightly varying experimental conditions, thereby testing the robustness and transferability of the learned constitutive model. This work presents the first experimental demonstration of unsupervised PANN model training and establishes a practical protocol for data-based hyperelastic characterization from volumetric measurements.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118592"},"PeriodicalIF":7.3,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731144","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}
引用次数: 0
Invariant-manifold-based model reduction for geometrically exact beam dynamics 基于不变流形的几何精确梁动力学模型约简
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-12 DOI: 10.1016/j.cma.2025.118643
Yifan Qi , Shilei Han , Minghe Shan , Mingwu Li , Qiang Tian
This paper presents an invariant-manifold (IM)-based reduction approach for geometrically exact beams undergoing finite rigid-body motions and elastic deformations. For a free beam or beam assembly, direct application of standard IM reduction is hindered by infinite system equilibria and secular growth in physical coordinates. To address this challenge, the floating frame of reference (FFR) method is employed to decompose the beam kinematics into a fictitious rigid-body motion and elastic deformation relative to the FFR. Furthermore, externally applied loads are recast as ordinary differential equations for an excitation vector. These steps render the system dynamics independent of the FFR motion tensor and the time variable, yielding an autonomous system with all parameters bounded in the vicinity of the state-space origin. Polynomial parameterization of the target invariant manifold is enabled using the FFR velocity, a set of dominant modal coordinates, and the excitation vector. A hierarchy of cohomological equations is formulated to compute the manifold mapping and the reduced dynamics. An efficient solution strategy is proposed to solve the high-dimensional cohomological equations, where higher-order problems are separated into resonant and non-resonant parts: the resonant part, associated with small-scale singular linear systems, is solved using least-squares minimization, while the non-resonant part is solved sequentially through a hierarchy of linear systems. Numerical results show that the proposed IM-based reduction yields highly accurate predictions of the nonlinear dynamic response of geometrically exact beams using only third–order polynomial expansions of the manifold and reduced dynamics. Moreover, the method achieves substantial computational savings, as dynamic simulation requires integrating only the low-dimensional reduced-order equations rather than the full state dynamics.
本文提出了一种基于不变流形(IM)的几何精确梁在有限刚体运动和弹性变形下的约简方法。对于一个自由梁或梁组合,直接应用标准的IM减少阻碍了无限系统平衡和长期增长的物理坐标。为了解决这一问题,采用浮动参照系(FFR)方法将梁的运动学分解为虚拟的刚体运动和相对于FFR的弹性变形。此外,外部施加的载荷被重铸为激励矢量的常微分方程。这些步骤使系统动力学独立于FFR运动张量和时间变量,从而产生一个所有参数都在状态空间原点附近的自治系统。使用FFR速度、一组主导模态坐标和激励向量实现目标不变流形的多项式参数化。建立了一套上同调方程来计算流形映射和简化动力学。提出了一种求解高维上同调方程的有效策略,其中高阶问题分为谐振部分和非谐振部分:谐振部分与小规模奇异线性系统相关,采用最小二乘最小化方法求解,而非谐振部分则通过线性系统的层次顺序求解。数值结果表明,所提出的基于im的约简方法仅使用流形和约简动力学的三阶多项式展开,就能高精度地预测几何精确梁的非线性动力响应。此外,由于动态模拟只需要积分低维降阶方程,而不需要积分全状态动力学,因此该方法节省了大量的计算量。
{"title":"Invariant-manifold-based model reduction for geometrically exact beam dynamics","authors":"Yifan Qi ,&nbsp;Shilei Han ,&nbsp;Minghe Shan ,&nbsp;Mingwu Li ,&nbsp;Qiang Tian","doi":"10.1016/j.cma.2025.118643","DOIUrl":"10.1016/j.cma.2025.118643","url":null,"abstract":"<div><div>This paper presents an invariant-manifold (IM)-based reduction approach for geometrically exact beams undergoing finite rigid-body motions and elastic deformations. For a free beam or beam assembly, direct application of standard IM reduction is hindered by infinite system equilibria and secular growth in physical coordinates. To address this challenge, the floating frame of reference (FFR) method is employed to decompose the beam kinematics into a fictitious rigid-body motion and elastic deformation relative to the FFR. Furthermore, externally applied loads are recast as ordinary differential equations for an excitation vector. These steps render the system dynamics independent of the FFR motion tensor and the time variable, yielding an autonomous system with all parameters bounded in the vicinity of the state-space origin. Polynomial parameterization of the target invariant manifold is enabled using the FFR velocity, a set of dominant modal coordinates, and the excitation vector. A hierarchy of cohomological equations is formulated to compute the manifold mapping and the reduced dynamics. An efficient solution strategy is proposed to solve the high-dimensional cohomological equations, where higher-order problems are separated into resonant and non-resonant parts: the resonant part, associated with small-scale singular linear systems, is solved using least-squares minimization, while the non-resonant part is solved sequentially through a hierarchy of linear systems. Numerical results show that the proposed IM-based reduction yields highly accurate predictions of the nonlinear dynamic response of geometrically exact beams using only third–order polynomial expansions of the manifold and reduced dynamics. Moreover, the method achieves substantial computational savings, as dynamic simulation requires integrating only the low-dimensional reduced-order equations rather than the full state dynamics.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118643"},"PeriodicalIF":7.3,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145732426","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}
引用次数: 0
Stabilization-free virtual element method for 2D third medium contact 二维第三介质接触的无稳定虚元法
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-12 DOI: 10.1016/j.cma.2025.118611
Bing-Bing Xu, Peter Wriggers
The third medium contact (TMC) has been proven to be an effective approach for simulating contact problems involving large deformations. Unlike traditional contact algorithms, this methodology introduces a third medium between two contacting bodies, thereby avoiding the complex treatment of the contact constraints. The approach has been successfully applied for different problems in the framework of the finite element method (FEM). As a generalization of the finite element method, the virtual element method (VEM) can handle arbitrary polygonal elements, providing greater flexibility for modeling third medium contact. However, due to the introduction of the projection operator, VEM requires additional stabilization terms to control the rank of the stiffness matrix. Moreover, the regularization term in the third medium contact formulation requires a second-order numerical scheme, which further complicates the application of classical virtual element schemes to such problems. In this work, the stabilization-free virtual element method (SFVEM) is adopted to TMC and applied to solve contact problems undergoing large deformations. Different from the traditional VEM, SFVEM does not require additional stabilization terms, which simplifies the construction of necessary regularization terms in third medium contact. Building upon the traditional second-order FEM framework, we present the specific format of SFVEM for solving third medium contact, including the construction of high-order projection operator and the tangent stiffness matrix. Numerical examples are provided to demonstrate the effectiveness and applicability of SFVEM in solving complex 2D contact problems with the TMC approach.
第三介质接触(TMC)已被证明是模拟大变形接触问题的有效方法。与传统的接触算法不同,该方法在两个接触体之间引入了第三种介质,从而避免了接触约束的复杂处理。该方法已成功地应用于有限元框架下的不同问题。虚拟元法作为有限元方法的推广,可以处理任意多边形单元,为第三介质接触建模提供了更大的灵活性。然而,由于引入了投影算子,VEM需要额外的稳定项来控制刚度矩阵的秩。此外,第三介质接触公式中的正则化项需要二阶数值格式,这进一步使经典虚元格式在此类问题中的应用复杂化。本文将无稳定虚元法(SFVEM)应用于TMC,并将其应用于大变形接触问题的求解。与传统VEM不同,SFVEM不需要附加稳定化项,简化了第三次介质接触中必要正则化项的构造。在传统二阶有限元框架的基础上,提出了求解三次介质接触的SFVEM的具体格式,包括构造高阶投影算子和切向刚度矩阵。通过数值算例验证了该方法在求解复杂二维接触问题时的有效性和适用性。
{"title":"Stabilization-free virtual element method for 2D third medium contact","authors":"Bing-Bing Xu,&nbsp;Peter Wriggers","doi":"10.1016/j.cma.2025.118611","DOIUrl":"10.1016/j.cma.2025.118611","url":null,"abstract":"<div><div>The third medium contact (TMC) has been proven to be an effective approach for simulating contact problems involving large deformations. Unlike traditional contact algorithms, this methodology introduces a third medium between two contacting bodies, thereby avoiding the complex treatment of the contact constraints. The approach has been successfully applied for different problems in the framework of the finite element method (FEM). As a generalization of the finite element method, the virtual element method (VEM) can handle arbitrary polygonal elements, providing greater flexibility for modeling third medium contact. However, due to the introduction of the projection operator, VEM requires additional stabilization terms to control the rank of the stiffness matrix. Moreover, the regularization term in the third medium contact formulation requires a second-order numerical scheme, which further complicates the application of classical virtual element schemes to such problems. In this work, the stabilization-free virtual element method (SFVEM) is adopted to TMC and applied to solve contact problems undergoing large deformations. Different from the traditional VEM, SFVEM does not require additional stabilization terms, which simplifies the construction of necessary regularization terms in third medium contact. Building upon the traditional second-order FEM framework, we present the specific format of SFVEM for solving third medium contact, including the construction of high-order projection operator and the tangent stiffness matrix. Numerical examples are provided to demonstrate the effectiveness and applicability of SFVEM in solving complex 2D contact problems with the TMC approach.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118611"},"PeriodicalIF":7.3,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731145","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}
引用次数: 0
An implicit RKPM formulation for large deformations in elastoplastic micropolar media: Framework, stabilization and assessment 弹塑性微极介质中大变形的隐式RKPM公式:框架、稳定和评估
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-11 DOI: 10.1016/j.cma.2025.118582
T. Mader , A. Dummer , P. Gamnitzer , H. Gimperlein , M. Neuner
Lagrangian finite element methods (FEM) can suffer from excessive mesh distortion or entanglements in simulations involving large deformations, which become critical in the modeling of quasi-brittle materials characterized by failure mechanisms involving highly localized deformations such as shear bands. To overcome these limitations, we present an implicit numerical framework based on the Reproducing Kernel Particle Method (RKPM) for the micropolar continuum. The latter provides an effective framework to mitigate pathological mesh dependencies in shear-dominated failure and allows to consider size effects due to its intrinsic length scale. The RK shape functions are formulated in a semi-Lagrangian description based on B-Spline kernel functions. The weak form of the governing equations of the micropolar continuum is then derived and discretized within the RKPM framework. For implicit time integration, algorithmic tangent operators are derived to obtain quadratic convergence within the Newton-Raphson solution procedure. Several stabilization techniques are adopted for the proposed framework. First, stabilized nodal integration (SNI) is employed. In this context, a novel quasi-conforming update scheme for the smoothing domains is proposed, leveraging the full deformation gradient. Second, Naturally Stabilized Nodal Integration (NSNI) is applied by expanding the internal virtual work contributions for both linear and angular momentum into linear Taylor series. As an alternative to NSNI, a sub-domain integration (SDI) is presented, where each smoothing domain is partitioned into smaller sub-domains to improve integration accuracy and stability. Finally, to achieve variationally consistent integration (VCI), a Petrov-Galerkin correction is applied to the gradients of the test functions. The performance of the method is demonstrated by 2D benchmark simulations of Cook’s membrane using an elastic Neo-Hookean micropolar material. In addition, simulations of a plane strain compression test considering perfect plasticity using a Drucker-Prager-based micropolar model are performed. The capabilities of the formulation are assessed by comparing the results with those obtained with the FEM and the Material Point Method (MPM), highlighting its potential for modeling large deformations in quasi-brittle materials.
拉格朗日有限元方法(FEM)在涉及大变形的模拟中可能遭受过度的网格畸变或纠缠,这对于以涉及高度局部变形(如剪切带)的破坏机制为特征的准脆性材料的建模至关重要。为了克服这些限制,我们提出了一个基于再现核粒子法(RKPM)的隐式数值框架。后者提供了一个有效的框架,以减轻剪切主导破坏中的病态网格依赖,并允许考虑由于其固有长度尺度而产生的尺寸效应。RK形状函数是基于b样条核函数的半拉格朗日描述。然后导出微极连续统控制方程的弱形式,并在RKPM框架内进行离散化。对于隐式时间积分,导出了切线算子算法,在Newton-Raphson解过程中获得二次收敛。该框架采用了几种稳定技术。首先,采用稳定节点集成(SNI)。在此背景下,提出了一种利用全变形梯度的光滑域准一致性更新方案。其次,通过将线性动量和角动量的内部虚功贡献扩展为线性泰勒级数,应用自然稳定节点积分(NSNI)。作为NSNI的替代方案,提出了子域集成(SDI),其中每个平滑域被划分为更小的子域,以提高集成精度和稳定性。最后,为了实现变相一致积分(VCI),对测试函数的梯度进行Petrov-Galerkin校正。采用弹性Neo-Hookean微极性材料对Cook膜进行了二维基准模拟,证明了该方法的性能。此外,采用基于drucker - prager的微极模型对考虑完美塑性的平面应变压缩试验进行了模拟。通过将结果与FEM和材料点法(MPM)的结果进行比较,评估了该公式的能力,突出了其在准脆性材料中模拟大变形的潜力。
{"title":"An implicit RKPM formulation for large deformations in elastoplastic micropolar media: Framework, stabilization and assessment","authors":"T. Mader ,&nbsp;A. Dummer ,&nbsp;P. Gamnitzer ,&nbsp;H. Gimperlein ,&nbsp;M. Neuner","doi":"10.1016/j.cma.2025.118582","DOIUrl":"10.1016/j.cma.2025.118582","url":null,"abstract":"<div><div>Lagrangian finite element methods (FEM) can suffer from excessive mesh distortion or entanglements in simulations involving large deformations, which become critical in the modeling of quasi-brittle materials characterized by failure mechanisms involving highly localized deformations such as shear bands. To overcome these limitations, we present an implicit numerical framework based on the Reproducing Kernel Particle Method (RKPM) for the micropolar continuum. The latter provides an effective framework to mitigate pathological mesh dependencies in shear-dominated failure and allows to consider size effects due to its intrinsic length scale. The RK shape functions are formulated in a semi-Lagrangian description based on B-Spline kernel functions. The weak form of the governing equations of the micropolar continuum is then derived and discretized within the RKPM framework. For implicit time integration, algorithmic tangent operators are derived to obtain quadratic convergence within the Newton-Raphson solution procedure. Several stabilization techniques are adopted for the proposed framework. First, stabilized nodal integration (SNI) is employed. In this context, a novel quasi-conforming update scheme for the smoothing domains is proposed, leveraging the full deformation gradient. Second, Naturally Stabilized Nodal Integration (NSNI) is applied by expanding the internal virtual work contributions for both linear and angular momentum into linear Taylor series. As an alternative to NSNI, a sub-domain integration (SDI) is presented, where each smoothing domain is partitioned into smaller sub-domains to improve integration accuracy and stability. Finally, to achieve variationally consistent integration (VCI), a Petrov-Galerkin correction is applied to the gradients of the test functions. The performance of the method is demonstrated by 2D benchmark simulations of Cook’s membrane using an elastic Neo-Hookean micropolar material. In addition, simulations of a plane strain compression test considering perfect plasticity using a Drucker-Prager-based micropolar model are performed. The capabilities of the formulation are assessed by comparing the results with those obtained with the FEM and the Material Point Method (MPM), highlighting its potential for modeling large deformations in quasi-brittle materials.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118582"},"PeriodicalIF":7.3,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731766","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}
引用次数: 0
Design of stiff elasto-plastic structures 刚性弹塑性结构设计
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-11 DOI: 10.1016/j.cma.2025.118622
Gunnar Granlund, Mathias Wallin
This study investigates “stiffness” optimization of path-dependent elasto-plastic structures by comparing the traditional secant stiffness formulation with a novel tangent stiffness formulation. Design updates are generated using the gradient-based Method of Moving Asymptotes (MMA), and the material behavior is modeled and limited to small-strain elasto-plasticity and isotropic hardening. Numerical examples compare and analyze the two stiffness definitions, revealing that they produce different designs and structural responses. Cross validation and comparison reveals that the secant stiffness optimized designs may have very low end tangent stiffness, and that the end tangent stiffness optimized designs may have very low secant stiffness, showing that the two formulations can contradict each other for elasto-plastic structures.
通过比较传统的割线刚度公式和新的切线刚度公式,研究了路径相关弹塑性结构的“刚度”优化问题。设计更新使用基于梯度的移动渐近线方法(MMA)生成,材料行为建模并仅限于小应变弹塑性和各向同性硬化。数值算例对两种刚度定义进行了比较和分析,揭示了它们产生不同的设计和结构响应。交叉验证和比较表明,割线刚度优化设计可能具有非常低的端切刚度,而端切刚度优化设计可能具有非常低的割线刚度,这表明两种公式对于弹塑性结构可能相互矛盾。
{"title":"Design of stiff elasto-plastic structures","authors":"Gunnar Granlund,&nbsp;Mathias Wallin","doi":"10.1016/j.cma.2025.118622","DOIUrl":"10.1016/j.cma.2025.118622","url":null,"abstract":"<div><div>This study investigates “stiffness” optimization of path-dependent elasto-plastic structures by comparing the traditional secant stiffness formulation with a novel tangent stiffness formulation. Design updates are generated using the gradient-based Method of Moving Asymptotes (MMA), and the material behavior is modeled and limited to small-strain elasto-plasticity and isotropic hardening. Numerical examples compare and analyze the two stiffness definitions, revealing that they produce different designs and structural responses. Cross validation and comparison reveals that the secant stiffness optimized designs may have very low end tangent stiffness, and that the end tangent stiffness optimized designs may have very low secant stiffness, showing that the two formulations can contradict each other for elasto-plastic structures.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118622"},"PeriodicalIF":7.3,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731148","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}
引用次数: 0
An efficient level set-based mesh adaptation for the particle finite element method 基于水平集的粒子有限元网格自适应方法
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-11 DOI: 10.1016/j.cma.2025.118644
Martin Lacroix, Eduardo Fernández, Simon Février, Luc Papeleux, Romain Boman, Jean-Philippe Ponthot
The Particle Finite Element Method (PFEM) is a discretization technique that combines the flexibility of particle-based methods with the precision of finite elements, using a Lagrangian approach to naturally track evolving interfaces and automatic remeshing to prevent mesh distortion. Historically, the PFEM has relied on the extraction of an α-shape from a Delaunay triangulation of the cloud of nodes forming fluid domain during the remeshing process. This approach helps to maintain good quality elements throughout the simulation, but introduces shortcomings that demand geometrical treatments tailored to each problem. In order to improve the remeshing process in PFEM, Falla et al. (2023) proposed a 2D mesh adaptation technique based on the edge splitting, showing promising results in terms of mass conservation and mesh quality. In parallel, Fernández et al. (2023) proposed the use of a level set (LS) function instead of the α-shape criterion, demonstrating improved mass conservation and free surface smoothness compared to standard approaches. While both innovative and foundational, these approaches have been limited to 2D applications, and the computation of the LS adds significant execution time to the PFEM. In this work, we propose a new remeshing algorithm, building upon the advances achieved by Falla and Fernández, designed to deliver good performance while extending the capability to handle 3D scenarios effectively. The interest of the LS lies in its ability to consider the overall fluid volume rather than focusing on the shape of individual elements as in the classical α-shape. Consequently, the LS allows for a better control over the connecting elements that are created during the fluid/fluid or fluid/solid contact, which helps to reduce spurious mass creation when merging free surfaces. The methodology is presented and validated using free surface flow problems in 2D and 3D. Finally, an overview of computation times is provided.
粒子有限元法(PFEM)是一种离散化技术,它结合了基于粒子的方法的灵活性和有限元的精度,使用拉格朗日方法自然地跟踪不断变化的界面和自动重划分以防止网格变形。历史上,PFEM依赖于重网格过程中形成流体域的节点云的Delaunay三角剖分提取α-形状。这种方法有助于在整个模拟过程中保持高质量的元素,但也引入了需要针对每个问题进行几何处理的缺点。为了改进PFEM中的网格重划分过程,Falla et al.(2023)提出了一种基于边缘分裂的二维网格自适应技术,在质量守恒和网格质量方面都取得了良好的效果。与此同时,Fernández等人(2023)提出使用水平集(LS)函数代替α-形状准则,与标准方法相比,证明了质量守恒和自由表面光滑性的改进。虽然这些方法都是创新和基础的,但它们仅限于2D应用,并且LS的计算增加了PFEM的执行时间。在这项工作中,我们提出了一种新的重网格算法,建立在Falla和Fernández取得的进步的基础上,旨在提供良好的性能,同时扩展有效处理3D场景的能力。LS的有趣之处在于它能够考虑整体流体体积,而不是像经典的α-形状那样专注于单个元素的形状。因此,LS可以更好地控制流体/流体或流体/固体接触过程中产生的连接元件,这有助于减少合并自由表面时产生的虚假质量。在二维和三维自由表面流动问题中提出并验证了该方法。最后,概述了计算时间。
{"title":"An efficient level set-based mesh adaptation for the particle finite element method","authors":"Martin Lacroix,&nbsp;Eduardo Fernández,&nbsp;Simon Février,&nbsp;Luc Papeleux,&nbsp;Romain Boman,&nbsp;Jean-Philippe Ponthot","doi":"10.1016/j.cma.2025.118644","DOIUrl":"10.1016/j.cma.2025.118644","url":null,"abstract":"<div><div>The Particle Finite Element Method (PFEM) is a discretization technique that combines the flexibility of particle-based methods with the precision of finite elements, using a Lagrangian approach to naturally track evolving interfaces and automatic remeshing to prevent mesh distortion. Historically, the PFEM has relied on the extraction of an <em>α</em>-shape from a Delaunay triangulation of the cloud of nodes forming fluid domain during the remeshing process. This approach helps to maintain good quality elements throughout the simulation, but introduces shortcomings that demand geometrical treatments tailored to each problem. In order to improve the remeshing process in PFEM, Falla et al. (2023) proposed a 2D mesh adaptation technique based on the edge splitting, showing promising results in terms of mass conservation and mesh quality. In parallel, Fernández et al. (2023) proposed the use of a level set (LS) function instead of the <em>α</em>-shape criterion, demonstrating improved mass conservation and free surface smoothness compared to standard approaches. While both innovative and foundational, these approaches have been limited to 2D applications, and the computation of the LS adds significant execution time to the PFEM. In this work, we propose a new remeshing algorithm, building upon the advances achieved by Falla and Fernández, designed to deliver good performance while extending the capability to handle 3D scenarios effectively. The interest of the LS lies in its ability to consider the overall fluid volume rather than focusing on the shape of individual elements as in the classical <em>α</em>-shape. Consequently, the LS allows for a better control over the connecting elements that are created during the fluid/fluid or fluid/solid contact, which helps to reduce spurious mass creation when merging free surfaces. The methodology is presented and validated using free surface flow problems in 2D and 3D. Finally, an overview of computation times is provided.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118644"},"PeriodicalIF":7.3,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731765","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}
引用次数: 0
Efficient bound-preserving numerical schemes for a phase-Field model of tumor growth with extracellular matrix degradation 具有细胞外基质降解的肿瘤生长相场模型的有效保界数值格式
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-10 DOI: 10.1016/j.cma.2025.118621
Qiumei Huang , Zhonghua Qiao , Cheng Wang , Huiting Yang
In cancer research, the role of the extracellular matrix (ECM) and its associated matrix-degrading enzyme (MDE) has been a significant area of focus. This study presents a numerical algorithm designed to simulate a previously established tumor model that incorporates various biological factors, including tumor cells, viable cells, necrotic cells, and the dynamics of MDE and ECM. The model consists of a system that includes a phase-field equation, two reaction-diffusion equations, and two ordinary differential equations (ODEs). We employ the fast exponential time differencing Runge-Kutta (ETDRK) method with stabilizing terms to solve this system, resulting in a decoupled, explicit, linear numerical algorithm. The objective of this algorithm is to preserve the physical properties of the model variables, including the maximum bound principle (MBP) for nutrient concentration and MDE volume fraction, as well as bound preserving for ECM density and tumor volume fraction. We perform simulations of 2D and 3D tumor models and discuss how different biological components impact growth dynamics. These simulations may help predict tumor evolution trends, offer insights for related biological and medical research, potentially reduce the number and cost of experiments, and improve research efficiency.
在癌症研究中,细胞外基质(ECM)及其相关基质降解酶(MDE)的作用一直是一个重要的关注领域。本研究提出了一种数值算法,旨在模拟先前建立的肿瘤模型,该模型包含各种生物因素,包括肿瘤细胞、活细胞、坏死细胞以及MDE和ECM的动态。该模型由一个相场方程、两个反应扩散方程和两个常微分方程组成。采用带稳定项的快速指数差分龙格-库塔(ETDRK)方法求解该系统,得到解耦、显式、线性的数值算法。该算法的目标是保留模型变量的物理性质,包括营养物质浓度和MDE体积分数的最大结合原理(MBP),以及ECM密度和肿瘤体积分数的结合保存。我们进行了二维和三维肿瘤模型的模拟,并讨论了不同的生物成分如何影响生长动力学。这些模拟可能有助于预测肿瘤的进化趋势,为相关的生物学和医学研究提供见解,有可能减少实验的数量和成本,并提高研究效率。
{"title":"Efficient bound-preserving numerical schemes for a phase-Field model of tumor growth with extracellular matrix degradation","authors":"Qiumei Huang ,&nbsp;Zhonghua Qiao ,&nbsp;Cheng Wang ,&nbsp;Huiting Yang","doi":"10.1016/j.cma.2025.118621","DOIUrl":"10.1016/j.cma.2025.118621","url":null,"abstract":"<div><div>In cancer research, the role of the extracellular matrix (ECM) and its associated matrix-degrading enzyme (MDE) has been a significant area of focus. This study presents a numerical algorithm designed to simulate a previously established tumor model that incorporates various biological factors, including tumor cells, viable cells, necrotic cells, and the dynamics of MDE and ECM. The model consists of a system that includes a phase-field equation, two reaction-diffusion equations, and two ordinary differential equations (ODEs). We employ the fast exponential time differencing Runge-Kutta (ETDRK) method with stabilizing terms to solve this system, resulting in a decoupled, explicit, linear numerical algorithm. The objective of this algorithm is to preserve the physical properties of the model variables, including the maximum bound principle (MBP) for nutrient concentration and MDE volume fraction, as well as bound preserving for ECM density and tumor volume fraction. We perform simulations of 2D and 3D tumor models and discuss how different biological components impact growth dynamics. These simulations may help predict tumor evolution trends, offer insights for related biological and medical research, potentially reduce the number and cost of experiments, and improve research efficiency.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118621"},"PeriodicalIF":7.3,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145732428","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}
引用次数: 0
On the role of locking mitigation in phase-field modelling of ductile fracture 塑性断裂相场模拟中锁紧减缓的作用
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2025-12-10 DOI: 10.1016/j.cma.2025.118607
A. Quintanas-Corominas , P. Olivares-Rodríguez , J. Reinoso , I.G. García
The development of nonlocal fracture models based on the phase-field fracture approach has been a subject of intensive research in recent years. This work focuses on the analysis of shear band localisation in ductile fracture through the application of well-established numerical techniques to mitigate volumetric locking. In particular, the corresponding performance of the B-Bar and the enhanced assumed strain formulations is examined, alongside a detailed energy-based analysis of ductile fracture models employing damage activation thresholds within the phase-field fracture framework. The standard von Mises (J2) plasticity model is used to describe the mechanical response at the material point level, which is coupled with the phase-field. Several representative benchmark problems are addressed, with special attention given to the discrepancies observed in shear band formation and the resulting failure patterns when locking alleviation techniques are applied. The present results demonstrate the efficiency and applicability of these formulations for predicting shear band localisation.
基于相场裂缝方法的非局部裂缝模型的发展是近年来研究的热点。这项工作的重点是通过应用成熟的数值技术来减轻体积锁定,分析韧性断裂中的剪切带局部化。特别地,研究了B-Bar的相应性能和增强的假设应变公式,以及在相场断裂框架内采用损伤激活阈值的韧性断裂模型的详细能量分析。采用标准的von Mises (J2)塑性模型来描述与相场耦合的材料点水平的力学响应。讨论了几个代表性的基准问题,特别注意在应用锁定缓解技术时观察到的剪切带形成和由此产生的破坏模式的差异。目前的结果证明了这些公式在预测剪切带局部化方面的有效性和适用性。
{"title":"On the role of locking mitigation in phase-field modelling of ductile fracture","authors":"A. Quintanas-Corominas ,&nbsp;P. Olivares-Rodríguez ,&nbsp;J. Reinoso ,&nbsp;I.G. García","doi":"10.1016/j.cma.2025.118607","DOIUrl":"10.1016/j.cma.2025.118607","url":null,"abstract":"<div><div>The development of nonlocal fracture models based on the phase-field fracture approach has been a subject of intensive research in recent years. This work focuses on the analysis of shear band localisation in ductile fracture through the application of well-established numerical techniques to mitigate volumetric locking. In particular, the corresponding performance of the B-Bar and the enhanced assumed strain formulations is examined, alongside a detailed energy-based analysis of ductile fracture models employing damage activation thresholds within the phase-field fracture framework. The standard von Mises (J2) plasticity model is used to describe the mechanical response at the material point level, which is coupled with the phase-field. Several representative benchmark problems are addressed, with special attention given to the discrepancies observed in shear band formation and the resulting failure patterns when locking alleviation techniques are applied. The present results demonstrate the efficiency and applicability of these formulations for predicting shear band localisation.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"450 ","pages":"Article 118607"},"PeriodicalIF":7.3,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731151","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}
引用次数: 0
期刊
Computer Methods in Applied Mechanics and Engineering
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1