International Journal for Numerical Methods in Engineering最新文献

{"title":"A Thermo-Flow-Mechanics-Fracture Model Coupling a Phase-Field Interface Approach and Thermo-Fluid-Structure Interaction","authors":"Sanghyun Lee,&nbsp;Henry von Wahl,&nbsp;Thomas Wick","doi":"10.1002/nme.7646","DOIUrl":"https://doi.org/10.1002/nme.7646","url":null,"abstract":"<p>This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7646","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143121085","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":"Accelerated Non-linear Stability Analysis Based on Predictions From Data-Based Surrogate Models","authors":"Anika Strauß,&nbsp;Jonas Kneifl,&nbsp;Anton Tkachuk,&nbsp;Jörg Fehr,&nbsp;Manfred Bischoff","doi":"10.1002/nme.7649","DOIUrl":"https://doi.org/10.1002/nme.7649","url":null,"abstract":"<p>In many applications in computer-aided engineering, like parametric studies, structural optimization, or virtual material design, a large number of almost similar models must be simulated. Although the individual scenarios may differ only marginally in both space and time, the same amount of effort is invested in each new simulation, without taking into account the experience and knowledge gained in previous simulations. Therefore, we have developed a method that combines data-based Model Order Reduction (MOR) and reanalysis, exploiting knowledge from previous simulation runs to accelerate computations in multi-query contexts. While MOR allows reducing model fidelity in space and time without significantly deteriorating accuracy, reanalysis uses results from previous computations as a predictor or preconditioner. In particular, this method enables acceleration of the exact computation of critical points, such as limit and bifurcation points, by the method of extended systems for systems that depend on a set of design parameters, such as material or geometric properties. Such critical points are of utmost engineering significance due to the special characteristics of the structural behavior in their vicinity. Conventional reanalysis methods, like the fold line analysis, can be used to accelerate the computation of critical points of almost similar systems but are limited in their applicability. For the fold line analysis, only small parameter variations are possible as the algorithm may not converge to the correct solution or fail to converge elsewise. Moreover, this method is only suited to finding the first critical points of limit point problems. In contrast to that, our developed data-based “reduced model reanalysis” method overcomes these problems. Thus, a larger parameter space can be covered. The efficiency of this method is demonstrated for a couple of numerical examples, including standard and isogeometric finite element models.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7649","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143120729","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":"Bar Axial Force Correction Method for the Elastic Performance of Truss Materials With Periodic Units","authors":"Mengfei Shang,&nbsp;Panxu Sun,&nbsp;Dongwei Wang,&nbsp;Shuxia Wang","doi":"10.1002/nme.7642","DOIUrl":"https://doi.org/10.1002/nme.7642","url":null,"abstract":"<div>\u0000 \u0000 <p>A bar axial force correction method is proposed on the basis of the representative volume element (RVE) method in this paper, which can predict the elastic performance of truss materials with periodic units. Based on the bar axial force correction method, the problem that the RVE method that node reaction forces on the common boundary may be unbalanced can be solved. For truss materials with periodic units, the effective elastic parameters and nodal displacement responses based on the proposed method are completely consistent those for the asymptotic homogenization (AH) method. Meanwhile, the proposed method is simpler than the AH method. Furthermore, compared with the axial compression experiment, the elastic modulus error of the proposed method is 3.81%. Thus, the correctness of this method is proved. Finally, the comparisons of different homogenization methods are analyzed in the numerical examples. The numerical examples show that the calculation results of the proposed method are the same as those of the AH method. Compared with the AH method, the maximum error of effective elastic parameters and nodal displacements for the RVE method are 22.51% and 13.65%, respectively. Compared with the AH method, the maximum error of effective elastic parameters and nodal displacements for the two-node (TN) method are 14.43% and 10.60%, respectively. Therefore, the bar axial force correction method has a wider range of applications and higher computational accuracy.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119322","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":"Regularized X-FEM Modeling of Arbitrary 3D Interacting Crack Networks","authors":"Endel V. Iarve,&nbsp;Eric Zhou,&nbsp;M. Keith Ballard,&nbsp;Zhenjia Gao,&nbsp;Hari K. Adluru,&nbsp;David Mollenhauer","doi":"10.1002/nme.7653","DOIUrl":"https://doi.org/10.1002/nme.7653","url":null,"abstract":"<div>\u0000 \u0000 <p>An extension to the Regularized eXtended Finite Element Method (RX-FEM) is proposed that allows arbitrary 3D cracks through a novel hierarchical enrichment algorithm. The technique avoids geometric consideration of how a crack cuts elements or intersects other cracks. Each crack is described separately in terms of its sign distance function and regularized step function, which are only recorded for nodes in the region where the gradient of the regularized step function is non-zero. The algorithm creates a set of superimposed nodes, referred to as node twins, and elements, referred to as twinned elements, for each crack and determines the connectivity of the element twins using the involved crack indices. The displacement jump is calculated between each pair of element twins corresponding to the same crack, and a cohesive zone model (CZM) is formulated for each pair of twins to model crack opening. Following the theory for the novel method, several examples are presented that illustrate capabilities of the new method that the traditional RX-FEM formulation lacked. A quasi-2D example of offset crack propagation is considered and successfully compared with published results. Additionally, two 3D examples involving perpendicularly intersecting cracks are considered, illustrating the intersection of two crack fronts and correct partitioning of the domain into eight fragments due to three crossing cracks, respectively.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119494","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":"Surface Structured Quadrilateral Mesh Generation Based on Topology Consistent-Preserved Patch Segmentation","authors":"Haoxuan Zhang,&nbsp;Haisheng Li,&nbsp;Xiaoqun Wu,&nbsp;Nan Li","doi":"10.1002/nme.7644","DOIUrl":"https://doi.org/10.1002/nme.7644","url":null,"abstract":"<div>\u0000 \u0000 <p>Surface-structured mesh generation is an important part of the Computational Fluid Dynamics (CFD) preprocessing stage. The traditional method cannot automatically divide the 3D surface topology in complex structures. Thus, we propose a surface-structured quadrilateral mesh generation method based on topology-consistent-preserved patch segmentation. The core idea is to segment the complex 3D model into several simple parts according to mesh quality and map each part to the 2D parametric domain based on the conformal parameterization method. Then, we utilize pattern-based topology partitioning to divide the parametric domain into multiple quadrilateral subdomains, facilitating the generation of 2D structured quadrilateral meshes. By using the inverse mapping algorithm based on barycentric weights, the generated 2D structured mesh is inversely mapped back to the 3D space. Finally, we splice each part accurately according to the structured mesh distribution. Experimental results show that our proposed method can generate higher-quality structured quadrilateral meshes than previous methods without losing the mesh topology of the original model.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119051","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":"Time-Domain Spectral BFS Plate Element With Lobatto Basis for Wave Propagation Analysis","authors":"Hela Ambati,&nbsp;Sascha Eisenträger,&nbsp;Santosh Kapuria","doi":"10.1002/nme.7617","DOIUrl":"https://doi.org/10.1002/nme.7617","url":null,"abstract":"<div>\u0000 \u0000 <p>A computationally efficient spectral Kirchhoff plate element is presented for time-domain analysis of wave propagation at high frequencies in thin isotropic plates. It employs a <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {C}^1 $$</annotation>\u0000 </semantics></math>-continuous spectral interpolation based on the modified bi-Hermite polynomials using the Gauss–Lobatto–Legendre (GLL) points as a basis with selective collocation of rotational and twisting degrees of freedom (DOFs) at element edge and corner nodes. The lowest order version of the proposed element reduces to the classical Bogner–Fox–Schmit (BFS) element for Kirchhoff plates. The GLL basis allows diagonalisation of the mass matrix using the nodal quadrature technique, which enhances the computational efficiency. The numerical properties of the proposed element are comprehensively evaluated, including the conditioning of the system matrices. Moreover, the effect of employing different numerical integration schemes and nodal sets is examined in both static and free vibration analyses. The effectiveness of the proposed element in wave propagation problems is evaluated by comparing its performance to the converged solutions achieved using the BFS element with a very fine mesh. Results demonstrate that the current element, without and even with mass matrix diagonalisation delivers exceptional accuracy while also exhibiting faster convergence and enhanced computational efficiency than the existing Kirchhoff plate elements.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143363005","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 Electro-Elastic Coupling Model for Piezoelectric Composites Based on the Voronoi Cell Finite Element Method","authors":"Huan Li,&nbsp;Nan Yang,&nbsp;Ran Guo","doi":"10.1002/nme.7631","DOIUrl":"https://doi.org/10.1002/nme.7631","url":null,"abstract":"<div>\u0000 \u0000 <p>The Voronoi cell finite element method (VCFEM) has successfully characterized the linear elastic behavior of the composites. This study is dedicated to develop an electro-elastic coupling VCFEM model mimicking the fully-coupled electro-elastic behavior of piezoelectric composites. For fiber-reinforced piezoelectric composites considering interfacial cracks, the interface traction reciprocity, the interface charge density reciprocity on bonded interfaces and the interface traction-free, the interface charge density-free on debonded interfaces are comprised in the new assumed stress and electric displacement hybrid variational functional. The new variational functional is derived on the base of the element multifield energy functionals. Independent stress/electric displacement fields are respectively assumed within the two-phase material domain. Several numerical examples considering perfectly-bonded interface and partially cracked interface were used to demonstrate the accuracy of the proposed method by comparing the piezoelectric Voronoi element model results with those obtained by ABAQUS. Then this model is used to study the effect of several microscopic details, such as the property ratio of fiber to matrix, volume fraction, interfacial crack length and polarization direction on macroscopic equivalent physical and mechanical properties, as well as local stress/electric displacement fields. It is clear that the proposed model is suitable for analyzing piezoelectric composites containing many microstructures with bonded interface or debonded interface.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142861666","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 Improved Spectral Element Differential Method in Solving Nonlinear Thermoelastic Coupling Problems With Discontinuous Interfaces","authors":"Jianning Zhao,&nbsp;Dong Wei,&nbsp;Yuxi Wang,&nbsp;Donghuan Liu","doi":"10.1002/nme.7645","DOIUrl":"https://doi.org/10.1002/nme.7645","url":null,"abstract":"<div>\u0000 \u0000 <p>In this paper, the spectral element differential method (SEDM) is improved to solve the nonlinear thermoelastic coupling problems with interface thermal resistance and interface gap in composite structures. The utilization of both Lobatto and Chebyshev node sets in SEDM significantly enhances solution efficiency by replacing integration with direct differential. Moreover, for strongly nonlinear problems caused by thermal radiation, unknown terms are incorporated into the stiffness matrix, and the relaxation iteration technique is also employed, the convergence has been improved compared to traditional methods. Importantly, the element format of the SEDM for 3D problems with discontinuous interfaces is given specifically in this paper, and element-by-element loop assembly of stiffness matrices is realized. Numerical examples confirm the effectiveness of the present method in efficiently and accurately solving 2D and 3D problems with discontinuous interfaces. The present method not only achieves faster convergence than the traditional finite element method, but also attains higher accuracy with fewer degrees of freedom and shorter computational time. Compared to the spectral element method (SEM), the proposed method significantly reduces the computation time of the stiffness matrix. Furthermore, by employing the coupled SEM-SEDM approach, computational efficiency is enhanced while maintaining high precision in sensitive regions.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142861667","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 Implicit Integration Algorithms With Identical Second-Order Accuracy and Flexible Dissipation Control for First-Order Transient Problems","authors":"Jinze Li,&nbsp;Naigang Cui,&nbsp;Yijun Zhu,&nbsp;Kaiping Yu","doi":"10.1002/nme.7639","DOIUrl":"https://doi.org/10.1002/nme.7639","url":null,"abstract":"<div>\u0000 \u0000 <p>For the solutions of first-order transient problems with optimal efficiency, conventional single-step implicit methodologies, unaided by additional post-processing techniques, encounter limitations in concurrently attaining identical second-order accuracy and controllable numerical dissipation. Addressing this challenge, the present study contributes not only a comprehensive analytical framework for formulating implicit integration algorithms but also leverages the auxiliary variable and the composite sub-step technique to propose two distinct types of implicit integration algorithms. Each type is characterized by self-initiation, unconditional stability, identical second-order accuracy, controllable numerical dissipation, and zero-order overshoots. Recognizing that single-step implicit methods can be conceptualized as composite single-sub-step ones, both algorithm types achieve identical effective stiffness matrices, thereby reducing the computational effort for solving linear problems. The utilization of auxiliary variables endows the proposed single-step method with numerical attributes akin to established counterparts, while achieving identical second-order accuracy. Conversely, the adoption of the composite sub-step technique in the proposed two-sub-step methods surpasses the published algorithms by significantly reducing the error constants. This superiority persists when enforcing the same sub-step sizes and sub-step dissipation levels. Numerical simulations affirm that the proposed methods consistently outperform existing alternatives without incurring additional computational expenses.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362997","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Reduced Order Modelling of Fully Coupled Electro-Mechanical Systems Through Invariant Manifolds With Applications to Microstructures","authors":"Attilio Frangi,&nbsp;Alessio Colombo,&nbsp;Alessandra Vizzaccaro,&nbsp;Cyril Touzé","doi":"10.1002/nme.7641","DOIUrl":"https://doi.org/10.1002/nme.7641","url":null,"abstract":"<p>This article presents the first application of the direct parametrisation method for invariant manifolds to a fully coupled multiphysics problem involving the nonlinear vibrations of deformable structures subjected to an electrostatic field. The formulation proposed is intended for model order reduction of electrostatically actuated resonating Micro-Electro-Mechanical Systems (MEMS). The continuous problem is first rewritten in a manner that can be directly handled by the parametrisation method, which relies upon automated asymptotic expansions. A new mixed fully Lagrangian formulation is thus proposed, which contains only explicit polynomial nonlinearities, which is then discretised in the framework of finite element procedures. Validation is performed on the classical parallel plate configuration, where different formulations using either the general framework or an approximation of the electrostatic field due to the geometric configuration selected are compared. Reduced-order models along these formulations are also compared to full-order simulations operated with a time integration approach. Numerical results show a remarkable performance both in terms of accuracy and the wealth of nonlinear effects that can be accounted for. In particular, the transition from hardening to softening behaviour of the primary resonance while increasing the constant voltage component of the electric actuation is recovered. Secondary resonances leading to superharmonic and parametric resonances are also investigated with the reduced-order model.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7641","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362372","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}