首页 > 最新文献

Computer Methods in Applied Mechanics and Engineering最新文献

英文 中文
System stabilization with policy optimization on unstable latent manifolds 在不稳定潜流形上通过策略优化实现系统稳定
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117483
Steffen W.R. Werner , Benjamin Peherstorfer
Stability is a basic requirement when studying the behavior of dynamical systems. However, stabilizing dynamical systems via reinforcement learning is challenging because only little data can be collected over short time horizons before instabilities are triggered and data become meaningless. This work introduces a reinforcement learning approach that is formulated over latent manifolds of unstable dynamics so that stabilizing policies can be trained from few data samples. The unstable manifolds are minimal in the sense that they contain the lowest dimensional dynamics that are necessary for learning policies that guarantee stabilization. This is in stark contrast to generic latent manifolds that aim to approximate all—stable and unstable—system dynamics and thus are higher dimensional and often require higher amounts of data. Experiments demonstrate that the proposed approach stabilizes even complex physical systems from few data samples for which other methods that operate either directly in the system state space or on generic latent manifolds fail.
稳定性是研究动态系统行为的基本要求。然而,通过强化学习来稳定动态系统具有挑战性,因为在不稳定性触发和数据变得毫无意义之前,只能在短时间内收集到少量数据。这项工作引入了一种强化学习方法,该方法是在不稳定动力学的潜在流形上制定的,因此可以从少量数据样本中训练出稳定策略。不稳定流形是最小的,因为它们包含最低维度的动态,而这些动态是学习保证稳定的策略所必需的。这与一般的潜在流形形成了鲜明对比,后者旨在近似所有稳定和不稳定系统动态,因此维度更高,通常需要更多数据。实验证明,即使是复杂的物理系统,所提出的方法也能从少量数据样本中实现稳定,而其他直接在系统状态空间或通用潜流形上运行的方法却无法做到这一点。
{"title":"System stabilization with policy optimization on unstable latent manifolds","authors":"Steffen W.R. Werner ,&nbsp;Benjamin Peherstorfer","doi":"10.1016/j.cma.2024.117483","DOIUrl":"10.1016/j.cma.2024.117483","url":null,"abstract":"<div><div>Stability is a basic requirement when studying the behavior of dynamical systems. However, stabilizing dynamical systems via reinforcement learning is challenging because only little data can be collected over short time horizons before instabilities are triggered and data become meaningless. This work introduces a reinforcement learning approach that is formulated over latent manifolds of unstable dynamics so that stabilizing policies can be trained from few data samples. The unstable manifolds are minimal in the sense that they contain the lowest dimensional dynamics that are necessary for learning policies that guarantee stabilization. This is in stark contrast to generic latent manifolds that aim to approximate all—stable and unstable—system dynamics and thus are higher dimensional and often require higher amounts of data. Experiments demonstrate that the proposed approach stabilizes even complex physical systems from few data samples for which other methods that operate either directly in the system state space or on generic latent manifolds fail.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117483"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142552280","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 numerical algorithm to solve hydrodynamic lubrication problems with cavitation 解决带有气蚀的流体动力润滑问题的高效数值算法
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117470
P. Gómez-Molina , L. Sanz-Lorenzo , J. Carpio
In this paper we present an efficient numerical algorithm to solve stationary problems of hydrodynamic lubrication with cavitation in bearings using the method of characteristics in a finite element framework. The problem is based on the Elrod–Adams mathematical model for the lubricant fluid behavior. To achieve realistic pressure solutions, cavitation must be considered. This leads to a non-linear system of equations including a multivalued operator. In order to solve this problem, we propose a modified Newton algorithm that presents a high convergence speed, allowing one to use small tolerances to solve the problem with high accuracy and efficiency using a low number of iterations, in comparison with alternative numerical algorithms in the literature to solve the same problem. Numerical results are presented and analyzed in the work.
本文提出了一种高效的数值算法,可在有限元框架内使用特征法解决轴承中带有气蚀的流体动力润滑静态问题。该问题基于润滑油流体行为的 Elrod-Adams 数学模型。为了获得逼真的压力解决方案,必须考虑气蚀问题。这将导致一个包含多值算子的非线性方程组。为了解决这个问题,我们提出了一种改进的牛顿算法,与文献中解决相同问题的其他数值算法相比,该算法收敛速度快,允许使用较小的公差,以较低的迭代次数高精度、高效率地解决问题。工作中展示并分析了数值结果。
{"title":"An efficient numerical algorithm to solve hydrodynamic lubrication problems with cavitation","authors":"P. Gómez-Molina ,&nbsp;L. Sanz-Lorenzo ,&nbsp;J. Carpio","doi":"10.1016/j.cma.2024.117470","DOIUrl":"10.1016/j.cma.2024.117470","url":null,"abstract":"<div><div>In this paper we present an efficient numerical algorithm to solve stationary problems of hydrodynamic lubrication with cavitation in bearings using the method of characteristics in a finite element framework. The problem is based on the Elrod–Adams mathematical model for the lubricant fluid behavior. To achieve realistic pressure solutions, cavitation must be considered. This leads to a non-linear system of equations including a multivalued operator. In order to solve this problem, we propose a modified Newton algorithm that presents a high convergence speed, allowing one to use small tolerances to solve the problem with high accuracy and efficiency using a low number of iterations, in comparison with alternative numerical algorithms in the literature to solve the same problem. Numerical results are presented and analyzed in the work.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117470"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142550564","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
Second-order computational homogenization for bridging poromechanical scales under large deformations 在大变形条件下桥接孔力学尺度的二阶计算均质化
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117481
J.L.M. Thiesen , B. Klahr , T.A. Carniel , G.A. Holzapfel , P.J. Blanco , E.A. Fancello
We introduce a second-order computational homogenization procedure designed to address heterogeneous poromechanical media. Our approach relies on the method of multiscale virtual power, a variational multiscale method that extends the Hill–Mandel principle of macro-homogeneity. Constraints on displacement and pore pressure fields are managed using periodic and second-order minimally constrained fluctuating spaces. Numerical comparisons reveal that first-order models fail to accurately represent nonzero net fluid flow and volume changes at the micro-scale. In contrast, our second-order approach effectively captures nonuniform fluid flow across representative volume element boundaries, in agreement with results from direct numerical simulations. Our findings indicate that the classical first-order expansion of the pressure field is inadequate for poromechanical homogenization in cases involving micro-scale volume changes, such as swelling or contraction. The proposed second-order approach not only overcomes these limitations but also proves effective in cases where the principle of separation of scales is not strictly upheld.
我们介绍了一种二阶计算均质化程序,旨在解决异质孔隙力学介质问题。我们的方法依赖于多尺度虚拟功率方法,这是一种扩展了宏观均质性希尔-曼德尔原理的变分多尺度方法。对位移和孔隙压力场的约束采用周期和二阶最小约束波动空间进行管理。数值比较显示,一阶模型无法准确表示微尺度的非零净流体流量和体积变化。相比之下,我们的二阶方法能有效捕捉跨代表性体积元素边界的非均匀流体流动,这与直接数值模拟的结果一致。我们的研究结果表明,在涉及微尺度体积变化(如膨胀或收缩)的情况下,压力场的经典一阶扩展不足以实现孔力学均质化。所提出的二阶方法不仅克服了这些局限性,而且在不严格遵守尺度分离原则的情况下也证明是有效的。
{"title":"Second-order computational homogenization for bridging poromechanical scales under large deformations","authors":"J.L.M. Thiesen ,&nbsp;B. Klahr ,&nbsp;T.A. Carniel ,&nbsp;G.A. Holzapfel ,&nbsp;P.J. Blanco ,&nbsp;E.A. Fancello","doi":"10.1016/j.cma.2024.117481","DOIUrl":"10.1016/j.cma.2024.117481","url":null,"abstract":"<div><div>We introduce a second-order computational homogenization procedure designed to address heterogeneous poromechanical media. Our approach relies on the method of multiscale virtual power, a variational multiscale method that extends the Hill–Mandel principle of macro-homogeneity. Constraints on displacement and pore pressure fields are managed using periodic and second-order minimally constrained fluctuating spaces. Numerical comparisons reveal that first-order models fail to accurately represent nonzero net fluid flow and volume changes at the micro-scale. In contrast, our second-order approach effectively captures nonuniform fluid flow across representative volume element boundaries, in agreement with results from direct numerical simulations. Our findings indicate that the classical first-order expansion of the pressure field is inadequate for poromechanical homogenization in cases involving micro-scale volume changes, such as swelling or contraction. The proposed second-order approach not only overcomes these limitations but also proves effective in cases where the principle of separation of scales is not strictly upheld.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117481"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142550566","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
Direct serendipity finite elements on cuboidal hexahedra 立方体六面体上的直接偶然性有限元
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117500
Todd Arbogast , Chuning Wang
We construct direct serendipity finite elements on general cuboidal hexahedra, which are H1-conforming and optimally approximate to any order. The new finite elements are direct in that the shape functions are directly defined on the physical element. Moreover, they are serendipity by possessing a minimal number of degrees of freedom satisfying the conformity requirement. Their shape function spaces consist of polynomials plus (generally nonpolynomial) supplemental functions, where the polynomials are included for the approximation property and supplements are added to achieve H1-conformity. The finite elements are fully constructive. The shape function spaces of higher order r3 are developed first, and then the lower order spaces are constructed as subspaces of the third order space. Under a shape regularity assumption, and a mild restriction on the choice of supplemental functions, we develop the convergence properties of the new direct serendipity finite elements. Numerical results with different choices of supplements are compared on two mesh sequences, one regularly distorted and the other one randomly distorted. They all possess a convergence rate that aligns with the theory, while a slight difference lies in their performance.
我们在一般立方体六面体上构建了直接的偶然性有限元,这些有限元具有 H1 顺应性,并能以最佳方式逼近任何阶次。新的有限元是直接的,因为形状函数是直接定义在物理元素上的。此外,它们还具有偶然性,因为它们拥有满足一致性要求的最少自由度。它们的形状函数空间由多项式加上(一般为非多项式)补充函数组成,其中多项式是为了逼近特性而加入的,而补充函数则是为了实现 H1 一致性而加入的。有限元是完全构造的。首先建立高阶 r≥3 的形状函数空间,然后作为三阶空间的子空间构建低阶空间。在形状正则性假设和对补充函数选择的温和限制下,我们发展了新的直接偶然性有限元的收敛特性。在两个网格序列(一个是规则变形网格,另一个是随机变形网格)上比较了不同补充函数选择的数值结果。它们都具有与理论一致的收敛速度,但在性能上略有不同。
{"title":"Direct serendipity finite elements on cuboidal hexahedra","authors":"Todd Arbogast ,&nbsp;Chuning Wang","doi":"10.1016/j.cma.2024.117500","DOIUrl":"10.1016/j.cma.2024.117500","url":null,"abstract":"<div><div>We construct direct serendipity finite elements on general cuboidal hexahedra, which are <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming and optimally approximate to any order. The new finite elements are <em>direct</em> in that the shape functions are directly defined on the physical element. Moreover, they are <em>serendipity</em> by possessing a minimal number of degrees of freedom satisfying the conformity requirement. Their shape function spaces consist of polynomials plus (generally nonpolynomial) supplemental functions, where the polynomials are included for the approximation property and supplements are added to achieve <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conformity. The finite elements are fully constructive. The shape function spaces of higher order <span><math><mrow><mi>r</mi><mo>≥</mo><mn>3</mn></mrow></math></span> are developed first, and then the lower order spaces are constructed as subspaces of the third order space. Under a shape regularity assumption, and a mild restriction on the choice of supplemental functions, we develop the convergence properties of the new direct serendipity finite elements. Numerical results with different choices of supplements are compared on two mesh sequences, one regularly distorted and the other one randomly distorted. They all possess a convergence rate that aligns with the theory, while a slight difference lies in their performance.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117500"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142551904","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
An assemblable interlocking joint generation method for multi-material topology optimization using interfacial partial stress constraints and dimensional constraints 利用界面局部应力约束和尺寸约束进行多材料拓扑优化的可装配联锁接头生成方法
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117475
Yukun Feng , Takayuki Yamada
Multi-material topology optimization has become a promising method in structural design due to its excellent structural performance. However, existing research assumes that the multi-material structures are joined by welding, adhesive, or other methods that do not support reassembly and disassembly and are unsuitable for manufacturing, limiting the practical application of topology optimization. An interlocking joint is a type of connection between two parts where the shapes of the parts are designed to fit together precisely, the multi-material structure joined by interlocking joints can be easily reassembled repeatedly. To solve the joint problem of multi-material structure, this study proposes an assemblable interlocking joint generation method for multi-material topology optimization, the connection between material components is achieved through compression at the joint areas. To generate the interlocking joints, a novel interfacial partial stress constraint is proposed by converting a part of the interface bearing tensile stress into the interface bearing compressive stress. A novel filtering process is used to control the shape of the interlocking joints and the filtered tensile stresses are integrated by a P-norm function. To constrain the distribution area of material components and ensure structural manufacturability, dimensional constraints are applied. The sensitivity is based on the topological derivative and adjoint variable method. The proposed method was applied to several numerical examples including one manufactured prototype to demonstrate its effectiveness and contribution to the practical application of topology optimization.
多材料拓扑优化因其优异的结构性能已成为结构设计中一种前景广阔的方法。然而,现有研究假设多材料结构通过焊接、粘合或其他方法连接,这些方法不支持重新组装和拆卸,不适合制造,限制了拓扑优化的实际应用。联锁接头是两部分零件之间的一种连接方式,零件的形状被设计成可以精确地配合在一起,通过联锁接头连接的多材料结构可以很容易地重复组装。为了解决多材料结构的连接问题,本研究提出了一种用于多材料拓扑优化的可装配互锁连接生成方法,通过在连接区域的压缩实现材料部件之间的连接。为生成互锁接头,提出了一种新的界面部分应力约束,即将界面承受的部分拉应力转换为界面承受的压应力。新颖的过滤过程用于控制互锁接缝的形状,过滤后的拉应力通过 P-norm 函数进行积分。为了限制材料成分的分布面积并确保结构的可制造性,应用了尺寸约束。灵敏度基于拓扑导数和邻接变量法。所提出的方法已应用于多个数值实例,包括一个已制造的原型,以证明其有效性以及对拓扑优化实际应用的贡献。
{"title":"An assemblable interlocking joint generation method for multi-material topology optimization using interfacial partial stress constraints and dimensional constraints","authors":"Yukun Feng ,&nbsp;Takayuki Yamada","doi":"10.1016/j.cma.2024.117475","DOIUrl":"10.1016/j.cma.2024.117475","url":null,"abstract":"<div><div>Multi-material topology optimization has become a promising method in structural design due to its excellent structural performance. However, existing research assumes that the multi-material structures are joined by welding, adhesive, or other methods that do not support reassembly and disassembly and are unsuitable for manufacturing, limiting the practical application of topology optimization. An interlocking joint is a type of connection between two parts where the shapes of the parts are designed to fit together precisely, the multi-material structure joined by interlocking joints can be easily reassembled repeatedly. To solve the joint problem of multi-material structure, this study proposes an assemblable interlocking joint generation method for multi-material topology optimization, the connection between material components is achieved through compression at the joint areas. To generate the interlocking joints, a novel interfacial partial stress constraint is proposed by converting a part of the interface bearing tensile stress into the interface bearing compressive stress. A novel filtering process is used to control the shape of the interlocking joints and the filtered tensile stresses are integrated by a P-norm function. To constrain the distribution area of material components and ensure structural manufacturability, dimensional constraints are applied. The sensitivity is based on the topological derivative and adjoint variable method. The proposed method was applied to several numerical examples including one manufactured prototype to demonstrate its effectiveness and contribution to the practical application of topology optimization.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117475"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142550563","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A new floating node-based element formulation for modelling pressure-driven fracture 用于模拟压力驱动型断裂的基于浮动节点的新要素配方
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-30 DOI: 10.1016/j.cma.2024.117482
E.S. Kocaman , B.Y. Chen , S.T. Pinho
When simulating pressure-driven fracture with the Finite Element Method (FEM), significant difficulties can arise upon representing newly formed complex damage surfaces and their concurrent crack face loading. Application of this loading can also be required when additional physics is involved as in the case of hydraulic fracture where fluid physics inside a damage need to be considered. This paper presents a new Finite Element based practical numerical framework which can model pressure-driven fractures as they form on-the-fly without remeshing. The exact location of physical discontinuities passing through the element domain can be represented in the numerical model. The numerical framework can be implemented as a user-defined element and can be integrated into any FE package. A new element (called pressure element) is formulated with the capability to apply pressure and associated forces onto the crack surfaces in an adaptive manner. This element is verified using relevant examples from literature. The framework can also be configured for multi-physics problems where crack face loading is dictated by an additional physics. The element formulation is then extended for multi-physics problems involving fluid–solid interaction. The formulation provides the capability for multi-physics coupling adaptively as the crack propagates. The element is used to successfully simulate a test case from literature using different solution procedures (iterative and simultaneous). This element is also used to model failure in different pressure vessel problems to demonstrate its potential use in structural applications. A new higher-scale vessel element is developed which can represent different size, partitioning and failure states of composite vessel systems at element level. Composite vessel failure involving high number of pressurized cracks and delaminations as well as their interaction is modelled, and burst pressures are predicted for different vessel systems. The proposed numerical framework can be used towards designing more damage-tolerant vessels critical for the sustainable propulsion technologies.
用有限元法(FEM)模拟压力驱动断裂时,在表示新形成的复杂损伤面及其同时产生的裂纹面加载时会遇到很大困难。当涉及额外的物理因素时,也需要应用这种加载,例如在水力压裂中,需要考虑损伤内部的流体物理因素。本文介绍了一种新的基于有限元的实用数值框架,该框架可在压力驱动型断裂形成时对其进行建模,而无需重塑。通过元素域的物理不连续性的确切位置可以在数值模型中表示出来。该数值框架可作为用户定义的元素实施,并可集成到任何有限元软件包中。制定的新要素(称为压力要素)能够以自适应的方式在裂缝表面施加压力和相关力。利用文献中的相关实例对该元素进行了验证。该框架还可配置为多物理场问题,其中裂纹面加载由附加物理场决定。然后,该元素公式可扩展用于涉及流固相互作用的多物理场问题。随着裂缝的扩展,该公式可提供多物理场自适应耦合的能力。该元素使用不同的求解程序(迭代和同步)成功模拟了文献中的一个测试案例。该元素还用于模拟不同压力容器问题中的故障,以证明其在结构应用中的潜在用途。开发了一种新的更高尺度容器元素,可在元素级表示复合容器系统的不同尺寸、分区和失效状态。对涉及大量受压裂缝和分层以及它们之间相互作用的复合材料容器失效进行了建模,并预测了不同容器系统的爆破压力。所提出的数值框架可用于设计对可持续推进技术至关重要的更耐损伤容器。
{"title":"A new floating node-based element formulation for modelling pressure-driven fracture","authors":"E.S. Kocaman ,&nbsp;B.Y. Chen ,&nbsp;S.T. Pinho","doi":"10.1016/j.cma.2024.117482","DOIUrl":"10.1016/j.cma.2024.117482","url":null,"abstract":"<div><div>When simulating pressure-driven fracture with the Finite Element Method (FEM), significant difficulties can arise upon representing newly formed complex damage surfaces and their concurrent crack face loading. Application of this loading can also be required when additional physics is involved as in the case of hydraulic fracture where fluid physics inside a damage need to be considered. This paper presents a new Finite Element based practical numerical framework which can model pressure-driven fractures as they form on-the-fly without remeshing. The exact location of physical discontinuities passing through the element domain can be represented in the numerical model. The numerical framework can be implemented as a user-defined element and can be integrated into any FE package. A new element (called <em>pressure</em> element) is formulated with the capability to apply pressure and associated forces onto the crack surfaces in an adaptive manner. This element is verified using relevant examples from literature. The framework can also be configured for multi-physics problems where crack face loading is dictated by an additional physics. The element formulation is then extended for multi-physics problems involving fluid–solid interaction. The formulation provides the capability for multi-physics coupling adaptively as the crack propagates. The element is used to successfully simulate a test case from literature using different solution procedures (iterative and simultaneous). This element is also used to model failure in different pressure vessel problems to demonstrate its potential use in structural applications. A new higher-scale <em>vessel</em> element is developed which can represent different size, partitioning and failure states of composite vessel systems at element level. Composite vessel failure involving high number of pressurized cracks and delaminations as well as their interaction is modelled, and burst pressures are predicted for different vessel systems. The proposed numerical framework can be used towards designing more damage-tolerant vessels critical for the sustainable propulsion technologies.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117482"},"PeriodicalIF":6.9,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142550567","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 variational phase-field framework for thermal softening and dynamic ductile fracture 热软化和动态韧性断裂的变分相场框架
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-28 DOI: 10.1016/j.cma.2024.117452
David E. Torres , Tianchen Hu , Andrew J. Stershic , Timothy R. Shelton , John E. Dolbow
A variational phase field model for dynamic ductile fracture is presented. The model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.
本文介绍了动态韧性断裂的变分相场模型。该模型是针对受快速变形影响的弹塑性材料设计的,在快速变形过程中,热生成和材料软化的影响占主导地位。变分框架允许将塑性耗散和热软化一致纳入热方程。它采用了凝聚函数来降低高塑性流动时的断裂能量。为与该框架配合使用,开发了约翰逊-库克(Johnson-Cook)模型的变化一致形式。本文介绍了动态韧性断裂中各种基准问题的结果,以展示其能力。特别是展示了该模型在二维和三维问题中规范剪切带形成和后续损伤演变的能力。重要的是,这些现象可以通过底层物理自然地捕捉到,而不需要现象学标准,如剪切带形成的稳定性阈值。
{"title":"A variational phase-field framework for thermal softening and dynamic ductile fracture","authors":"David E. Torres ,&nbsp;Tianchen Hu ,&nbsp;Andrew J. Stershic ,&nbsp;Timothy R. Shelton ,&nbsp;John E. Dolbow","doi":"10.1016/j.cma.2024.117452","DOIUrl":"10.1016/j.cma.2024.117452","url":null,"abstract":"<div><div>A variational phase field model for dynamic ductile fracture is presented. The model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117452"},"PeriodicalIF":6.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142529236","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 study on the energy consistency in SPH surface tension modelling SPH 表面张力建模中的能量一致性研究
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-28 DOI: 10.1016/j.cma.2024.117473
S. Marrone , M. Antuono , A. Agresta , A. Colagrossi
In the present work the evolution of viscous drops oscillating under the action of surface tension is tackled. Thanks to its structure, the SPH scheme allows for an analysis of the energy balance that is rarely addressed to in the general Computational Fluid Dynamics literature for this kind of flows. A procedure for checking the consistency between the energy of the surface-tension force and the free-surface evolution is proposed. Such a procedure relies on well-known analytical relations for the surface tension and on the evaluation of the free-surface area through a level-set function. Several test cases, in both two and three dimensional frameworks, are considered for validation. The study is performed by selecting a specific SPH scheme with a specific single-phase surface tension model. In any case, the procedure proposed is general and extendable to other SPH surface tension models and SPH schemes.
本研究探讨了在表面张力作用下摆动的粘性液滴的演变过程。由于其结构,SPH 方案允许对能量平衡进行分析,而一般的计算流体力学文献很少涉及这类流动。本文提出了一种检查表面张力能量与自由表面演变之间一致性的程序。这种程序依赖于众所周知的表面张力分析关系,以及通过水平集函数对自由表面积的评估。在二维和三维框架下,考虑了几个测试案例进行验证。研究是通过选择特定的 SPH 方案和特定的单相表面张力模型进行的。无论如何,所提出的程序是通用的,可扩展到其他 SPH 表面张力模型和 SPH 方案。
{"title":"A study on the energy consistency in SPH surface tension modelling","authors":"S. Marrone ,&nbsp;M. Antuono ,&nbsp;A. Agresta ,&nbsp;A. Colagrossi","doi":"10.1016/j.cma.2024.117473","DOIUrl":"10.1016/j.cma.2024.117473","url":null,"abstract":"<div><div>In the present work the evolution of viscous drops oscillating under the action of surface tension is tackled. Thanks to its structure, the SPH scheme allows for an analysis of the energy balance that is rarely addressed to in the general Computational Fluid Dynamics literature for this kind of flows. A procedure for checking the consistency between the energy of the surface-tension force and the free-surface evolution is proposed. Such a procedure relies on well-known analytical relations for the surface tension and on the evaluation of the free-surface area through a level-set function. Several test cases, in both two and three dimensional frameworks, are considered for validation. The study is performed by selecting a specific SPH scheme with a specific single-phase surface tension model. In any case, the procedure proposed is general and extendable to other SPH surface tension models and SPH schemes.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117473"},"PeriodicalIF":6.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142529322","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 extra-dof-free generalized finite element method for incompressible Navier-Stokes equations 不可压缩纳维-斯托克斯方程的无差别广义有限元法
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-28 DOI: 10.1016/j.cma.2024.117466
Wenhai Sheng , Qinglin Duan
The generalized finite element method (GFEM) without extra degrees of freedom (dof) is extended to solve incompressible Navier-Stokes (N-S) equations. Unlike the existing extra-dof-free GFEM, we propose a new approach to construct the nodal enrichments based on the weighted least-squares. As a result, the essential boundary conditions can be imposed more accurately. The Characteristic-Based Split (CBS) scheme is used to suppress oscillations due to the standard Galerkin discretization of the convective terms, and the pressure is further stabilized by the finite increment calculus (FIC) formulation. Hence, equal velocity-pressure interpolation and the incremental version of the split scheme can be used without inducing spurious oscillations. The developed extra-dof-free GFEM is very flexible and can achieve high-order spatial accuracy and convergence rates by adopting high-order polynomial enrichments. In particular, better accuracy could be obtained with special enrichments reflecting a-priori knowledge about the solution. This is demonstrated by numerical results. Benchmark examples such as the Lid-Driven Cavity flow and the flow past a circular cylinder are also presented to further verify the effectiveness of the proposed extra-dof-free GFEM for incompressible flow.
无额外自由度(dof)的广义有限元法(GFEM)被扩展用于求解不可压缩纳维-斯托克斯(N-S)方程。与现有的无额外自由度 GFEM 不同,我们提出了一种基于加权最小二乘法构建节点富集的新方法。因此,可以更精确地施加基本边界条件。基于特性的分割(CBS)方案用于抑制由于对流项的标准 Galerkin 离散化而产生的振荡,而压力则通过有限增量微积分(FIC)公式得到进一步稳定。因此,使用等速等压插值和增量版分割方案不会引起虚假振荡。所开发的无差错外 GFEM 非常灵活,通过采用高阶多项式富集,可以获得高阶空间精度和收敛速率。特别是,通过反映解的先验知识的特殊富集,可以获得更好的精度。数值结果证明了这一点。此外,还介绍了一些基准实例,如Lid-Driven Cavity流动和圆筒流动,以进一步验证所提出的无损外GFEM在不可压缩流动中的有效性。
{"title":"An extra-dof-free generalized finite element method for incompressible Navier-Stokes equations","authors":"Wenhai Sheng ,&nbsp;Qinglin Duan","doi":"10.1016/j.cma.2024.117466","DOIUrl":"10.1016/j.cma.2024.117466","url":null,"abstract":"<div><div>The generalized finite element method (GFEM) without extra degrees of freedom (dof) is extended to solve incompressible Navier-Stokes (N-S) equations. Unlike the existing extra-dof-free GFEM, we propose a new approach to construct the nodal enrichments based on the weighted least-squares. As a result, the essential boundary conditions can be imposed more accurately. The Characteristic-Based Split (CBS) scheme is used to suppress oscillations due to the standard Galerkin discretization of the convective terms, and the pressure is further stabilized by the finite increment calculus (FIC) formulation. Hence, equal velocity-pressure interpolation and the incremental version of the split scheme can be used without inducing spurious oscillations. The developed extra-dof-free GFEM is very flexible and can achieve high-order spatial accuracy and convergence rates by adopting high-order polynomial enrichments. In particular, better accuracy could be obtained with special enrichments reflecting a-priori knowledge about the solution. This is demonstrated by numerical results. Benchmark examples such as the Lid-Driven Cavity flow and the flow past a circular cylinder are also presented to further verify the effectiveness of the proposed extra-dof-free GFEM for incompressible flow.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117466"},"PeriodicalIF":6.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142529238","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
Lattice Boltzmann for linear elastodynamics: Periodic problems and Dirichlet boundary conditions 线性弹性力学的玻尔兹曼网格:周期问题和迪里夏特边界条件
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-28 DOI: 10.1016/j.cma.2024.117469
Oliver Boolakee , Martin Geier , Laura De Lorenzis
We propose a new second-order accurate lattice Boltzmann formulation for linear elastodynamics that is stable for arbitrary combinations of material parameters under a CFL-like condition. The construction of the numerical scheme uses an equivalent first-order hyperbolic system of equations as an intermediate step, for which a vectorial lattice Boltzmann formulation is introduced. The only difference to conventional lattice Boltzmann formulations is the usage of vector-valued populations, so that all computational benefits of the algorithm are preserved. Using the asymptotic expansion technique and the notion of pre-stability structures we further establish second-order consistency as well as analytical stability estimates. Lastly, we introduce a second-order consistent initialization of the populations as well as a boundary formulation for Dirichlet boundary conditions on 2D rectangular domains. All theoretical derivations are numerically verified by convergence studies using manufactured solutions and long-term stability tests.
我们为线性弹性动力学提出了一种新的二阶精确晶格玻尔兹曼公式,在类似 CFL 的条件下,这种公式对于任意材料参数组合都是稳定的。数值方案的构建使用了等效的一阶双曲方程组作为中间步骤,并为此引入了矢量晶格玻尔兹曼公式。与传统玻尔兹曼网格公式的唯一区别是使用了矢量值种群,从而保留了算法的所有计算优势。利用渐近展开技术和前稳定性结构概念,我们进一步建立了二阶一致性和分析稳定性估计。最后,我们介绍了种群的二阶一致性初始化,以及二维矩形域上 Dirichlet 边界条件的边界表述。所有理论推导都通过使用人造解的收敛性研究和长期稳定性测试进行了数值验证。
{"title":"Lattice Boltzmann for linear elastodynamics: Periodic problems and Dirichlet boundary conditions","authors":"Oliver Boolakee ,&nbsp;Martin Geier ,&nbsp;Laura De Lorenzis","doi":"10.1016/j.cma.2024.117469","DOIUrl":"10.1016/j.cma.2024.117469","url":null,"abstract":"<div><div>We propose a new second-order accurate lattice Boltzmann formulation for linear elastodynamics that is stable for arbitrary combinations of material parameters under a CFL-like condition. The construction of the numerical scheme uses an equivalent first-order hyperbolic system of equations as an intermediate step, for which a vectorial lattice Boltzmann formulation is introduced. The only difference to conventional lattice Boltzmann formulations is the usage of vector-valued populations, so that all computational benefits of the algorithm are preserved. Using the asymptotic expansion technique and the notion of pre-stability structures we further establish second-order consistency as well as analytical stability estimates. Lastly, we introduce a second-order consistent initialization of the populations as well as a boundary formulation for Dirichlet boundary conditions on 2D rectangular domains. All theoretical derivations are numerically verified by convergence studies using manufactured solutions and long-term stability tests.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117469"},"PeriodicalIF":6.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142529237","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 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学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1