{"title":"Continuum and computational modeling of surface effects in flexoelectric materials","authors":"Mònica Dingle , Irene Arias , David Codony","doi":"10.1016/j.cma.2025.117971","DOIUrl":"10.1016/j.cma.2025.117971","url":null,"abstract":"<div><div>In recent times, with the rise of nanoscale technologies, miniaturization of devices has prompted the need to study electromechanical phenomena at small scales. Most studies focus on the phenomena occurring at the bulk portion of the material, such as flexoelectricity, but neglect the effects that arise from the surfaces of the samples. Given the fact that, at such scales, surface-to-volume ratio is inherently large, surface effects cannot be ignored if the full and accurate description of the material’s response wants to be provided. In this work, we present a model that successfully integrates flexoelectricity and the effects of surfaces, and we properly derive the governing equations and boundary conditions for the boundary value problem. We also present a numerical approach in order to computationally solve it, converging at high-order optimal rates. In addition, we present an analytical 1D Euler–Bernoulli electromechanical beam model. Numerically, we find the presence of boundary layers in the transversal electric field across the beam thickness, which are not accounted for in the analytical 1D model. Finally, we find numerical solutions for geometrically-polarized flexoelectric lattice metamaterials, which have large area-to-volume ratios, giving rise to very relevant surface effects. This work emphasizes the importance of accounting for surface effects in modeling and design of flexoelectric devices, including geometrically-polarized metamaterials.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117971"},"PeriodicalIF":6.9,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143838562","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}
{"title":"A staggered grid shell particle method for shell structural damage subjected to underwater explosion","authors":"Jiasheng Li , Yong Liang , Zhixin Zeng , Xiong Zhang","doi":"10.1016/j.cma.2025.117996","DOIUrl":"10.1016/j.cma.2025.117996","url":null,"abstract":"<div><div>A novel staggered grid shell particle method (SGSPM) is proposed in this paper to model the shell structural damage subjected to underwater explosion. The material point method (MPM) is used to model the fluid in underwater explosion, and the solid shell material point method (SSMPM) is adopted to model the shell structures. A staggered grid scheme is employed to eliminate the cell crossing noise and improve the accuracy of fluid simulation, and a conversion algorithm is proposed to handle the dynamic fracture of shell structures. In addition, a local multi-mesh contact method is introduced into the staggered grid scheme for modeling the fluid–structure interaction. Several numerical examples, including full hemispherical shell, penetration of a thin plate, large deformation of a plate subjected to underwater explosion, fragmentation of a plate and structural damage of a ship hull subjected to contact underwater explosion, are studied by the proposed SGSPM, and the numerical results agree well with the data in the literature and experiments.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117996"},"PeriodicalIF":6.9,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143834451","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}
Juan C. Alzate Cobo , Xiang-Long Peng , Bai-Xiang Xu , Oliver Weeger
{"title":"A finite swelling 3D beam model with axial and radial diffusion","authors":"Juan C. Alzate Cobo , Xiang-Long Peng , Bai-Xiang Xu , Oliver Weeger","doi":"10.1016/j.cma.2025.117983","DOIUrl":"10.1016/j.cma.2025.117983","url":null,"abstract":"<div><div>We present a geometrically exact 3D beam model that incorporates axial and radial swelling strains, both small and large, resulting from a rotationally symmetric, thermal or chemical diffusion. Isogeometric collocation is employed to discretize both the mechanical momentum balances and the axis-symmetric, steady-state 2D diffusion equation along the beam. The resulting coupled nonlinear problem for displacements, rotations, and temperatures or concentrations is solved using a staggered scheme. The approach is further extended to include beam-to-beam interfaces and is therefore well suited for the simulation of lattice structures. The model and its discretization are validated against 3D continuum models in various numerical examples and prove to be both accurate and numerically efficient. The novelty of the presented method is twofold. First, it relates beam theory, and consequently small elastic strains, with large swelling deformation stemming from anisotropic diffusion phenomena. Second, it also provides insight into the implementation of isogeometric collocation for solving diffusion equations subject to large deformations. Ultimately, this novel finite swelling beam model can present the starting point for the efficient modeling of lattice structures under diffusion conditions, such as microstructured Li-ion electrodes or thermoelectric semiconductors.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117983"},"PeriodicalIF":6.9,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143834450","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}
{"title":"Predicting change, not states: An alternate framework for neural PDE surrogates","authors":"Anthony Zhou , Amir Barati Farimani","doi":"10.1016/j.cma.2025.117990","DOIUrl":"10.1016/j.cma.2025.117990","url":null,"abstract":"<div><div>Neural surrogates for partial differential equations (PDEs) have become popular due to their potential to quickly simulate physics. With a few exceptions, neural surrogates generally treat the forward evolution of time-dependent PDEs as a black box by directly predicting the next state. While this is a natural and easy framework for applying neural surrogates, it can be an over-simplified and rigid framework for predicting physics. In this work, we evaluate an alternate framework in which neural solvers predict the temporal derivative and an ODE integrator forwards the solution in time, which has little overhead and is broadly applicable across model architectures and PDEs. We find that by simply changing the training target and introducing numerical integration during inference, neural surrogates can gain accuracy and stability in finely-discretized regimes. Predicting temporal derivatives also allows models to not be constrained to a specific temporal discretization, allowing for flexible time-stepping during inference or training on higher-resolution PDE data. Lastly, we investigate why this framework can be beneficial and in what situations does it work well.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117990"},"PeriodicalIF":6.9,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143838577","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}
Sichen Dong, Lei Li, Tianyu Yuan, Xiaotan Yu, Pan Wang, Fusen Jia
{"title":"A random interval coupling-based active learning Kriging with meta-model importance sampling method for hybrid reliability analysis under small failure probability","authors":"Sichen Dong, Lei Li, Tianyu Yuan, Xiaotan Yu, Pan Wang, Fusen Jia","doi":"10.1016/j.cma.2025.117992","DOIUrl":"10.1016/j.cma.2025.117992","url":null,"abstract":"<div><div>In this study, a novel active learning method is proposed and combined with Meta-IS-AK for hybrid reliability analysis under small failure probability. Considering the proportion of responses falling into the failure domain, the interval failure degree is introduced to describe the probability of misjudging the state for random samples. The novel active learning method (IAD) is proposed to select valuable samples for updating Kriging model, considering the interval failure degree and the sample clustering. Additionally, a corresponding convergence criterion based on the similarity of the indicator functions in importance sampling samples is proposed to further enhance efficiency. The accuracy and superiority of the proposed method are validated through seven illustrative examples, accompanied by detailed explanations.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117992"},"PeriodicalIF":6.9,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143834449","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}
{"title":"Geometrically exact beam finite element with generalized B-spline interpolation on the special Euclidean group SE(3)","authors":"Kunming Ren , Tingting Yuan , Jinyang Liu","doi":"10.1016/j.cma.2025.117979","DOIUrl":"10.1016/j.cma.2025.117979","url":null,"abstract":"<div><div>This work aims to address the challenge of achieving continuity in beam element interpolation by introducing a geometrically exact beam finite element based on generalized B-spline interpolation on the special Euclidean group <span><math><mrow><mi>S</mi><mi>E</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math></span>. The beam’s configuration is represented within the <span><math><mrow><mi>S</mi><mi>E</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math></span> framework and interpolated using a generalized B-spline approach, enabling high-order interpolation with enhanced <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> continuity. The interpolation is implemented through an extension of the classical De Boor algorithm to Lie groups via the pyramid algorithm. A systematic method is proposed for computing derivatives and linearizations essential for finite element formulations. Both static and dynamic equilibrium equations are derived on <span><math><mrow><mi>S</mi><mi>E</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math></span>, and the finite element formulations are established accordingly. Numerical examples validate the proposed element, confirming its correctness and adherence to critical properties, including objectivity, path-independence, and the absence of locking. The combination of <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> continuity and the high degree of the interpolation substantially enhances the convergence performance. In particular, the degree-2 beam element achieves improved accuracy with sixth-order convergence rates when the degrees of freedom are relatively low. These characteristics make the proposed element highly suitable for high-accuracy simulations of beam structures undergoing large deformations and rotations.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117979"},"PeriodicalIF":6.9,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143829216","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}
{"title":"Computational model and simulation of the solid-state sintering process as a thermal treatment for metal extrusion additive manufacturing: Microstructural and Multiphysics approach","authors":"Judice Cumbunga , Saïd Abboudi , Dominique Chamoret","doi":"10.1016/j.cma.2025.117978","DOIUrl":"10.1016/j.cma.2025.117978","url":null,"abstract":"<div><div>A numerical model has been developed to simulate microstructure evolution in stainless steel 316L components produced via Metal Extrusion Additive Manufacturing (MExAM). The model integrates key phenomena, including heat conduction, mechanical fields influenced by gravity, and phase-field equations, offering a robust framework to understand and control changes in thermomechanical properties during pressureless solid-state sintering. Advanced numerical techniques, such as the Finite Element Method (FEM) and the Physics-based Preconditioned Jacobian-free Newton-Krylov Method, were employed to solve the complex nonlinear system efficiently. Validation against literature data demonstrated the model’s accuracy and reliability, while tests across varying particle sizes highlighted its adaptability. Simulation results underscore the model’s potential for optimizing sintered materials by providing detailed insights into microstructural, thermal, and mechanical behavior.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117978"},"PeriodicalIF":6.9,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143834452","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}
Oded Ovadia , Vivek Oommen , Adar Kahana , Ahmad Peyvan , Eli Turkel , George Em Karniadakis
{"title":"Real-time inference and extrapolation with Time-Conditioned UNet: Applications in hypersonic flows, incompressible flows, and global temperature forecasting","authors":"Oded Ovadia , Vivek Oommen , Adar Kahana , Ahmad Peyvan , Eli Turkel , George Em Karniadakis","doi":"10.1016/j.cma.2025.117982","DOIUrl":"10.1016/j.cma.2025.117982","url":null,"abstract":"<div><div>Neural Operators are fast and accurate surrogates for nonlinear mappings between functional spaces within training domains. Extrapolation beyond the training domain remains a grand challenge across all application areas. We present Time-Conditioned UNet (TC-UNet) as an operator learning method to solve time-dependent PDEs continuously in time without any temporal discretization, including in extrapolation scenarios. TC-UNet incorporates the temporal evolution of the PDE into its architecture by combining a parameter conditioning approach with the attention mechanism from the Transformer architecture. After training, TC-UNet makes real-time inferences on an arbitrary temporal grid. We demonstrate its extrapolation capability on a climate problem by estimating the global temperature for several years and also for inviscid hypersonic flow around a double cone. We propose different training strategies involving temporal bundling and sub-sampling. We demonstrate performance improvements for several benchmarks, performing extrapolation for long time intervals and zero-shot super-resolution time.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117982"},"PeriodicalIF":6.9,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143834453","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}
Fei Ren , Pei-Zhi Zhuang , Xiaohui Chen , Hai-Sui Yu , He Yang
{"title":"Physics-Informed Extreme Learning Machine (PIELM) for Stefan problems","authors":"Fei Ren , Pei-Zhi Zhuang , Xiaohui Chen , Hai-Sui Yu , He Yang","doi":"10.1016/j.cma.2025.118015","DOIUrl":"10.1016/j.cma.2025.118015","url":null,"abstract":"<div><div>Stefan problems describe heat transfer through a material undergoing phase change, and solving these problems poses a real challenge due to the existence of a time-dependent moving boundary at the phase change interface. We propose an efficient and reliable physics-informed extreme learning machine (PIELM) framework for solving Stefan problems, which is achieved by replacing deep neural networks in the widely used physics-informed neural network (PINN) with extreme learning machines (ELM). We use a dual-network structure to approximate the latent solution and the moving boundary by two separate ELM networks, and in each ELM we incorporate physical laws of governing equations as well as initial and boundary conditions. Then, determining ELM layer weights is transformed from minimising loss into solving a system of equations. These equations are nonlinear because of the moving boundary, and we tackle them using an iterative least-squares procedure. The feasibility and validity of the proposed PIELM framework are demonstrated by carrying out six numerical case studies. Compared to conventional PINN frameworks, it shows that our PIELM framework can significantly improve the accuracy and efficiency for solving Stefan problems, reducing relative <em>L</em><sub>2</sub> errors from the orders of 10<sup>–3</sup>∼10<sup>–5</sup> to 10<sup>–6</sup>∼ 10<sup>–8</sup>.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 118015"},"PeriodicalIF":6.9,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143829214","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}
{"title":"A wrinkling model for general hyperelastic materials based on tension field theory","authors":"H.M. Verhelst , M. Möller , J.H. Den Besten","doi":"10.1016/j.cma.2025.117955","DOIUrl":"10.1016/j.cma.2025.117955","url":null,"abstract":"<div><div>Wrinkling is the phenomenon of out-of-plane deformation patterns in thin walled structures, as a result of a local compressive (internal) loads in combination with a large membrane stiffness and a small but non-zero bending stiffness. Numerical modelling typically involves thin shell formulations. As the mesh resolution depends on the wrinkle wave lengths, the analysis can become computationally expensive for shorter ones. Implicitly modelling the wrinkles using a modified kinematic or constitutive relationship based on a taut, slack or wrinkled state derived from a so-called tension field, a simplification is introduced in order to reduce computational efforts. However, this model was restricted to linear elastic material models in previous works. Aiming to develop an implicit isogeometric wrinkling model for large strain and hyperelastic material applications, a modified deformation gradient has been assumed, which can be used for any strain energy density formulation. The model is an extension of a previously published model for linear elastic material behaviour and is generalised to other types of discretisation as well. The extension for hyperelastic materials requires the derivative of the material tensor, which can be computed numerically or derived analytically. The presented model relies on a combination of dynamic relaxation and a Newton–Raphson solver, because of divergence in early Newton–Raphson iterations as a result of a changing tension field, which is not included in the stress tensor variation. Using four benchmarks, the model performance is evaluated. Convergence with the expected order for Newton–Raphson iterations has been observed, provided a fixed tension field. The model accurately approximates the mean surface of a wrinkled membrane with a reduced number of degrees of freedom in comparison to a shell solution.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"441 ","pages":"Article 117955"},"PeriodicalIF":6.9,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143829215","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}