首页 > 最新文献

Computer Methods in Applied Mechanics and Engineering最新文献

英文 中文
Wavelet-based enrichment for physics informed neural networks to approximate localized and heterogeneous solutions in solid mechanics 基于小波的物理富集告诉神经网络在固体力学中近似局域和非均质解
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-24 DOI: 10.1016/j.cma.2026.118768
Duc-Vinh Nguyen , Mohamed Jebahi , Francisco Chinesta
Recent research has highlighted the potential of physics-informed neural networks (PINNs) as an efficient methodology for approximating solutions of boundary value problems in solid mechanics. Nevertheless, their ability to accurately capture highly heterogeneous solutions of complex problems remains limited and requires further investigation. The present paper explores new strategies to address this challenge. In line with existing approaches based on local refinement of collocation (training) points, a weighted version of the loss function is first proposed to better balance the physical residuals across the entire computational domain. Although this modification improves overall performance, the approximation accuracy remains unsatisfactory. To overcome this limitation, an enriched version of PINN is developed to more effectively capture locally heterogeneous distributions of state variables. Specifically, wavelet-based enrichment functions are designed to approximate local high-frequency components of the full-field solution, thereby simplifying the task of the neural network, which is then required only to approximate the global smooth component of the solution. This approach achieves satisfactory accuracy even with relatively simple neural network architectures and few collocation points, as demonstrated through several benchmark problems. Therefore, the proposed enrichment concept represents a promising direction for further improving the performance of PINNs as solvers in computational mechanics, paving the way for their application to more complex problems.
最近的研究强调了物理信息神经网络(pinn)作为固体力学中边值问题近似解的有效方法的潜力。然而,它们准确捕获复杂问题的高度异构解决方案的能力仍然有限,需要进一步研究。本文探讨了应对这一挑战的新策略。在现有的基于搭配(训练)点局部细化的方法基础上,首先提出了损失函数的加权版本,以更好地平衡整个计算域的物理残差。虽然这种修改提高了整体性能,但近似精度仍然令人不满意。为了克服这一限制,开发了一个丰富的PINN版本,以更有效地捕获状态变量的局部异构分布。具体来说,基于小波的富集函数被设计为近似全场解的局部高频分量,从而简化了神经网络的任务,然后只需要近似解的全局光滑分量。这种方法即使在相对简单的神经网络架构和很少的搭配点下也能达到令人满意的精度,通过几个基准测试问题证明了这一点。因此,所提出的富集概念为进一步提高pin神经网络作为计算力学求解器的性能提供了一个有希望的方向,为其应用于更复杂的问题铺平了道路。
{"title":"Wavelet-based enrichment for physics informed neural networks to approximate localized and heterogeneous solutions in solid mechanics","authors":"Duc-Vinh Nguyen ,&nbsp;Mohamed Jebahi ,&nbsp;Francisco Chinesta","doi":"10.1016/j.cma.2026.118768","DOIUrl":"10.1016/j.cma.2026.118768","url":null,"abstract":"<div><div>Recent research has highlighted the potential of physics-informed neural networks (PINNs) as an efficient methodology for approximating solutions of boundary value problems in solid mechanics. Nevertheless, their ability to accurately capture highly heterogeneous solutions of complex problems remains limited and requires further investigation. The present paper explores new strategies to address this challenge. In line with existing approaches based on local refinement of collocation (training) points, a weighted version of the loss function is first proposed to better balance the physical residuals across the entire computational domain. Although this modification improves overall performance, the approximation accuracy remains unsatisfactory. To overcome this limitation, an enriched version of PINN is developed to more effectively capture locally heterogeneous distributions of state variables. Specifically, wavelet-based enrichment functions are designed to approximate local high-frequency components of the full-field solution, thereby simplifying the task of the neural network, which is then required only to approximate the global smooth component of the solution. This approach achieves satisfactory accuracy even with relatively simple neural network architectures and few collocation points, as demonstrated through several benchmark problems. Therefore, the proposed enrichment concept represents a promising direction for further improving the performance of PINNs as solvers in computational mechanics, paving the way for their application to more complex problems.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118768"},"PeriodicalIF":7.3,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023763","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
Physics informed surface autoencoders for thin shell analysis 物理通知表面自编码器薄壳分析
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-24 DOI: 10.1016/j.cma.2026.118764
Aswanth Thani , Adrian Buganza Tepole
We present a physics-informed surface autoencoder (PISA) framework for Kirchhoff-Love thin shell analysis. The method constructs global C1 surface parameterizations directly from unstructured point clouds for both single-patch surfaces homeomorphic to disks, and multi-patch parameterizations for closed genus-zero surfaces. In the multi-patch case, a classification network assigns probabilistic labels to points, and the autoencoder learns overlapping charts with smooth transitions, ensuring global C1 continuity. With the learned parameterizations, we introduce a decoder for the displacement field and compute differential geometric quantities such as the metric and second fundamental form in the reference and deformed surfaces. Then, we enforce equilibrium by minimizing the total potential energy. The approach is validated on classical shell benchmarks, including the Scordelis-Lo roof, pinched cylinder, and hemisphere under pressure. We showcase the flexibility of the framework with complex geometries such as the Stanford Bunny and dura mater. Compared with traditional spline-based parameterizations and existing machine learning approaches, PISA offers a pipeline for generating smooth surface maps for complex geometries and integrates the surface representation into the physics-informed solver. Importantly, the thin shell analysis pipeline proposed works directly with unstructured point cloud data. Thus, this PISA framework’s potential applications range from engineering structures to biological membranes such as heart valves, skin, and dura mater.
我们提出了一个用于Kirchhoff-Love薄壳分析的物理通知表面自动编码器(PISA)框架。该方法直接从非结构化点云构建全局C1曲面参数化,用于单片曲面与磁盘同胚,以及闭合零属曲面的多片参数化。在多补丁情况下,分类网络为点分配概率标签,自动编码器学习平滑过渡的重叠图,确保全局C1连续性。利用学习到的参数化,我们引入了位移场的解码器,并计算了参考曲面和变形曲面上的度量和第二基本形式等微分几何量。然后,我们通过最小化总势能来实现平衡。该方法在经典的壳体基准测试中得到了验证,包括Scordelis-Lo顶板、挤压圆柱体和压力下的半球。我们展示了框架的灵活性与复杂的几何形状,如斯坦福兔和硬脑膜。与传统的基于样条的参数化和现有的机器学习方法相比,PISA提供了一个为复杂几何图形生成光滑表面映射的管道,并将表面表示集成到物理信息求解器中。重要的是,所提出的薄壳分析管道直接适用于非结构化点云数据。因此,这个PISA框架的潜在应用范围从工程结构到生物膜,如心脏瓣膜、皮肤和硬脑膜。
{"title":"Physics informed surface autoencoders for thin shell analysis","authors":"Aswanth Thani ,&nbsp;Adrian Buganza Tepole","doi":"10.1016/j.cma.2026.118764","DOIUrl":"10.1016/j.cma.2026.118764","url":null,"abstract":"<div><div>We present a physics-informed surface autoencoder (PISA) framework for Kirchhoff-Love thin shell analysis. The method constructs global <em>C</em><sup>1</sup> surface parameterizations directly from unstructured point clouds for both single-patch surfaces homeomorphic to disks, and multi-patch parameterizations for closed genus-zero surfaces. In the multi-patch case, a classification network assigns probabilistic labels to points, and the autoencoder learns overlapping charts with smooth transitions, ensuring global <em>C</em><sup>1</sup> continuity. With the learned parameterizations, we introduce a decoder for the displacement field and compute differential geometric quantities such as the metric and second fundamental form in the reference and deformed surfaces. Then, we enforce equilibrium by minimizing the total potential energy. The approach is validated on classical shell benchmarks, including the Scordelis-Lo roof, pinched cylinder, and hemisphere under pressure. We showcase the flexibility of the framework with complex geometries such as the Stanford Bunny and dura mater. Compared with traditional spline-based parameterizations and existing machine learning approaches, PISA offers a pipeline for generating smooth surface maps for complex geometries and integrates the surface representation into the physics-informed solver. Importantly, the thin shell analysis pipeline proposed works directly with unstructured point cloud data. Thus, this PISA framework’s potential applications range from engineering structures to biological membranes such as heart valves, skin, and dura mater.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118764"},"PeriodicalIF":7.3,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023859","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
Optimal control of a hemivariational inequality of stationary convective Brinkman-Forchheimer extended Darcy equations with numerical approximation 平稳对流Brinkman-Forchheimer扩展Darcy方程半变分不等式的最优控制
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-24 DOI: 10.1016/j.cma.2026.118755
Wasim Akram, Manil T. Mohan
We investigate an optimal control problem governed by a stationary convective Brinkman-Forchheimer extended Darcy (CBFeD) model formulated as a hemivariational inequality in both two- and three-dimensional settings. This framework captures complex incompressible fluid flow through porous media by simultaneously accounting for convection, viscous damping, and nonlinear resistance effects, while naturally incorporating non-smooth frictional interactions through a subdifferential boundary condition. A key contribution of this work is a rigorous stability analysis of the CBFeD hemivariational inequality with respect to perturbations in both the external force density and the associated superpotential. Building on this analysis, we establish the existence of optimal controls when the external force density is treated as the control variable under admissible constraints. This result extends existing optimal control theories to a broader class of nonsmooth, nonlinear flow models in porous media. From a computational perspective, we propose a fully implementable numerical scheme for the resulting optimal control problem and prove its convergence. The method is based on finite element discretization and is applicable in both two and three dimensions, making it suitable for practical simulations. Numerical experiments are presented to illustrate the effectiveness of the proposed approach and to confirm the theoretical findings.
我们研究了一个由静态对流Brinkman-Forchheimer扩展Darcy (CBFeD)模型控制的最优控制问题,该模型在二维和三维环境中被表述为半分不等式。该框架通过同时考虑对流、粘性阻尼和非线性阻力效应来捕获复杂的不可压缩流体在多孔介质中的流动,同时通过次微分边界条件自然地纳入非光滑摩擦相互作用。这项工作的一个关键贡献是对CBFeD半变分不等式在外力密度和相关超势的扰动下的严格稳定性分析。在此基础上,建立了在允许约束条件下以外力密度为控制变量的最优控制的存在性。这一结果将现有的最优控制理论扩展到更广泛的非光滑、非线性多孔介质流动模型。从计算的角度,我们提出了一个完全可实现的最优控制问题的数值格式,并证明了其收敛性。该方法基于有限元离散化,适用于二维和三维,适合于实际仿真。数值实验证明了所提方法的有效性,并证实了理论结论。
{"title":"Optimal control of a hemivariational inequality of stationary convective Brinkman-Forchheimer extended Darcy equations with numerical approximation","authors":"Wasim Akram,&nbsp;Manil T. Mohan","doi":"10.1016/j.cma.2026.118755","DOIUrl":"10.1016/j.cma.2026.118755","url":null,"abstract":"<div><div>We investigate an optimal control problem governed by a stationary convective Brinkman-Forchheimer extended Darcy (CBFeD) model formulated as a hemivariational inequality in both two- and three-dimensional settings. This framework captures complex incompressible fluid flow through porous media by simultaneously accounting for convection, viscous damping, and nonlinear resistance effects, while naturally incorporating non-smooth frictional interactions through a subdifferential boundary condition. A key contribution of this work is a rigorous stability analysis of the CBFeD hemivariational inequality with respect to perturbations in both the external force density and the associated superpotential. Building on this analysis, we establish the existence of optimal controls when the external force density is treated as the control variable under admissible constraints. This result extends existing optimal control theories to a broader class of nonsmooth, nonlinear flow models in porous media. From a computational perspective, we propose a fully implementable numerical scheme for the resulting optimal control problem and prove its convergence. The method is based on finite element discretization and is applicable in both two and three dimensions, making it suitable for practical simulations. Numerical experiments are presented to illustrate the effectiveness of the proposed approach and to confirm the theoretical findings.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118755"},"PeriodicalIF":7.3,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146048313","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
A new family of explicit generalized single-step single-stage integration methods for structural/multibody dynamics with improved stability for viscous damping 一种新的结构/多体动力学单步单级显式广义积分方法,提高了粘性阻尼的稳定性
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-24 DOI: 10.1016/j.cma.2026.118732
Dean J. Maxam, Pranav Chengala Madhusoodana, Kumar K. Tamma
This article details the analysis and synthesis of a novel family of explicit single-step integration methods for the second-order dynamics of structural and multibody systems. The family is derived from the Generalized Single-Step Single-Solve framework, pertaining to the class of linear multistep methods. The new explicit methods achieve second-order accuracy with optimal starting error, controllable numerical dissipation, and an option for explicit or implicit treatment of damping; the latter yields a stability limit which scales optimally with modal damping ratio, unlike prior methods with only incidental gains. The new family is compared with existing explicit methods on the basis of numerical accuracy and stability. Its superior performance for linear and nonlinear systems is demonstrated by numerical examples.
本文详细分析和综合了一类新的结构和多体系统二阶动力学的显式单步积分方法。该族由广义单步单解框架衍生而来,属于线性多步方法类。新的显式方法实现了二阶精度,具有最优的启动误差,可控的数值耗散,并可选择显式或隐式处理阻尼;后者产生一个稳定性极限,它与模态阻尼比最佳缩放,不像以前的方法只有附带增益。在数值精度和稳定性方面,与现有的显式方法进行了比较。数值算例证明了该方法在线性和非线性系统中的优越性能。
{"title":"A new family of explicit generalized single-step single-stage integration methods for structural/multibody dynamics with improved stability for viscous damping","authors":"Dean J. Maxam,&nbsp;Pranav Chengala Madhusoodana,&nbsp;Kumar K. Tamma","doi":"10.1016/j.cma.2026.118732","DOIUrl":"10.1016/j.cma.2026.118732","url":null,"abstract":"<div><div>This article details the analysis and synthesis of a novel family of explicit single-step integration methods for the second-order dynamics of structural and multibody systems. The family is derived from the Generalized Single-Step Single-Solve framework, pertaining to the class of linear multistep methods. The new explicit methods achieve second-order accuracy with optimal starting error, controllable numerical dissipation, and an option for explicit or implicit treatment of damping; the latter yields a stability limit which scales optimally with modal damping ratio, unlike prior methods with only incidental gains. The new family is compared with existing explicit methods on the basis of numerical accuracy and stability. Its superior performance for linear and nonlinear systems is demonstrated by numerical examples.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118732"},"PeriodicalIF":7.3,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146048404","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
Eigenstate based homogenization 基于特征态的均匀化
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-24 DOI: 10.1016/j.cma.2025.118718
Jacob Fish, Junhe Cui
The goal of this work is to develop a homogenization approach that achieves the computational efficiency of unresolved (reduced-order) schemes while retaining the accuracy of fully resolved computational homogenization. In the proposed Eigenstate-Based Homogenization (EBH), the fine-scale solution space is represented in terms of eigenstates—such as plastic multipliers or damage parameters—resulting in far fewer unknowns than in finite element discretizations. The representative volume element (RVE) response is expressed as a linear combination of precomputed eigenmodes, together with partition-level state variables updated to satisfy local consistency, enabling accurate recovery of both fine- and coarse-scale fields.
The nonlinear evolution of the state variables is solved using a Controllable Fixed-Point (CFP) iteration. For plasticity, the spectral radius of the iteration matrix scales linearly with the plastic-multiplier increment; for regularized damage, convergence depends on the load increment, maximum damage, and viscous regularization, ranging from sublinear to superlinear. The resulting tangent operator of the nonlinear system of equation of eigenstates has a condition number several orders of magnitude smaller than that of finite-element direct homogenization, yielding substantial computational savings. To further reduce computational cost, we study a truncated variant of the fixed-point method, termed the Truncated Fixed-Point (TFP) scheme, which leverages the rapid initial error decrease of CFP to obtain low-tolerance solutions with far fewer iterations. Its CPU time is comparable to well-known reduced-order homogenization approaches employing one partition per phase with calibrated properties to fit coarse-scale behavior, yet its accuracy approaches that of fully converged CFP solutions.
Numerical studies show that EBH reproduces both fine-scale fields and macroscopic responses in close agreement with fully resolved finite-element simulations, but at a fraction of the computational cost.
本工作的目标是开发一种均质化方法,在保持完全解析计算均质化的准确性的同时,实现未解析(降阶)格式的计算效率。在提出的基于特征态的均质化(EBH)中,精细尺度解空间用特征态(如塑性乘数或损伤参数)表示,导致的未知数比有限元离散化少得多。代表性体积元(RVE)响应表示为预先计算的特征模态的线性组合,以及更新以满足局部一致性的分区级状态变量,从而能够精确恢复细尺度和粗尺度场。采用可控不动点迭代法求解状态变量的非线性演化问题。对于塑性,迭代矩阵的谱半径随塑性乘法器的增量呈线性增长;对于正则化损伤,收敛取决于载荷增量、最大损伤和粘性正则化,范围从亚线性到超线性。所得到的非线性特征态方程系统的正切算子的条件数比有限元直接均匀化的条件数小几个数量级,从而大大节省了计算量。为了进一步减少计算成本,我们研究了一种截断的定点方法,称为截断定点(TFP)方案,该方案利用CFP的快速初始误差减小来获得低公差解,迭代次数少得多。它的CPU时间与众所周知的降阶均匀化方法相当,该方法采用每相位一个分区,具有校准特性以适应粗尺度行为,但其精度接近完全收敛的CFP解决方案。数值研究表明,EBH重现的细尺度场和宏观响应与完全解析的有限元模拟非常接近,但计算成本只是其中的一小部分。
{"title":"Eigenstate based homogenization","authors":"Jacob Fish,&nbsp;Junhe Cui","doi":"10.1016/j.cma.2025.118718","DOIUrl":"10.1016/j.cma.2025.118718","url":null,"abstract":"<div><div>The goal of this work is to develop a homogenization approach that achieves the computational efficiency of unresolved (reduced-order) schemes while retaining the accuracy of fully resolved computational homogenization. In the proposed Eigenstate-Based Homogenization (EBH), the fine-scale solution space is represented in terms of eigenstates—such as plastic multipliers or damage parameters—resulting in far fewer unknowns than in finite element discretizations. The representative volume element (RVE) response is expressed as a linear combination of precomputed eigenmodes, together with partition-level state variables updated to satisfy local consistency, enabling accurate recovery of both fine- and coarse-scale fields.</div><div>The nonlinear evolution of the state variables is solved using a Controllable Fixed-Point (CFP) iteration. For plasticity, the spectral radius of the iteration matrix scales linearly with the plastic-multiplier increment; for regularized damage, convergence depends on the load increment, maximum damage, and viscous regularization, ranging from sublinear to superlinear. The resulting tangent operator of the nonlinear system of equation of eigenstates has a condition number several orders of magnitude smaller than that of finite-element direct homogenization, yielding substantial computational savings. To further reduce computational cost, we study a truncated variant of the fixed-point method, termed the Truncated Fixed-Point (TFP) scheme, which leverages the rapid initial error decrease of CFP to obtain low-tolerance solutions with far fewer iterations. Its CPU time is comparable to well-known reduced-order homogenization approaches employing one partition per phase with calibrated properties to fit coarse-scale behavior, yet its accuracy approaches that of fully converged CFP solutions.</div><div>Numerical studies show that EBH reproduces both fine-scale fields and macroscopic responses in close agreement with fully resolved finite-element simulations, but at a fraction of the computational cost.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118718"},"PeriodicalIF":7.3,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023762","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 prospects of interpolatory spline bases for accurate mass lumping strategies in isogeometric analysis 等几何分析中精确质量集总策略的插值样条基展望
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-23 DOI: 10.1016/j.cma.2026.118762
Yannis Voet , Espen Sande
While interpolatory bases such as the Lagrange basis form the cornerstone of classical finite element methods, they have been replaced in the more general finite element setting of isogeometric analysis in favor of other desirable properties. Yet, interpolation is a key property for devising accurate mass lumping strategies that are ubiquitous in explicit dynamic analyses of structures. In this article, we explore the possibility of restoring interpolation for spline bases within isogeometric analysis for the purpose of mass lumping. Although reminiscent of the spectral element method, this technique comes with its lot of surprises and challenges, which are critically assessed.
虽然诸如拉格朗日基之类的插值基构成了经典有限元方法的基石,但在更一般的等几何分析的有限元设置中,它们已被其他理想性质所取代。然而,插值是设计精确的质量集总策略的关键属性,在结构的显式动力分析中无处不在。在这篇文章中,我们探讨了在等几何分析中为质量集总目的恢复样条基插值的可能性。虽然让人想起光谱元素方法,但这项技术带来了许多惊喜和挑战,这些都是经过严格评估的。
{"title":"On the prospects of interpolatory spline bases for accurate mass lumping strategies in isogeometric analysis","authors":"Yannis Voet ,&nbsp;Espen Sande","doi":"10.1016/j.cma.2026.118762","DOIUrl":"10.1016/j.cma.2026.118762","url":null,"abstract":"<div><div>While interpolatory bases such as the Lagrange basis form the cornerstone of classical finite element methods, they have been replaced in the more general finite element setting of isogeometric analysis in favor of other desirable properties. Yet, interpolation is a key property for devising accurate mass lumping strategies that are ubiquitous in explicit dynamic analyses of structures. In this article, we explore the possibility of restoring interpolation for spline bases within isogeometric analysis for the purpose of mass lumping. Although reminiscent of the spectral element method, this technique comes with its lot of surprises and challenges, which are critically assessed.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118762"},"PeriodicalIF":7.3,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023846","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
Topology-aware stress analysis in shell structures 壳结构的拓扑感知应力分析
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-23 DOI: 10.1016/j.cma.2026.118770
Junpeng Wang , Yingjian Liu , Jun Wu , Rüdiger Westermann
We present a stable and accurate algorithm for tracing principal stress lines (PSLs) in shell structures, applicable to both first- and second-order triangular and quadrilateral elements. The algorithm operates directly in the isoparametric space of the elements, leveraging their inherent shape functions to account for curved geometry without resorting to artificial subdivision. This approach enables, for the first time, a consistent stress topology analysis for shell elements, including a rigorous treatment of stress degeneracies. Our PSL seeding strategy integrates stress topology with the curved shell surface, ensuring a uniform and consistent PSL distribution. We evaluate the algorithm’s performance through a series of numerical experiments, demonstrating its utility for advanced stress analysis. Furthermore, we demonstrate the generation of a globally consistent, space-filling PSL structure, which is directly applicable to downstream tasks such as lightweight structural design. To support practical use, we provide a publicly available MATLAB implementation. The implementation features a unified file interface that supports diverse mesh types and is compatible with standard finite element method (FEM) output, offering a versatile tool for stress investigation and design evaluation in computational mechanics. The code is available at https://github.com/PSLer/PSLshell.
本文提出了一种稳定、准确的壳结构主应力线追踪算法,适用于一阶和二阶三角形和四边形单元。该算法直接在元素的等参空间中操作,利用其固有的形状函数来解释弯曲的几何形状,而无需诉诸人工细分。该方法首次实现了对壳单元的一致应力拓扑分析,包括对应力退化的严格处理。我们的PSL播种策略将应力拓扑与弯曲的壳表面相结合,确保PSL分布均匀一致。我们通过一系列的数值实验来评估算法的性能,证明了它在高级应力分析中的实用性。此外,我们还展示了一种全球一致的、填充空间的PSL结构的生成,它直接适用于轻量化结构设计等下游任务。为了支持实际使用,我们提供了一个公开可用的MATLAB实现。该实现具有统一的文件接口,支持多种网格类型,并与标准有限元方法(FEM)输出兼容,为计算力学中的应力调查和设计评估提供了一个通用的工具。代码可在https://github.com/PSLer/PSLshell上获得。
{"title":"Topology-aware stress analysis in shell structures","authors":"Junpeng Wang ,&nbsp;Yingjian Liu ,&nbsp;Jun Wu ,&nbsp;Rüdiger Westermann","doi":"10.1016/j.cma.2026.118770","DOIUrl":"10.1016/j.cma.2026.118770","url":null,"abstract":"<div><div>We present a stable and accurate algorithm for tracing principal stress lines (PSLs) in shell structures, applicable to both first- and second-order triangular and quadrilateral elements. The algorithm operates directly in the isoparametric space of the elements, leveraging their inherent shape functions to account for curved geometry without resorting to artificial subdivision. This approach enables, for the first time, a consistent stress topology analysis for shell elements, including a rigorous treatment of stress degeneracies. Our PSL seeding strategy integrates stress topology with the curved shell surface, ensuring a uniform and consistent PSL distribution. We evaluate the algorithm’s performance through a series of numerical experiments, demonstrating its utility for advanced stress analysis. Furthermore, we demonstrate the generation of a globally consistent, space-filling PSL structure, which is directly applicable to downstream tasks such as lightweight structural design. To support practical use, we provide a publicly available MATLAB implementation. The implementation features a unified file interface that supports diverse mesh types and is compatible with standard finite element method (FEM) output, offering a versatile tool for stress investigation and design evaluation in computational mechanics. The code is available at <span><span>https://github.com/PSLer/PSLshell</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118770"},"PeriodicalIF":7.3,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023847","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
Isogeometric block BDDC/FETI-DP preconditioners for the three-field Biot’s consolidation model 三场Biot固结模型的等几何块BDDC/FETI-DP预调节器
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-23 DOI: 10.1016/j.cma.2026.118757
Hanyu Chu , Luca Franco Pavarino , Stefano Zampini
This paper presents a block dual-primal preconditioner for the three-field mixed isogeometric discretization of the stationary Biot’s consolidation model, formulated in terms of displacement, pressure, and total pressure. After decomposing the computational domain into subdomains and eliminating the displacement variables and the interior components of pressure and total pressure within each subdomain, we reduce the problem to a symmetric positive definite system for the subdomain interface unknowns and the Lagrange multiplier. The reduced system is solved by a preconditioned conjugate gradient method with a block preconditioner, based on a Balancing Domain Decomposition by Constraints (BDDC) method with deluxe scaling for the interface block and a FETI-DP preconditioner for the Lagrange multiplier block. We prove that the algorithm is scalable with respect to the number of subdomains and achieves a quasi-optimal convergence rate bound that is polylogarithmic in the ratio of subdomain to element sizes and robust with respect to the model parameters. Numerical experiments confirm the efficiency of the proposed preconditioner, even in the presence of discontinuous Lamé parameters, and illustrate its robustness with respect to the spline polynomial degree, regularity, and domain deformation.
本文提出了固定Biot固结模型的三场混合等几何离散化的块双原预条件,用位移、压力和总压力表示。将计算域分解为子域,消去每个子域内的位移变量、压力和总压力的内部分量,将问题简化为子域界面未知数和拉格朗日乘子的对称正定系统。基于基于约束的平衡域分解(BDDC)方法和基于拉格朗日乘法器块的FETI-DP预调节器,采用带块预调节器的预条件共轭梯度法求解了简化后的系统。我们证明了该算法在子域数量上是可扩展的,并实现了一个拟最优收敛速度界,该界在子域与元素大小的比值上是多对数的,并且对模型参数具有鲁棒性。数值实验证实了所提出的预调节器的有效性,即使在存在不连续lam参数的情况下,也证明了它在样条多项式度、规则性和域变形方面的鲁棒性。
{"title":"Isogeometric block BDDC/FETI-DP preconditioners for the three-field Biot’s consolidation model","authors":"Hanyu Chu ,&nbsp;Luca Franco Pavarino ,&nbsp;Stefano Zampini","doi":"10.1016/j.cma.2026.118757","DOIUrl":"10.1016/j.cma.2026.118757","url":null,"abstract":"<div><div>This paper presents a block dual-primal preconditioner for the three-field mixed isogeometric discretization of the stationary Biot’s consolidation model, formulated in terms of displacement, pressure, and total pressure. After decomposing the computational domain into subdomains and eliminating the displacement variables and the interior components of pressure and total pressure within each subdomain, we reduce the problem to a symmetric positive definite system for the subdomain interface unknowns and the Lagrange multiplier. The reduced system is solved by a preconditioned conjugate gradient method with a block preconditioner, based on a Balancing Domain Decomposition by Constraints (BDDC) method with deluxe scaling for the interface block and a FETI-DP preconditioner for the Lagrange multiplier block. We prove that the algorithm is scalable with respect to the number of subdomains and achieves a quasi-optimal convergence rate bound that is polylogarithmic in the ratio of subdomain to element sizes and robust with respect to the model parameters. Numerical experiments confirm the efficiency of the proposed preconditioner, even in the presence of discontinuous Lamé parameters, and illustrate its robustness with respect to the spline polynomial degree, regularity, and domain deformation.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118757"},"PeriodicalIF":7.3,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023860","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
The inference of Fokker-Planck equations via transport maps 通过输运图的福克-普朗克方程的推断
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-22 DOI: 10.1016/j.cma.2026.118760
Saem Han , Krishna Garikipati
We present a framework, which, from the trajectories detailing the spatiotemporal dynamics of a population, simultaneously reconstructs a transport map as well as the Fokker-Planck equation governing the coarse-grained probability distribution. Leveraging the Knothe-Rosenblatt rearrangement, we model the transport map from a fixed reference distribution to the target distribution, and derive the velocity fields of the flows from the trajectory of transport maps. Exploiting the velocity fields, we circumvent spatial gradients to infer the Fokker-Planck equation’s potential and diffusivity. The sparsity of trajectories injects uncertainty, which we treat in a Bayesian setting using variational inference. The approach is applied to inferring the Fokker-Planck dynamics in spaces of up to five dimensions, demonstrating both accurate identification of the system and efficiency with respect to data size.
我们提出了一个框架,该框架从详细描述人口时空动态的轨迹中,同时重建了运输图以及控制粗粒度概率分布的福克-普朗克方程。利用Knothe-Rosenblatt重排法,我们建立了从固定参考分布到目标分布的输运图模型,并从输运图的轨迹推导出了流的速度场。利用速度场,我们绕过空间梯度来推断Fokker-Planck方程的势和扩散率。轨迹的稀疏性注入了不确定性,我们使用变分推理在贝叶斯设置中处理。该方法被应用于在多达五个维度的空间中推断福克-普朗克动力学,证明了系统的准确识别和数据大小方面的效率。
{"title":"The inference of Fokker-Planck equations via transport maps","authors":"Saem Han ,&nbsp;Krishna Garikipati","doi":"10.1016/j.cma.2026.118760","DOIUrl":"10.1016/j.cma.2026.118760","url":null,"abstract":"<div><div>We present a framework, which, from the trajectories detailing the spatiotemporal dynamics of a population, simultaneously reconstructs a transport map as well as the Fokker-Planck equation governing the coarse-grained probability distribution. Leveraging the Knothe-Rosenblatt rearrangement, we model the transport map from a fixed reference distribution to the target distribution, and derive the velocity fields of the flows from the trajectory of transport maps. Exploiting the velocity fields, we circumvent spatial gradients to infer the Fokker-Planck equation’s potential and diffusivity. The sparsity of trajectories injects uncertainty, which we treat in a Bayesian setting using variational inference. The approach is applied to inferring the Fokker-Planck dynamics in spaces of up to five dimensions, demonstrating both accurate identification of the system and efficiency with respect to data size.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118760"},"PeriodicalIF":7.3,"publicationDate":"2026-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146006490","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
Reinforcement learning closures for underresolved partial differential equations using synthetic data 使用合成数据的欠解偏微分方程的强化学习闭包
IF 7.3 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-22 DOI: 10.1016/j.cma.2026.118767
Lothar Heimbach , Sebastian Kaltenbach , Petr Karnakov , Francis J. Alexander , Petros Koumoutsakos
Partial Differential Equations (PDEs) describe phenomena ranging from turbulence and epidemics to quantum mechanics and financial markets. Despite recent advances in computational science, solving such PDEs for real-world applications remains prohibitively expensive because of the necessity of resolving a broad range of spatiotemporal scales. In turn, practitioners often rely on coarse-grained approximations of the original PDEs, trading off accuracy for reduced computational resources. To mitigate the loss of detail inherent in such approximations, closure models are employed to represent unresolved spatiotemporal interactions. We present a framework for developing closure models for PDEs using synthetic data acquired through the method of manufactured solutions. These data are used in conjunction with reinforcement learning to provide closures for coarse-grained PDEs. We illustrate the efficacy of our method using the one-dimensional and two-dimensional Burgers’ equations and the two-dimensional advection equation. Moreover, we demonstrate that closure models trained for inhomogeneous PDEs can be effectively generalized to homogeneous PDEs. The results demonstrate the potential for developing accurate and computationally efficient closure models for systems with scarce data.
偏微分方程(PDEs)描述从湍流和流行病到量子力学和金融市场的各种现象。尽管计算科学在最近取得了进展,但由于需要解决大范围的时空尺度,因此在实际应用中解决这种偏微分方程仍然非常昂贵。反过来,从业者通常依赖于原始偏微分方程的粗粒度近似值,以降低计算资源的准确性为代价。为了减轻这种近似中固有的细节损失,闭合模型被用来表示未解决的时空相互作用。我们提出了一个框架,利用通过制造解决方案的方法获得的合成数据来开发pde的闭合模型。这些数据与强化学习一起使用,为粗粒度的pde提供闭包。我们用一维和二维Burgers方程以及二维平流方程来说明我们的方法的有效性。此外,我们证明了针对非齐次偏微分方程训练的闭包模型可以有效地推广到齐次偏微分方程。研究结果表明,为具有稀缺数据的系统开发精确且计算效率高的闭包模型具有潜力。
{"title":"Reinforcement learning closures for underresolved partial differential equations using synthetic data","authors":"Lothar Heimbach ,&nbsp;Sebastian Kaltenbach ,&nbsp;Petr Karnakov ,&nbsp;Francis J. Alexander ,&nbsp;Petros Koumoutsakos","doi":"10.1016/j.cma.2026.118767","DOIUrl":"10.1016/j.cma.2026.118767","url":null,"abstract":"<div><div>Partial Differential Equations (PDEs) describe phenomena ranging from turbulence and epidemics to quantum mechanics and financial markets. Despite recent advances in computational science, solving such PDEs for real-world applications remains prohibitively expensive because of the necessity of resolving a broad range of spatiotemporal scales. In turn, practitioners often rely on coarse-grained approximations of the original PDEs, trading off accuracy for reduced computational resources. To mitigate the loss of detail inherent in such approximations, closure models are employed to represent unresolved spatiotemporal interactions. We present a framework for developing closure models for PDEs using synthetic data acquired through the method of manufactured solutions. These data are used in conjunction with reinforcement learning to provide closures for coarse-grained PDEs. We illustrate the efficacy of our method using the one-dimensional and two-dimensional Burgers’ equations and the two-dimensional advection equation. Moreover, we demonstrate that closure models trained for inhomogeneous PDEs can be effectively generalized to homogeneous PDEs. The results demonstrate the potential for developing accurate and computationally efficient closure models for systems with scarce data.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"452 ","pages":"Article 118767"},"PeriodicalIF":7.3,"publicationDate":"2026-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023844","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