{"title":"An unsymmetric 4-node quadrilateral Reissner–Mindlin plate element using radial–polynomial interpolation for linear and nonlinear analyses","authors":"Yan-Liang Ju, Ying-Qing Huang, Hai-Bo Chen","doi":"10.1016/j.finel.2025.104408","DOIUrl":"10.1016/j.finel.2025.104408","url":null,"abstract":"<div><div>Traditional Reissner–Mindlin plate elements often encounter considerable challenges, particularly their sensitivity to mesh distortion and limited accuracy in stress prediction. This study introduces a novel unsymmetric quadrilateral plate element formulation, extending its application to doubly-curved shells and geometrically nonlinear analysis. The formulation differentiates between test and trial functions to independently construct virtual and real displacement fields, respectively. The virtual displacement field is constructed using standard isoparametric interpolation combined with the mixed interpolation of tensorial components method, effectively suppressing shear locking whilst maintaining inter-element displacement continuity. This construction technique also retains the advantages of isoparametric elements in applying boundary conditions and calculating equivalent nodal external force. Meanwhile, the real displacement field is generated through a radial–polynomial interpolation strategy, using element supports to substantially improve inter-element stress continuity. Numerical examples confirm that the new element successfully eliminates shear locking, exhibits high resistance to mesh distortion, achieves high stress accuracy and ensures excellent inter-element stress continuity.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104408"},"PeriodicalIF":3.5,"publicationDate":"2025-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144687405","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Paul Larousse , David Dureisseix , Anthony Gravouil , Jean Di Stasio
{"title":"A thermodynamically based explicit transient dynamics framework for large transformation contact problems","authors":"Paul Larousse , David Dureisseix , Anthony Gravouil , Jean Di Stasio","doi":"10.1016/j.finel.2025.104411","DOIUrl":"10.1016/j.finel.2025.104411","url":null,"abstract":"<div><div>An explicit framework to solve fast dynamic problems with large transformation and rigid-deformable contact, involving non-regular and non-linear behaviors is under concern. Based on previous works, a framework combining thermodynamically-based behaviors and the so-called explicit symplectic time integrator CD-Lagrange owning good energy properties is developed. In this article, the interface behavior is modeled with a cohesive zone model, the RCCM delayed damage model, and for the deformable body, a hyper-elastic Saint-Venant–Kirchhoff model coupled with viscous effects is chosen. The modular proposed framework shows that switching between a large transformation or a small perturbation problem (and a wide range of non-linear laws both on the interface or in the bulk) is non-intrusive in terms of numerical implementation and dedicated spatial interface integration in the framework of finite elements. In this work, illustrations and feasibility are exemplified for a simplified unmolding industrial case.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104411"},"PeriodicalIF":3.5,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144680254","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A novel mesh-reduction technique for deriving closed-form solutions of framed structures using a single finite element per structural member","authors":"Juan Camilo Molina-Villegas, Cristian Posso","doi":"10.1016/j.finel.2025.104406","DOIUrl":"10.1016/j.finel.2025.104406","url":null,"abstract":"<div><div>This paper presents a novel mesh reduction technique designed to enhance the accuracy of Finite Element Method (FEM) results for the analysis of framed structures idealized using rod and beam theories. The proposed method facilitates the derivation of closed-form solutions, including exact nodal displacements, reaction forces, displacement and internal forces response fields within structural elements, while employing a single finite element per structural member. Inspired by the Green’s Functions Stiffness Method, this approach is based on the decomposition of the element’s response into homogeneous and particular or fixed components, both calculated exclusively using analytical shape functions and the stiffness matrix coefficients, which are key inputs in standard FEM formulations for framed elements. This technique can be seamlessly integrated into the post-processing stage of traditional FEM implementations, offering improved computational efficiency and precision of results.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104406"},"PeriodicalIF":3.5,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144672256","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Cheap and stable quadrature on polyhedral elements","authors":"Alvise Sommariva, Marco Vianello","doi":"10.1016/j.finel.2025.104409","DOIUrl":"10.1016/j.finel.2025.104409","url":null,"abstract":"<div><div>We discuss a cheap tetrahedra-free approach to the numerical integration of polynomials on polyhedral elements, based on hyperinterpolation in a bounding box and Chebyshev moment computation via the divergence theorem. No conditioning issues arise, since no matrix factorization or inversion is needed. The resulting quadrature formula is theoretically stable even in the presence of some negative weights.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104409"},"PeriodicalIF":3.5,"publicationDate":"2025-07-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144664703","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A. Francisca Carvalho Alves , Bernardo P. Ferreira , F.M. Andrade Pires
{"title":"Multi-scale optimization of PC/ABS polymer blends: Microstructural design for superior toughness, strength, and weight efficiency","authors":"A. Francisca Carvalho Alves , Bernardo P. Ferreira , F.M. Andrade Pires","doi":"10.1016/j.finel.2025.104407","DOIUrl":"10.1016/j.finel.2025.104407","url":null,"abstract":"<div><div>The PC/ABS polymer blend is widely used in automotive and consumer electronics due to its balanced combination of thermal, mechanical, and processing properties. Its behavior depends on deformation mechanisms such as rubber particle cavitation and debonding at the PC/ABS interface, which vary with loading conditions and morphology. Modeling and optimizing the PC/ABS microstructure is a complex challenge. This work proposes a multi-scale framework to model and optimize different PC/ABS blends, based on: (i) efficient generation of representative volume elements, (ii) accurate constitutive models for the blend phases, and (iii) an unsupervised optimization process for microstructural design. The optimization considers ABS content in the blend, rubber fraction in ABS, and ABS particle orientation to maximize toughness and strength while minimizing cost and weight — a challenging task due to the negative correlation between toughness and strength. To handle the high-dimensional solution space, Multi-Criteria Decision Making techniques are employed to select optimal solutions within the Pareto front. Two case studies are explored: (i) a lightweight, high-toughness application and (ii) a high-strength application. Additionally, the framework is tested in a functionally graded material optimization problem, where a Cook’s membrane is discretized into PC/ABS layers, with the ABS fraction in each layer adjusted to simultaneously minimize maximum displacement and structural weight. The numerical results validate the proposed design framework for PC/ABS while demonstrating its flexibility for other materials and structural applications.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104407"},"PeriodicalIF":3.5,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633894","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A novel class of Hessian recovery-based numerical methods for solving biharmonic equations and their applications in phase field modeling","authors":"Minqiang Xu , Lei Zhang , Boying Wu , Kai Liu","doi":"10.1016/j.finel.2025.104405","DOIUrl":"10.1016/j.finel.2025.104405","url":null,"abstract":"<div><div>In this paper, we introduce unified Hessian recovery-based <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>0</mn></mrow></msup></math></span> finite element methods (HRB–FEM) and finite volume methods (HRB–FVM) for 2D biharmonic equations. Within the framework of Petrov–Galerkin methods, we propose a novel <span><math><mrow><msup><mrow><mi>H</mi></mrow><mrow><mn>3</mn></mrow></msup><mo>−</mo><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></mrow></math></span> formulation. Initially, we employ the Hessian recovery operator to discretize the Laplacian operator, subsequently integrating it into both the standard <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>0</mn></mrow></msup></math></span> Lagrange finite element framework and finite volume framework. Through tailored treatments of Neumann-type boundary conditions aimed at reducing computational overhead, we extend our Hessian recovery-based FEM to address phase field equations. Numerical experiments confirm optimal order of convergence under <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> and <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> norms, demonstrating rates of <span><math><mrow><mi>O</mi><mrow><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></mrow></math></span> and <span><math><mrow><mi>O</mi><mrow><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mi>k</mi></mrow></msup><mo>)</mo></mrow></mrow></math></span> respectively for both proposed methods. Furthermore, a series of benchmark tests highlight the robustness of our approach and its ability to faithfully capture the physical characteristics during prolonged simulations of phase field equations.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104405"},"PeriodicalIF":3.5,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144614377","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Reduced-integration hexahedral finite element for static and vibration analysis of micropolar continuum","authors":"Yu Yao , Xin Zhao , Linghao Chen , Yujie Gu , Tianqi Zhou","doi":"10.1016/j.finel.2025.104412","DOIUrl":"10.1016/j.finel.2025.104412","url":null,"abstract":"<div><div>The micropolar elasticity finite element method is widely used for analyzing advanced materials with complex microstructures, but existing implementations often suffer from computational inefficiency due to full integration schemes and shear locking in bending scenarios. This study proposes a high-performance, reduced-integration, first-order hexahedral micropolar element to address these limitations. The formulation combines standard Lagrange interpolation with uniform strain and curvature fields, ensuring patch test satisfaction and accuracy in skewed configurations. An artificial stiffness method is introduced to suppress displacement and rotational hourglass instabilities. Rigorous numerical validations, including force and displacement patch tests, cantilever beam bending, and free vibration analysis, demonstrate the superior accuracy and computational efficiency of the element. Furthermore, applications to star-shaped lattices and 3D chiral metamaterials highlight its effectiveness in capturing microstructure-dependent mechanical behaviors, such as unexpected bending deformation and tension-twist coupling. The proposed element significantly enhances computational efficiency in homogenization simulation, providing a robust and practical tool for the simulation-driven design of advanced mechanical metamaterials with complex deformation mechanisms.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"251 ","pages":"Article 104412"},"PeriodicalIF":3.5,"publicationDate":"2025-07-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144611919","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Three-dimensional simulation of finite-strain debonding using immersed meshes","authors":"Andrew B. Groeneveld , Pinlei Chen","doi":"10.1016/j.finel.2025.104404","DOIUrl":"10.1016/j.finel.2025.104404","url":null,"abstract":"<div><div>We propose a method for modeling interfacial damage and debonding under quasi-static loads using immersed meshes in 3D at finite strains. This is an extension of our previous work on an immersed variational multiscale discontinuous Galerkin (VMDG) method in 2D. The variational approach remains the same, but transitioning from 2D to 3D introduces significant complications in the computational geometry aspects. The immersed VMDG method is a stabilized interface formulation derived using variational multiscale (VMS) ideas to apply discontinuous Galerkin (DG) treatment to the interface while employing a continuous Galerkin (CG) approximation elsewhere. Key benefits of VMDG are the variationally derived stabilization terms that evolve during deformation and are free of user-defined parameters. Also, the transition from perfect bond to damage behavior at the interface is handled naturally by incorporating an interfacial gap variable governed by a yield criterion and a flow rule. To support 3D simulations, we introduce algorithms for integrating cut elements, forming interface segments, and computing the VMDG stabilization tensor. Cut-element integration is performed using voxel-based moment-fitting integration to avoid the robustness issues associated with using mesh Booleans and tetrahedral integration cells. A simplification of the stabilization tensor is also proposed to reduce the computational cost while retaining the variational character of the stabilization. Several numerical examples are presented to demonstrate the robustness, efficiency, and range of applicability of the method.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104404"},"PeriodicalIF":3.5,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144536086","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Runge–Kutta discontinuous Galerkin method based on flux vector splitting for hyperbolic conservation laws","authors":"Zhengrong Xie","doi":"10.1016/j.finel.2025.104398","DOIUrl":"10.1016/j.finel.2025.104398","url":null,"abstract":"<div><div>The flux vector splitting (FVS) method has firstly been incorporated into the Runge–Kutta Discontinuous Galerkin (RKDG) framework for reconstructing the numerical fluxes required for the spatial semi-discrete formulation, setting it apart from the conventional RKDG approaches that typically utilize the Lax–Friedrichs flux scheme or classical Riemann solvers such as HLLC. The control equations are initially reformulated into a flux-split form. Subsequently, a variational approach is applied to this flux-split form, from which a DG spatial semi-discrete scheme based on FVS is derived. Then, FVS-RKDG is implemented in two-dimensional case by splitting the normal flux on cell interfaces instead of splitting dimension by dimension in the x and y directions Finally, the concept of “flux vector splitting based on Jacobian eigenvalue decomposition” has been applied to the conservative linear scalar transport equations and the nonlinear Burgers’ equation. This approach has led to the rederivation of the classical Lax–Friedrichs flux scheme and the provision of a Steger–Warming flux scheme for scalar cases.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104398"},"PeriodicalIF":3.5,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144518096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Two-grid domain decomposition methods for the coupled Dual-Porosity-Navier-Stokes system with Beavers-Joseph interface condition","authors":"Chongxin Zhang, Guangzhi Du, Xinxin Sun","doi":"10.1016/j.finel.2025.104403","DOIUrl":"10.1016/j.finel.2025.104403","url":null,"abstract":"<div><div>In this paper, two kinds of two-grid domain decomposition methods for the coupled Dual-Porosity-Navier-Stokes system are proposed and analyzed by integrating the established robin-type domain decomposition approach with a two-grid strategy. Initially, we apply the established robin-type domain decomposition approach on a coarse grid to address the coupled problem. Subsequently, on a fine grid, we employ two distinct approaches: first, to solve the matrix and microfracture subproblems, followed by the Navier–Stokes subproblem. Both approaches fundamentally approximate the interface term using the coarse-grid solution. The proposed algorithms integrate the two-grid approach with the established domain decomposition method, capitalizing on the strengths of both techniques while addressing their respective limitations. Comprehensive theoretical analysis is established, and four in-depth numerical investigations are conducted to assess the efficiency, accuracy, and robustness of the proposed algorithms by comparing them with the domain decomposition method.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104403"},"PeriodicalIF":3.5,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144501860","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}