Finite Elements in Analysis and Design最新文献

{"title":"Computation of the elastodynamic response of finite doubly periodic cylinders by the wave finite element method","authors":"D. Duhamel","doi":"10.1016/j.finel.2026.104512","DOIUrl":"10.1016/j.finel.2026.104512","url":null,"abstract":"<div><div>Periodic media are widely studied for their industrial applications and unique ability to block wave propagation within certain frequency bands. For one-dimensional periodic systems, the Wave Finite Element (WFE) method efficiently computes dispersion relations and dynamic responses. Research on two-dimensional periodic structures extended this approach to periodicity along two directions, using finite element and reduction techniques such as Craig–Bampton and Bloch mode projection, including cases with damping and anisotropy. Beyond planar geometries, curved and helical periodic structures have been modeled with WFE and semi-analytical finite element methods to capture complex cyclic or screw symmetries for computing dispersion relations. Cylindrical configurations have also been explored, from simple vibration studies to wave propagation in layered or ribbed cylinders and metamaterial shells mainly for dispersion analysis or studies of infinite structures. As real structures are bounded, the present work focuses on finite elastic cylinders with double periodicity, using WFE to compute their dynamic response. Based on finite element matrices of a substructure, circumferential wavenumbers are imposed to obtain axial modes and responses as linear combination of modes. Numerical examples illustrate the method’s effectiveness for modeling finite, doubly periodic cylindrical systems such as homogeneous structures, structures with holes and finally structures with resonators. The low computing time of the present approach allows the consideration of structures with a large number of substructures.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104512"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145929284","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":"Isogeometric implementation of a single-variable shear deformable plate theory","authors":"Loc V. Tran ,&nbsp;Thang N. Dao ,&nbsp;Vuong Nguyen Van Do","doi":"10.1016/j.finel.2025.104502","DOIUrl":"10.1016/j.finel.2025.104502","url":null,"abstract":"<div><div>The shear deformation plate theory is applicable to moderate and thick plates by accounting for transverse shear effects. In an effort to optimize the number of unknowns, this study proposes a new form for shear-locking free Reissner-Mindlin plate theory that employs only a single variable–the bending deflection. Consequently, the governing equation is expressed as a fourth-order partial differential equation, retaining the same form as the classical Kirchhoff–Love theory, while fully considering shear deformation. Based on that, analytical solutions for transverse displacement are derived for rectangular plates with arbitrary slenderness ratios. Additionally, the weak form of the plate problem is derived and includes the second- and third-order derivatives. To address these higher-order continuity requirements, a conforming Galerkin method based on an isogeometric analysis (IGA) is adopted. In particular, the basis functions-based IGA with order <em>p</em> ≥ 3, naturally satisfy the <em>C</em>²-continuity requirement mandated by the proposed model. Moreover, these basis functions facilitate a straightforward enforcement of natural boundary conditions, such as prescribed slopes and curvatures, that are inherent in the present plate formulation. Numerical examples demonstrate that the proposed model, despite adopting a single unknown, provides highly accurate results for thin and thick plates and achieves high convergence rates for all quantities of interest.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104502"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145904202","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 displacement-based Timoshenko finite element of multi-cracked circular arch","authors":"Ilaria Fiore,&nbsp;Francesco Cannizzaro","doi":"10.1016/j.finel.2025.104510","DOIUrl":"10.1016/j.finel.2025.104510","url":null,"abstract":"<div><div>This study presents a Displacement-Based (DB) finite element formulation of the multi-cracked planar circular arch within the Timoshenko theory. The cracks presence is simulated with the localised flexibility approach, considering the tangential, radial and rotational kinematic jumps at the damaged sections. The kinematic discontinuities are modelled with Dirac's deltas in the generalised governing equations defined over a unique domain for a generic number of damaged sections, which represents a significant computational advantage over classic approaches requiring enforcement of continuity conditions at the discontinuous sections. The governing equations of the statics of the multi-cracked Timoshenko circular arch are integrated, leading to the closed-form response as a function of six boundary conditions only. Then, the static discontinuous shape functions are inferred and exploited to obtain the stiffness and consistent mass matrices of the multi-cracked circular arch, considering just three degrees of freedom for each end of the member irrespective of the amount of cracks along the axis of the element. The obtained circular multi-cracked element, which is completely locking-free, is implemented in a finite element environment that encompasses an analogous finite element introduced for cracked straight members, thus allowing the study of general planar damaged frames both in the static and dynamic contexts.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104510"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146014797","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":"Parametric reduced order models for graded lattice structures","authors":"Max Nezdyur ,&nbsp;Lynn Munday ,&nbsp;Wilkins Aquino","doi":"10.1016/j.finel.2026.104521","DOIUrl":"10.1016/j.finel.2026.104521","url":null,"abstract":"<div><div>Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. To address these challenges, this paper introduces a novel Element Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method’s modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to nuclear engineering.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104521"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145962065","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":"Modelling of 3D woven textile reinforced cement composites behaviour accounting for through-thickness reinforcement","authors":"Christian Toderascu ,&nbsp;Badadjida Wintiba ,&nbsp;Karim Ehab Moustafa Kamel ,&nbsp;Thierry J. Massart ,&nbsp;Tine Tysmans","doi":"10.1016/j.finel.2025.104509","DOIUrl":"10.1016/j.finel.2025.104509","url":null,"abstract":"<div><div>Textile Reinforced Cement (TRC) composites provide slender concrete material solutions. Using 3D textile reinforcements significantly improves the bending performance in the post-cracking stage compared to 2D textiles. Yet, no computational model takes into account explicitly the woven connections in 3D TRC. This contribution develops a novel strategy for generating the complex geometry of a mesoscale through-thickness Representative Volume Element (RVE) of 3D TRC shells, thanks to advanced geometric algorithms such as Rotation Minimising Frames (RMF), subsequently meshed and used in finite element (FE) analysis. The RVE with realistic reinforcement geometry enables the numerical evaluation of different fine scale processes contributing to the composite material performance, in particular the effect of the woven reinforcement architecture. RVE simulations under two types of bending loading conditions, using computational homogenisation procedures, illustrate that the proposed approach enables the investigation of the average macroscopic bending properties of 3D TRC based on the fine scale morphology of their reinforcement, together with an assessment of local strain fields.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104509"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145903371","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":"Implicit numerical integration of multisurface plasticity for both plane stress and three-dimensional stress conditions","authors":"Rafael Abreu ,&nbsp;Cristian Mejia ,&nbsp;Deane Roehl","doi":"10.1016/j.finel.2025.104500","DOIUrl":"10.1016/j.finel.2025.104500","url":null,"abstract":"<div><div>Formulating robust integration algorithms for elastoplastic models is crucial for efficient and accurate numerical simulations of materials such as concrete, rocks, and soil. While traditional elastoplastic models typically employ single yield surfaces, more intricate behaviors can be captured using multiple yield surfaces. On the other hand, implementing these models within a finite element framework requires sophisticated numerical methods, particularly implicit integration schemes based on the backward Euler method, to ensure accuracy and stability. In this context, this paper introduces a novel Newton-Raphson-based implicit integration algorithm for multisurface plasticity models, accommodating both plane stress and three-dimensional stress conditions. The proposed algorithm employs well-established smooth complementary functions to handle multisurface plasticity without needing a scheme to identify active surfaces. In addition, this algorithm addresses plane stress plasticity by modifying calculations based on plane strain conditions. The study includes an assessment of the computational efficiency of three complementary functions, considering various finite element problems. The robustness of the algorithm is demonstrated through a series of numerical experiments, highlighting its potential for challenging engineering applications.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104500"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145929283","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 third medium approach for thermo-mechanical contact based on low order ansatz spaces","authors":"P. Wriggers","doi":"10.1016/j.finel.2026.104522","DOIUrl":"10.1016/j.finel.2026.104522","url":null,"abstract":"<div><div>The third medium contact approach has been successfully employed in structural applications and extended to various optimization problems. This discretization technique replaces classical contact formulations and algorithms by introducing a compliant interfacial layer – referred to as the third medium – between the contacting bodies. Unlike traditional contact methods, this formulation naturally accommodates finite deformations at the interface. As the two bodies approach each other, the third medium undergoes compression and effectively acts as a deformable barrier, preventing interpenetration and transmitting contact forces in a smooth and numerically stable manner. In thermo-mechanical problems, heat conduction must be incorporated into the model, which typically requires specialized interface laws when using classical contact formulations. These laws aim to capture the complex thermal behavior at the contact interface, including discontinuities and varying conductance. In contrast, the third medium approach offers a significant advantage: the thermo-mechanical formulation inherently accounts for the interface behavior without the need for additional interface conditions. This includes the gradual heat transfer through the surrounding gas when the bodies are near each other, as well as the localized heat conduction that occurs upon physical contact. As a result, the third medium naturally captures both non-contact and contact-phase thermal conduction within a unified framework. In this paper, we propose a new thermo-mechanical model based on a continuum formulation for finite strains and show by means of examples the behavior of the associated finite element formulation based on linear ansatz functions.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104522"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145995152","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":"hp-adaptive finite element simulation of a static anti-plane shear crack in a nonlinear strain-limiting elastic solid","authors":"S.M. Mallikarjunaiah ,&nbsp;Pavithra Venkatachalapthy","doi":"10.1016/j.finel.2026.104520","DOIUrl":"10.1016/j.finel.2026.104520","url":null,"abstract":"<div><div>An <span><math><mrow><mi>h</mi><mi>p</mi></mrow></math></span>-adaptive continuous Galerkin finite element method is developed to analyze a static anti-plane shear crack embedded in a nonlinear, strain-limiting elastic body. The geometrically linear material is described by a constitutive law relating stress and strain that is algebraically nonlinear. Such a formulation is advantageous as it regularizes the stress and strain fields in the neighborhood of a crack-tip, thereby circumventing the non-physical strain singularities inherent to linear elastic fracture mechanics. In this investigation, the constitutive relation utilized is <em>uniformly bounded</em>, <em>monotone</em>, <em>coercive</em>, and <em>Lipschitz continuous</em>, ensuring the well-posedness of the mathematical model. The governing equation, derived from the balance of linear momentum coupled with the nonlinear constitutive relationship, is formulated as a second-order quasi-linear elliptic partial differential equation. For a body with an edge crack, this governing equation is augmented with a classical traction-free boundary condition on the crack faces. An <span><math><mrow><mi>h</mi><mi>p</mi></mrow></math></span>-adaptive finite element scheme is proposed for the numerical approximation of the resulting boundary value problem. The adaptive strategy is driven by a dual-component error estimation scheme: mesh refinement (<span><math><mi>h</mi></math></span>-adaptivity) is guided by a residual-based a posteriori error indicator of the <em>Kelly type</em>, while the local polynomial degree (<span><math><mi>p</mi></math></span>-adaptivity) is adjusted based on an estimator of the local solution regularity. The performance, accuracy, and convergence characteristics of the proposed method are demonstrated through numerical experiments. The structure of the regularized crack-tip fields is examined for various modeling parameters. Furthermore, the presented framework establishes a robust foundation for extension to more complex and computationally demanding problems, including quasi-static and dynamic crack propagation in brittle materials. A comparative analysis reveals that while standard <span><math><mi>h</mi></math></span>-refinement yields only algebraic convergence due to the crack-tip singularity, the proposed <span><math><mrow><mi>h</mi><mi>p</mi></mrow></math></span>-scheme recovers exponential convergence, reducing the degrees of freedom required for high-precision solutions by nearly two orders of magnitude.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104520"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145957280","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 fast and robust third medium contact approach using the neighbored element method","authors":"Max von Zabiensky,&nbsp;Dustin R. Jantos,&nbsp;Philipp Junker","doi":"10.1016/j.finel.2025.104489","DOIUrl":"10.1016/j.finel.2025.104489","url":null,"abstract":"<div><div>In several engineering applications, self-contact is a major effect or even desired feature, e.g., energy absorption in severe structural deformations, gripper in (soft-)robotics, and the optimization of those. Usual contact approaches require extensive calculation effort, but a fast and robust simulation approach could lead to remarkable improvements in product design. The third medium contact is a very promising approach, which can fulfill the requirements mentioned above. This method was also already applied in topology optimization approaches. However, the most concurrent regularization approaches require at least quadratic shape functions for hexahedron finite elements and at least cubic shape functions for tetrahedrons. In recent works, a new regularization technique was developed: herein, the gradients of deformation measures are computed by additional nodal degrees of freedom in the finite elements, i.e., as additional fields of unknowns. This is accompanied by a reduction of the order of the shape function. In this work, we derive a special discretization of the additional unknowns which enables self-contact based on the neighbored element method. This approach leads to a fast algorithm for the computation of the regularized set of partial differential equations: the displacements are the only nodal degrees of freedom and can be computed with linear finite element shape functions while the additional unknowns are discretized in the quadrature points of the mesh and solved in a staggered manner. The results are critically analyzed and illustrated for several two-dimensional problems.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104489"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146047949","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":"Formal asymptotic derivation of one-dimensional models for thin-walled beams based on shell theory","authors":"Jun-Sik Kim ,&nbsp;Tuan Anh Bui ,&nbsp;Junyoung Park","doi":"10.1016/j.finel.2026.104511","DOIUrl":"10.1016/j.finel.2026.104511","url":null,"abstract":"<div><div>Thin-walled beams with complex cross-sections require efficient reduced-order models to overcome the high computational cost of three-dimensional finite element analyses. Among existing approaches, the Formal Asymptotic Method (FAM) provides an efficient framework for slender beams through a two-step procedure: a cross-sectional analysis to determine warping-related quantities, followed by the construction of an equivalent one-dimensional beam model within a Timoshenko-type framework. The resulting macroscopic model involves only six degrees of freedom per node and is therefore very compact.</div><div>In conventional FAM formulations, the cross-sectional analysis is based on three-dimensional solid elasticity and requires discretization using two-dimensional finite elements, which can be computationally expensive, especially for thin-walled and multilayer composite structures. In this paper, a new one-dimensional beam formulation is proposed by integrating shell theory into the FAM framework. The cross-section is represented by its reference surface and discretized using two-node line elements, leading to a significant reduction in the computational cost of the sectional analysis while preserving the essential deformation characteristics of thin to moderately thick walls.</div><div>The accuracy and efficiency of the proposed formulation are demonstrated through numerical examples, including rectangular box beams and a wind turbine blade. Comparisons with other reduced-order models and with three-dimensional finite element results obtained using Abaqus show that the proposed approach accurately predicts global displacements at a significantly lower computational cost.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"255 ","pages":"Article 104511"},"PeriodicalIF":3.5,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145957281","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}