Finite Elements in Analysis and Design最新文献

{"title":"A Hybrid finite element implementation of two-potential constitutive model of dielectric elastomers","authors":"Kamalendu Ghosh ,&nbsp;Bhavesh Shrimali","doi":"10.1016/j.finel.2025.104348","DOIUrl":"10.1016/j.finel.2025.104348","url":null,"abstract":"<div><div>There has been an increasing interest in the constitutive modeling of dielectric elastomers due to their potential in enabling new technologies such as soft robotics, actuators and haptic devices. Under realistic time-dependent loadings, dielectric elastomers are inherently dissipative. They dissipate energy both through viscous deformation and through friction in their electric polarization process. However, a majority of constitutive models and their corresponding Finite Element (FE) implementations consider only mechanical dissipation. The main reason for this bias is that the mechanical relaxation time of dielectric elastomers is much larger than their electric relaxation time. However, accounting for electric dissipation, in addition to mechanical dissipation, is crucial when dealing with applied alternating electric fields. A fully coupled 3-D constitutive model for isotropic and incompressible dielectric elastomers was proposed by Ghosh and Lopez-Pamies (2021). In this paper, we critically investigate the numerical scheme proposed in this paper to solve the initial boundary value problem (IBVP) that describes the time-dependent behavior of dielectric elastomers. We find that the scheme in Ghosh and Lopez-Pamies (2021), employing a fifth-order explicit Runge–Kutta time discretization, may lead to excessively small or nonphysical time steps for IBVPs simulating the behavior of real-world elastomers. This is because of the stark contrast in the relaxation times of mechanical dissipation and electric polarization. To this end, we first present a stable implicit time-integration algorithm that overcomes the unrealistic time-step constraints imposed by the fifth-order explicit Runge–Kutta algorithm in Ghosh and Lopez-Pamies (2021). We then deploy the algorithm along with a conforming FE discretization to solve the IBVP. We present implementations of the mixed-FE formulation of the governing equations for dielectric elastomers in <span>FEniCSx</span>. We also show that the numerical scheme is robust, accurate, capable of handling finite deformations, the incompressibility constraint of the rubber, and general time-dependent loading conditions. In the last part, the FE code is deployed to validate the experimental findings describing the electromechanical behavior of VHB 4910 (from 3M) under a complex time-dependent electromechanical load as studied in Hossain et al. (2015).</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"247 ","pages":"Article 104348"},"PeriodicalIF":3.5,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143759997","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":"Explicit dynamics and buckling simulations with 7-p shell elements and enhanced assumed strain","authors":"Anh-Khoa Chau,&nbsp;Michael Brun,&nbsp;Pascal Ventura,&nbsp;Hamid Zahrouni,&nbsp;Michel Potier-Ferry","doi":"10.1016/j.finel.2025.104346","DOIUrl":"10.1016/j.finel.2025.104346","url":null,"abstract":"<div><div>Explicit strategies for shell dynamics are presented using 7-parameter shell elements and the Central Difference scheme. The formulation of the 7-parameter shell element is based on the widely used Enhanced Assumed Strain (EAS), allowing the use of a 3D constitutive law in the shell element without the need to condense the transverse normal stress component in the material law. The 7-parameter shell elements were mainly employed in the context of non-linear quasi-static loading and implicit dynamics for reproducing buckling phenomena. Explicit dynamics is the focus of this work, which proposes different strategies to handle the EAS field. In addition, kinematic constraints at the intersections between shell components are prescribed in terms of velocity. The critical time step size for thin-shell structures modeled with 7-parameter shell elements is increased thanks to the Selective Mass Scaling technique (SMS). The relevance of the proposed approaches is based on the ability to conserve momenta and energy and reproduce complex dynamic buckling phenomena. Numerical applications include classical benchmark tests for assessing the relevance of momentum-energy conserving time integration schemes: the free fly of a toss rule and the three intersecting plates. Buckling phenomena are also investigated for a roof cylindrical shell and a closed cylinder under follower external pressure. The proposed non-linear explicit dynamic strategies are attractive due to their enhanced capability to conserve momenta and energy and efficient prediction of dynamic buckling phenomena.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"247 ","pages":"Article 104346"},"PeriodicalIF":3.5,"publicationDate":"2025-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143714995","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Four-point contact slewing bearing dynamics. Guidelines for FE modelling and mechanistic model correlation","authors":"Martin Eizmendi,&nbsp;Iker Heras,&nbsp;Mikel Abasolo,&nbsp;Josu Aguirrebeitia","doi":"10.1016/j.finel.2025.104347","DOIUrl":"10.1016/j.finel.2025.104347","url":null,"abstract":"<div><div>The vibrational response of mechanical systems including four-point contact slewing bearings is heavily influenced by the stiffness and damping properties of the bearing joint itself. As an initial approach to study the dynamic response of these components, in this work several aspects are addressed. With a view toward dynamic modelling, first, a FE-based modification of the force-deflection Hertz formula is proposed to simulate more accurately the ball-raceway contact, including the effect of the conformity, and thus providing a more accurate formula that can be used also for solving static load distribution problems. Then, novel guidelines regarding the dynamic FE modelling of these components are provided for two complexity levels: an accurate one and a simplified one. Simultaneously, a mechanistic model to simulate the dynamic response of these bearings under axial loads is proposed. This latter model is validated against an ad-hoc developed FE model, demonstrating its efficiency. Finally, the energy dissipation due to material hysteresis is implemented into the mechanistic model through different damping models, and their performance is compared with results provided by a FE accurate model with a view toward implementing them in FE simplified modelling techniques for future experimental correlation.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"247 ","pages":"Article 104347"},"PeriodicalIF":3.5,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679916","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 mixed finite element method based on the volume coordinate system for stress analysis of plates","authors":"Jintao Zhou,&nbsp;Guanghui Qing","doi":"10.1016/j.finel.2025.104332","DOIUrl":"10.1016/j.finel.2025.104332","url":null,"abstract":"<div><div>Traditional bilinear isoparametric coordinate systems exhibit sensitivity to mesh distortion due to their fully high-order polynomials being only equivalent to first-order polynomials in Cartesian coordinate systems when confronted with mesh distortion. This paper combines the concept of 3- and 6-component volume coordinate systems (VCS) with the generalized mixed element method to develop a novel bivariate method called the non-conforming generalized mixed element method in the volume coordinate system (NGMVC). Established a bivariate field analysis mode in the VCS. Taking VCS as local coordinates significantly improves the morbidity relationship between local and Cartesian coordinate systems in conventional isoparametric elements during mesh distortion. Also avoids the calculation of the Jacobian inverse in the element strain matrix, greatly simplifies mathematical expressions, and lowers computational costs while ensuring that elements are insensitive to mesh distortion. On the other hand, in the analyzing procedure, the NGMVC describes the finite element model more objectively and rationally by considering both stress and displacement boundary conditions simultaneously. Thus, addressing the limitation of traditional displacement methods lacking consideration of stress boundary conditions. Based on these, considering the objective fact that the in-plane stress in composite laminated structures may not be continuous between layers, the non-conforming generalized partial mixed method in the volume coordinates (NGPMVC) was established by separating in-plane stress and out-of-plane stress in the mixed element formulation. The proposed method was verified through benchmark problems for laminates. The numerical results demonstrate that the method is not sensitive to mesh distortion and has a good ability to capture each stress component for different mesh densities, aspect ratios, and geometric structures.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"247 ","pages":"Article 104332"},"PeriodicalIF":3.5,"publicationDate":"2025-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143521129","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":"Simple finite element algorithm for solving antiplane problems with Gurtin–Murdoch material surfaces","authors":"María A. Herrera-Garrido ,&nbsp;Sofia G. Mogilevskaya ,&nbsp;Vladislav Mantič","doi":"10.1016/j.finel.2025.104318","DOIUrl":"10.1016/j.finel.2025.104318","url":null,"abstract":"<div><div>The finite element algorithm is developed to solve antiplane problems involving elastic domains whose boundaries or their parts are coated with thin and relatively stiff layers. These layers are modeled by the vanishing thickness Gurtin–Murdoch material surfaces that could be open or closed, and smooth or non-smooth. The governing equations for the problems are derived using variational arguments. The domains are discretized using triangular finite elements. In general, standard linear elements are used to approximate displacements in the domain. However, to capture the singular behavior of the elastic fields near the tips of the open Gurtin–Murdoch surfaces, a novel blended singular element is devised. Numerical examples are presented to demonstrate the accuracy and robustness of the algorithm developed.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104318"},"PeriodicalIF":3.5,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143487296","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":"Numerical modelling of shear cutting in complex phase high strength steel sheets: A comprehensive study using the Particle Finite Element Method","authors":"Olle Sandin ,&nbsp;Patrick Larour ,&nbsp;Juan Manuel Rodríguez ,&nbsp;Sergi Parareda ,&nbsp;Samuel Hammarberg ,&nbsp;Jörgen Kajberg ,&nbsp;Daniel Casellas","doi":"10.1016/j.finel.2025.104331","DOIUrl":"10.1016/j.finel.2025.104331","url":null,"abstract":"<div><div>The study examines the shear cutting process of Advanced High Strength Steel using the Particle Finite Element Method. Shear cutting, a crucial process in sheet metal forming, often leads to microcracks and plastic deformation that degrades the material performance in subsequent applications, such as cold forming, crashworthiness, and fatigue resistance. This work utilises the Particle Finite Element Method as an alternative to conventional Finite Element Methods to address the challenges of large deformation solid mechanics, offering high predictive accuracy in localised shearing deformation and fracture. The model was validated against experimental data from sheet punching tests, with evaluations at both macroscopic and mesoscopic levels, including cut edge profiles and microstructural deformation within the shear-affected zone. The Particle Finite Element Method approach demonstrated a high level of accuracy in predicting cut edge shape and shear-induced damage across various cutting conditions. As an unconventional numerical technique, usage of the Particle Finite Element Method advances modelling of large deformations solid mechanics and providing a robust tool for optimising manufacturing processes of materials sensitive to sheared edge damage.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104331"},"PeriodicalIF":3.5,"publicationDate":"2025-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143479354","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Application of zonal Reduced-Order-Modeling to tire rolling simulation","authors":"D. Danan ,&nbsp;R. Meunier ,&nbsp;T. Dairay ,&nbsp;T. Homolle ,&nbsp;M. Yagoubi","doi":"10.1016/j.finel.2025.104330","DOIUrl":"10.1016/j.finel.2025.104330","url":null,"abstract":"<div><div>Physic-based simulation remains a key enabler for real-world ever-growing complex industrial systems especially when crucial decisions are needed. While classical approaches have proven their accuracy and robustness over the years and come with a rich mathematical foundation, they suffer from several limitations depending of the underlying physics and use cases. For instance, especially concerning the resolution of Partial Differential Equations (PDEs) in 3 dimensions (3D), classical approaches are known to be computationally expensive. However, it turns out that simple pure data-driven approaches, while allegedly much more efficient from a computational point of view, do not necessarily hold up well regarding physical considerations. In this work, our aim is to investigate the tradeoff between accuracy and computational cost to design efficient and robust physical simulation methods under industrial constraints. In particular, as it is not easy to generate a large dataset through numerical simulations for such a problem, our aim is to design an approach addressing the data scarcity issue. To do so, we propose to hybridize a standard Finite Element Method (FEM) physics-based solver with a zonal Reduced Order Model (ROM) approach to simulate a rolling tire.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104330"},"PeriodicalIF":3.5,"publicationDate":"2025-02-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143471709","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":"An enhanced single Gaussian point continuum finite element formulation using automatic differentiation","authors":"Njomza Pacolli ,&nbsp;Ahmad Awad ,&nbsp;Jannick Kehls ,&nbsp;Bjorn Sauren ,&nbsp;Sven Klinkel ,&nbsp;Stefanie Reese ,&nbsp;Hagen Holthusen","doi":"10.1016/j.finel.2025.104329","DOIUrl":"10.1016/j.finel.2025.104329","url":null,"abstract":"<div><div>This contribution presents an improved low-order 3D finite element formulation with hourglass stabilization using automatic differentiation (AD). Here, the former Q1STc formulation is enhanced by an approximation-free computation of the inverse Jacobian. To this end, AD tools automate the computation and allow a direct evaluation of the inverse Jacobian, bypassing the need for a Taylor series expansion. Thus, the enhanced version, Q1STc+, is introduced. Numerical examples are conducted to compare the performance of both element formulations for finite strain applications, with particular focus on distorted meshes. Moreover, the performance of the new element formulation for an elasto-plastic material is investigated. To validate the obtained results, a volumetric locking-free element based on scaled boundary parametrization is used. Both the implementation of the element routine Q1STc+ and the corresponding material subroutine are made accessible to the public at <span><span>https://doi.org/10.5281/zenodo.14259791</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104329"},"PeriodicalIF":3.5,"publicationDate":"2025-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143444282","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Robust multi-physical-material topology optimization with thermal-self-weight uncertain loads","authors":"Minh-Ngoc Nguyen ,&nbsp;Joowon Kang ,&nbsp;Soomi Shin ,&nbsp;Dongkyu Lee","doi":"10.1016/j.finel.2025.104319","DOIUrl":"10.1016/j.finel.2025.104319","url":null,"abstract":"<div><div>Most topology optimization techniques for enhanced designs rely on the premise of deterministic loads. Nevertheless, in actuality, variables such as placements, weights, and orientations of applied loads may inadvertently fluctuate. Deterministic load-based designs may exhibit suboptimal structural performance in the presence of loading uncertainties. Uncertain aspects must be considered in topological optimization to provide robust outcomes. This work introduces an innovative robust multi-physics topology optimization technique for the design of multi-materials in response to unforeseen load variations. A combination of thermo-mechanical and self-weight loads, along with loading uncertainties, is provided based on the extended SIMP technique to achieve resilient designs. The optimized structures can be concurrently refined by minimizing the weighted sum of predicted compliance and standard deviation. The impact of self-weight and heat loads is examined through various cases to validate the proposed strategy.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104319"},"PeriodicalIF":3.5,"publicationDate":"2025-02-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143395080","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":"An assumed enhanced strain finite element formulation for modeling hydraulic fracture growth in a thermoporoelastic medium","authors":"Fushen Liu","doi":"10.1016/j.finel.2025.104320","DOIUrl":"10.1016/j.finel.2025.104320","url":null,"abstract":"<div><div>This paper presents an assumed enhanced strain finite element framework for simulating hydraulic fracture propagation in saturated thermoporoelastic media, considering the influence of thermal effects. The proposed approach combines classical thermoporoelasticity theory with a cohesive fracture model to describe the coupled behaviors of fluid flow, rock deformation and fracture propagation. Within this framework, fractures are represented using constant strain triangular elements enriched with constant displacement jumps. The mechanical response of fractures is governed by a trilinear cohesive law, and fracture initiation and propagation are both determined by using standard Newton’s method while maintaining global equilibrium. The numerical framework is verified through a series of examples, including cases without fractures, cases with rigid and deformable fractures, and hydraulic fracture propagation with thermal effects. The results show that thermal stress primarily affects the region near the injection point but has limited impact on fracture length evolution and fluid pressure distribution within the fracture. In contrast, temperature-dependent viscosity can significantly influence hydraulic fracture propagation. This work can be beneficial to our understanding of hydraulic fracture modeling in thermoporoelastic media and provide a potential useful numerical tool for simulating hydraulic fracturing processes with consideration of thermal effects.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"246 ","pages":"Article 104320"},"PeriodicalIF":3.5,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143387536","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}