International Journal for Numerical Methods in Engineering最新文献

{"title":"A Novel Strain-Continuous Finite Element Formulation for Geometrically Exact Thin-Walled Beam on Lie Algebra","authors":"Ziheng Huang, Ju Chen, Shixing Liu, Yongxin Guo","doi":"10.1002/nme.70024","DOIUrl":"https://doi.org/10.1002/nme.70024","url":null,"abstract":"<div>\u0000 \u0000 <p>Most existing geometrically exact thin-walled beam formulations on Lie group SE(3)<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>×</mo>\u0000 </mrow>\u0000 <annotation>$$ times $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℝ</mi>\u0000 </mrow>\u0000 <annotation>$$ mathbb{R} $$</annotation>\u0000 </semantics></math> considering warping are <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <annotation>$$ C $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mo> </mo>\u0000 <mrow>\u0000 <mn>0</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {}^0 $$</annotation>\u0000 </semantics></math>-continuous. In this study, a novel strain-continuous element for geometrically exact thin-walled beam on Lie algebra <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>s</mi>\u0000 </mrow>\u0000 <annotation>$$ mathfrak{s} $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>e</mi>\u0000 </mrow>\u0000 <annotation>$$ mathfrak{e} $$</annotation>\u0000 </semantics></math>(3)<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>×</mo>\u0000 </mrow>\u0000 <annotation>$$ times $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℝ</mi>\u0000 </mrow>\u0000 <annotation>$$ mathbb{R} $$</annotation>\u0000 </semantics></math> is originally proposed, in which the torsion-related Wagner effect and warping are considered in the constitutive relations. The proposed beam element is not only locking-free intrinsically, but also is <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <annotation>$$ C $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mo> </mo>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {}^1 $$</annotation>\u0000 </semantics></math>-continuous by using the proposed geometrical Hermite interpolation on Lie algebra. This interpolation based on Magnus expansion can simplify the description of interpolated element stress tensor and strain energy. Then, the sim","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852803","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 Deep Learning-Based Surrogate Model for Seismic Data Assimilation in Fault Activation Modeling","authors":"Caterina Millevoi,&nbsp;Claudia Zoccarato,&nbsp;Massimiliano Ferronato","doi":"10.1002/nme.70040","DOIUrl":"https://doi.org/10.1002/nme.70040","url":null,"abstract":"<p>Assessing the safety and environmental impacts of subsurface resource exploitation and management is critical and requires robust geomechanical modeling. However, uncertainties stemming from model assumptions, intrinsic variability of governing parameters, and data errors challenge the reliability of predictions. In the absence of direct measurements, inverse modeling and stochastic data assimilation methods can offer reliable solutions, but in complex and large-scale settings, the computational expense can become prohibitive. To address these challenges, this paper presents a deep learning-based surrogate model (SurMoDeL) designed for seismic data assimilation in fault activation modeling. The surrogate model leverages neural networks to provide simplified yet accurate representations of complex geophysical systems, enabling faster simulations and analyses essential for uncertainty quantification. The work proposes two different methods to integrate an understanding of fault behavior into the model, thereby enhancing the accuracy of its predictions. The application of the proxy model to integrate seismic data through effective data assimilation techniques efficiently constrains the uncertain parameters, thus bridging the gap between theoretical models and real-world observations.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70040","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852802","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":"A Chemo-Damage-Mechanical Coupled Phase-Field Model for Three-Dimensional Hydrogen-Assisted Dynamic Cracking","authors":"Hui Li,&nbsp;Shanyong Wang","doi":"10.1002/nme.70038","DOIUrl":"https://doi.org/10.1002/nme.70038","url":null,"abstract":"<p>This study develops a chemo-damage-mechanical coupled phase-field method for modeling two-dimensional and/or three-dimensional hydrogen-assisted transient dynamic cracking in metallic materials. In this method, hydrogen diffusion in solids is described by the evolution of bulk hydrogen concentration governed by the diffusion equation with an extended Fick's law. The hydrogen concentration and the inertial force of solids are incorporated into the governing equations of a non-standard quasi-brittle phase-field model (known as the phase-field regularized cohesive zone model, PF-CZM) capable of modeling complicated multiple crack nucleation, initiation, and propagation with insensitivity to mesh size and length scale. The resultant displacement-damage-hydrogen concentration coupled three-field equation is derived by the finite element Galerkin method and solved using a staggered Newton–Raphson iteration algorithm with an unconditionally stable implicit Newmark integration scheme. The new method was verified by four benchmark examples with hydrogen-free numerical/experimental results for comparison, including quasi-static/dynamic fracture of a notched plate under uniaxial tension, a dynamic crack branching experiment under constant traction, the Kalthoff-Winkler impact fracture experiment, and dynamic fragmentation of a cylinder under internal pressures, with the effects of loading velocity, initial bulk hydrogen concentration, and diffusion time on crack propagation investigated in detail. It is found that the present method is capable of modeling complex 2D and 3D hydrogen-assisted dynamic crack propagation and bifurcation under impact or internal pressure loadings, and thus holds the potential to be used for the structural design of hydrogen storage.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70038","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852865","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":"Investigation of Using High-Order Discontinuous Galerkin Methods in Adjoint Gradient-Based 3D Aerodynamic Shape Optimization","authors":"Yiwei Feng,&nbsp;Lili Lv,&nbsp;Tiegang Liu,&nbsp;Kun Wang,&nbsp;Bangcheng Ai","doi":"10.1002/nme.70033","DOIUrl":"https://doi.org/10.1002/nme.70033","url":null,"abstract":"<div>\u0000 \u0000 <p>Based on recent work by He et al. (<i>Computer Physics Communications</i> <b>286</b> (2023): 108660), this work develops a robust and efficient workflow for 3D aerodynamic shape optimization (ASO) based on the discontinuous Galerkin method (DGM) with solution remapping technique. It is theoretically found and numerically validated through 2D cases that the DGMs can provide more precise adjoint-enabled gradients even on a coarse mesh as compared with the FVMs under coequal computational costs. A number of 2D and 3D intricate cases are tested to illustrate the potential advantages of the DGM-based ASO workflow in terms of computational efficiency and final optimization results. The results demonstrate that the solution remapping technique can save around <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>50</mn>\u0000 <mo>%</mo>\u0000 </mrow>\u0000 <annotation>$$ 50% $$</annotation>\u0000 </semantics></math> of the computational time. Moreover, the DGM-based workflow extends the capability to explore the design space, leading to aerodynamic shapes with superior performance.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852805","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":"Adaptive Coupling of Peridynamic and Classical Continuum Mechanical Models Driven by Broken Bond/Strength Criteria for Structural Dynamic Failure","authors":"JiuYi Li,&nbsp;ShanKun Liu,&nbsp;Fei Han,&nbsp;Yong Mei,&nbsp;YunHou Sun,&nbsp;FengJun Zhou","doi":"10.1002/nme.70021","DOIUrl":"https://doi.org/10.1002/nme.70021","url":null,"abstract":"<div>\u0000 \u0000 <p>Peridynamics (PD) is widely used to simulate structural failure. However, PD models are time consuming. To improve the computational efficiency, we developed an adaptive coupling model between PD and classical continuum mechanics (PD-CCM) based on the Morphing method, driven by the broken bond or strength criteria. We derived the dynamic equation of the coupled models from the Lagrangian equation and then the discretized finite element formulation. An adaptive coupling strategy was introduced by determining the key position using the broken bond or strength criteria. The PD subdomain was expanded by altering the value of the Morphing function around the key position. Additionally, the PD subdomain was meshed by discrete elements (DEs) (i.e., nodes were not shared between elements), allowing the crack to propagate freely along the boundary of the DE. The remaining subdomains were meshed by continuous elements (CEs). Following the PD subdomain expansion, the CEs were converted into DEs, and new nodes were inserted. The displacement vector and mass matrix were reconfigured to ensure calculation consistency throughout the solving process. Furthermore, the relationship between the expansion radius of the PD subdomain and the speed of crack propagation was also discussed. Finally, the effectiveness, efficiency, and accuracy of the proposed model were verified via three two-dimensional numerical examples.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143793300","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":"Novel Study on Strain Modes-Based Interval Damage Identification Methodology Utilizing Orthogonal Polynomials and Collocation Theories","authors":"Lei Wang,&nbsp;Lihan Cheng,&nbsp;Qinghe Shi","doi":"10.1002/nme.70032","DOIUrl":"https://doi.org/10.1002/nme.70032","url":null,"abstract":"<div>\u0000 \u0000 <p>The sensitivity of strain modes to local stiffness changes within a structure underscores their potential as robust indicators of damage, enhancing the efficacy of damage identification processes. This study establishes sensitivity matrices of natural frequencies and strain modes to damage parameters, laying the groundwork for a novel fusion index that integrates both metrics to assess structural damage extent. In order to quantify the impact of uncertainty information on the results of damage identification processes, a non-probabilistic structural damage identification method rooted in the collocation methodology is proposed in this study. In consideration of computational efficiency, a two-step damage identification strategy encompassing localization and quantification is proposed. Initially, damage localization is achieved through the dynamic fingerprints, followed by the quantification of the uncertainty of damage extent. The proposed methodology is validated through a detailed numerical example, illustrating that the fusion index outperforms individual indices in terms of accuracy and computational efficiency. The non-probabilistic structural damage identification method based on collocation methodology can identify the damage extent and uncertainty interval even under the influence of uncertain factors.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143786701","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 X-FFT Solver for Two-Dimensional Thermal Homogenization Problems","authors":"Flavia Gehrig,&nbsp;Matti Schneider","doi":"10.1002/nme.70022","DOIUrl":"https://doi.org/10.1002/nme.70022","url":null,"abstract":"<p>We introduce an approach to computational homogenization which unites the accuracy of interface-conforming finite elements (FEs) and the computational efficiency of methods based on the fast Fourier transform (FFT) for two-dimensional thermal conductivity problems. FFT-based computational homogenization methods have been shown to solve multiscale problems in solid mechanics effectively. However, the obtained local solution fields lack accuracy in the vicinity of material interfaces, and simple fixes typically interfere with the numerical efficiency of the solver. In the work at hand, we identify the extended finite element method (X-FEM) with modified absolute enrichment as a suitable candidate for an accurate discretization and design an associated fast Lippmann-Schwinger solver. We implement the concept for two-dimensional thermal conductivity and demonstrate the advantages of the approach with dedicated computational experiments.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70022","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143762221","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":"Plastic-Damage Thin-Layer Element Model for Seismic Failure Analysis of Concrete Dams","authors":"Yi-Xiang Qiu,&nbsp;Tian-Yu Zhou,&nbsp;Jin-Ting Wang,&nbsp;Jian-Wen Pan,&nbsp;Chu-Han Zhang","doi":"10.1002/nme.70030","DOIUrl":"https://doi.org/10.1002/nme.70030","url":null,"abstract":"<div>\u0000 \u0000 <p>Finite element methods based on the small deformation assumption have been widely used in the seismic response analysis of concrete dams. However, these methods are not effective in simulating the entire process of concrete dams transitioning from small deformation damage cracking to larger deformation collapse under extreme earthquake events. This paper proposes a plastic-damage thin-layer element model for analyzing large deformation failure of concrete dams under strong earthquakes. Using the equivalence principles, the smeared crack model is equivalently transformed into a layered separation format, which can consider both small and large deformations. Numerical verification through single thin-layer elements and Petersson's three-point bending beam demonstrates that the proposed model can comprehensively consider plastic deformation, damage, and stiffness changes in concrete. Furthermore, the seismic failure process of Koyna gravity dam is simulated using the proposed plastic-damage thin-layer element method. The results show that this model can realistically simulate the entire process from initial cracking to ultimate failure in concrete dams.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143749857","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":"Normalized Field Product Approach: A Parameter-Free Density Evaluation Method for Close-To-Binary Solutions in Topology Optimization With Embedded Length Scale","authors":"Nikhil Singh,&nbsp;Prabhat Kumar,&nbsp;Anupam Saxena","doi":"10.1002/nme.7673","DOIUrl":"https://doi.org/10.1002/nme.7673","url":null,"abstract":"<div>\u0000 \u0000 <p>This article provides a normalized field product approach for topology optimization to achieve close-to-binary optimal designs. The method uses a parameter-free density measure that enforces a specified minimum length scale on the solid phase, ensuring smooth and transition-free topologies. The density evaluation does not rely on weight functions; however, the associated density functions are required to confined between 0 and 1. The method combines the SIMP scheme with the introduced density function for material stiffness interpolation. The success and efficacy of the approach are demonstrated through the design of both two- and three-dimensional designs, including stiff structures and compliant mechanisms. The structure's compliance is minimized for the former, whereas the latter involves optimizing a multicriteria objective. The presented numerical examples consider different volume fractions, length scales, and density functions. The proposed method is also seamlessly extended with advanced elements for solving 3D problems. The optimized designs obtained are close to binary without any user intervention while satisfying the desired feature size on the solid phase.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143749858","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":"Deep Material Networks for Fiber Suspensions With Infinite Material Contrast","authors":"Benedikt Sterr,&nbsp;Sebastian Gajek,&nbsp;Andrew Hrymak,&nbsp;Matti Schneider,&nbsp;Thomas Böhlke","doi":"10.1002/nme.70014","DOIUrl":"https://doi.org/10.1002/nme.70014","url":null,"abstract":"<p>We extend the laminate based framework of direct deep material networks (DMNs) to treat suspensions of rigid fibers in a non-Newtonian solvent. To do so, we derive two-phase homogenization blocks that are capable of treating incompressible fluid phases and infinite material contrast. In particular, we leverage existing results for linear elastic laminates to identify closed form expressions for the linear homogenization functions of two-phase layered emulsions. To treat infinite material contrast, we rely on the repeated layering of two-phase layered emulsions in the form of coated layered materials. We derive necessary and sufficient conditions which ensure that the effective properties of coated layered materials with incompressible phases are non-singular, even if one of the phases is rigid. With the derived homogenization blocks and non-singularity conditions at hand, we present a novel DMN architecture, which we name the flexible DMN (FDMN) architecture. We build and train FDMNs to predict the effective stress response of shear-thinning fiber suspensions with a Cross-type matrix material. For 31 fiber orientation states, six load cases, and over a wide range of shear rates relevant to engineering processes, the FDMNs achieve validation errors below 4.31% when compared to direct numerical simulations with fast-Fourier-transform based computational techniques. Compared to a conventional machine learning approach introduced previously by the consortium of authors, FDMNs offer better accuracy at an increased computational cost for the considered material and flow scenarios.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70014","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143749856","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}